牛頓法求解非線性方程組matlab源程序(共3頁)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上牛頓法求解非線性方程組matlab源程序Newton-Raphson 求解非線性方程組matlab源程序matlab程序如下：functionhomP,iter,err=newton(f,JF,7.8e-001;4.9e-001; 3.7e-001,0.01,0.001,1000);(P);(iter);(err);function Y=f(x,y,z)Y=x2+y2+z2-1;2*x2+y2-4*z;3*x2-4*y+z2;function y=JF(x,y,z)f1=x2+y2+z2-1;f2=2*x2+y2-4*z;f3=3*x2-4*y+z2;df1x=(sy

2、m(f1),x);df1y=(sym(f1),y);df1z=(sym(f1),z);df2x=(sym(f2),x);df2y=(sym(f2),y);df2z=(sym(f2),z);df3x=(sym(f3),x);df3y=(sym(f3),y);df3z=(sym(f3),z);j=df1x,df1y,df1z;df2x,df2y,df2z;df3x,df3y,df3z;y=(j);function P,iter,err=newton(F,JF,P,tolp,tolfp,)%輸入P為初始猜測值，輸出P則為近似解%JF為相應(yīng)的Jacobian矩陣%tolp為P的允許誤差%tolfp為f

3、(P)的允許誤差%max：循環(huán)次數(shù)Y=(F,P(1),P(2),P(3);for k=1:=(JF,P(1),P(2),P(3);Q=P-()*Y;Z=(F,Q(1),Q(2),Q(3);err=(Q-P);P=Q;Y=Z;iter=k;if (errtolp)|(Y)tolfp|abs(Y)0.0001)breakendendfunctionhomework4P,iter,err=newton(f,JF,7.8e-001;4.9e-001; 3.7e-001,0.01,0.001,1000);(P);(iter);(err);function Y=f(x,y,z)Y=x2+y2+z2-1;2

4、*x2+y2-4*z;3*x2-4*y+z2;function y=JF(x,y,z)f1=x2+y2+z2-1;f2=2*x2+y2-4*z;f3=3*x2-4*y+z2;df1x=(sym(f1),x);df1y=(sym(f1),y);df1z=(sym(f1),z);df2x=(sym(f2),x);df2y=(sym(f2),y);df2z=(sym(f2),z);df3x=(sym(f3),x);df3y=(sym(f3),y);df3z=(sym(f3),z);j=df1x,df1y,df1z;df2x,df2y,df2z;df3x,df3y,df3z;y=(j);function P,iter,err=newton(F,JF,P,tolp,tolfp,)%輸入P為初始猜測值，輸出P則為近似解%JF為相應(yīng)的Jacobian矩陣%tolp為P的允許誤差%tolfp為f(P)的允許誤差%max：循環(huán)次數(shù)Y=(F,P(1),P(2),P(3);for k=1:=(JF,P(1),P(2),P(3);Q=P-()*Y;Z=(F,Q(1),Q(2),Q(3);err