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（その１ 東大寺→奈良の大仏） »
「学問・資格」カテゴリの記事
- 高周波・RFニュース 2026年4月15日 Microwave Journalはアンプと発振器特集、Signal Integrity Journalは100GHz越えのインターコネクトのAIを使うHFSSモデル化、ローデ・シュワルツが潜水艦通信をUDT2026で発表、Xiaomi Poco X8 Pro分解動画、atisの5Gポリシーレポート(2026.04.15)
- 高周波・RFニュース 2026年4月14日 IEEE Microwave Magazineは高周波エンジニア向け量子コンピュータ入門、Antenna and Propagation Magazineはニューラルネット電磁界シミュレーションなど、第106回ARFTG論文公開、QorvoのSバンドスイッチトフィルターバンクなど(2026.04.14)
- RF Weekly Digest (Gemini 3.1 Pro・Google AI Studio BuildによるAIで高周波・RF情報の週刊まとめアプリ）2026/4/5-4/12(2026.04.12)
- GLM-5.1(Ollamaから利用）でPythonのscikit-rfを使ってTouchstoneフォーマットのSパラメータファイルを読んでdB, 位相, スミスチャート, TDRを表示するGUIアプリを作ってもらった。5分など長く考えた後、Gemma 4:31bよりさらに出来が良く、思った通りのものができた。(2026.04.09)
« 京都初詣シリーズ2013（その8 猿田彦神社） | トップページ | 奈良初詣シリーズ2013（その１ 東大寺→奈良の大仏） »
コメント