時間依存シュレーディンガー方程式（波束の散乱）をExcel VBAでルンゲクッタ8次で計算してGIFアニメに。
だいぶ前に書いた話をGIFアニメにしてみる。シュレーディンガー方程式も空間微分を差分化すると普通の常微分方程式になるので、ルンゲクッタ8次のDormand-Prince法で計算できる。
具体的には空間の拡散項は2次の差分
(φi+1 - 2*φi + φi-1) / (dx*dx)
で近似。中央に井戸型ポテンシャルを置いて、左から波束が入ってくるときの計算。
これは”計算物理”にでてたもので、もともとはシッフの例題にあるものだそうです。
量子力学の本で私が勉強したシッフとかメシアとか、朝永さんとかは最近の人は読まないんだろうなあ、、、JJサクライさんのは今でもOK?
トンネル効果の計算もしてみた。それはこちら。
https://sci.tea-nifty.com/blog/2013/01/excelvba8gif-07.html
« Sパラメータのスミスチャートを描くExcelシートを作った。アップロードしときます。 | トップページ | 布亀の救急箱のCMがなぜか大晦日にたくさん流れていますが、あの赤ちゃんは迅くん。 »
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2021年2月16日 IEEE Microwave Magazineの特集はオールデジタルのRFID、Microwave JournalはＥバンド ミリ波通信に衛星や気球を使う話、アメリカの半導体企業がバイデンに投資を迫る、(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)
« Sパラメータのスミスチャートを描くExcelシートを作った。アップロードしときます。 | トップページ | 布亀の救急箱のCMがなぜか大晦日にたくさん流れていますが、あの赤ちゃんは迅くん。 »
コメント