Newton迭代法Matlab程序求函数极小值点.docVIP

Newton迭代法Matlab程序求函数极小值点 clear allclc%Newton迭代法求解极小值点%090311%=====================================%定义函数disp 函数 f(x) 为:syms x1 x2f=(x1-2)^4+(x1-2)^2*x2^2+(x2+1)^2disp 初始点的值：x0=[1;1]%=====================================%求函数的梯度和海色阵disp 函数f的梯度：g=jacobian(f,[x1;x2])disp 函数f的Hesse矩阵：G=jacobian([g(1);g(2)],[x1,x2])%=====================================%定义迭代的最大次数n=10;%=====================================%一些初始值的计算g0=subs(g,{x1,x2},{x0(1),x0(2)});G0=subs(G,{x1,x2},{x0(1),x0(2)});f0=subs(f,{x1,x2},{x0(1),x0(2)});%=====================================%迭代点集合 x和函数值F的初始化x=zeros(2,n);F=zeros(1,n);%运用Newton方程解出下一近似值x(:,1)=x0-inv(G0)*g0; %注：用点乘有误A=x(:,1);F(1)=subs(f,{x1,x2},{A(1),A(2)});%=====================================%定义误差初始值为10deta=10;i=1;%循环用求出的近似解迭代Newton方程%求出下一个近似解，并在规定的误差范围内while deta=1e-10i10??? A=x(:,i);??? gi=subs(g,{x1,x2},{A(1),A(2)});??? Gi=subs(G,{x1,x2},{A(1),A(2)});??? i=i+1; ??? x(:,i)=x(:,i-1)-inv(Gi)*gi;??? A=x(:,i);??? F(i)=subs(f,{x1,x2},{A(1),A(2)});??? deta=F(i)-F(i-1);endk=(1:n+1);F=[f0 F];x=[x0,x];disp =====================================================disp 迭代的各步结果如下：disp ????? k???????????? x(k)???????? F{x(k)}[k????????? x???????????? F]%=====================================