非线性方程组迭代法的matlab实现详细介绍
我正在用matlab来求解非线性方程组,想用下面的方法来算,但可惜我还没有弄出来,希望各位matlab高手的大哥、大姐们帮帮我,在此不胜感激!!!方法一:牛顿迭代法方法二...
我正在用matlab来求解非线性方程组,想用下面的方法来算,但可惜我还没有弄出来,希望各位matlab高手的大哥、大姐们帮帮我,在此不胜感激!!!
方法一:牛顿迭代法
方法二:牛顿下山法
方法三:修正牛顿法
由于用数学符号无法写在这里,计算的迭代格式(网上有)和例子(只要是非线性方程组就可以)就不写了,就由你们了。 展开
方法一:牛顿迭代法
方法二:牛顿下山法
方法三:修正牛顿法
由于用数学符号无法写在这里,计算的迭代格式(网上有)和例子(只要是非线性方程组就可以)就不写了,就由你们了。 展开
3个回答
展开全部
我之前用过迭代法解超越方程,不知道你那个好不好编。
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
展开全部
牛顿迭代法:
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
另外两个在此基础上稍作改动就可以了。
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
另外两个在此基础上稍作改动就可以了。
本回答被提问者采纳
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询