Excelで5次方程式をDKA法で計算（VBA複素数ライブラリの例）
だいぶ前に作ったExcel VBAで複素数/FFTが使えるライブラリだけど、あまりにも使い方の説明をしてなくて質問サイトでも使い方がわからん、とか質問されてたorz
クリスマスプレゼントとしてなにか簡単な例がないかな、と思って5次方程式でも解いてみようかと。
a5*z^5 + a4*z^4 + a3*z^3 + a2* z^2 + a1*z + a0 =0
（係数も全部複素数でOK）
は解の公式はないけれど、DKA法で数値的に計算できる。
これがそのExcelシート↓
係数を複素数で入力して、計算ボタンを押すと答えも5個、複素数で出てくる。
x=1,2,3,4,5が答えになるような式でチェックしてみた。こうやって収束していくんだな。
めんどくさいんで5次方程式にしたけど、何次でも拡張できるように作ってるので適宜修正してもらえば。
« おおきく振りかぶって Vol. 13を買った。 | トップページ | 高校生がべき和の公式で賞をもらったというのを聞いてほほえましく思う。 »
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年9月17日 IEEE microwave magazine10月号はマイクロ波顕微鏡特集、 Microwave Journalはミリ波パッケージ。Vayyarの4Dイメージングレーダの仕組み、MIPIのA-PHY、など。(2020.09.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年8月17日 Microwave Journal 8月号はドローンで衛星通信をテスト、220GHzの半導体テストシステム。DesignConで導体の粗さモデルと因果律について解説。Qualcommとワシントンポストの5Gクラス、Samsung Galaxy Watch3分解(2020.08.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年8月10日 Microwave Magazine 9月号はシリコンフォトニクスをマイクロ波・ミリ波通信に使う、コロナと5G、IEEE RFICでimecが低消費電力60GHzレーダを発表、Qorvoの5G向けフィルタ解説、など。(2020.08.11)
- Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使う（その16) 時間依存のシュレーディンガー方程式のトンネル効果を計算してGIFアニメに。(2020.07.20)
- Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使う（その15) 蔵本モデルを計算してGIFアニメに。 (2020.07.18)
« おおきく振りかぶって Vol. 13を買った。 | トップページ | 高校生がべき和の公式で賞をもらったというのを聞いてほほえましく思う。 »
コメント