matlab代码--非线性方程组求解.pdfVIP

非线性方程组求解

1.mulStablePoint用不动点迭代法求非线性方程组的一个根

function[r,n]=mulStablePoint(F,x0,eps)

%非线性方程组：f

%初始解：a

%解的精度：eps

%求得的一组解：r

%迭代步数：n

ifnargin==2

eps=1.0e-6;

end

x0=transpose(x0);

n=1;

tol=1;

whiletoleps

r=subs(F,findsym(F),x0);%迭代公式

tol=norm(r-x0);%注意矩阵的误差求法，norm为矩阵的欧几里

德范数

n=n+1;

x0=r;

if(n100000)%迭代步数控制

disp(迭代步数太多，可能不收敛！);

return;

end

end

2.mulNewton用牛顿法法求非线性方程组的一个根

function[r,n]=mulNewton(F,x0,eps)

i