20080511~
13と7と11の倍数の論理積は13と7と11の積の倍数である。
和ァ・・・
× [PR]上記の広告は3ヶ月以上新規記事投稿のないブログに表示されています。新しい記事を書く事で広告が消えます。
いらっしゃいませこんにちは、期待値観測するには確率の問題に超絶疎い、量子きのこです
好奇心はあるのにメンタル不足でソレ系の職業に就けなかったので 何か出来ないかと思ったわけじゃないですが 結果的に物理や数学とかの面白さを伝えられたらいいなー なブログになりました。 このブログをきっかけに、一人でも多くの人が楽しんで物理や数学を学んでくれるとサイワイです。ヽ゚ー゚ノ サイエンスファンタジーなアニメの感想も書きます。 あまり伝わらないようなことを毎日ぼちぼちやっています。 伝わらない・人を選ぶ内容なので、ついカッとなって今は公開しています。 「ブログ村」というブログ紹介サイトにいつもお世話になっております。 ヨカタ!と思ったら押してネ! ※同姓同名の同一人物「量子きのこ」が、pixivにも生息してます。 よかったら覗いてやってください PR
量子力学の本を返したのでこのチャンスを活かして、調和振動子のほうの昇降演算子の、
空で書けるようになった部分を書きます。 まず、調和振動子ポテンシャルのハミルトニアンを無次元化します。 ハミルトニアンの次元はエネルギーなので、 同じくエネルギーの次元を持つ の2倍でハミルトニアン を割ります。 hにバーがついたのはプランク定数を2πで割ったディラック定数 kはばね定数、mは質量、xは位置、pは運動量、ωは角振動数で の関係があるので、無次元化したハミルトニアンは以下のようになります。 そしてこれをガウス素数っぽく半分に因数分解します。 この片方、どっちか忘れましたがプラスがついたほうをaと置くと マイナスがついたほうはaの複素共役(というかたぶんエルミート共役)a†と書くことができます。 †をどっちのaにつける慣習かは忘れましたが 運動量pにプラスiがついたほうを下降演算子、マイナスiがついたほうを上昇演算子と呼び 波動関数に上昇演算子を左からかけると、量子数1個分だけ固有値がアップした波動関数を 波動関数に下降演算子を左からかけると、量子数1個分だけ固有値がダウンした波動関数を 算出することができます。 ではこの芋づるの端っこはといいますと 一番下の波動関数に下降演算子をかけると恒等的にゼロになることが知られているため、それを利用して微分方程式を解くことで算出できます。 つまり具体的にやってみますと 物理量xとpを微分演算子にして このように解くことができます。(Cは任意の定数です)
前回までは、粘性抵抗と慣性抵抗両方が加わった落下運動の運動方程式の解について
粘性抵抗と慣性抵抗片方ずつに係数を制限して、元々どちらかしか考慮していない運動方程式の解と同じになるかを検証していた。また、それぞれの終端速度についても議論していた。 今回はいよいよ、粘性抵抗と慣性抵抗両方を含んだ状態の終端速度について考えてみる。 微分方程式の解は以下のようになるのだったが このとき、時間t→∞の極限で終端速度vはどうなるだろうか。 cothの中身が大きくなると1に収束するので 分母分子に をかけて これは、元の微分方程式である mv'=mg-c1v-c2v^2の加速度(力)v'をv'=0としてvについて解いた2次方程式の解そのものである。 そのうえ、数回前に としたDがまさに判別式であったこともわかるだろう。 また、すべての変数、関数、パラメータが正の実数と仮定すると 2次方程式の判別式Dのルートは必ず正の実数であり Dのルート√(c1^2+4mgc2)は必ずc1より大きいため、v=-c1-√(c1^2+4mgc2)/(2c2)の解はありえないこともわかる。
前回はc1≠0、c2=0の条件で解いたが、今度は逆にc1=0、c2≠0の条件で解いてみよう。
となるが、元の運動方程式 も解いてみると 部分分数分解を用いて (A1+A2)v=0 からA1=-A2 初期条件t=0でv=0を与えるとG2=-1となって ここで、式を整理してvの式にする 両辺にv+√(mg/c2)をかけると ここで、右辺の分母分子にをかけると これは双曲線関数ハイパボリックタンジェントなので を得、見事に両者は一致する。 また、終端速度については t→∞でv=√(mg/c2)になるが、 これについても元の運動方程式 mv'=mg-c2v^2のv'=0の条件で解くことで、単なる代数方程式として 終端速度v=√(mg/c2)を簡単に得ることができる。 前回、粘性抵抗と慣性抵抗両方を受ける落下運動について運動方程式を解いた。
質量をm、重力加速度をg、速度v、時間t
速度の1乗に比例する空気抵抗の係数をc1、速度の2乗に比例する空気抵抗の係数をc2とすると 運動方程式は以下のように書ける。 速度は高さy方向の下向きをプラスに定義する。 ここで、v^2の係数を1にしたいので、両辺を-c2で割ると ここで、式の簡単化のためG=m/c2、B=c1/c2、A=mg/c2とすると dvの分母に右辺を持ってくると以下のようになる。 両辺を積分してvの式にすると、速度vについての運動が解けるのだが、vの2乗が邪魔でうまく積分できないので、部分分数分解を用いる。 という恒等式を作り、A1、A2、Dがいくつになるのか定めることにする。 まず、 になるようにDを定めたい。 なので、 であることがわかる。なぜいきなりDという文字を仮定したのかはあとでわかるが、これは2次方程式の判別式Dの意味である。 次にA1とA2を定める。 という恒等式からA1とA2を定めたいので v(A1+A2)=0 からA1=-A2が導かれるので B/2(A1+A2)もゼロである。 D(A1-A2)=1にA1=-A2を大入して、 を得る。 そうすると元の微分方程式は なので積分ができ となる。G2は任意の積分定数。 対数同士の引き算は対数の中の割り算として引っ込めることができるため 両辺に2Dを掛け算して (ここでG2はまだ任意の定数なので、2DG2を新たなG2として上書きする) 両辺の指数をとって (ここでもexp(G2)を新たなG2として上書きしている) 左辺の分母分子に2を掛け算してやると 初期条件 t=0でv=0を与えて任意の定数G2を定める。 微分方程式の解は となる。 これをvの式に直す。 両辺に(B+2D)と(2v+B+2D)を掛け算する。 展開すると 全部左辺に移行し、vを含むか含まないかで整理すると ここで1/(2v)の式に着目してみることにする。 第1項は約分できるので 問題は第2項だが の分母分子にexp(x)を掛け算すると となって、これは双曲線関数ハイパボリックタンジェントtanhの逆数ハイパボリックコタンジェントcothであることがわかる。 ここでようやく2vの式に戻すと シークバーで見れるバージョンはpixivにあります。 べき乗関数の宝庫である熱力学の状態図、 両対数方眼紙でべき乗関数は直線になりますが、 変数が3つ以上あるので、とりあえずエントロピーを除いて、3D対数ブロックにしたらどうなるか試してみました。 下にDL用Excelファイル(マクロなし)へのリンクを貼りましたので、自粛期間中に暇だったら遊んでみてください。(いつもながら対象年齢高めですみません) シート保護の関係でサンプルgif動画と異なる表示がされています。 DLしたら、「編集を可能にする」にしてからご使用ください。 3DカルノーサイクルExcelファイル(ダウンロード)
有限深さ井戸型ポテンシャルは以前解いてブログに載せたのですが
そのときは力技で解いていました。 もし、偶数番目の量子数で波動関数が奇関数、奇数番目の量子数で偶関数になるとわかっているなら、どれくらい簡略化されるのか これについてまだブログには載せていませんでした。 幅2aの井戸型ポテンシャルの、中のポテンシャルはV=0、外のポテンシャルはV=V0とし、 井戸の中心をx=0とします。(-a≦x≦a) x<0の場合は、偶関数なら波動関数f(-x)=f(x)とし、奇関数ならf(-x)=-f(x)とするわけです。 境界条件はx=±aでf(±a)とその微分f'(±a)が連続であることですが x<0を考慮しなくて済むので、 4本の連立方程式が2本まで簡略化できます。 つまり、4次ではなく2次の行列方程式を解けばいいだけになるので、かなり楽です。 波数をk、積分定数をA,Bとして(iは虚数単位) f(x)=Aexp(ikx)+Bexp(-ikx) という波動関数の境界条件を解けばいいのですが、 左右対称か左右反対称どちらかの解だけになるので、波動関数fは三角関数 井戸の中はf(x)=Asinkxかf(x)=Bcoskxのどちらかだけになりますし、 井戸の外はf(x)=Cexp(-k1x)になります。 量子数が奇数の場合、波動関数は偶関数なので f(a)=Bcoska=Cexp(-k1a) f'(a)=-kBsinka=-k1Cexp(-k1a) 量子数が偶数の場合、波動関数は奇関数なので f(a)=Asinka=Cexp(-k1a) f'(a)=kAcoska=-k1Cexp(-k1a) これを解けばいいだけになります。 太線から下は、無限深さ井戸型ポテンシャルへの近似で、検算です。 k1=∞になるので n番目の量子数がそのまま量子数nとなります。 量子数nが偶数だと波動関数が奇関数になり 量子数nが奇数だと波動関数が偶関数になることが確認できました
次のページ>>
|
カレンダー
カテゴリー
最新CM
[12/30 buy steroids credit card]
[09/26 Rositawok]
[03/24 hydraTep]
[03/18 Thomaniveigo]
[03/17 Robertaverm]
最新記事
(01/01)
(02/11)
(05/30)
(05/28)
(05/28)
(05/27)
(08/04)
(10/24)
(06/08)
(05/22)
(01/13)
(11/04)
(11/02)
(10/28)
(10/27)
最新TB
プロフィール
HN:
量子きのこ
年齢:
43
HP:
性別:
男性
誕生日:
1981/04/04
職業:
WinDOS.N臣T
趣味:
妄想・計算・測定・アニメ
自己紹介:
日記タイトルの頭についてるアルファベットは日記の番号です
26進数を右から読みます 例:H→7番目、XP→15(P)×26+23(X)=413番目。 A=0とする仕様につき一番右の桁はAにできませんのでご了承くださいズコー
ブログ内検索
アーカイブ
最古記事
(05/11)
(05/11)
(05/13)
(05/13)
(05/13)
(05/13)
(05/13)
(05/13)
(05/14)
(05/14)
(05/14)
(05/14)
(05/16)
(05/16)
(05/16)
アクセス解析
|