MATLAB在优化模型中的应用.docVIP

1、线性规划问题的MATLAB辅助计算及工程应用实例 假设线性规划M题的数学模型为： min f1 x s.t. A^xb Aeq * x = beq lb x ub t.述模型中，/，x，/?，/?叫，//?和M为向量，A，Aeq为矩阵。 用于线性规划的MATLAB函数主要是linprog. Linprog函数的调用格式为： 1*，卜刈=11叩1*（^（匕八山，人69山叫，11），111）），求解上述的线性规划问题，如果没有不等式约朿则 令A=[]，b=[];如果没有等式约束则令Aeq=[]，beq=[];最终得到解x和x处的目标闲数位fval. 例： min - 2%, 一 x2 + 3%3 - 5x4 s.t. + 2x2 + 4x3 -%4 6 2x, +3x2 -x3 +%4 12 Xj + x3 + x4 4 x,,x2,x3,%4 0 解MATLAB求解程序为： f=[-2,-l,3,-5]f; A=[l，2,4，-l;2,3，-l，l;l，0，l，l]; b=「6，12,4L lb=[O,O,O,Or; [x，fval]=linprog（f，A，b，[]，[]，lb） 结果为： x = 0.0000 2.6667 0.0000 4.0000 fval = -22.6667 1、 生产决策问题 某工厂生产甲、乙两种产品，己知生产lkg产品甲需要用原料A5kg，原料B6kg;生产 lkg产品乙耑要用原料A3kg，原料B7kg,原料C5kg.若lkg产品甲和乙的销售价格分别为6 万元和5万元，三种原料的限用量分别为100kg，160kg，180kg.试确定应生产这两种产品各多 少千克才能使总销售价格最高？ 解令生产产品甲的数量为 ,生产产品乙的数S:为x2,有题意可建立下而的模型： max 6x, + 5%2 s.t. 5%, + 3%2 100 6x, +7x2 160 5x, 180 x,, 100 该模型要求使目标函数最大化，则按照MATLAB的要求进行转换，将目标函数最小化，即 min -6%, 一 5x2 s.t. 5x, +3x2 100 6x, +7x2 160 5x2180 xpx2 100 MATLAB程序清单为： f=[-6，-5].; A=[5,3；6J;0,5]; b=l 100,160,180J*; lb=[0,0r; [x，fval]=linprog（f，A，b,[]，[]，lb） 结果输出为： x = 12.9412 11.7647 fval = -136.4706 说明：生产产品甲、乙的数量分别为12.94kg，11.76kg吋，创造最高总售价为136.47万元。 2、 工程项目投资问题： 某公司有一批资金欲投资到5个工程项H中，各工程项目的净收益（投入资金的百分比） 见下表： 工程项目的净收益表 工程项目 A B C D E 收益（％） 10 12 15 12 8 由于一些原因，公司决定用于项目A的投资不大于其他各项目投资之和，而用于项目B和 项目D的投资要大于项目C和项目E的投资。试确定投资分配方案，使该公司收益达到最 大。 解：设6，久，；^，；14，15分别表示用于项目八上，：，0和￡的投资百分数，由于各项目的投资百 分数之和等于100%，所以 + X，+ + x5 = 1 根据题设可以建立下面的数学模型，即： max 0. lx, + 0.12x2 + 0.15x3 + 0.12x4 + 0.08x5 S.t. — X-j — — 一 0 %2 + x4 - x3 - %5 0 x x{ + x2 + x3 + x4 + x5 = 1 x,,x2,x3,x4,x5 0 将该模型转化为标准的线性规划模型，即： mins.t.MATLABA求解程序为:-0.1%,-0.12x2 -0.15x3-0.12x min s.t. MATLABA求解程序为: -x2 - %3 - x4 - x5 0 -- x4 + x3 + x5 0 r, + ;v2 + x3 + ;v4 + x5 = 1 ;px2,x3,x4,x5 0 f=[-0.1,-0.12,-0.15,-0.12,-0.08]*; A=[ l，-l，-l，-l，-l;0，-l，l,-l，l]; b=[0,0]’; Aeq=[l，l，l，l，l]; beq=[l]; lb=f0,0,0,0,OP; [x,fval]=linprog(f,A,b,Aeq,beq,lb) 结果输山为： 0.0000 0.2500 0.5000 0.2500 0.0000 fval = -0.1350 说明：5个项目的投资便分数分别为0, 25%, 50%, 25%和0时可以使该公司获得最大的 收益，最大收益为13.5%.