« 「パプリカ」を読んだ。 | トップページ | 「Self-Reference ENGINE」を読んだ。 »

2010年7月18日 (日)

スティフなvan der Pol方程式を陰的ルンゲクッタで計算(Excel VBA)

van der Pol方程式

y1' = y2

y2' = ((1-y1^2)*y2 - y1) /ε

は、εが1程度なら普通の陽的ルンゲクッタでOK。でもε= 1.e-6みたいなときは非常に硬い(Stiffな)方程式になってしまう。そこで陰的ルンゲクッタで計算。スキームはラダウIIA。

ただし、ステップ幅は可変にしないとめちゃくちゃ効率悪い。で、計算結果はこちら。

Vanderpolstiff

そりゃy2が0から1.e6まで急激に変わるんだから解きにくいはずだわ。

« 「パプリカ」を読んだ。 | トップページ | 「Self-Reference ENGINE」を読んだ。 »

学問・資格」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

トラックバック

この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/48735906

この記事へのトラックバック一覧です: スティフなvan der Pol方程式を陰的ルンゲクッタで計算(Excel VBA):

« 「パプリカ」を読んだ。 | トップページ | 「Self-Reference ENGINE」を読んだ。 »

最近のコメント

2018年11月
        1 2 3
4 5 6 7 8 9 10
11 12 13 14 15 16 17
18 19 20 21 22 23 24
25 26 27 28 29 30  
フォト
無料ブログはココログ