高木関数が出るラプラス方程式をGIFアニメで。
以前にも書いたけれど、”無限・カオス・ゆらぎ”という本を読んで若いころ影響をうけました。
その中でいちばんおもしろかったのは、差分化の方法によっては元の微分方程式とは似ても似つかないカオス解が出ると言うところで、その関連として”変なラプラシアン”というのが書いてありました。
(こちらの文献でも同様な話が読めます。
http://repository.kulib.kyoto-u.ac.jp/dspace/bitstream/2433/92310/1/KJ00004774034.pdf
)
ラプラス方程式は、ΔF(x)=0ですが、一次元だと単にd^2 F / dx^2 =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アニメにしたのがこちら。
kが大きくなるにつれて高木関数に。初めてこれを見た時は面白かったなあ。
GIFアニメにはGIAMを使わせてもらいました。
http://homepage3.nifty.com/furumizo/giamd.htm
わざわざ画像データをファイルとして落とす必要すらなく、エクセルでグラフを選択して張り付ければいいだけできわめて簡単。
« 映画”007 スカイフォール”を観てきた。 | トップページ | 永久にビールが循環するマシンができた?! »
「学問・資格」カテゴリの記事
- TensorFlow(Keras)でモデルをsaveで保存してload_modelで読み込むときに突然エラー(utf8で読めないとかディレクトリが存在しないとか)が出始めた。なんで?といろいろやると、単に日本語がフォルダ名に使われているときにエラーになるだけだった…(Windowsネイティブの場合)(2023.05.28)
- 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(4) 10進31桁まで計算できることを生かして4段4次のルンゲクッタ法でローレンツ方程式を計算して図示してみる。テキストベース(アスキーアート)で!(2023.05.23)
- 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(3) ロジスティック写像の分岐図を31桁まで10進計算ができることを生かして描いてみる。テキストベースで!今回は2次元配列の練習。(2023.05.19)
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース(5/18) IEEE Microwave Magazineはグラフェンのマイクロ波応用、Microwave Journalはテラヘルツ波の発生、Next G allianceの新6Gレポート、Hexa-Xの新動画、ESAが軌道上に6Gラボ、De-embeddingって何?など。(2023.05.18)
- カシオの高精度計算サイトkeisan.casio.jpが復旧してログインや自作式の作成ができるようになった。ということで作った自作式のアクセスベスト10とお勧めを紹介。(2023.05.17)
コメント