拡散方程式の解をkeisan.casio.jpにUP!
Wolfram|AlphaやMicrosoft Math3.0に浮気していたけれど、久しぶりに高精度計算サイトkeisan.casio.jpに自作式をUPした。これです↓
一次元拡散方程式 ∂φ/∂t= D∂^2φ/∂x^2の解を計算します。
初期条件は高さA,幅εの矩形とします。
計算結果はこんな感じ。
どうやって計算しているかというと、初期条件が矩形なら、誤差関数erf(x)を用いて
φ(x,t) = A/2 (erf((x+ε/2)/√(4Dt)) - erf((x-ε/2)/√(4Dt)) )
と厳密に解ける。ここで
erf(x) = (2/√π)∫[0,x] exp(-y^2) dy
です。誤差関数はkeisan.casio.jpにちゃんと入っているので簡単。
« インスタントラーメン発明記念館に行ってきた。 | トップページ | 一文橋のフジタカの看板もウインクする。 »
「学問・資格」カテゴリの記事
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2021年2月16日 IEEE Microwave Magazineの特集はオールデジタルのRFID、Microwave JournalはEバンド ミリ波通信に衛星や気球を使う話、アメリカの半導体企業がバイデンに投資を迫る、(2021.02.17)
- カオスを生じる電気回路、Chua’s circuitをLTspiceで回路シミュレーションしてみる。(2021.02.19)
- Labyrinth Chaos(迷宮カオス)を生むThomas-Rössler方程式のパラメータbを色々変えて、Python+Scipyでルンゲクッタ8次のDOP853(Dormand&Prince)を使って計算してGIFアニメ(2021.02.16)
- フィッツヒュー・南雲 (FitzHugh-Nagumo) 方程式をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算。(2021.02.23)
- 「水晶振動子の等価回路計算」をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP! インピーダンスの大きさと位相がグラフ化できる。(2021.02.12)
コメント