owner古典的な波動方程式の解き方ポストシェアはてブ送るPocket古典的波動方程式は、量子力学でおなじみのシュレディンガー方程式を理解するうえでの基礎になります。今回は、一次元の波動方程式を解いていきます。分離定数K<0のときの変位u(x,t)振動する弦の波動方程式まず、両端が固定された弦を振動させたときの弦の振る舞いを考えます。上図のように、左端を原点として、時刻\(t\)における位置\(x\)の変位を\(u(x,t)\)とします。端点間の長さは\(l\)とします。このとき、\(u(x,t)\)は以下の方程式を満足することが知られています。これを古典的波動方程式といいます。\tag{1}$$このとき、\(v\)は弦を伝わる波の速度になります。今回は、この方程式を既知のものとして解き進めていきます。上図の弦は両端が固定端になっているため、時刻\(t\)によらず変位は\(0\)になります。このように、\(u(x,t)\)が満たすべき条件のことを境界条件と呼びます。式に表すと以下のようになります。\tag{2}$$波動方程式は変数分離法を使って解くことができる波動方程式は、偏微分方程式の解法の1つである変数分離法によって解くことができます。変数分離法では、2変数関数である\(u(x,t)\)が2つの1変数関数\(X(x)\)と\(T(t)\)の積であると仮定します。\tag{3}$$変数分離法で仮定した(3)式を波動方程式(1)に代入します。$$このとき、方程式の微分の部分が1変数関数の微分、すなわち常微分に変わったので微分記号が\(\partial{}\)から\(d\)になっています。\tag{4}$$上の(4)式では、左辺の変数が\(x\)のみで、右辺の変数は\(t\)のみになっています。xとtは互いに独立な変数ですから、(4)式が任意の\(x\)および\(t\)で成り立つためには(左辺)=(右辺)=(定数)が必要です。\tag{6}$$となります。このとき、定数としておいた\(K\)を分離定数と呼び、実際の値は波動方程式を解く過程で求まります。\tag{8}$$上の(7)(8)式はいずれも1変数関数の微分方程式であるため、常微分方程式になります。変数とその微分が1乗の形でしか現れていないことと、方程式の各項の係数が定数であることから、(7)(8)式は定係数の線形微分方程式でもあります。このような方程式は簡単に解くことができます。具体的には、\(K\)の値が0か正か負かで場合分けします。分離定数K=0のときまずは\(K=0\)の場合を仮定します。$$ただし、\(a\)と\(b\)は積分定数です。確かに、上の2式を2階微分すれいずれも\(0\)になります。$$$$u(l,t)=X(l)T(t)=(a_1l+b_1)(a_2t+b_2)=0$$ここで、弦の変位\(u(x,t)\)は時間変化するため、つねに\(T(t)=0\)となることはありません。上2式がtの値に関わらず成り立つためには\(X(0)=X(l)=0\)、すなわち\(a_1=b_1=0\)であることが必要です。\tag{9}$$になります。これは、弦の変位が時間と位置に関係なく\(0\)になってしまうため物理的に意味のない解です。このような解を無意味な解と呼びます。分離定数K>0のとき次に、分離定数\(K>0\)のときです。(7)(8)式を以下に示します。\tag{9}$$の形で表される方程式であることがわかります。このとき、\(k\)は実数定数であるとします。形をとることがわかっています。\(α\)は定数で、方程式を解く過程で決まります。\mathrm{e}^{αx}=0$$\(\mathrm{e}^{αx}=0\)のときは\(y(x)=0\)で無意味な解となるため、\(α^{\:2}-k^2=0\)です。\tag{10}$$になります。\(c_1\)と\(c_2\)はともに定数です。定係数の線形微分方程式の解き方について詳しく知りたい方は常微分方程式の教科書を参照してください。直感的な理解としては、得られた2つの解の線形和もまた方程式の解であり、2階微分を含む方程式を解くためには積分を2回実行する必要があることから、(10)式のように定数が2つ現れることになります。あとは(7)(8)式に応じて\(k\)の値を決定すれば良いですし、境界条件などを利用すれば定数cの値も求まります。分離定数K<0のとき最後に、分離定数\(K<0\)のときです。解き方の流れは\(K>0\)のときと同じです。再度、(7)(8)式を以下に示します。\tag{11}$$の形で表すことができます。\(k\)は実数定数です。\mathrm{e}^{αx}=0$$\(\mathrm{e}^{αx}=0\)のときは\(y(x)=0\)で無意味な解となるため、\(α^{\:2}+k^2=0\)です。\tag{12}$$このとき、(\(c_1\)と\(c_2\)はともに定数です。\mathrm{i}k\)に変わっただけなのがわかります。\(\mathrm{e}\)を底とする指数関数の指数部分が虚数になった際には、以下に示すオイラーの式と呼ばれる式を用いて式を書き換えることができます。(\mathrm{i}c_1-\mathrm{i}c_2)\sin{kx}\end{align}$$になりました。このとき、\(c_1+c_2\)および\(\mathrm{i}c_1-\mathrm{i}c_2\)は定数であるので、それぞれ\(c_3\)と\(c_4\)に置き換えて見やすくします。\tag{13}$$以上、\(K\)の範囲で場合分けして方程式を解くことによって、\(u(x,t)=X(x)T(t)\)の\(X(t)\)と\(T(t)\)がどのような式になるのかを文字定数を用いて表すことができました。次の項からは、\(X(x)\)と\(T(t)\)をそれぞれ具体的に求めて\(u(x,t)\)を決定します。波動の変位u(x,t)を求めるここまでは、古典的な波動方程式を、変数分法を用いながら分離定数\(K\)の範囲で場合分けをすることで解き進めてきました。ここからは、求められた関数\(u(x,t)\)の文字定数の部分を決定していきます。以下は、(7)(8)式の再掲です。\tag{8}$$\(K=0\)のときは無意味な解になるため考えません。分離定数K>0のときの変位u(x,t)\(K>0\)のとき、\(K=β^{\:2}\)を満たす実数\(β\)をおきます。先の議論で、\(K>0\)のときの\(X(x)\)と\(T(t)\)が(10)式の形で表されることがわかりました。l}=0$$この連立方程式を\(c_1\)と\(c_2\)について解くと\(c_1=c_2=0\)になりますが、\(X(x)=0\)になっていしまいます。よって、\(u(x,t)=0\)となり、分離定数\(K>0\)のときも方程式の解は無意味です。残るは\(K<0\)のときです。分離定数K<0のときの変位u(x,t)最後に、分離定数\(K<0\)のときです。流れは\(K>0\)のときと同じです。先の議論で、\(K<0\)のときの\(X(x)\)と\(T(t)\)が(13)式の形で表されることがわかりました。$$\(c_4=0\)のとき、\(X(x)=0\)より\(u(x,t)=0\)ですから、このときの解は無意味なものになります。\tag{16}$$(8)式も解きます。\(T(t)\)の解も\(X(t)\)と同じく(13)式の形で表されますが、定数\(k\)の重複を避けるために\(k_1\)とします。v^2\)の関係を得ます。$$このとき、係数\(c\)についても重複を避けるため、\(X(t)\)の係数とは別の文字\(c_5,c_6\)を用いました。$$ここで、簡素化のために\(c_4c_5=B\)、\(c_4c_6=C\)とおきました。\tag{18}$$と表し直します。このことから、変位\(u\)は複数の解をもつことになるわけですが、そもそもの始まりである波動方程式(1)が線形微分方程式であるため、一般解は(18)式の無限間数列の和になるわけです。\tag{19}$$これが古典的な波動方程式の解です。pp.43-52物理化学ポストシェアはてブ送るPocketfeedly時間に依存しないシュレディンガー方程式を導出しようハロアルカンの性質と求核置換反応RECOMMENDこちらの記事も人気です。