声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2992|回复: 18

[综合讨论] 如何确定牛顿迭代的初始值

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

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

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

x
现在下预用牛顿迭代法求解方程的根,而根的范围是在复数域,但是没有初值,请问有确定初始值的好办法吗?谢谢!!!
回复
分享到:

使用道具 举报

发表于 2008-1-2 16:26 | 显示全部楼层

回复 #1 kuang333 的帖子

没用过,不过在下想是不是可以慢慢试呢?
最后看试出来,哪个合适就哪个
发表于 2008-1-2 16:26 | 显示全部楼层
问题很久前就有讨论

http://forum.vibunion.com/forum/viewthread.php?tid=23149&extra=page%3D1

再次发帖注意下版规

[ 本帖最后由 花如月 于 2008-1-2 16:27 编辑 ]

评分

1

查看全部评分

 楼主| 发表于 2008-1-2 16:44 | 显示全部楼层
我已经看了之前相关的帖子了,可是现在有个关键的问题是我想求复数解,这个范围不好定啊,相应的初值也就不好定,类似于这方面的问题各位之前有过研究吗?请指教!
发表于 2008-1-2 16:54 | 显示全部楼层

回复 #4 kuang333 的帖子

复数解你的迭代公式怎么写?
发表于 2008-1-2 18:35 | 显示全部楼层
方程不是很复杂的话,我可以帮你求

[ 本帖最后由 eight 于 2008-1-2 18:44 编辑 ]
 楼主| 发表于 2008-1-2 19:02 | 显示全部楼层
迭代公式还是常规的牛顿迭代公式,只不过将解的范围扩大到了复数范围。
式子还算比较复杂,如果有思路的话也可直接与我联系,我的QQ是271659505,谢谢!
发表于 2008-1-2 19:05 | 显示全部楼层

回复 #7 kuang333 的帖子

建议你把方程用图片的形式,附件上来
 楼主| 发表于 2008-1-2 19:34 | 显示全部楼层
式子如下所示: clip_image002.jpg

u为变量,就是要找u的初值。
发表于 2008-1-3 09:29 | 显示全部楼层
暂时列出这59个根吧
当然了,不一定全
跟:
-4.41304421045121 - 4.79581948323569i
-4.05137399488236 - 4.44726391081319i
-3.65609750055567 - 4.06697754955189i
  0.06214854384538 - 0.07479458777251i
  4.79581424766862 - 4.41303778916129i
-4.41304421045120 + 4.79581948323570i
-10.91413907715508 -11.15304584946110i
15.92327850897341 -15.73394769476522i
15.11178210723974 -14.91569548273951i
-11.15304543630551 -10.91413863986255i
  4.44725733757458 - 4.05136574802695i
  4.06696894228589 - 3.65608638184522i
  3.64389740262993 - 3.21585839823413i
-3.21587453472969 + 3.64390940396964i
-3.65609750055567 + 4.06697754955189i
-4.05137399488237 + 4.44726391081320i
-14.91569565449129 +15.11178227312519i
-15.73394784113157 +15.92327865074456i
-17.43567246798397 -17.61249181957545i
-17.34549764310308 -17.52293353833762i
-17.61249171483751 -17.43567236037770i
-17.52293343198426 -17.34549753381211i
-10.76987205948511 -11.01072219397109i
16.02181434424667 -15.83326833055867i
15.82412667476269 -15.63399866530976i
15.21559630039253 -15.02040425578323i
15.00724661649096 -14.81024978061342i
-11.01072176455814 -10.76987160445121i
-0.06214854384984 - 0.07479458712177i
-0.06214854384538 + 0.07479458777251i
  0.00811316383591 + 0.00094522677081i
15.11178210723976 +14.91569548273954i
-14.81024995605276 +15.00724678587076i
-15.02040442397468 +15.21559646290222i
-15.63399881449650 +15.82412681921823i
-15.83326847419232 +16.02181448341573i
14.91569565449119 +15.11178227312503i
21.62137054153009 +21.46796162595834i
21.46796168364436 +21.62137059811860i
28.07465423856983 +27.94728955024177i
27.94728957640216 +28.07465426440884i
-0.06214854383774 + 0.07479458777246i
  0.06214854400273 + 0.07479458786650i
  8.77646806345800 + 8.49963704166806i
  8.49963796629286 + 8.77646891742103i
15.00724661649089 +14.81024978061350i
14.81024995605285 +15.00724678587086i
34.59481665402060 +34.48546175693473i
34.48546177086147 +34.59481666782734i
  8.59439994432930 + 8.31409365327278i
  8.31409463785461 + 8.59440084845787i
