1、解n元非线性方程组的牛顿法的MATLAB主程序
(相关资料图)
2、function[ci,D,danfan,xddf,hanfan,Xk,Yk]=newtonzun(X,tol,ftol,gxmax)
3、Y=Z(X);
4、fori=1:gxmax
5、dY=JZ(X);D=det(dY);Xk=X-(dY\Y")";
6、hanfan=norm(Y);danfan=norm(Xk-X);
7、xddf=danfan/(norm(Xk)+eps);X=Xk;Y=Z(X);ci=i;
8、ifD~=0
9、ci=i;Xk=X-(dY\Y")";Yk=Y;[ci,D,danfan,xddf,hanfan,X,Y];
10、else
11、disp("请注意!迭代方程组的系数行列式的值等于零.")
12、end
13、if(hanfan 14、[ci,D,danfan,xddf,hanfan,X,Y]; 15、return; 16、end 17、end 18、ifi>gxmax 19、disp("请注意:迭代次数超过给定的最大值gxmax,请重新输入初始值.") 20、return; 21、end 22、既然非常专业 为何不尝试自己写一个? 本文到此讲解完毕了,希望对大家有帮助。 关键词:
质检
推荐