KdV方程式をExcel VBA＆擬スペクトル法＆高次ルンゲクッタ法で計算
この前、新PCを購入したので今まで時間がかかるかと思って敬遠していたような計算をしてみよう。まずは偏微分方程式を空間部分は擬スペクトル法、時間部分はルンゲクッタ法（ここでは6段5次のCash-Karp法）という手法を使ってExcel VBAでプログラムを書く。題材は一番ポピュラーなKdV方程式にしてみる。
で早速計算した結果がこちら。
オリジナルのザブスキーとクルスカルの数値計算は蛙飛び（の工夫したもの）で、それより今回の結果はやっぱりパルスがシャープに計算できてる。でも逆にあの時代にああいった計算をしたのはすごいな。前にうちのブログでも3次精度の風上差分で計算したときはパルス鈍りまくりだったし。
時間発展を追っていったものはこちら。これもきれいに計算できてる。
« 「どんどん橋、落ちた」を読んだ。 | トップページ | 蔵本・シバシンスキー方程式をExcel VBA&擬スペクトル法&高次ルンゲクッタ法で計算 »
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年1月16日 5G、レーダで使われるミリ波フェイズドアレイアンテナ各社まとめ、HuaweiはSAWまで取り込む、iPhone 11 ProがFBIにロック解除される、など。(2020.01.16)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年1月15日 iPhone 5G mmWのAiP、TSMCとともにASEにも発注と報道、Microwave Journal 1月号の特集はレーダ＆アンテナ、ResonantのWhitepaperなど。(2020.01.15)
- ”楕円の円周と一方の径からもう一方の径を求める。”をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUP.(2020.01.13)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年1月13日 IMECがbeyond 5GのRF半導体技術（III-V/CMOS)について解説、ミンチークオはiPhone 5GモデルはSub6GHzもミリ波も同時に出ると語る、など。(2020.01.13)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年1月4日 Qorvoが6GHz対応のWiFiを支持＆5G RF for dummies 2nd edを公開、TSMCがiPhone mmW向けにInFO AiPの注文を受ける、Akouticsが軍用レーダ向けBAW出荷など。(2020.01.04)
« 「どんどん橋、落ちた」を読んだ。 | トップページ | 蔵本・シバシンスキー方程式をExcel VBA&擬スペクトル法&高次ルンゲクッタ法で計算 »
コメント