最优化方法及其matlab程序设计 马昌凤 课后答案(21)

最优化方法及其matlab程序设计 马昌凤版 课后答案 杭电课件

%gval=norm(gfun(x));%%目标函数(I)t=0.5

functiony=Fk(x)

y(1)=x(1)2+x(2)2+x(3)2 1;y(2)=x(1)+x(2)+x(3) 1;y(3)=x(1)2+x(2)2+(x(3) 2)2 1;y(4)=x(1)+x(2) x(3)+1;

y(5)=x(1)3+3 x(2)2+(5 x(3) x(1)+1)2 36 0.5;y=y(:);

%%%%Jacobi阵%%%%%%functionJF=JFk(x)JF=[2*x(1),2*x(2),2*x(3);1,1,1;

2*x(1),2*x(2),2*(x(3)-2);1,1-1;

3 x(1)2 2 (5 x(3) x(1)+1),6 x(2),10 (5 x(3) x(1)+1)];>>x0=[1,1,1]′;[x,val,k]=lmm(′Fk′,′JFk′,x0)x=

0.339361063668441-0.2001835788046710.714384339944574val=

0.486062168183995

21