ローレンツ方程式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)を使う(その4) ローレンツ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”)にしないとだめっぽい。 | トップページ | 和泉中央駅近くの中央公園をぶらぶら。メアンダパターンのようなうねった道があった。 »
「学問・資格」カテゴリの記事
- 高周波・RFニュース 2026年3月11日 STMicroelectronicsが新UWBチップ発表、Ericssonが主導のVICTOR6G発足、Silicon LabsのBluetooth SoCがBANFのタイヤモニタリングシステムに採用、京セラが新しい差動クロック水晶発振器を発表など(2026.03.11)
- 高周波・RFニュース 2026年3月10日 IEEE Microwave MagazineはHF-VHF-UHF特集、Pythonの高周波ライブラリscikit-rfがv1.11.0に、Samsung Galaxy S26 Ultra分解動画、フジクラが4000心SWR/WTC製品化、Perasoの60GHzモジュールが軍用ドローン識別に採用など(2026.03.10)
- RF Weekly Digest (Gemini 3.1 Pro・Google AI Studio BuildによるAIで高周波・RF情報の週刊まとめアプリ)2026/3/1-3/8(2026.03.08)
- MATLAB OnlineでAntenna ToolboxのantennaDesigner機能を使って様々なアンテナ(ホーン、フラクタル(スノーフレーク)パッチ、Vivaldi、誘電体共振器)のSパラメータ、指向性を計算する。(2026.03.11)
- MATLAB OnlineでAntenna ToolboxのantennaArrayDesigner機能を使って一行もスクリプトを書かずにパッチアンテナアレイを設計してSパラメータ、指向性などを計算する。1素子とちがって計算にはかなり時間がかかるのでとりあえず2素子のアレイで。(2026.03.09)
「日記・コラム・つぶやき」カテゴリの記事
- 高周波・RFニュース 2026年3月11日 STMicroelectronicsが新UWBチップ発表、Ericssonが主導のVICTOR6G発足、Silicon LabsのBluetooth SoCがBANFのタイヤモニタリングシステムに採用、京セラが新しい差動クロック水晶発振器を発表など(2026.03.11)
- 高周波・RFニュース 2026年3月10日 IEEE Microwave MagazineはHF-VHF-UHF特集、Pythonの高周波ライブラリscikit-rfがv1.11.0に、Samsung Galaxy S26 Ultra分解動画、フジクラが4000心SWR/WTC製品化、Perasoの60GHzモジュールが軍用ドローン識別に採用など(2026.03.10)
- RF Weekly Digest (Gemini 3.1 Pro・Google AI Studio BuildによるAIで高周波・RF情報の週刊まとめアプリ)2026/3/1-3/8(2026.03.08)
- MATLAB OnlineでAntenna ToolboxのantennaDesigner機能を使って様々なアンテナ(ホーン、フラクタル(スノーフレーク)パッチ、Vivaldi、誘電体共振器)のSパラメータ、指向性を計算する。(2026.03.11)
- MATLAB OnlineでAntenna ToolboxのantennaArrayDesigner機能を使って一行もスクリプトを書かずにパッチアンテナアレイを設計してSパラメータ、指向性などを計算する。1素子とちがって計算にはかなり時間がかかるのでとりあえず2素子のアレイで。(2026.03.09)
« Python+Scipyでルンゲクッタ8次のDOP853(Dormand&Prince)を使う(その17) 4つの羽をもつFour-wing attractorを描く。ところでMatplotlib3.4.0からfig.gca(projection=”3d”)が警告が出てfig.add_subplot(projection=”3d”)にしないとだめっぽい。 | トップページ | 和泉中央駅近くの中央公園をぶらぶら。メアンダパターンのようなうねった道があった。 »






コメント