リッカチ方程式をドルマン・プリンスDOP853で計算(Excel VBA)
だいぶ前にリッカチ/リカッチ論争にけりがついたリッカチ(Riccati)方程式ですが、
具体的に解こうと思ってすっかり忘れていた。
まずは、狭義のリッカチ方程式、
y' + a * y^2 = b * x^m
を計算しよう。mが特別な値、
m=-2, 4k/(1-2k) (ただしkは整数)
のときは厳密に初等的に計算できる。(たぶん、こういう話が「リッカチのひ・み・つ」に載っているんだろうと思うが、何せ読んでない。。。ベッセル関数を使えば一般解が求められるが、それもめんどくさい。。。)
でめんどくさいので、まずは数値的に計算。
« マラゴニー効果って何?(トイレスタンプクリーナー)、というかバブル君効果って何? | トップページ | 神足神社でお参り »
「学問・資格」カテゴリの記事
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2021年2月16日 IEEE Microwave Magazineの特集はオールデジタルのRFID、Microwave JournalはEバンド ミリ波通信に衛星や気球を使う話、アメリカの半導体企業がバイデンに投資を迫る、(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)
« マラゴニー効果って何?(トイレスタンプクリーナー)、というかバブル君効果って何? | トップページ | 神足神社でお参り »
コメント