高次空間差分を手計算で(1階微分)
ちょっと必要になって差分法で用いる空間微分の高次のものを手計算でやってみた。まずは1階微分。
精度の低いものから順番に:
(u(x+h) - u(x) ) / h ; O(h)
(u(x+h) - u(x-h)) / (2*h) ; O(h^2)
(u(x-2h) -8u(x-h) + 8u(x+h) - u(x+2h))/(12*h) ;O(h^4)
ここまではたまに見る。次は
(-u(x-3h)/9 + u(x-2h)-5u(x-h)+5u(x+h) -u(x+2h)+u(x+3h)/9)/(20h/3)
がO(h^6)。
手計算だから間違っている可能性も大なのでお気をつけて。
« 高木関数が出るラプラス方程式 | トップページ | 高次空間差分を手計算で(2階微分) »
「学問・資格」カテゴリの記事
- 高周波回路シミュレータQucsStudioを使ってみる(その5)マルチポート(4ポートだとs4pみたいな)のSパラメータデータをTouchstoneファイルに出力しようとするとExport to SnPでは2ポートになってしまう。Project→Import Dataからポート数を手で入れないと。(2022.05.24)
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2022年5月19日 IEEE Microwave Magazineで磁性体を使わない非可逆デバイス、Microwave Journalは5Gスモールセル、THz、Movanoがミリ波血圧&血糖値計発表、三菱の3Dプリンタ衛星アンテナ、など。(2022.05.19)
- シン・ウルトラマンをIMAXで観てきた。面白かった!物理学者の有岡君がラグランジアン,AdS,プランクブレーンとかしゃべってる!Randall–Sundrumも。マグカップは標準理論+重力で、物理監修は予想通り橋本幸士さんでした!ブラックホールの構造もタイムリー!(2022.05.14)
- 個人の正解確率pとして多数決をとったときの正解確率が、pが0.5に近づくとどのくらいの人数で0.9を超えるか計算してみた。p=0.50005だと日本の人口くらいの多数決。(神とさざなみの密室読んで正当性確率が気になったので)(2022.05.10)
- LTspiceで高周波解析(2)伝送線路のマイクロストリップライン素子を作ったのでそれをスタブとして使ってスタブBPF(バンドパスフィルタ)のSパラメータを計算してみる。QucsStudioのツールを使ってフィルタ合成したもの。(2022.05.09)
コメント