非线性规划.docVIP

用MATLAB优化工具箱解线性规划 MATLAB软件求解线性规划的命令如下： x=linprog(c,A,b) 用于求解模型： Min z=cX (1) s.t. AXb A是不等式约束的系数矩阵，b是相应的常数列向量，若没有不等式约束，则均用[ ]代替，即令A=[ ],b=[ ]. x=linprog(c,A,b,Aeq,beq) 用于求解模型： Min z=cX (2) Aeq是等式约束的系数矩阵，beq是相应的常数列向量，若没有不等式约束，则均用[ ]代替. x=linprog(c,A,b,Aeq,beq,vlb,vub) 用于求解模型： Min z=cX (3) 若没有等式约束： 则令Aeq=[ ]，beq=[ ]. 如果某个变量无下界，则用-inf表示；如果某个变量无上界，则用inf表示；若决策变量x无下界和上界，则vlb和vub均用[ ]代替. x=linprog(c,A,b,Aeq,beq,vlb,vub,) 也用于求解模型（3），其中表示线性规划的初始解，这种设计仅对中规模算法有效，通常可以缺省.. [x,fval]=linprog(...) 返回最优解x及x处的目标函数值fval. [x,fval,exitflag,output,lambda]=linprog(c,A,b,Aeq,beq,vlb,vub) 也用于将最优目标值返回变量fval中. 输出项exitflag的值描述程序运行的情况：当exitflag的值为1时，说明程序收敛于解x；当exitflag的值为0时，说明迭代次数达到了设定的最大次数；当exitflag的值小于0时，说明的情况比较多，可在MATLAB的help中通过搜索函数linprog查看. 输出项output中的iterations表示程序的迭代次数；algorithm表示程序所用的算法，cgiterations给出共轭梯度迭代次数（仅对大规模算法）. 输出项lambda是解x处的Lagrange乘子，其中lower是对应于vlb的，upper是对应于vub的，ineqlin是对应于不等式约束的，eqlin是对应于等式约束的。在经济问题中，Lagrange乘子的含义是影子价格. [x,fval,exitflag,output,lambda]=linprog(c,A,b,Aeq,beq,vlb,vub,,options) 输入中多了两项，是初始点，单纯型法不需要初始点. options是选择项，应用options为函数设置选项，其调用格式是： options=optimset(‘param1’,value1, ‘param2’,value2,……) 可以用options选择使用单纯型算法，具体格式： options=optimset(‘LargeScale’, ‘off’, ‘Simplex’, ‘on’) 注. linprog对应看MATLAB5.3中的lp. lp目前在6.3版中仍然使用，但在以后的新版本中将废除. 例1 任务分配问题：某车间有甲乙两台机床，可用于加工三种工件. 假定这两台车床的可用台时数分别为800和900，三种工件的数量分别为400、600和500，且已知用三种不同车床加工单位数量不同工件所需的台时数和加工费用如下表. 问怎样分配车床的加工任务，才能既满足加工工件的要求，又使加工费用最低？ 车床 类型 单位工件所需加工台时数 单位工件的加工费用 可用 台时数 工件1 工件2 工件3 工件1 工件2 工件3 甲 0.4 1.1 1.0 13 9 10 800 乙 0.5 1.2 1.3 11 12 8 900 解 设在甲车床上加工工件1、2、3的数量分别是、、，在乙车床上加工工件1、2、3的数量分别是、、. 可建立以下线性规划模型： 改写为 编写M文件xxgh3.m如下： f=[13 9 10 11 12 8]; A=[0.4 1.1 1 0 0 0 0 0 0 0.5 1.2 1.3]; b=[800;900]; Aeq=[1 0 0 1 0 0 0 1 0 0 1 0 0 0 1 0 0 1]; beq=[400 600 500]; vlb=zeros(6,1); vub=[]; [x,fval]=linprog(f,A,b,Aeq,