试验8约束优化.DOCVIP

实验9 约束优化 实验目的： 1）掌握用matlab优化工具箱解线性规划和非线性规划的方法 2）练习建立实际问题的线性规划和非线性规划模型 实验内容: 第9题：某厂向用户提供发动机，合同规定，第一，二，三季度分别交货40，60，80台，每季生产费用为，交货后有剩余，可下季交货，每台每季储存费c元，每季最大生产能力100台，第一季度没有存货，设a=50,b=0.2,c=4,问如何安排生产使总费用最小，并讨论a,b,c的影响，作出合理解释。 ? ? 10）炼油厂将A、B、C三种原料加工成甲乙丙三种汽油。一桶原油加工成汽油的费用为4元，每天至多能加工汽油14，000桶。原油的买入价、买入量、辛烷值、硫含量，及汽油的卖出价、需求量、辛烷值、硫含量由下表给出。问如何安排生产计划，在满足需求的条件下使利润最大？ 一般来说，作广告可以增加销售，估计一天向一种汽油投入一元广告费，可以使这种汽油日销量增加10桶。问如何安排生产计划和广告计划使利润最大？ 原油类别 买入价（元/桶） 买入量（桶/天） 辛烷值（％） 硫含量（％） A 45 ≤5000 12 0.5 B 35 ≤5000 6 2.0 C 25 ≤5000 8 3.0 ? 汽油类别 卖出价（元/桶） 需求量（桶/天） 辛烷值（％） 硫含量（％） 甲 70 3000 ≥10 ≤1.0 乙 60 2000 ≥ ≤ 丙 50 1000 ≥6 ≤1.0 ? ? 11）一基金管理人的工作是，每天将现有的美元、英镑、马克、日元四种货币按当天的汇率兑换，使在满足需要的条件下，按美元计算的价值最高。设某天的汇率、现有货币和当天需求如下： ? 美元 英镑 马克 日元 现有量（） 需求量（） 美元 1 .58928 1.743 138.3 8 6 英镑 1.697 1 2.9579 234.7 1 3 马克 .57372 .33808 1 79.346 8 1 日元 .007233 .004426 .0126 1 0 10 问该天基金管理人应如何操作（“按美元计算的价值”指兑入、兑出汇率的平均值，入1英镑相当于（1.697+(1/.589.28)）/2=1.696993美元）。 ? ? 问题解答: 9): 本题是一个有约束条件的二次规划问题。决策变量是工厂每季度生产的台数，分别记为，目标函数是总费用（包括生产费用和贮存费），记为y。约束条件是生产合同，生产能力的限制。问题可抽出为求解： s.t. 将这个二次规划写成标准形式可得： ，， ，，， 用matlab优化工具箱计算如下： a=50;b=0.2;c=4; H=diag(2*b*ones(1,3));C=[a+2*c,a+c,a]; A1=[-1,0,0;-1,-1,0];b1=[-40,-100]; A2=[1 1 1];b2=180; v1=[0 0 0];v2=[100 100 100]; [x,faval,exitflag,output,lambada]=quadprog(H,C,A1,b1,A2,b2,v1,v2,[]) y=x*H*x/2+C*x-140*c 计算结果： x = 50.0000 60.0000 70.0000 y =11280 进一步分析，讨论参数a，b，c对生产计划的影响： 1）、固定b，c不变，a变化（分别取a＝20、60），仍运行上述程序，结果为： a y 20 50 60 70 5880 60 50 60 70 13080 2）、固定a，c不变，b变化（分别取b＝0、0.15、0.5），结果为： b y 0 40 60 80 9000 0.15 46.7 60 73.3 10727 0.5 56 60 64 14544 3）、固定a，b不变，c变化（分别取c＝0，2，6），结果为： c y 0 60 60 60 11160 2 55 60 65 11230 6 45 60 75 11310 结果分析： 由于生产总量是恒定的，而 ，故a的变化不会影响生产计划；b是x的二次项的系数，它反映了生产费用。当b比较大时，生产费用占主导地位，应趋于相等；而当b较小时，贮存费占主导地位，此时应使每季度的贮存量较少。c反映了贮存费。当c较大时，贮存费占主导地位，此时应使贮存量尽量少；而当c较小时，生产费用占主导地位，应趋于相等。 ? ? ? ? 10）： 问题分析求解： A种原油加工的汽油桶数分别为 B种原油加工的汽油桶数分别为 C种原油加工的汽油桶数分别为 依题意知，目标函数为总利润，记为f，约束条件为买入量、需求量的限制，加工能力的限制以及辛烷值、硫含量的要求。可得： S.T. 每天至多加工汽油总量为14 000桶 ；