安监局网站做模拟,邓州网站推广,自己怎么做搬家网站,苏州建网站提clear allclc%Newton迭代法求解极小值点%090311%%定义函数disp 函数 f(x) 为:syms x1 x2f(x1-2)^4(x1-2)^2*x2^2(x21)^2disp 初始点的值#xff1a;x0[1;1]%%求函数的梯度和海色阵disp 函数f的梯度#xff1a;gjacobian(f,[x1;x2])disp 函数f的Hesse矩阵#xff1a;Gjacobia…clear allclc%Newton迭代法求解极小值点%090311%%定义函数disp 函数 f(x) 为:syms x1 x2f(x1-2)^4(x1-2)^2*x2^2(x21)^2disp 初始点的值x0[1;1]%%求函数的梯度和海色阵disp 函数f的梯度gjacobian(f,[x1;x2])disp 函数f的Hesse矩阵Gjacobian([g(1);g(2)],[x1,x2])%%定义迭代的最大次数n10;%%一些初始值的计算g0subs(g,{x1,x2},{x0(1),x0(2)});G0subs(G,{x1,x2},{x0(1),x0(2)});f0subs(f,{x1,x2},{x0(1),x0(2)});%%迭代点集合 x和函数值F的初始化xzeros(2,n);Fzeros(1,n);%运用Newton方程解出下一近似值x(:,1)x0-inv(G0)*g0; %注用点乘有误Ax(:,1);F(1)subs(f,{x1,x2},{A(1),A(2)});%%定义误差初始值为10deta10;i1;%循环用求出的近似解迭代Newton方程%求出下一个近似解并在规定的误差范围内while deta1e-10i10Ax(:,i);gisubs(g,{x1,x2},{A(1),A(2)});Gisubs(G,{x1,x2},{A(1),A(2)});ii1;x(:,i)x(:,i-1)-inv(Gi)*gi;Ax(:,i);F(i)subs(f,{x1,x2},{A(1),A(2)});detaF(i)-F(i-1);endk(1:n1);F[f0 F];x[x0,x];disp disp 迭代的各步结果如下disp k x(k) F{x(k)}[k x F]%
