KdV方程式をExcel VBAで計算してGIFアニメにしてみた。
最近、GIFアニメに凝っている。方程式の解が目に見える形にできるのがいいなあ、と。
今回はKdV方程式
∂u/∂t +u*∂u/∂x+δ^3 ∂^3u/∂x^3 = 0
を計算。
方法は、空間微分をO(h^6)という高次の差分で近似、かつ時間微分をルンゲクッタ8次のスキームであるDormand-Prince法で計算している。
（内容としてはオリジナルのザブスキーとクルスカルのものと同じ）。
ソリトンが複数衝突を繰り返す様子がよくわかる。
« 京都初詣シリーズ2013（その8 猿田彦神社） | トップページ | 奈良初詣シリーズ2013（その１ 東大寺→奈良の大仏） »
「学問・資格」カテゴリの記事
- 温度と湿度を高温・高湿にして加速試験をするときに相対湿度を使う「温湿度加速モデル(Hallberg-Peck model)」をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP!(2021.09.18)
- スマホRF駆け出しエンジニアの豆知識: ⑩パワーアンプの近くにある電源系のICはエンベロープトラッカー。ETモードとAPTモードがある。(2021.09.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2021年9月15日 iPhone13の内部構造、IEEE Microwave magazineの特集はマイクロ波と機械学習、Microwave Journalは車載へのUWB応用、Huaweiは2030年に6Gを狙う、など。(2021.09.16)
- 高周波回路シミュレータQucsStudioを使ってみる（その４） 電磁界シミュレーション（FDTD）機能を使ってKeysight ADSのMomentumの例題のフィルタ（スタブLPF)を計算。(2021.09.13)
- コラッツ予想を複素数に拡張したものをPythonで図示してみる。(2021.09.07)
« 京都初詣シリーズ2013（その8 猿田彦神社） | トップページ | 奈良初詣シリーズ2013（その１ 東大寺→奈良の大仏） »
コメント