matlab实验8 约束优化.doc

热动71 马千里 970669 实验八 约束优化 实验目的 掌握用MATLAB 优化工具箱解线性规划和非线性规划的方法； 练习建立实际问题的线性规划和非线性规划模型。 实验内容 3.对初值为(-1，1)，增加以下条件求解非线性规划： ，； ，，； ，，； ，；用分析梯度计算。 给出最优解、最优值，比较收敛速度，能得到什么结论。 解：按MATLAB的constr函数的要求编制程序如下 a. ，； function [f,g]=fun1(x) f=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1); g(1)=x(1)*x(2)-x(1)-x(2)+1.5; g(2)=-x(1)*x(2)-10; opt(1)=1; x=constr(‘fun1’,[-1 1],opt) 得出最优解x=[-9.5474 1.0474] 最优值 f=0.0236 迭代次数n=29 ，，； v1=[0 0]; opt(1)=1; x=constr(‘fun1’,[-1 1],opt,v1) 得出最优解x=[0 1.500] 最优值 f=0.500 迭代次数n=10 ，，； function [f,g]=fun2(x) f=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1); g(1)=x(1)*x(2)-x(1)-x(2)+1.5; g(2)=-x(1)*x(2)-10; g(3)=x(1)+x(2); g(4)=-x(1)-x(2); opt(1)=1; x=constr(‘fun2’,[-1 1],opt) 得出最优解x=[-1.2247 1.2247] 最优值 f=8.500 迭代次数n=13 ，；用分析梯度计算。 function [df,dg]=grad(x) df(1)=exp(x(1))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1)+exp(x(1))*(8*x(1)+4*x(2)); df(2)=exp(x(1))*(4*x(2)+4*x(1)+2); dg=[x(2)-1,-x(2);x(1)-1,-x(1)]; opt(1)=1; x=constr(‘fun1’,[-1 1],opt,[],[],’grad’) 得出最优解x=[-9.5474 1.0474] 最优值 f=0.0236 迭代次数n=11 可以看到，给出分析梯度时收敛较快。 10.炼油厂将A, B, C, 三种原油加工成甲，乙丙三种汽油。一桶原油加工成一桶汽油的费用为4元，每天至多能加工汽油14000桶。原油的买入价、买入量、辛烷值、硫含量，及汽油的卖出价、需求量、辛烷值、硫含量由下表给出。 问如何安排生产计划，在满足需求的条件下使利润最大？ 一般来说，作广告可以增加销售，估计一天向一种汽油投入一元广告费，可使这种汽油日销量增加10桶，问如何安排生产和广告计划使利润最大？ 原油类别 买入价(元/桶) 买入量(桶/天) 辛烷值(%) 硫含量(%) A B C 45 35 25 ≤5000 ≤5000 ≤5000 12 6 8 0.5 2.0 3.0 汽油类别 卖出价(元/桶) 需求量(桶/天) 辛烷值(%) 硫含量(%) 甲 乙 丙 70 60 50 3000 2000 1000 ≥10 ≥8 ≥6 ≤1.0 ≤2.0 ≤1.0 解：首先考虑不投入广告费的情况。设9个变量如下 生产甲类汽油消耗x1a桶A类原油，x1b桶B类原油，x1c桶C类原油； 生产乙类汽油消耗x2a桶A类原油，x2b桶B类原油，x2c桶C类原油； 生产丙类汽油消耗x3a桶A类原油，x3b桶B类原油，x3c桶C类原油； 按已知条件列出约束不等式： x1a+x2a+x3a≤5000 x1b+x2b+x3b≤5000 x1c+x2c+x3c≤5000 x1a+x2a+x3a+x1b+x2b+x3b+x1c+x2c+x3c≤14000 x1a+x1b+x1c≥3000 x2a+x2b+x2c≥2000 x3a+x3b+x3c≥1000 2*x1a-4*x1b-2x1c≥0 4*x2a-2x2b≥0 6*x3a+2*x3c≥0 -0.5*x1a+1*x1b+2*x1c≤0 -1.5*x2a+1*x2c≤0 -0.5*x3a+1*x3b+2*x3c≤0 最大利润： max f=3000*70+2000*60+1000*50-(x1a+x2a+x3a)*49-(x1b+x2b+x3b)*39 -(x1c+x2c+x3c)*29 用MATLA