最优化计算方法.pptVIP

牛顿法求较精确的近似值牛顿法见书p.56x=[x1;x2];F=diff(z,x1);G=diff(z,x2);Dz=[F;G];%thegradientvectorofz即求方程组Dz=0的解。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)])xmaxfigure,ezcontourf(z,[0.1100.110])holdonplot3(xmax(1),xmax(2),zmax,mo,LineWidth,2,...MarkerEdgeColor,k,MarkerFaceColor,[.491.63],...MarkerSize,12);title(Countourplotandoptimalvalue);例3.4一个家庭农场有625英亩的土地可用来种植农作物。这个家庭可考虑种植的农作物有玉米、小麦、燕麦。预计有1000英亩-英尺的灌溉用水，农场工人每周可以投入的工作时间为300小时。其他数据如下表。为获得最大收益，每种作物应各种植多少？条件（每英亩）作物玉米小麦燕麦灌溉用水（英亩-英尺）3.01.01.5劳力（人-小时/周）0.80.20.3收益（美元）400200250变量：x1,x2,x3=种植玉米、小麦、燕麦的亩数w=需要的灌溉用水（英亩-英尺）l=需要的劳力（人-小时/周）t=种植作物的总英亩数y=总收益（美元）假设：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的最大值建模方法—线性规划01线性规划简介见书p.5902可以用lindo/lingo软件求解MAX400X1+200X2+250X3SUBJECTTOX1+X2+1.5X3=10000.8X1+0.2X2+0.3X3=300X1+X2+X3=625END01reducedcost值表示当该非基变量增加一个单位时（其他非基变量保持不变）目标函数减少的量(对max型问题)02也可理解为：为了使该非基变量变成基变量，目标函数中对应系数应增加的量增加1英亩-英尺灌溉水量对最优解的影响MAX400X1+200X2+250X3SUBJECTTOX1+X2+1.5X3=10010.8X1+0.2X2+0.3X3=300X1+X2+X3=625END农作物每英亩收益会随气候及市场变化MAX450X1+200X2+250X3SUBJECTTOX1+X2+1.5X3=10000.8X1+0.2X2+0.3X3=300X1+X2+X3=625ENDMAX400X1+200X2+260X3SUBJECTTOX1+X2+1.5X3=10000.8X1+0.2X2+0.3X3=300X1+X2+X3=625END这种玉米新品种需要较少的灌溉用水2.5英亩-英尺（而不是3.0）。MAX400X1+200X2+250X3SUBJECTTOX1+X2+1.5X3=1000X1+0.2X2+0.3X3=300X1+X2+X3=625END一英亩大麦需要1.5英亩-英尺的水和0.25人-小时的劳力，预期可获得200美元的收益。用一个新的决策