运筹学II练习题.docx
运筹学Il练习题1试判定下述非线性规划是否为凸规划：firt/(X)=x12+8x12-x2O-x_%2+2=0x1,x2OMin/(X)=2x12+x22+x32-x1X2x12+x2245+=10X1,x2,x30max/(X)=x1+x2s,t.Mzh(X)=xi2÷x22+84。)=内2_/02(X)=-xi-x22+2=0X1,x20/(X),g(X),g2(X)的海赛矩阵的行列式:IhI=x2x打仆)x12，8(x)x2xlO(X)xix2ae%(x)xx1%(x)%2>00gz(X)IIix：1.I(X)x2xlMX)x.xOO72=o»g2(X)O-2x22知/()为严格凸函数，g(x),为凸函数，g2(x)为凹函数，所以不是一个凸规划问题。Min/(X)=2xl2+x22+X32-xlX2g；(X)=X2+疗4<z>1(X)=-(x12+x22)+40g2(X)=5x12+x3=10x1,2,x30同上有/(x),&(x)&(x)的海赛矩阵的行列式4网=-10-120>0Igb-2-2010,是凹函数，g2=是凸函数，不是凸规划问题。(3)min(-f(X)=-x-x2g(X)=IT202(X)=x10g3(X)=x20HT(X)=WpX)二>0,Hy(X)=Hr(X)=t2”20说明/(X)是凸函数，gI(X)、g2(X)、g3(X)是凹函数。因此，本模型是一个凸规划。2试用斐波那契法求函数/(x)=x2-3x+2在区间。10上的极小点，要求缩短后的区间长度不大于原区间长度的8沆(1.5)产lb=12.5,"=6;0=0,b0=10;F5t=b()(。040)=3.846;F6r=÷-(Z?0-0)=6.154;F6/(11)=5.254;/()=21.409;/(rl)<fS”H=O;bl=6.154"2'=3.846;F4t2=b-(bi-a)=2.308;f(t2)=0.403;/(r2)</(z2,)=5.254;/a2=0/2=3.846/3'=2.308;t3=h2-S2-a2)=1.538;F4f(t3)=-0.248<f(t3,)=0.403;/a3=0;Z?3=2.308;4=1.538;t4=b3-S3-。3)=0.769;F3/(r4)=0.284>/(r4,)=-0.248;-,a4=0.769;M=2.30815=1.538;/5'=+(M-¢/4)=1.538F23用分数法求/二/一f+2在区间-1,3上的近似极小点，要求缩短后的区向长度不大于原区间长的8%。(0.538)4试用最速下降法求函数/(X)=-(-2)2-2石的极大点，先以X(°)=(o,o)/为初始点进行计算，求出极大点；再以X=(0,1)7为初始点进行两次迭代。最后比较从上述两个不同初始点出发的寻优过程。(2,0)解求"X)=-&-2)2-2石的极大点，即求g(X)=-2)2+2君的极小点。(1)取初始点X=(0,0)，取精度£=0.1Vg(X)=2(xl-2)94x29Vg(0)=(4,0)TRg(X叫=(J(4+()j=16>£/J26H(XHoj)Vg(x®)/"(0)Vg(0)(40一/20Y-4"(4(j(,_162322XJ,Fg(X叫啕-犷用V(x)=2(2-2),4×7=(0,0)r即X即为极小点。/.1为/(X)的极大点。(2)取初始点(°)=(0,l)7取精度£=0.1,同上方法进行两次迭代有：两次步长4斗4.两次迭代结果比较：对于目标函数的等值线为椭圆的问题来说，椭圆的圆心即为最小值，负梯度方向指向圆心，但初值点与圆心在同一水平直线上时，收敛很快，即尽量使搜索路径呈现较少的直角锯齿状。3IT+-205求/(X)=+万尺玉工22%的极小点，取Xo=24o(1,1)X2V,(X)=3xi-x2-2x2-x,tVf(X0)=-126=-Vf(X0)=12-6A)=-12-6一不可(Xo)1.-1263-1X1=X0+>=-24+Ai2V(X1)=,1217121717126-6,=17121726173817.17.-126二W(X)TYX)_为一守(XJ守(乙)一126P1=-Vf(Xi)+0P0=-61277+12-6=28990210289289由于可(妈)=0o,故羽二口1即极小点，计算经两步终止。6试用牛顿法求解(0,0)Maxf(X)=取初始点X=(4,0),解求Mar"X)=的极大点，即求g(X)=H+¥+2的极小点。1/2O01/2X=1/20所以极大点为20g(X)=2玉2x2'Vg(X°)=8(Hm=V(X,)=00x'=(o,o)/(x*)=gf(X)=-XAX的极小点，此处“x)=JXzX=(x12+2xix2+E)=xl+x2,x1+2x27现从(°)=(l,l)7开始可(X)=(2,3)，M)=T/(X(。)=(-2,_3),V/(X(0)yPA)=(P(O)AP。唱1334于是v()A.±-lT=R,-43342>k34343434J(3434JP,二Jf(XYT(X)_"4341.W-.0守甲W(X)(2,3)f21342M)=-守(X)+Z°)v(x)rP(PJAP(1-AY.li更丫34,34jl342,342Jrp1K11I342,342Jkl2)342,342J-3×104-2×65二34，二34-39×(-104)+26×6513344故104'X®=X+4尸=P_,_2+但姿=P)(3434J1365(0_34r_故得到极小值点X(2)=(0,0)78考虑下面的非线性规划：maxf(X)=ln(xl+l)+x22x1+x23s.一x1,x20验证它为凸规划，并用K-T条件求解。(0.3)解原问题可写为min/(X)=-ln(xl+l)-x2s.t.gl(X)=3-2x1-220g2(X)=x103(X)=x2O计算目标和约束函数的海赛阵7o0oHf(X)=(x1÷D220,%(X)=(X)=(X)=oo0ooj1.故此问题是凸规划。W(X)=热T(X)=-2-,v2(X)=o,V3(X)=oK-T条件表达式为1C-=-21+2< i÷1< =-1÷3< l(3-2x1-x2)=0r2=°/,2W32°若%=0,则无解，于是差工0,有3-2%一马二0令玉=o,%=o,则有l-21+2=0< 1%=03-x2=0解得内=0,X2=3,%=力=1.%=O，显然03是可行点，从而是极小点。9试写出下述非线性规则的Kuhn-Tucker条件并进行求解：Mr(x)=(x-3)21 x5清华版，7章例110求解二次规划minf(X)=(2xjz+2xj)-8x1-1Ox263x.2x,20x10,x220见天大版例3-1611试解二次规划Min/(X)=2x12-4xix2+-6x1-3x2xl+x234x1+x29x1,x20解将上述二次规划改写为Min/(X)=g(4%：一8x1x2+-6xl-3x23NX?No9-4x1-x20x1,x2O可知目标函数为严格凸函数，此外ci=-6,C2=-3,4=4,c22=2c2=c2i=_4,b=3,au=-1,a2=-1b2=9,a2=-4,¾=-1由于q和C2小于零，故引入的人工变量Zl和Z2前面取负号，这样得到线性规划问题如下ming(z)=z1+z2-%一4以+必一4王+4%一4=-6一%-+y2+4x1-4x2-z2=-3一王一%-%3+3=0-4xl-x2-x4+9=Ox1,x2,x3,x4,ypy2,y3,y4,z1,z20解此线性规划问题得,202120,x；=0,E=320z；=O,z；=0,.21=O/(XA唱2-4××÷4×2-6×-3×21_4412020<20>202012试用SUMT外点法求解(1,2)MaX/(X)=X1(x2-2)÷(xi-1)30(x1-1)3-(x2-2)0xi,x2O解原问题转化为min/(X)=-x1-(x2-2)-(x,-1)30<-(xi-1)3+(x2-2)0xiOx2O构造惩罚函数P(x,M)=-x1+Mmin+minmin(0,-(x1-1)3+(x2-2)j+min(0,x1)+min(0,x2)1=-l+2Mmixl+2Mmin仪,(x1-1)3-(x2-2)j(3(x1-I)2)+2Mmin(O,x1)P=2Mminx2-(0,(x2-2)+(xi-1)j+2f(min0,-(x1-1)3+(x2-2)jj+2Mmin(0,x2)解得最优解为X*=(1,2)713一工人管2台机器，每台机器发生故障前的运转时间为具有均值为1/2小时的负指数分布，修理时间也属负指数分布，均值为1/3小时。(1)画出转速图。(2)列出平衡方程式求出状态概率P。，P1,P2。(3) 求故障机器数的均值1.o(4) 一台机器每次停机时间均值W5o3P二2台/小时,=3台/小时4248Pl=-Po,P2=一×-Po=Po333948Po+Pl+P2=PoTPoPo=1p。=39C12Pl=298P2=22%（台）=09666029=28（分钟）26（
