应用MATLAB求解规划问题.pdfVIP

MATLAB6.0数学手册 176 第 5章 优化问题 5.1 线性规划问题 线性规划问题是目标函数和约束条件均为线性函数的问题，MATLAB6.0解决的线性规 划问题的标准形式为： min nRxxf  sub.to： bxA  beqxAeq  ubxlb  其中 f、x、b、beq、lb、ub为向量，A、Aeq为矩阵。 其它形式的线性规划问题都可经过适当变换化为此标准形式。 在 MATLAB6.0 版中，线性规划问题（Linear Programming）已用函数 linprog 取代了 MATLAB5.x版中的 lp函数。当然，由于版本的向下兼容性，一般说来，低版本中的函数在 6.0版中仍可使用。 函数 linprog 格式 x = linprog(f,A,b) %求 min f *x sub.to bxA  线性规划的最优解。 x = linprog(f,A,b,Aeq,beq) %等式约束 beqxAeq  ，若没有不等式约束 bxA  ，则 A=[ ]，b=[ ]。 x = linprog(f,A,b,Aeq,beq,lb,ub) %指定 x的范围 ubxlb  ，若没有等式约束 beqxAeq  ，则 Aeq=[ ]，beq=[ ] x = linprog(f,A,b,Aeq,beq,lb,ub,x0) %设置初值 x0 x = linprog(f,A,b,Aeq,beq,lb,ub,x0,options) % options为指定的优化参数 [x,fval] = linprog(…) % 返回目标函数最优值，即 fval= f *x。 [x,lambda,exitflag] = linprog(…) % lambda为解 x的 Lagrange乘子。 [x, lambda,fval,exitflag] = linprog(…) % exitflag为终止迭代的错误条件。 [x,fval, lambda,exitflag,output] = linprog(…) % output为关于优化的一些信息 说明 若 exitflag0表示函数收敛于解 x，exitflag=0表示超过函数估值或迭代的最大数 字，exitflag0表示函数不收敛于解 x；若 lambda=lower 表示下界 lb，lambda=upper表示上 界 ub，lambda=ineqlin表示不等式约束，lambda=eqlin表示等式约束，lambda中的非 0元素 表示对应的约束是有效约束；output=iterations表示迭代次数，output=algorithm表示使用的 运算规则，output=cgiterations表示 PCG迭代次数。 例 5-1 求下面的优化问题 min 321 x6x4x5  sub.to 20xxx 321  42x4x2x3 321  第 5章 优化问题 177 30x2x3 21  321 x0,x0,x0  解： f = [-5; -4; -6]; A = [1 -1 1;3 2 4;3 2 0]; b = [20; 42; 30]; lb = zeros(3,1); [x,fval,exitflag,output,lambda] = linprog(f,A,b,[],[],lb) 结果为： x = %最优解 0.0000 15.0000 3.0000 fval = %最优值 -78.0000 exitflag = %收敛 1 output = iterations: 6 %迭代次数 cgiterations: 0 algorithm: lipsol %所使用规则 lambda = ineqlin: [3x1 double] eqlin: [0x1 double] upper: [3x1 double] lower: [3x1 double] lambda.ineqlin ans = 0.0000 1.5000 0.5000 lambda.lower ans = 1.0000 0.0000 0.0000 表明：不等约束条件 2和 3以及第 1个下界是有效的 5.2 foptions函数 对于优化控制，MATLAB提供了 18个参数，这些参数的具体意义为： options(1)-参数显示控制（默认值为 0）。等于 1时显示一些结果。 options(2)-优化点 x的精度控制(默认值为 1e-4)。 options(3)-优化函数 F的精度控制(默认值为 1e-4)。 options(4)-违反约束的结束标准(默认值为 1e-6)。 options(5)-算法选择，不常用。 MATLAB6.0数学手册 178 options