シュレディンガーじゃないほうの波動方程式の数値解をExcelで

くすぶっていてもしゃーない、とりあえず中間報告だ！


シュレディンガー波動方程式、
iħ∂u/∂t=[(-iħ∇)^2/(2m)+V]u

じゃなくて純粋な波動方程式
∂^2u/∂t^2=c^2∇^2u
のほうのやつ、とりあえずようやく数値的に解けた！

Excelで！

といってもまだまだ一次元で
こんな感じのやつです。


ホントはこいつ、循環参照でやりたいんですわ。

でもなんか不安で、何日もくすぶってるんですわ
なんでこんなことでくすぶってんのって自分でも思うんですけど、なんか不安なんですよね。


上の表で下の表をvlookupで列参照してるんです(横が距離xで縦が時間tってことになってます)。
下の表自体はすぐに収束して静的なんすよ。この時点では循環参照まったく使ってません。

下の表の1行目と2行目が位置と速度の初期条件で、右端と左端がそれぞれ境界条件なんすわ
今回は左端が波源でsinを入れといて、右端で固定タン反射させてます。
「右端=0」の境界条件で固定タン反射、「右端＝すぐ左」で自由タン反射なのはわかったんですが
そのあいだの、たとえばインピーダンスマッチングを行なって反射しない場合の境界条件がまだよくわからないんですよね


それと、時間に1次元配列使ってるのがもったいなくて・・・
3行目以降の白セルの中身は同じ演算で2個上と1つ上、それと斜め右と斜め左しか参照してないんで
ぶっちゃけ3行目と4行目をお互いに循環参照しあえば出来ちゃいそうな気もするんですよ

そしたら縦のセルが空くんで、満を持して横：x、縦：yの2次元の波にしたいんですよねぇ

いやーこれもZ状参照が許してくれるかどうかわかんないんですが。

課題が多いんですけどやるべきことはほぼわかってるし、地道にやればいいだけなんですけどね
どうも調子が出ません。


出ないので説明も兼ねて中間報告なんすよ


「シュレディンガータイプばっかり目立っちまいやがって・・・！
俺のほうが本家だっつーのによ」

にほんブログ村