高次ニュートン法(Householder法)をExcelで試してみる。
最近、twitterでニュートン法より高次の反復計算があるのを知った。それとはちょっと違うけど、同様に高次のニュートン法が↓に記載されていたので試してみよう。
http://www.sztaki.hu/~bozoki/oktatas/nemlinearis/SebahGourdon-Newton.pdf
Householder法は
xn+1= xn + (p+1) *( (1/f) (p)/ (1/f) (p+1) )xn
です。(1/f) (p)はfの逆数のp階微分。
ではp=0のときは普通のニュートン法になる。2次収束
xn+1 = xn - f '/f
p=1のときは3次収束で
xn+1 = xn - f / (f ' - f ''*f /(2* f ' ) )
p=2のときは4次収束で
xn+1 = xn - f * (f ' ^2 - f ''*f / 2)/(f ' ^3 - f * f ' * f '' + f '''*f^2 /6)
ここまでは頑張って手計算。もう無理。なのでMaximaでやってみよう。
p=3のときは5次収束で、
xn+1=xn+(4*f(xn)*(f(xn)^2*(f(xn) ''')-6*f(xn)*(f(xn) ')*(f(xn) '')+6*(f(xn) ')^3))/(f(xn)^3*(f(xn) '''')-8*f(xn)^2*(f(xn) ')*(f(xn) ''')-6*f(xn)^2*(f(xn) '')^2+36*f(xn)*(f(xn) ')^2*(f(xn) '')-24*(f(xn) ')^4)
でした。
ではExcelで計算して比較してみる。√10の計算を初期値1としてグラフにしてみた。
予想通り収束は速いなあ。でもそのtwitterでも言われていたけど、結局高次の方法は乗算数など多くなるので有利ともいえないという。でもおもしろいな。
« 「ハッピーエンドにさよならを」を読んだ。 | トップページ | 「法月綸太郎の新冒険」を読んだ。 »
「学問・資格」カテゴリの記事
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2021年2月16日 IEEE Microwave Magazineの特集はオールデジタルのRFID、Microwave JournalはEバンド ミリ波通信に衛星や気球を使う話、アメリカの半導体企業がバイデンに投資を迫る、(2021.02.17)
- カオスを生じる電気回路、Chua’s circuitをLTspiceで回路シミュレーションしてみる。(2021.02.19)
- Labyrinth Chaos(迷宮カオス)を生むThomas-Rössler方程式のパラメータbを色々変えて、Python+Scipyでルンゲクッタ8次のDOP853(Dormand&Prince)を使って計算してGIFアニメ(2021.02.16)
- フィッツヒュー・南雲 (FitzHugh-Nagumo) 方程式をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算。(2021.02.23)
- 「水晶振動子の等価回路計算」をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP! インピーダンスの大きさと位相がグラフ化できる。(2021.02.12)
コメント