实验三非线性规划.docVIP

数学建模试验报告（三） 姓名 学号 班级 0903 问题：（非线性规划） 某厂向用户提供发动机，合同规定，第一、二、三季度末分别交货40台、60台、80台．每季度的生产费用为 （元），其中x是该季生产的台数．若交货后有剩余，可用于下季度交货，但需支付存储费，每台每季度c元．已知工厂每季度最大生产能力为100台，第一季度开始时无存货，设a=50、b=0.2、c=4，问工厂应如何安排生产计划，才能既满足合同又使总费用最低．讨论a、b、c变化对计划的影响，并作出合理的解释．. 问题的分析和假设： 若每季度的生产费用为 f(x) = ax + bx^2（元） 设三季度分别生产x , y , 180-x-y台。 且应满足40≤x≤100，100≤x+y≤180，0≤y≤100，x,y∈N+（正整数）=50、b=0.2、c=4 则第一季度生产费用T1=50 x + 0.2x^2 剩余产品存储到下一季度的费用K1=4（x-40） 同理T2=50y + 0.2y^2 K2=4(x+y-100) T3=50(180-x-y) + 0.2(180-x-y )^2 总费用F=T1+T2+T3+K1+K2=9000+0.2(x^2+ y^2)+0.2(180-x-y) ^2+4(2x+y-140) 令Fx=0Fy=0 即 0.4x-0.4（180-x-y）+8=0 0.4y-0.4（180-x-y）+4=0 解得x=50 y=60 易验证该点处 Fxx≥0 Fyy≥0 即为F的极小值点。在通过和边界值的比较知其是定义域上的最小值点。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 求解的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 faval = 11840 exitflag = 1 output = iterations: 1 algorithm: medium-scale: active-set firstorderopt: [] cgiterations: [] message: Optimization terminated. lambada = lower: [3x1 double] upper: [3x1 double] eqlin: -78 ineqlin: [2x1 double] y = 11280 计算结果与问题分析讨论： 问题分析： 费用总量最低生产方案是：三个季度分别生产50、60、70台a,b,c对生产方案的影响： a增大或减小对生产方案完全没有影响（无论a为多少，方案都是50、60、70）。 b逐渐增大，则三个季度的生产量趋近交付总量的平均值，即同趋于60台（第一季度生产量增加，第二季度不变，第三季度减少）。 c逐渐增大，三季度的生产量分别趋近于每季度的交付量，即分别趋于40、60、80（第一季度生产量减少，第二季度不变，第三季度增加）。