ラプラス方程式をある差分化で数値計算すると高木関数になる例（F((2i+1)/^(2k+1))=0.5*(F(i/2^k)+F((i+1/2^k)) + Ck）
だいぶ前に書いたものの再掲：
昔、”無限・カオス・ゆらぎ”という本を読んで非常に面白かったのを覚えている。
その中でいちばんおもしろかったのは、差分化の方法によっては
元の微分方程式とは似ても似つかないカオス解が出ると言うところで、
その一つはロジスティック方程式の中心差分。
さらにもう一つ、”変なラプラシアン”というのが書いてあった。
（こちらの文献でも同様な話が読めます。
）
ラプラス方程式は、ΔF(x)=0ですが、
一次元だと単にd2 F / dx2 =0。
単純な差分化だとF(x) = 0.5 * (F(x+Δx) + F(x-Δx)) ということになる。これをメッシュをちょっと工夫して
F((2i+1)/2k+1)=0.5*(F(i/2k)+F((i+1）/2k)) + Ck
というように書きなおしてみる。
つまりあるkの点列が得られたとき、
次のk+1ではある点のデータは前のkの値で、
その点と隣り合ったデータを使って求めると見る（マルチグリッド）
でCk=1/2^(k+1)と選ぶ。Ckはk→∞では消えるはずの項。
（山口さんは、これはポアソン方程式じゃなくて変なラプラシアンだと言いきってますが。）
でkをだんだん増やして行った時のExcelのグラフをGIFアニメにしたのがこちら。
しかもこれはおなじくCkをk→∞で消える項1/4^kと選ぶと、
なだらかになるという。。。
単純な数値計算でも恐ろしいことになる例は多い。。。
« 京都・一乗寺の高安で唐揚げ定食を食す。巨大な唐揚げが3個。ラーメンも美味しい。 | トップページ | ラーメン横綱で麻辣担々麺を食す。 »
「学問・資格」カテゴリの記事
- 高周波・RFニュース 2026年1月16日 Microwave JournalにDirect RF特集記事、NXPが5G RF事業を縮小、RCRTechのテスト・計測についてのレポート、MediaTekがDimensity 9500sと8000を発表、TDKの車載チップバリスタ、I-PEXのロック付きRF同軸コネクタ(2026.01.16)
- 高周波・RFニュース 2026年1月15日 IEEE Microwave MagazineはCMOS LNAレビューなど、Journal of Microwavesはアボカドのマイクロ波イメージングなど、LitePointが3つのワイヤレストレンドを調査、Ericsson ConsumerLabのCSP差別化レポート、Motorola Moto G Power 2026 分解(2026.01.15)
「日記・コラム・つぶやき」カテゴリの記事
- 高周波・RFニュース 2026年1月16日 Microwave JournalにDirect RF特集記事、NXPが5G RF事業を縮小、RCRTechのテスト・計測についてのレポート、MediaTekがDimensity 9500sと8000を発表、TDKの車載チップバリスタ、I-PEXのロック付きRF同軸コネクタ(2026.01.16)
- 高周波・RFニュース 2026年1月15日 IEEE Microwave MagazineはCMOS LNAレビューなど、Journal of Microwavesはアボカドのマイクロ波イメージングなど、LitePointが3つのワイヤレストレンドを調査、Ericsson ConsumerLabのCSP差別化レポート、Motorola Moto G Power 2026 分解(2026.01.15)
トラックバック
この記事へのトラックバック一覧です： ラプラス方程式をある差分化で数値計算すると高木関数になる例（F((2i+1)/^(2k+1))=0.5*(F(i/2^k)+F((i+1/2^k)) + Ck）:
« 京都・一乗寺の高安で唐揚げ定食を食す。巨大な唐揚げが3個。ラーメンも美味しい。 | トップページ | ラーメン横綱で麻辣担々麺を食す。 »
コメント