28.01861216326493 +27.89106348173092i
27.89106350804990 +28.01861218925924i
34.50384804262111 +34.39428068040120i
34.54936231339855 +34.43990136884590i
34.43990138282860 +34.54936232726024i
21.39474741599969 +21.54851998768447i
34.39428069443948 +34.50384805653769i
21.54851993051972 +21.39474735772002i

带入回去误差:
     1.818989403545857e-012 +4.547473508864641e-013i
     2.273736754432321e-012 -9.094947017729282e-013i
     4.547473508864641e-013 -3.410605131648481e-013i
     8.073282778517550e-021 -1.799945012915388e-021i
    -1.364242052659392e-012 -1.818989403545857e-012i
                          0 +1.364242052659392e-012i
     5.820766091346741e-011 +1.164153218269348e-010i
    -9.313225746154785e-010 -4.656612873077393e-010i
    -6.984919309616089e-010 -6.984919309616089e-010i
     2.328306436538696e-010 -1.164153218269348e-010i
    -1.364242052659392e-012 -9.094947017729282e-013i
    -1.724400817693095e-007 -1.298690222029109e-007i
    -1.136868377216160e-013                        
     4.547473508864641e-013 +1.136868377216160e-013i
                          0 -3.410605131648481e-013i
    -4.547473508864641e-013 -1.818989403545857e-012i
    -4.656612873077393e-010 -9.313225746154785e-010i
                          0 -4.656612873077393e-010i
                          0 -1.862645149230957e-009i
                          0 +1.862645149230957e-009i
     4.656612873077393e-009 -4.190951585769653e-009i
     1.862645149230957e-009 -1.396983861923218e-009i
                          0 +1.746229827404022e-010i
    -6.984919309616089e-010 -9.313225746154785e-010i
     4.656612873077393e-010 +4.656612873077393e-010i
     1.164153218269348e-009 +1.396983861923218e-009i
    -4.656612873077393e-010 -6.984919309616089e-010i
    -5.820766091346741e-011 -5.820766091346741e-011i
     3.933891038505578e-015 +4.870370757634161e-016i
     6.617444900424221e-023 -1.058791184067875e-022i
     2.574639240131468e-008 -5.505186161540935e-008i
    -2.328306436538696e-010 +6.984919309616089e-010i
     9.313225746154785e-010 +9.313225746154785e-010i
    -4.656612873077393e-010                        
    -2.328306436538696e-010 -2.328306436538696e-010i
     2.328306436538696e-010 +2.328306436538696e-010i
     1.164153218269348e-009 -4.656612873077393e-010i
     3.725290298461914e-009 -3.725290298461914e-009i
     1.862645149230957e-009                        
     3.725290298461914e-009 -2.235174179077148e-008i
    -1.490116119384766e-008                        
     5.735630662773291e-018 +4.621833291459620e-017i
    -6.798476950901698e-016 +8.854966604495589e-016i
     8.621573215350509e-007 -1.803273335099220e-006i
    -3.675464540719986e-006 +1.123597030527890e-006i
    -1.862645149230957e-009 +1.164153218269348e-009i
    -1.629814505577087e-009 +9.313225746154785e-010i
    -2.980232238769531e-008                        
                          0 -4.470348358154297e-008i
    -3.779859980568290e-007 +1.975859049707651e-007i
     5.998299457132816e-008 -1.844164216890931e-007i
     1.117587089538574e-008 +7.450580596923828e-009i
    -2.980232238769531e-008 +3.725290298461914e-009i
     2.980232238769531e-008 +2.980232238769531e-008i
     1.490116119384766e-008 -1.490116119384766e-008i
     2.980232238769531e-008 -2.980232238769531e-008i
    -3.725290298461914e-009                        
                          0 -1.490116119384766e-008i
    -1.862645149230957e-009 -1.862645149230957e-009i

[ 本帖最后由 eight 于 2008-1-3 11:45 编辑 ]

评分

1

查看全部评分

发表于 2008-1-3 09:33 | 显示全部楼层

回复 #10 zjj 的帖子

把程序也发上来吧:lol
发表于 2008-1-3 09:49 | 显示全部楼层
发表于 2008-1-3 10:42 | 显示全部楼层

回复 #12 zjj 的帖子

可不可以把他的思路说下?
 楼主| 发表于 2008-1-3 14:18 | 显示全部楼层
强悍,同上所问,如果程序不能发上来的话,那就把思路说一下吧,谢谢啦!:@)
 楼主| 发表于 2008-1-3 14:34 | 显示全部楼层
各位,实在不好意思,昨天发的式子有个错误,刚才看的时候才发现,现将式子更正如下,红笔处为错误,改为“+”号, clip_image002.jpg
大家再看看。谢谢啦!!!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-11 20:57 , Processed in 0.097278 second(s), 32 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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