移流方程式を高次ルンゲクッタで計算してみる（その１）
∂u/∂t=c∂u/∂xを計算するのに時間方向は5次のルンゲクッタ、空間は4次差分を使ってみよう。4次差分は
u'(x) ≒ (u[i-2] - 8u[i-1] + 8u[i+1] - u[i+2]) / (12*Δx）で計算。初期条件は矩形波。すると、
やっぱり波形がくずれていく。これは風上差分を使わないといくら高次までとってもだめということですな。予想通り。
« 拡散方程式を高次ルンゲクッタで計算してみる（その2）。 | トップページ | 移流方程式を高次ルンゲクッタで計算してみる（その2） »
「学問・資格」カテゴリの記事
- 高周波・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のantennaArrayDesigner機能を使って一行もスクリプトを書かずにパッチアンテナアレイを設計してSパラメータ、指向性などを計算する。1素子とちがって計算にはかなり時間がかかるのでとりあえず2素子のアレイで。(2026.03.09)
« 拡散方程式を高次ルンゲクッタで計算してみる（その2）。 | トップページ | 移流方程式を高次ルンゲクッタで計算してみる（その2） »
コメント