ExcelでKdV方程式を解いてみる。
ちょっと前にザブスキーとクルスカルのオリジナルの計算法を使って解いてみた。
https://sci.tea-nifty.com/blog/2008/05/kdvexcel_99e6.html
最初が正弦波の時は、
のようになって論文（Physical Revew Letters,vol.15 No.6(1965))と同じ結果に。
では3次精度風上差分を使うとどうなるか？
山の高さが低くなってるなあ。やっぱりちょっとなまってるのかな。でも工夫せずにそのまま使えそうなのはこっち？
« 亀in桂川(IXY 20IS) | トップページ | POs-CaのCMソングは中塚武のKiss&Ride »
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2021年2月16日 IEEE Microwave Magazineの特集はオールデジタルのRFID、Microwave JournalはＥバンド ミリ波通信に衛星や気球を使う話、アメリカの半導体企業がバイデンに投資を迫る、(2021.02.17)
- カオスを生じる電気回路、Chua’s circuitをLTspiceで回路シミュレーションしてみる。(2021.02.19)
- Labyrinth Chaos(迷宮カオス）を生むThomas-Rössler方程式のパラメータbを色々変えて、Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使って計算してGIFアニメ(2021.02.16)
- フィッツヒュー・南雲 (FitzHugh-Nagumo) 方程式をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算。(2021.02.23)
- 「水晶振動子の等価回路計算」をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP！ インピーダンスの大きさと位相がグラフ化できる。(2021.02.12)
コメント