优化设计作业.docVIP

优化设计及matlab程序设计 一、优化设计题目 15、目标函数： 要求：分别采用1）阻尼牛顿法；2）共轭梯度法；3）变尺度法 等3种方法进行求解 二、牛顿阻尼法 1、阻尼牛顿法流程图 是 否 2.阻尼牛顿法命令流 X0=input(请输入迭代初向量X0=\n); n=1; eps=input(请输入精度eps=\n); syms x1 x2 x3; f=x1^2+2*x2^2+2*x3^2+0.5*x1*x2-2*x1*x3+2*x2*x3+3*x1+2*x2+2; df=[diff(f,x1);diff(f,x2);diff(f,x3)];%df表示梯度表达式 G=[diff(df(1),x1) diff(df(1),x2) diff(df(1),x3); diff(df(2),x1) diff(df(2),x2) diff(df(2),x3); diff(df(3),x1) diff(df(3),x2) diff(df(3),x3)]%G表示海森矩阵表达式 while(1) df0=subs(df,[x1;x2;x3],X0); H=subs(G,[x1;x2;x3],X0);%H表示海森矩阵 fprintf(第%d次迭代:\n,n) fprintf(输出海森矩阵\n) disp(H) df02=inv(H); %df02表示海森矩阵的逆矩阵 fprintf(输出海森矩阵的逆矩阵\n); disp(df02) d0=-df02*df0; syms a; X1=X0+d0*a;%a表示最优步长 f1=subs(f,[x1;x2;x3],X1); da=eval(diff(f1,a)); digits(5) a=vpa(solve(da,a)); fprintf(输出最优步长因子a%d:\n,n) disp(a) X1=X0+d0*a; eps1=norm(X1-X0); if(eps1=eps) X0=X1; n=n+1; else break; end end fprintf(\n 运用阻尼牛顿法经过%d迭代次后得到极小点X及极小值f为:\n,n); X=X1 f=subs(f,[x1;x2;x3],X) 3.阻尼牛顿法命令流运行结果 请输入迭代初向量X0= [1;1;1] 请输入精度eps=0.02 G = [ 2, 1/2, -2] [ 1/2, 4, 2] [ -2, 2, 4] 第1次迭代: 输出海森矩阵 2.0000 0.5000 -2.0000 0.5000 4.0000 2.0000 -2.0000 2.0000 4.0000 输出海森矩阵的逆矩阵 4.0000 -2.0000 3.0000 -2.0000 1.3333 -1.6667 3.0000 -1.6667 2.5833 输出最优步长因子a1:1.0 第2次迭代: 输出海森矩阵 2.0000 0.5000 -2.0000 0.5000 4.0000 2.0000 -2.0000 2.0000 4.0000 输出海森矩阵的逆矩阵 4.0000 -2.0000 3.0000 -2.0000 1.3333 -1.6667 3.0000 -1.6667 2.5833 输出最优步长因子a2:0.0 运用阻尼牛顿法经过2迭代次后得到极小点X及极小值f为: X = -8.0 3.3333 -5.6667 f =-6.6667 三、共轭梯度法 1.共轭梯度法流程图 2.共轭梯度法命令流 X0=input(请输入初始迭代点X0(以列向量形式表示）：X0=); e=input(请输入允许的误差：e=); syms x1 x2 x3 t=1; n=3;