模糊线性规划算法

药品最优加工方案 某药品加工生产甲、乙两种药品，甲种药品的利润是3万元/kg ，乙种药品的利润是4万元/kg，生产1kg甲种药品需要原料A略少于4kg，需要原料B约12kg。生产乙种药品需要原料A略多于20kg，需要原料B约6.4kg。现有原料A约4600kg，原料B约4800kg。如何安排甲、乙两种药品的产量使利润最大？ （假设原料A最少有4500kg，伸缩指标为100kg；原料B最少有4700kg，伸缩指标为100kg）

%求解模糊线性规划题目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;%目标函数最优值

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         932  2016-04-27 15:16  2_FLP.m