Newtons.m函数
function[x_star,index,it]=Newtons(fun,x,ep,it_max)
%求解非线性方程组的牛顿法,其中,fun(x)为需要求根的函数
%第一个分量是函数值,第二个分量是导数值
%x为初始向量(列向量)
% ep为精度,当 | |x(k)-x(k-1)| |<ep时,终止计算,缺省值为1e-5
% it_max为最大迭代次数,缺省值为100
% x_star为当迭代成功时,输出方程的根
% 当迭代失败时,输出最后的迭代值
% index为指标向量,当index=1时,表明迭代成功
% 当index=0时,表明迭代失败(迭代次数>=in+max)
% it为迭代次数
if nargin<4 it_max=100;end
if nargin<3 ep=1e-5;end
index=0;k=1;
while k<it_max
x1=x;[f,J]=feval(fun,x);
if abs(det(J))<ep break;end
x=x-J\f;
if norm(x-x1)<ep
index=1;break;
end
k=k+1;
end
x_star=x;it=k;
funs.m函数(输入任意一个非线性方程组进行测试)
function[f,J]=funs(x)
f=[x(1)^2+x(2)^2-5;(x(1)+1)*x(2)-(3*x(1)+1)];
J=[2*x(1) 2*x(2);x(2)-3 x(1)+1];
结果: