kuang333 发表于 2008-1-2 15:29

如何确定牛顿迭代的初始值

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

sigma665 发表于 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 编辑 ]

kuang333 发表于 2008-1-2 16:44

我已经看了之前相关的帖子了,可是现在有个关键的问题是我想求复数解,这个范围不好定啊,相应的初值也就不好定,类似于这方面的问题各位之前有过研究吗?请指教!

花如月 发表于 2008-1-2 16:54

回复 #4 kuang333 的帖子

复数解你的迭代公式怎么写?

zjj 发表于 2008-1-2 18:35

方程不是很复杂的话,我可以帮你求

[ 本帖最后由 eight 于 2008-1-2 18:44 编辑 ]

kuang333 发表于 2008-1-2 19:02

迭代公式还是常规的牛顿迭代公式,只不过将解的范围扩大到了复数范围。
式子还算比较复杂,如果有思路的话也可直接与我联系,我的QQ是271659505,谢谢!

sigma665 发表于 2008-1-2 19:05

回复 #7 kuang333 的帖子

建议你把方程用图片的形式,附件上来

kuang333 发表于 2008-1-2 19:34

式子如下所示:

u为变量,就是要找u的初值。

zjj 发表于 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 编辑 ]

sigma665 发表于 2008-1-3 09:33

回复 #10 zjj 的帖子

把程序也发上来吧:lol

zjj 发表于 2008-1-3 09:49

原帖由 sigma665 于 2008-1-3 09:33 发表 http://www.chinavib.com/forum/images/common/back.gif
把程序也发上来吧:lol
不好意思
这个程序比较大,是别人编写的
我具有使用权,不拥有所有权,不能发上来。

sigma665 发表于 2008-1-3 10:42

回复 #12 zjj 的帖子

可不可以把他的思路说下?

kuang333 发表于 2008-1-3 14:18

强悍,同上所问,如果程序不能发上来的话,那就把思路说一下吧,谢谢啦!:@)

kuang333 发表于 2008-1-3 14:34

各位,实在不好意思,昨天发的式子有个错误,刚才看的时候才发现,现将式子更正如下,红笔处为错误,改为“+”号,
大家再看看。谢谢啦!!!
页: [1] 2
查看完整版本: 如何确定牛顿迭代的初始值