最优化问题的MATLAB求解（课堂PPT）.pptVIP

主要参考书： 高等应用数学问题的MATLAB求解 薛定宇 陈阳泉 著 清华大学出版社 精通MATLAB科学计算 王正林等编著 电子工业出版社 科学计算引论－基于MATLAB的数值分析 ［美］Shoichiro Nakamura 电子工业出版社 MATLAB与科学计算 王谟然 编著 电子工业出版社 谢谢！ 3、模型：min z=fX VLB≤X≤VUB 命令：[1] x=linprog（f，A，b，Aeq,beq, VLB，VUB） [2] x=linprog（f，A，b，Aeq,beq, VLB，VUB, X0） 注意：[1] 若没有等式约束: , 则令Aeq=[ ], beq=[ ]. [2]其中X0表示初始点 4、命令：[x,fval]=linprog(…) 返回最优解ｘ及ｘ处的目标函数值fval. 解: 编写M文件如下： c=[-7 -5]; A=[3 2; 4 6; 0 7]; b=[90;200;210]; Aeq=[]; beq=[]; vlb=[0,0]; vub=[inf,inf]; [x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub) 问题2解答 例3：求 编写M-脚本文件如下： f=[-2 -1 -4 -3 -1]; % 线性的目标函数的系数向量（行、列均可） A=[0 2 1 4 2;3 4 5 -1 -1];B=[54;62]; % 线性不等式约束的系数矩阵 C=[]; D=[]; % 线性等式约束的系数矩阵。 % 因题中无等式约束，所以该处用空矩阵占位。 xm=[0;0;3.32;0.678;2.57]; xM=[]; x0=[]; % 这三个向量写成行向量列向量均可 [x,fopt,flag,c]=linprog(f,A,B,C,D,xm,xM,x0) 执行结果： Optimization terminated. x = 19.7850 0.0000 3.3200 11.3850 2.5700 fopt = -89.5750 flag = 1 % 表明求解成功 c = iterations: 6 algorithm: large-scale: interior point cgiterations: 0 message: Optimization terminated. 练：若 x 要求满足 试写出x的下、上界向量xm和xM 。 xm=[-1;-inf;-inf;1;-inf]; xM=[inf;2;5;2;inf]; 解： Exiting: One or more of the residuals, duality gap, or total relative error has stalled: the dual appears to be infeasible (and the primal unbounded). (The primal residual TolFun=1.00e-008.) x = 1.0e+007 * 0.0000 0.0000 -0.7408 4.0745 -7.7785 fopt = -1.4816e+007 flag = -3 % 求解失败 练：求 编写M-脚本文件如下： f=[-3/4 150 -1/50 6]; % 求f (x)的最大值，即求-f (x)的最小值。 A=[1/4 -60 -1/50 9;1/2 -90 -1/50 3];B=[0;0]; % 线性不等式约束的系数矩阵 C=[]; D=[]; % 线性等式约束的系数矩阵。 xm=[-5;-5;-5;-5]; xM=[inf;inf;1;