24小时热门版块排行榜

返回列表

mjluanwu

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 1
在线: 45分钟
虫号: 3901524
注册: 2015-06-01

[交流] MATLAB牛顿法解非线性方程

function root=NewtonRootReal(f,a,b,eps)
f1=real(subs(sym(f),findsym(sym(f)),a));
f2=real(subs(sym(f),findsym(sym(f)),b));
if(f1==0)
root=f1;
end
if(f2==0)
root=f2;
end
if(f1*f2>0)
disp('>0')
return;
else
tol=1;
fun=diff(sym(f));
fa=real(subs(sym(f),findsym(sym(f)),a));
fb=real(subs(sym(f),findsym(sym(f)),b));
dfa=real(subs(sym(fun),findsym(sym(fun)),a));
dfb=real(subs(sym(fun),findsym(sym(fun)),b));
if(dfa>dfb)
      root=a-fa/dfa;
else
      root=b-fb/dfb;
end
while(tol>eps)
      r1=root;
      fx=real(subs(sym(f),findsym(sym(f)),r1));
      dfx=real(subs(sym(fun),findsym(sym(fun)),r1));
      root=r1-fx/dfx;
      tol=abs(root-r1);
end
end

函数f既有实数部分，又有虚数部分，我把牛顿法修改了只求函数f实数部分为零的解。基于原始的NewtonRoot法，增加了real以求得实数部分。其他未曾修改。求大神看一下可以这样修改么？应该没有问题吧

回复此楼

» 收录本帖的淘帖专辑推荐

matlab编程绘图

matlab学习精华

» 猜你喜欢

1楼 2015-09-06 15:02:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 mjluanwu 的主题更新

返回列表