%求解模糊線性規劃題目2
%化為幾個定系數的線性規劃，間接求出模糊規化最優解
%自變量（設計變量）x1x2
clear?all
%先求普通線性規劃：
%max：?3x1+4x2
%4x1+20x2<=4500;
%12x1+6.4x2<=4700;
%x1x2>=0;
%轉化為min：-3x1-4x2求解
g=[-3-4];%目標函數的系數
A=[4?20;12?6.4];%約束函數系數矩陣
b=[4500?4700];%約束函數邊界
lb=[0;0];%設計變量下邊界
[xfval0exitflagoutput]=linprog（gAb[][]lb）;%if?exitflag=1?找到最優解
f0=-fval0;%取反，得到man?3x1+4x2
d1=100;d2=100;%伸縮性指標
%轉化至普通線性規劃：
%max：?3x1+4x2
%4x1+20x2<=4600;
%12x1+6.4x2<=4800;
%x1x2>=0;
b1=[4500?4700]+[d1?d2];
[x1fval1exitflag1output1]=linprog（gAb1[][]lb）;%if?exitflag=1?找到最優解
f1=-fval1;%變為man?3x1+4x2
d0=f1-f0;%目標伸縮性指標
g2=[0?0?-1];
A2=[A（1:）/d1?1;A（2:）/d2?1;g/d0?1];
b2=[b（1）/d1+1?b（2）/d2+1?-f0/d0];
lb2=[0;0;0];
[x2fval2exitflag2outpu2t]=linprog（g2A2b2[][]lb2）;
lambda=-fval2;
f2=-[g?0]*x2;%目標函數最優值