スティフなvan der Pol方程式を陰的ルンゲクッタで計算(Excel VBA)
van der Pol方程式
y1' = y2
y2' = ((1-y1^2)*y2 - y1) /ε
は、εが1程度なら普通の陽的ルンゲクッタでOK。でもε= 1.e-6みたいなときは非常に硬い(Stiffな）方程式になってしまう。そこで陰的ルンゲクッタで計算。スキームはラダウIIA。
ただし、ステップ幅は可変にしないとめちゃくちゃ効率悪い。で、計算結果はこちら。
そりゃy2が0から1.e6まで急激に変わるんだから解きにくいはずだわ。
« 「パプリカ」を読んだ。 | トップページ | 「Self-Reference ENGINE」を読んだ。 »
「学問・資格」カテゴリの記事
- 高周波・RFニュース 2025年1月13日 IEEE Microwave Magazineの特集はニューラルネットワークとマイクロ波、Siversがミリ波ビームフォーマー開発を受注、バイデン・ハリス政権が ワイヤレス革命に1億 1,700 万ドル、HoneywellとNXPが航空機技術で提携(2025.01.13)
- UnityでVisual C#用の数値計算ライブラリMath.NET numericsを使う(1) まずはNuGetForUnityを使ってインストール。2Dゲーム画面に連立方程式を解いた結果を表示。(2025.01.14)
- 高周波・RFニュース 2025年1月9日 CES2025に合わせて各社プレスリリース、特にQualcomm、NVIDIA、INTELが大量。SEMCOのC0G MLCC (1210 inch, 22nF, 1000V)解説、TIのAI搭載60GHz車内レーダ、MarvellのCPO、Qorvoの車載UWB SoC、TDKのセンサがAI白杖に採用(2025.01.09)
- NHK パンサー尾形さんの笑わない数学 微分・積分 スペシャルがもうすぐ始まる。これから見てリアルタイムでポストしたのでそのスレッドを残す。(2024.12.29)
コメント