数学建模课件线性规划.pptVIP

牛顿法代码实现dFdx1=diff(F,x1);dFdx2=diff(F,x2);dGdx1=diff(G,x1);dGdx2=diff(G,x2);D2z=[dFdx1dFdx2;dGdx1dGdx2];%JacobianofDz(sameasHessianofD2z)x0=[5;5];%initialguessN=10;%numberofiterationsfori=1:NDz0=subs(Dz,[x1,x2],[x0(1),x0(2)]);D2z0=subs(D2z,[x1,x2],[x0(1),x0(2)]);xnew=x0-inv(D2z0)*Dz0;x0=xnew;endxmax=xnewzmax=subs(z,[x1,x2],[xmax(1),xmax(2)])第29页，共77页，星期日，2025年，2月5日xmaxfigure,ezcontourf(z,[0.1100.110])holdonplot3(xmax(1),xmax(2),zmax,mo,LineWidth,2,...MarkerEdgeColor,k,MarkerFaceColor,[.491.63],...MarkerSize,12);title(Countourplotandoptimalvalue);第30页，共77页，星期日，2025年，2月5日3.3线性规划例3.4一个家庭农场有625英亩的土地可用来种植农作物。这个家庭可考虑种植的农作物有玉米、小麦、燕麦。预计有1000英亩-英尺的灌溉用水，农场工人每周可以投入的工作时间为300小时。其他数据如下表。为获得最大收益，每种作物应各种植多少？第31页，共77页，星期日，2025年，2月5日农场问题的有关数据条件（每英亩）作物玉米小麦燕麦灌溉用水（英亩-英尺）3.01.01.5劳力（人-小时/周）0.80.20.3收益（美元）400200250第32页，共77页，星期日，2025年，2月5日变量：x1,x2,x3=种植玉米、小麦、燕麦的亩数w=需要的灌溉用水（英亩-英尺）l=需要的劳力（人-小时/周）t=种植作物的总英亩数y=总收益（美元）第33页，共77页，星期日，2025年，2月5日假设：w=3.0x1+1.0x2+1.5x3=1000l=0.8x1+0.2x2+0.3x3=300t=x1+x2+x3=625y=400x1+200x2+250x3x1,x2,x3=0目标：求y的最大值第34页，共77页，星期日，2025年，2月5日建模方法—线性规划线性规划简介见书p.59可以用lindo/lingo软件求解第35页，共77页，星期日，2025年，2月5日模型求解MAX400X1+200X2+250X3SUBJECTTO3X1+X2+1.5X3=10000.8X1+0.2X2+0.3X3=300X1+X2+X3=625END第36页，共77页，星期日，2025年，2月5日reducedcost值表示当该非基变量增加一个单位时（其他非基变量保持不变）目标函数减少的量(对max型问题)也可理解为：为了使该非基变量变成基变量，目标函数中对应系数应增加的量第37页，共77页，星期日，2025年，2月5日灵敏性分析增加1英亩-英尺灌溉水量对最优解的影响MAX400X1+200X2+250X3SUBJECTTO3X1+X2+1.5X3=10010.8X1+0.2X2+0.3X3=300X1+X2+X3=625END第38页，共77页，星期日，2025年，2月5日玉米收益的少量提高对最优解的影响农作物每英亩收益会随气候及市场变化MAX450X1+200X2+250X3SUBJECTTO3X1+X2+1.5X3=10000.8X1+0.2X2+0.3X3=300X1+X2+X3=625END第39页，共77页，星期日，2025年，2月5日燕麦收益的少量提高对最优解的影响MAX400X1+20