Question

非线性方程组迭代法的matlab实现详细介绍

我正在用matlab来求解非线性方程组，想用下面的方法来算，但可惜我还没有弄出来，希望各位matlab高手的大哥、大姐们帮帮我，在此不胜感激！！！ 方法一：牛顿迭代法 方法二：牛顿下山法 方法三：修正牛顿法 由于用数学符号无法写在这里，计算的迭代格式（网上有）和例子（只要是非线性方程组就可以）就不写了，就由你们了。

Answer 1

牛顿迭代法：
function[x0,n]=newton(fx,dfx,x0,tol,N)
%
牛顿迭代法
%
第一个参数
fx
是关于变量x的所要求的函数表达式.
%
第二个参数
dfx
是fx的一阶导数.
%
x0
是迭代初始值.
%
tol
是迭代误差限.
%
N
最大迭代次数.
x=x0;f0=eval(fx);df0=eval(dfx);
n=0;
disp('[
n
xn
xn+1
delta
]');
while
n<=N
x1=x0-f0/df0;
x=x1;f1=eval(fx);df1=eval(dfx);
delta=abs(x0-x1);
%
X=[n,x0,x1,delta];
disp(X);
%用于显示中间结果
if
delta<tol
fprintf('迭代计算成功')
return
else
n=n+1;
x0=x1;f0=f1;df0=df1;
end
end
if
n==N+1
fprintf('迭代计算失败
')
end
另外两个在此基础上稍作改动就可以了。