%阻尼Newton算法.docVIP

%阻尼Newton算法 function y=myznnd(x0,e) syms x1 x2; f1=x1^2+2*x2^2-2*x1*x2-4*x1; x0=[1,1]; e=0.001; gf1=jacobian(f1,[x1,x2]); gf2=jacobian(gf1,[x1,x2]); n=0; while 1 d1=subs(subs(gf1,x0(1)),x0(2)); d2=subs(subs(gf2,x0(1)),x0(2)); syms k; t=x0-d2\d1*k; d=diff(subs(f1,{x1,x2},t),k); k0=solve(d,k); x=x0; a=x0-d2\d1*k0; x0=a; n=n+1; if norm(d1)=e; break; end end x0,f1=subs(f1,{xi,yi},x0) function f=grad1(x0,e) x0=[2,2]; e=0.001; n=2; k=0; for k=0:1:n x=x0; syms xi yi a ; f=xi^2-xi*yi+yi^2+2*xi-4*yi; fx=diff(f,xi); fy=diff(f,yi); fx=subs(fx,{xi,yi},x0); fy=subs(fy,{xi,yi},x0); fi=vpa([fx,fy]); s=-fi; x=x+a*s; f=subs(f,{xi,yi},x); f1=diff(f); f1=solve(f1); if f1~=0 ai=vpa(f1); else break end x=subs(x,a,ai); f=xi^2-xi*yi+yi^2+2*xi-4*yi; fxi=diff(f,xi); fyi=diff(f,yi); fxi=subs(fxi,{xi,yi},x); fyi=subs(fyi,{xi,yi},x); fii=[fxi,fyi]; while (norm(fii)e) x=x+a*s; f=subs(f,{xi,yi},x); f1=diff(f); f1=solve(f1); if f1~=0 ai=vpa(f1); else break x,f=subs(f,{xi,yi},x) end x=subs(x,a,ai); f=xi^2-xi*yi+yi^2+2*xi-4*yi; fxi=diff(f,xi); fyi=diff(f,yi); x=subs(x,a,ai); f=xi^2-xi*yi+yi^2+2*xi-4*yi; fxi=diff(f,xi); fyi=diff(f,yi); fxi=subs(fxi,{xi,yi},x); fyi=subs(fyi,{xi,yi},x); fii=[fxi,fyi]; if kn-1 d=(fxi^2+fyi^2)/(fx^2+fy^2); s1=-fii+d*s; s=s1; else x0=x; k=0; break end end end x,f=subs(f,{xi,yi},x)