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を買った。 | トップページ | 高校生がべき和の公式で賞をもらったというのを聞いてほほえましく思う。 »
「学問・資格」カテゴリの記事
- Wordleを100回やってみて正解率の分布を見た。1回目で当たる確率をpとして、2回目からはk倍ずつ確率が上がっていく(kp,k²p,k³p…）とすると、私の結果はp=0.0089, k=2.4になった。(2022.05.27)
- 高周波回路シミュレータQucsStudioを使ってみる（その５）マルチポート（４ポートだとs4pみたいな）のSパラメータデータをTouchstoneファイルに出力しようとするとExport to SnPでは2ポートになってしまう。Project→Import Dataからポート数を手で入れないと。(2022.05.24)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2022年5月19日 IEEE Microwave Magazineで磁性体を使わない非可逆デバイス、Microwave Journalは5Gスモールセル、THz、Movanoがミリ波血圧＆血糖値計発表、三菱の3Dプリンタ衛星アンテナ、など。(2022.05.19)
- シン・ウルトラマンをIMAXで観てきた。面白かった！物理学者の有岡君がラグランジアン,AdS,プランクブレーンとかしゃべってる！Randall–Sundrumも。マグカップは標準理論＋重力で、物理監修は予想通り橋本幸士さんでした！ブラックホールの構造もタイムリー！(2022.05.14)
- 個人の正解確率pとして多数決をとったときの正解確率が、pが0.5に近づくとどのくらいの人数で0.9を超えるか計算してみた。p=0.50005だと日本の人口くらいの多数決。（神とさざなみの密室読んで正当性確率が気になったので）(2022.05.10)
« おおきく振りかぶって Vol. 13を買った。 | トップページ | 高校生がべき和の公式で賞をもらったというのを聞いてほほえましく思う。 »
コメント