一个matlab高手十年编程经验的吐血推荐之线性规划.doc

%目标函数所对应的设计变量的系数，为求极小，故取原目标函数系数的相反数 c=[-1;-1]; %不等式约束 A=[-4 2;4 2;0 -2]; b=[-1;11;-1]; %设计变量的边界约束，无上界约束 lb=[0;0]; %均要求为整数变量 M=[1;2]; %判断是否整数的误差限 Tol=1e-8; %求解原问题松弛线性规划的最优解x和目标函数在x处的值fval [x,fval]=linprog(c,A,b,[],[],lb,[]) %求最优解整数解x1和目标函数在x处的值fval1，结果为原问题最优值的相反数 [x1,fval1]=intprog(c,A,b,[],[],lb,[],M,Tol) Optimization terminated. x = 1.5000 2.5000 fval = -4.0000 x1 = 2 1 fval1 = -3.0000 %目标函数所对应的设计变量的系数 c=ones(1,6); %不等式约束 A=[ 0 0 0 -1 -1 -2; 0 -1 -2 0 -1 0; -4 -2 -1 -2 0 0]; b=[-20;-40;-80]; %设计变量的边界约束 lb=zeros(1,6); %所有变量均为整数变量，故将所有序号组成向量M M=[1;2;3;4;5;6]; %判定为整数的误差限 Tol=1e-8; %求最优解x和目标函数值fval，并返回状态指示 [x,fval]=intprog(c,A,b,[],[],lb,[],M,Tol) x = 5 0 20 20 0 0 fval = 45.0000 %目标函数所对应的设计变量的系数 c=[1;3;1;1;3;2]; %不等式约束 A=[ 0 0 0 -1 -1 -2 0 -1 -2 0 -1 0 -4 -2 -1 -2 0 0]; b=[-20;-40;-80]; %设计变量的边界约束 lb=zeros(1,6); %所有变量均为整数变量，故将所有序号组成向量M M=[1;2;3;4;5;6]; %判定为整数的误差限 Tol=1e-8; %求最优解x和目标函数值fval，并返回状态指示 [x,fval]=intprog(c,A,b,[],[],lb,[],M,Tol) %求出此时所用钢管的数目 s=sum(x) x = 5 0 20 20 0 0 fval = 45.0000 s = 45 %目标函数所对应的设计变量的系数 f=ones(1,6); %不等式约束 A=[-1 -1 0 0 0 0; 0 -1 -1 0 0 0; 0 0 -1 -1 0 0; 0 0 0 -1 -1 0; 0 0 0 0 -1 -1; -1 0 0 0 0 -1]; b=[-120;-80;-60;-30;-50;-100]; %边界约束，由于无上界，故设置ub=[Inf;Inf;Inf;Inf;Inf;Inf]; lb=[0;0;0;0;0;0]; ub=[Inf;Inf;Inf;Inf;Inf;Inf]; %所有变量均为整数变量，故将所有序号组成向量M M=[1;2;3;4;5;6]; %判定为整数的误差限 Tol=1e-8; %求最优解x和目标函数值fval，并返回状态指示 [x,fval,exitflag]=intprog(f,A,b,[],[],lb,ub,M,Tol) x = 50 70 25 35 0 50 fval = 230.0000 exitflag = 1 %目标函数所对应的设计变量的系数 c=[-20;-6;-8;-9]; %不等式约束 A=[10 6 5 2; 7 2 2 4; 2 1 1 10; 0 -1 -1 1]; b=[19;11;12;0]; %用于intprog函数的设计变量的边界约束，由于为0-1规划，故设置上界为全1 lb=[0;0;0;0]; ub=[1;1;1;1]; %均要求为整数变量 M=1:4; %判断是否整数的误差限 Tol=1e-8; %调用intprog函数求最优解整数解x和目标函数在x处的值fval [x,fval,exitflag]=intprog(c,A,b,[],[],lb,ub,M,Tol) %调用bintprog函数求最优解整数解x1和目标函数在x1处的值fval1 [x1,fva