声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3487|回复: 4

[综合讨论] 为何用牛顿迭代法对非线性方程组求解时总是出错

[复制链接]
发表于 2008-5-8 15:42 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我在做东西时,需要用到牛顿迭代法对非线性方程组求解,可是我的程序总是出错的,请高手帮我修改一下啊,下面附上我写的程序啊,谢谢啊
  function [r,n]=mulNewton3(F,x0,eps)
if nargin==1
    eps=1.0e-4;
end
J=jacobian(sym(F));
fa=subs(sym(F),x0);
fb=subs(sym(J),x0);
r=x0-fa/fb;
n=1;
tol=1;
while tol>eps
    x0=r;
    r=x0-fa/fb;
    tol=norm(r-x0);
    n=n+1;
    if(n>100000)
        disp('******');
        return;
    end
end

[ 本帖最后由 eight 于 2008-5-8 22:16 编辑 ]
回复
分享到:

使用道具 举报

发表于 2008-5-8 18:14 | 显示全部楼层

回复 楼主 的帖子

另外请附上你的错误提示
发表于 2008-5-8 22:17 | 显示全部楼层
原帖由 wlinwlin 于 2008-5-8 15:42 发表
我在做东西时,需要用到牛顿迭代法对非线性方程组求解,可是我的程序总是出错的,请高手帮我修改一下啊,下面附上我写的程序啊,谢谢啊
  function [r,n]=mulNewton3(F,x0,eps)
if nargin==1
    eps=1.0e-4;
end
J= ...
新人发帖前请先阅读置顶的会员守则,警告一次。养成一个习惯:无论上哪个论坛,先看置顶帖
发表于 2012-11-15 22:26 | 显示全部楼层
怎么看啊
发表于 2012-11-28 15:49 | 显示全部楼层
最好附上主程序
让程序能够运行
这样大家方便找问题
否者光看代码费事费力
很难找出问题所在

点评

赞成: 5.0
赞成: 5
  发表于 2012-11-28 21:19
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 02:48 , Processed in 0.059308 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表