标题: 请问牛顿迭代法如何求方程的重根?谢谢。
性别:男-离线 阿尔法孝直
(雀力日进)

闽国公
遂安军节度使
★★★★★★

Rank: 19Rank: 19Rank: 19Rank: 19
柱国(正二品) 轩辕春秋年度最佳(游戏人生区)
组别 节度使
级别 卫将军
好贴 2
功绩 1796
帖子 6034
编号 19070
注册 2004-10-16
家族 轩辕雀党


发表于 2009-8-13 18:26 资料 个人空间 短消息 只看该作者 QQ
请问牛顿迭代法如何求方程的重根?谢谢。

并求出该重根是几重根。

最近在写用牛顿法解方程的程序,但牛顿法貌似只能求单根。

[ 本帖最后由 阿尔法孝直 于 2009-8-13 18:28 编辑 ]


顶部
性别:未知-离线 lzy0274

Rank: 3Rank: 3Rank: 3
组别 士兵
级别 忠义校尉
功绩 2
帖子 244
编号 56773
注册 2006-1-4


发表于 2009-8-14 21:54 资料 短消息 只看该作者
需要对牛顿迭代法进行适当的改变。
迭代格式采用:X(k)=X(k-1)-Mf(X(k-1))/f'(x(k-1));
其中,M为重根的阶次。


顶部
性别:男-离线 阿尔法孝直
(雀力日进)

闽国公
遂安军节度使
★★★★★★

Rank: 19Rank: 19Rank: 19Rank: 19
柱国(正二品) 轩辕春秋年度最佳(游戏人生区)
组别 节度使
级别 卫将军
好贴 2
功绩 1796
帖子 6034
编号 19070
注册 2004-10-16
家族 轩辕雀党


发表于 2009-8-14 22:00 资料 个人空间 短消息 只看该作者 QQ
我是说先把所有根(无论单或重)都求出来,然后一个一个判断是几重根。
顶部
性别:未知-离线 lzy0274

Rank: 3Rank: 3Rank: 3
组别 士兵
级别 忠义校尉
功绩 2
帖子 244
编号 56773
注册 2006-1-4


发表于 2009-8-15 11:00 资料 短消息 只看该作者
感觉有难度,用牛顿法实现这个应该有困难。
顶部
性别:未知-离线 KYOKO
(★御姐控★)

唐国公
荆南节度使
★★

Rank: 22Rank: 22Rank: 22Rank: 22
柱国(正二品)
组别 节度使
级别 大将军
功绩 1456
帖子 65615
编号 32
注册 2003-8-19
来自 BWL


发表于 2009-8-15 12:48 资料 个人空间 短消息 只看该作者
假设知道一个方程所有的根,又其中一个是二重根,如何最简单的方法判断哪个根是重根?

(X-A)(X-B)。。。(X-α)=0同原方程做判断?
顶部
性别:男-离线 阿尔法孝直
(雀力日进)

闽国公
遂安军节度使
★★★★★★

Rank: 19Rank: 19Rank: 19Rank: 19
柱国(正二品) 轩辕春秋年度最佳(游戏人生区)
组别 节度使
级别 卫将军
好贴 2
功绩 1796
帖子 6034
编号 19070
注册 2004-10-16
家族 轩辕雀党


发表于 2009-8-17 01:00 资料 个人空间 短消息 只看该作者 QQ
可以对逐个根不断求导,知道导数不为零为止。例如对某个方程F(X)=0的一个根X0求导:假设F'(X0)=0,F''(X0)=0,……,F(n-1阶导)(X0)=0,而F(n阶导)(X0)≠0,那么X0就是方程F(X)=0的n重根。
顶部
性别:女-离线 天宫公主
(司徒家的颖颖)

虞国公主

Rank: 12Rank: 12Rank: 12
组别 限制发言用户
级别 大将军
好贴 6
功绩 517
帖子 11552
编号 1037
注册 2004-10-25
来自 天津
家族 司徒实业


发表于 2009-8-17 09:03 资料 主页 短消息 只看该作者 QQ
假设 F(x) = 0 的解是一个有限集,任何一个根的次数也都是有限的*,且 F 足够光滑(具体条件还没想好),以下步骤可以得出 F(x) = 0 的所有解:

记,N(F, x_0) 为初始值为 x_0 时,用牛顿迭代法对方程 F(x) = 0 找到的那个根。

令 f_0 = F,随便挑一个 x_0,N(f_0, x_0) 是 F(x) = 0 的一个根。

i = 1;
while f_i (x) = 0 仍然有根 {
   
    令 f_i = f_{i-1} (x)/(x - N(f_{i-1}, x_0));
    求 N(f_i, x_0);/* 它仍然是 F(x) = 0 的一个根。
          如果没有重复根,N(f_i, x_0)
          必然不同于以前找到的根。如
          果有重复根,f_{i+1} 必然能让
          f_i 在 N(f_i, x_0) 这个重复根
          上降次。*/
    i ++;

}

* 例如,如果 f(x) = e^(1/x),g(x) = f(x) 的 analytic continuation,以上程序对解 g(x) = 0 就会死机。

[ 本帖最后由 天宫公主 于 2009-8-17 09:19 编辑 ]

本帖最近评分记录
青木风亮 2009-8-17 13:04 +300 好帖奖励
顶部
性别:男-离线 阿尔法孝直
(雀力日进)

闽国公
遂安军节度使
★★★★★★

Rank: 19Rank: 19Rank: 19Rank: 19
柱国(正二品) 轩辕春秋年度最佳(游戏人生区)
组别 节度使
级别 卫将军
好贴 2
功绩 1796
帖子 6034
编号 19070
注册 2004-10-16
家族 轩辕雀党


发表于 2009-8-17 14:57 资料 个人空间 短消息 只看该作者 QQ
嗯,明白了,谢谢天宫。
顶部
性别:未知-离线 KYOKO
(★御姐控★)

唐国公
荆南节度使
★★

Rank: 22Rank: 22Rank: 22Rank: 22
柱国(正二品)
组别 节度使
级别 大将军
功绩 1456
帖子 65615
编号 32
注册 2003-8-19
来自 BWL


发表于 2009-8-17 23:03 资料 个人空间 短消息 只看该作者
公主的那个只用了牛顿迭代法?

我看光用大学里学的高数都不行啊
顶部
性别:男-离线 阿尔法孝直
(雀力日进)

闽国公
遂安军节度使
★★★★★★

Rank: 19Rank: 19Rank: 19Rank: 19
柱国(正二品) 轩辕春秋年度最佳(游戏人生区)
组别 节度使
级别 卫将军
好贴 2
功绩 1796
帖子 6034
编号 19070
注册 2004-10-16
家族 轩辕雀党


发表于 2009-8-18 01:12 资料 个人空间 短消息 只看该作者 QQ
反正我是看懂了,每求出一根,就把当前函数除以(x-当前根)再解,求出一个除一次,就像因式分解那样。
顶部

正在浏览此帖的会员 - 4 人在线 - 0 位会员(0 隐身), 4 位游客




当前时区 GMT+8, 现在时间是 2024-11-28 23:40
京ICP备2023018092号 轩辕春秋 2003-2023 www.xycq.org.cn

Powered by Discuz! 5.0.0 2001-2006 Comsenz Inc.
Processed in 0.009577 second(s), 9 queries , Gzip enabled

清除 Cookies - 联系我们 - 轩辕春秋 - Archiver - WAP