高次ニュートン法(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ニュース 2025年6月25日 Ericssonモビリティレポート6月号は5G FWAのマネタイズについて、QorvoがオートロックへのUWB応用について解説、NATOのダイアナチャレンジは「先進通信技術」と「電磁環境の競合」、TDKがQEIよりRF電源事業譲受(2025.06.25)
- 高周波・RFニュース 2025年6月20日 GSMA MWC25上海開催、3GPPのCCWでのプレゼン資料、Kyocera AVXのNB-NTN向けIoTデバイスのホワイトペーパー、TDKがAI・スマートグラスのSoftEye買収、Mini-CircuitsのAEC-Q200 qualifiled LTCC製品、Fibocomのスマートリビング向けFG390(2025.06.20)
- 高周波・RFニュース 2025年6月19日 QorvoがSバンドレーダ用のBAW switched filter bank発表、Ericssonがミッションクリティカル用途のアンテナ発表、SEMCOが125℃保証の0201インチX7T 1.0㎌ 6.3V MLCC発表、iFixitがトルクスプラスねじについて解説(2025.06.19)
- 高周波・RFニュース 2025年6月18日 Qorvoが5Gインフラ向けBAWフィルタとプリドライバアンプ発表、KeysightとNTTらが300GHz帯で280Gbpsを達成する信号発生システム発表、TDKが自動車用パワー・オーバー・コアクス・インダクター発表、NordicがNeuton AI買収(2025.06.18)
- 高周波・RFニュース2025年6月17日 everythingRF magazineはIMS2025特別号、MITの6Gに向け光でディープラーニングを行うチップ論文、NGMNが6Gに向けたキーメッセージを出版、Litepoint、Spirent、ViaviのTest and Measurementのトレンドレポート(2025.06.17)
コメント