ローレンツ方程式3種類（通常のLorenzモデル、 Lorenz84モデル、Lorenz96モデル）をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算（その18）
このサイトを見ていて、
https://www.dynamicmath.xyz/strange-attractors/
Lorenz84というモデルがあることを知った（って上のサイトは誤植でLorenz83と書いてありますが）。
こんな式：
http://analog-ontology.blogspot.com/2015/02/lorenz-84-system.html
そういや36変数のLorenz96モデルも計算したことがある。
Python+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)を使う（その４） ローレンツ96モデル 36変数の方程式
dx_i/dt=(x_i+1 - x_i-2)*x_i-1 -x_i +F
という式です。
で通常のLorenz方程式は
Python+Scipyでルンゲクッタ8次のDOP853(Dormand&Prince)を使う（その2） やはり最初はローレンツ方程式。
で計算したので、全部まとめてDOP853で計算して図示してみる。
通常：
Lorenz84
Lorenz96
もう少し長生きしてもらえれば今年のノーベル物理学賞も受賞のはず、ということで返す返すも残念。。。
« Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使う(その17) 4つの羽をもつFour-wing attractorを描く。ところでMatplotlib3.4.0からfig.gca(projection=”3d”)が警告が出てfig.add_subplot(projection=”3d”)にしないとだめっぽい。 | トップページ | 和泉中央駅近くの中央公園をぶらぶら。メアンダパターンのようなうねった道があった。 »
「学問・資格」カテゴリの記事
- アメリカの空港で5G C-bandをを使うなという話で揉めてますが、どう考えても電波高度計(4.2-4.4GHz）が無茶苦茶で、衛星通信(3.7-4.2GHz)を邪魔しているのにC-Band(3.7–3.98GHz)が入ってくるとクレームとは、、、ジャイアニズムか！(2022.01.22)
- 不定方程式a*x+b*y=gcd(a,b) （ディオファントス方程式、ベズーの等式）の拡張ユークリッド互除法を使った計算をカシオの高精度計算サイトkeisan.casio.jpにUP!2022年共通テストの数学I・数学Aの問題4に出てても結果に自信ないので機械に頼る。。。(2022.01.19)
- 2022年度共通テストの数学I・数学Aの問題、例の防衛省のイージス・アシュアの縮尺違い地図で角度間違ってる話をそのまま持ってきてる！秋田県男鹿市の本山か。。。(2022.01.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2022年1月18日 IEEE Microwave Magazineは女性高周波研究者特集、Microwave Journalはレーダで月面が見える。ドコモの誘電体線路アンテナ、3.45GHzのオークション結果、hBNを使ったスイッチ、Scikit-rfの使い方など。(2022.01.18)
- 犬や猫などのペットにマイクロチップを埋め込むことが義務化されるが、どんな無線通信技術を使っているか調べた。ISO11784/11785に準拠した134.2kHzを使うRFIDで、フェライトロッドアンテナ、キャパシタ、チップ内蔵。15桁の個体識別番号が読み取れる。(2022.01.10)
「日記・コラム・つぶやき」カテゴリの記事
- 新型コロナウイルス、日本の陽性者数＆ワクチン接種者数総計をプロット＆中国、韓国、アメリカ、ドイツ、フランス、イギリスの陽性者数もプロット（1/23更新）。オミクロン株で日本の垂直増加がやばい。私の周りでも陽性者が出始めた。(2022.01.23)
- アメリカの空港で5G C-bandをを使うなという話で揉めてますが、どう考えても電波高度計(4.2-4.4GHz）が無茶苦茶で、衛星通信(3.7-4.2GHz)を邪魔しているのにC-Band(3.7–3.98GHz)が入ってくるとクレームとは、、、ジャイアニズムか！(2022.01.22)
- 不定方程式a*x+b*y=gcd(a,b) （ディオファントス方程式、ベズーの等式）の拡張ユークリッド互除法を使った計算をカシオの高精度計算サイトkeisan.casio.jpにUP!2022年共通テストの数学I・数学Aの問題4に出てても結果に自信ないので機械に頼る。。。(2022.01.19)
- 2022年度共通テストの数学I・数学Aの問題、例の防衛省のイージス・アシュアの縮尺違い地図で角度間違ってる話をそのまま持ってきてる！秋田県男鹿市の本山か。。。(2022.01.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2022年1月18日 IEEE Microwave Magazineは女性高周波研究者特集、Microwave Journalはレーダで月面が見える。ドコモの誘電体線路アンテナ、3.45GHzのオークション結果、hBNを使ったスイッチ、Scikit-rfの使い方など。(2022.01.18)
« Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使う(その17) 4つの羽をもつFour-wing attractorを描く。ところでMatplotlib3.4.0からfig.gca(projection=”3d”)が警告が出てfig.add_subplot(projection=”3d”)にしないとだめっぽい。 | トップページ | 和泉中央駅近くの中央公園をぶらぶら。メアンダパターンのようなうねった道があった。 »
コメント