第九章最优化方法教程分析.docVIP

第九章 最优化方法 本章主要介绍线性规划、0-1规划、非线性规划等问题的MATLAB求解。 9.1 线性规划（Linear Programming，简写为LP）问题 线性规划问题MATLAB解决的线性规划问题的标准形式为：min 其中为向量，为矩阵。 其它形式的线性规划问题都可经过适当变换化为此标准形式。 在M线性规划问题函数linprog，其使用格式为： [x, fval, exitflag, output, lambda] = linprog(f, A, b, Aeq, beq, lb, ub) 输入部分：其中各符号对应线性规划问题标准形式中的向量和矩阵，如果约束条件中有缺少，则其相应位置用空矩阵[]代替。 输出部分：其中x为最优解，用列向量表示；fval为最优值；exitflag为退出标志，若exitflag表示函数exitflag=0表示超过的迭代最大，exitflag=-2，表示若output表明算法和迭代情况；lambda表示存储情况。 例1 下面的问题 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) 注意：由于该问题没有等式约束，所以输入格式中相应的位置用[]代替，变量没有上限约束，所以ub也用[]代替，但由于其在最后，可以不写。 输出结果如下Optimization terminated. x = % 最优解 0.0000 15.0000 3.0000 fval = %最优值 -78.0000 exitflag = %函数收敛于解 1 output = iterations: 6 algorithm: large-scale: interior point cgiterations: 0 message: Optimization terminated. lambda = ineqlin: [3x1 double] eqlin: [0x1 double] upper: [3x1 double] lower: [3x1 double] 例2 一家家具公司生产桌子和椅子，用于生产的全部劳动力共计450个工时，原料是400个单位的木材每张桌子使用15个工时的劳动力，20个单位的木材售价为80元每张椅子使用10个工时用5个单位售价45元问为达到最大效益应如何安排生产个，椅子个，建立如下模型： f = [-80,-45]; A = [20, 5; 15, 10]; b = [400, 450]; lb = zeros(2,1); [x, fval, exitflag] = linprog(f,A,b,[],[],lb) 结果如下： Optimization terminated. x = 14.0000 24.0000 fval = -2.2000e+003 exitflag = 1 注意：由于linprog是求目标函数的最小值，如求目标函数的最大值，可先求出的最小值fval，则-fval就是的最大值。本例只输出最优解、最优值和退出标志，可见生产14个桌子，24个椅子，可获得最大利润2200元。 9.2 0-1规划 0-1规划是一种特殊形式的整数规划，它的决策变量仅取值0或1.一般用0表示放弃，1表示选取，故0-1规划可以用来处理选址问题、指派问题、装箱问题、项目评价、资金分配、生产计划安排等问题。 在M规划问题函数prog，其针对下述0-1规划： min 其中为向量，为矩阵。 [x, fval, exitflag, output] = bintprog(f, A, b, Aeq, beq) 输入部分：其中各符号对应0-1规划问题标准形式中的向量和矩阵，如果约束条件中有缺少，则其相应位置用空矩阵[]代替。 输出部分：其中x为最优解，用列向量表示；fval为最优值；exitflag为退出标志，若exitflag表示函数exitflag=-1表示超过的迭代最大，exitflag=-2，表示若output包含使用算法、迭代次数、搜索过的节点数、算法执行时间、算法终止原因。 例3 求解下述0-1规划问题。 利用bintprog函数求解如下： f=-[1;2;2;-6;-4]; A=[3,2,-1,4,2; 2,4, -2,-1,-2]; b=[5;5] ; [x,fval,exit,