ラプラス方程式をある差分化で数値計算すると高木関数になる例（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・マイクロ波・ミリ波・5G）関連ニュース2020年1月16日 5G、レーダで使われるミリ波フェイズドアレイアンテナ各社まとめ、HuaweiはSAWまで取り込む、iPhone 11 ProがFBIにロック解除される、など。(2020.01.16)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年1月15日 iPhone 5G mmWのAiP、TSMCとともにASEにも発注と報道、Microwave Journal 1月号の特集はレーダ＆アンテナ、ResonantのWhitepaperなど。(2020.01.15)
- ”楕円の円周と一方の径からもう一方の径を求める。”をカシオの高精度計算サイトkeisan.casio.jpに自作式としてUP.(2020.01.13)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年1月13日 IMECがbeyond 5GのRF半導体技術（III-V/CMOS)について解説、ミンチークオはiPhone 5GモデルはSub6GHzもミリ波も同時に出ると語る、など。(2020.01.13)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2020年1月4日 Qorvoが6GHz対応のWiFiを支持＆5G RF for dummies 2nd edを公開、TSMCがiPhone mmW向けにInFO AiPの注文を受ける、Akouticsが軍用レーダ向けBAW出荷など。(2020.01.04)
「日記・コラム・つぶやき」カテゴリの記事
- 1999-2020年までのマグロの初競り価格のプロット(2020.01.05)
- 2020 あけましておめでとうございます。Happy New Yearの方程式(w + log |N| = log|ap^2H + e^aR| + log|y|) を解いてみよう。Scratchでもやってみた。(2020.01.01)
- GoogleとNORADのサンタトラッカー2019スタート。20:00時点でGoogleサンタはニュージーランドに、NORADサンタは北極から出発！(2019.12.25)
- Scratchで描く、クリスマスツリー2種（セルオートマトン、反復関数系）とクリスマスの数学（クリスマスの方程式）(2019.12.24)
- 自動販売機の左上にアンテナがついてるのに初めて気付いた。(2019.12.14)
トラックバック
この記事へのトラックバック一覧です： ラプラス方程式をある差分化で数値計算すると高木関数になる例（F((2i+1)/^(2k+1))=0.5*(F(i/2^k)+F((i+1/2^k)) + Ck）:
« 京都・一乗寺の高安で唐揚げ定食を食す。巨大な唐揚げが3個。ラーメンも美味しい。 | トップページ | ラーメン横綱で麻辣担々麺を食す。 »
コメント