最优化 刘志斌 课后习题3 5参考答案.doc

练习题三 1、用0.618法求解问题 3??2tt?min1(t)? t?0??。要求最后区间精度 的单谷区间为，的近似最优解，已知0.5?],3[0)(t答：t=0.8115；最小值-0.0886.（调用golds.m函数）（见例题讲解5） 2、求无约束非线性规划问题 222=min )x,xf(x,?4x2xx?x?2131312的最优解 解一：由极值存在的必要条件求出稳定点： ?f?f?f??xx?2x?2?8?2x?1x?0x?0 得，，则由，，，0x?f? 132132x???xx321再用充分条件进行检验： 222222f?f?f??ff??f， ，，，，00??08??2?2? 222xx?x??x?xx??xx???x313221213200????2T008f??(1,0,0)*?x，最优值为即-1为正定矩阵得极小点为。 ????200??解二：目标函数改写成 222)xx,f(x,?4?x?11)x(x?= min 312231易知最优解为（1,0,0），最优值为-1。 3、用最速下降法求解无约束非线性规划问题。 22?2xx?x?)?x?x2x(minfX 211122T0T)0(x)0X,?X?(x,。其中 ，给定初始点21?f(x)???? )x(?x?21?4x????121?)(x??f 解一：目标函数的梯度)f(x???1?2x?2x)(?fx????21?? )x(???2． 1?1????(0)(1)(0)(1)(0)方向作一维寻出发，沿令搜索方向再从dX?f(?f(X?)Xd)???????11???????0?1??????(1)(0)???? ，最优步长为优，令步长变量为，则有??X??d??????1?10??????22(0)(1)2??????????? 故))??(?(f(x)?fX??2(d?)?(??)?2(??2)11?1?0??????'(1)(1)(1)(0)?????1?点之后，可得令 求出 0?2?2(?)X??X??Xd???????111101??????1?1????(1)(1)(2) 与上类似地，进行第二次迭代：令?f)?(Xd)????f(X????1?1?????? ，最优步长为令步长变量为，则有2?1?11???????(2)(1)?? ???dX????????111???????故22(1)(2)2???????????)?(1?1)??5d?)((?1)?(?1)?2(??1)?2(2?1)(??(xf()?fX1)?211??0.8??????11')2(2)(1)(?????0(2)?10??? 可得令 ??X?X??d?????? 22221.1155??????0.2??(2)(2) 此时所达到的精度 )X??f(0.2828)X??f(???0.2???1????25?X1,)?(f ，本题最优解?X??1.5?? 解二：利用matlab程序求解 首先建立目标函数及其梯度函数的M文件 function f=fun(x) f=x(1)-x(2)+2*x(1)*x(1)+2*x(1)*x(2)+x(2)*x(2); function g=gfun(x) g=[1+4*x(1)+2*x(2),-1+2*x(1) +2* x(2) ]; 调用grad.m文件 x0=[0,0]; [x,val,k]=grad('fun','gfun',x0) 结果x=[ -1.0000 ,1.5000] val= -1.2500 k=33 。次的到最优解x=[ -1.0000 ,1.5000]；最优值val= -1.2500即迭代33 题。、试用Newton法求解第34 阵解一：计算目标函数的梯度和Hesse)x?f(???? )(x?xx?21?4????121?x)??f( 目标函数的梯度)f(x??xx?2?1?2)x?f(????21?? )x?(??20.50.5?42????2?1 ，其逆矩阵为?G)??G?f(X????0.5122?????0.5?0.5??TTT??????(0)(1)(0)?11,1.5??1,?G??f(X1)?0,0X??X ???0.51??(1)。 计算0(fX?)??1????)??1,25f(X 本题最优解，?X??1.5??解二：除了第3题建立两个M文件外，还需建立Hesse矩阵的M文件 利用matlab程序求解 首先建立目标函数及其梯度函数的M文件 function f=fun(x) f=x(1)-x(2)+2*x(1)*x(1)+2*x(1)*x(2)+x(2)*x(2); func