FPUT問題(Fermi-Pasta-Ulam-Tsingou)の2次の非線形の振動子の動きを Python＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使って計算
昨日はカテナリになるのをみたけれど、ほぼ同じようなプログラムでFPUT問題（過去はFPU問題と言われていた）が計算できる。
やっているのはこちらと同じく、
d2xi/ dt2 = (xi+1 + xi-1 - 2xi) + α[(xi+1 - xi)2- (xi - xi-1)2]
の計算。通常はエネルギーのモードごとの図を描きますが、
実際の動きをGIFアニメにした。
こちら。
これだけ見てても再帰現象はわからんな、、、
« 鎖をそのまま置くと懸垂曲線（カテナリ―）になるというのをPython＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使って計算 | トップページ | 松屋で海老のチリソース定食ポテトサラダセットを食す。海老が予想以上に多く、甘み多めの甘辛のエビチリでした。 »
「学問・資格」カテゴリの記事
- 高周波・RFニュース 2025年1月23日 5G Americasの新ホワイトペーパー「AI時代のセルラーネットワークの信頼性とセキュリティ」、KyoceraAVXの新薄膜フィルタ、TDKの車載／一般用C0G特性1,250V 3225サイズMLCC、Semtechの5G LPWAモジュール(2025.01.23)
- 高周波・RFニュース 2025年1月22日 everythingRFマガジンにMarkiの宇宙向けミリ波部品の記事、NordicのRF52810を使った太陽電池で動き暗闇でも3週間持つアセットトトラッカー、KnowlessのMRIの技術解説記事、Broadcomの3.5Dパッケージング解説(2025.01.22)
- UnityでVisual C#用の数値計算ライブラリMath.NET numericsを使う(3) 3D画面に補間（Interpolate) を行って表示する。リニア、３次スプライン、有理関数などいろいろ使える。(2025.01.23)
« 鎖をそのまま置くと懸垂曲線（カテナリ―）になるというのをPython＋Scipyでルンゲクッタ8次のDOP853(Dormand＆Prince)を使って計算 | トップページ | 松屋で海老のチリソース定食ポテトサラダセットを食す。海老が予想以上に多く、甘み多めの甘辛のエビチリでした。 »
コメント