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

kuang333

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

sigma665

没用过,不过在下想是不是可以慢慢试呢?
最后看试出来,哪个合适就哪个

花如月

问题很久前就有讨论

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

再次发帖注意下版规

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

kuang333

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

花如月

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

zjj

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

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

kuang333

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

sigma665

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

kuang333

式子如下所示：

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

zjj

暂时列出这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

把程序也发上来吧:lol

zjj

原帖由 sigma665 于 2008-1-3 09:33 发表

                               
登录/注册后可看大图

把程序也发上来吧:lol

不好意思
这个程序比较大，是别人编写的
我具有使用权，不拥有所有权，不能发上来。

sigma665

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

kuang333

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

kuang333

各位，实在不好意思，昨天发的式子有个错误，刚才看的时候才发现，现将式子更正如下，红笔处为错误，改为“+”号，
大家再看看。谢谢啦！！！