最优化-刘志斌-课后习题3-5参考答案要点word版本.docVIP

最优化-刘志斌-课后习题3-5参考答案要点 精品文档 精品文档 收集于网络，如有侵权请联系管理员删除 收集于网络，如有侵权请联系管理员删除 精品文档 收集于网络，如有侵权请联系管理员删除 练习题三 1、用0.618法求解问题 的近似最优解，已知的单谷区间为，要求最后区间精度。 答：t=0.8115；最小值-0.0886.（调用golds.m函数）（见例题讲解5） 2、求无约束非线性规划问题 min = 的最优解 解一：由极值存在的必要条件求出稳定点： ，，，则由得，， 再用充分条件进行检验： ，，，，， 即为正定矩阵得极小点为，最优值为-1。 解二：目标函数改写成 min = 易知最优解为（1,0,0），最优值为-1。 3、用最速下降法求解无约束非线性规划问题。 其中，给定初始点。 解一：目标函数的梯度 令搜索方向再从出发，沿方向作一维寻优，令步长变量为，最优步长为，则有 故 令可得 求出点之后，与上类似地，进行第二次迭代： 令 令步长变量为，最优步长为，则有 故令可得 此时所达到的精度 本题最优解， 解二：利用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 即迭代33次的到最优解x=[ -1.0000 ,1.5000]；最优值val= -1.2500。 4、试用Newton法求解第3题。 解一：计算目标函数的梯度和Hesse阵 目标函数的梯度 ，其逆矩阵为 计算。 本题最优解， 解二：除了第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); function g=gfun(x) g=[1+4*x(1)+2*x(2),-1+2*x(1) +2* x(2) ]; function h=hess(x) g=[4 2;2 2 ]; 调用newton.m文件 x0=[0,0]; [x,val,k]=newton(fun,gfun,hess,x0) 结果 x=[ -1.0000 ,1.5000] val= -1.2500 k=1 5、用Fletcher—Reeves法求解问题 其中，要求选取初始点。 解一： ， 第一次迭代：令， 即， 第二次迭代： ，， 第三次迭代： ……（建议同学们自己做下去，注意判别） 解二：利用matlab程序求解 首先建立目标函数及其梯度函数的M文件 function f=fun(x) f=x(1)^2+25* x(2)*x(2); function g=gfun(x) g=[2*x(1), 50* x(2) ]; 调用frcg.m文件 x0=[2,2]’;epsilon=1e-6; [x,val,k]=frcg(fun,gfun,x0, epsilon) 结果 x = 1.0e-006 *[ 0.2651, 0.0088] val =7.2182e-014 k = 61 6、试用外点法（二次罚函数方法）求解非线性规划问题  其中 解：设计罚函数 采用Matlab编程计算，结果x=[1 0]；最优结果为1。（调用waidianfa.m） 7、用内点法（内点障碍罚函数法）求解非线性规划问题： 解：容易看出此问题最优解为x=[1 0]；最优值为8. 给出罚函数为 令； 从而当时， （建议同学自己编程序计算） 8、用乘子法求解下列问题 解：建立乘子法的增广目标函数： 令： 解上述关于x的二元一次方程组得到稳定点 当乘子取2时，或发参数趋于无穷时，得到即最优解。 （建议同学自己编程序计算） 练习题四 1、石油输送管道铺设最优方案的选择问题：考察网络图4-6，设A为出发地，F为目的地，B，C，D，E分别为四个必须建立油泵加压站的地区。图中的线段表示管道可铺设的位置，线段旁的数字表示铺设这些管线所需的费用。问如何铺设管道才能使总费用最小？ 图4- SEQ 图4- \* ARABIC 6 解：? 第五阶段：E1—F 4；E2—F 3； 第四阶段：D1—E1?—F? 7；D2—E2—F ??5；D3—E1—F?? 5； 第