1偏微分方程式を解く
崎間＠物理のかぎプロジェクト

2003-02-09

量子力学の主役「シュレディンガー方程式」をはじめ，物理では偏微分方程式を使うことがよくあります．この方程式の解き方をなんとかマスターしようということで，少々大変ですが偏微分方程式を解いてみましょう．

1 次元波動方程式の例題
いまから解いてみるのは 1 次元波動方程式
2∂ 2 u(x, t)
2 ∂ u(x, t)=c,∂t2∂x2(0 < x < 1)(1)です．ここで c は波の速さで，正の定数です．偏微分方程式を完全に解くには境界条件と初期条件が必要ですから，境界条件 ：u(0, t) = u(1, t) = 0(2)∂u(x,0)
初期条件 ：= 0,∂tu(x, 0) = x(1 − x)(3)を課しておきます．波動方程式という名前は凄そうですが，ただ単に波の運動をあらわす式です．u(x, t)が波の関数で，これに境界条件と初期条件を付け加えることでこの波がどんな運動をしているのかを知ることができます．つまり「解く」というのは 式 (1) を u(x, t) = の形にするという事です．
ところで，波の運動というのがイマイチわかりません． u(x, t) を波の高さだと考えるとどうでしょう
か．u(x, t) = という形になっていたら， x と t を代入することで u(x, t) の値が決まります．つまり波がいつ（時間 t ），どこで（位置 x ）どんな高さなのか（波の高さ u ）が分かるようになりますね．
そうすると，たとえば初期状態から x = 3.5 メートルの場所で 2.4 秒後の波の高さは -1.3 メートルで，x = −10.2 メートルの場所で 5.9 秒後の波の高さは 2.1 メートルだよ，とすぐに計算できます．これらを全部ひっくるめてあるのが波動方程式で，様々な条件（境界条件，初期条件）での波を導き出す操作が波2偏微分方程式を解く

動方程式を解く，ということになります．

まずは変数分離x と t の 2 つの変数がある偏微分方程式では難しいので，変数を分離して 2 つの常微分方程式に分けます．変数を分離するには， u(x, t) の解としてu(x, t) = f (x)g(t)を代入するのが定石です．実際に波動方程式に代入すると
2∂ 2 f (x)g(t)2 ∂ f (x)g(t)=c∂t2∂x2となります．左辺は t で偏微分してあるので x の関数は関係ないですね．だから x の関数である f (x) を偏微分から出しておきます．

f (x)

2
∂ 2 g(t)2 ∂ f (x)g(t)
=c
∂t2
∂x2

同様に右辺は g(t) を偏微分から出しておきます．

f (x)∂ 2 g(t)∂ 2 f (x)2=cg(t)∂t2∂x2

左辺は t の関数，右辺は x の関数にまとめたいので，両辺を f (x)g(t) で割ります．ついでに c2 も移動
すると

1 ∂ 2 g(t)1 ∂ 2 f (x)=c2 g(t) ∂t2f (x) ∂x2となります．さて，この式は何を意味しるんでしょうか？ 左側は t の関数，右側は x の関数です．両辺がイコールで結ばれているので「 t だけ変化しても， x だけ変化しても両辺は等しくなる」ということをいっています．
そんなことが起きるのは両辺が定数のときだけです．定数なら x だけ変わろうが t だけ変わろうが， xと t が同時に変わろうが関係ありません．そこで，この両辺を1 ∂ 2 g(t)1 ∂ 2 f (x)==λc2 g(t) ∂t2f (x) ∂x2というふうに定数 λ と等しいとします（この定数は変数分離定数と呼ばれます）．するとこの式は，つぎのように簡単に 2 つに分けることができます．

1 ∂ 2 g(t)=λc2 g(t) ∂t2(4)1 ∂ 2 f (x)=λf (x) ∂x2(5)— 物理のかぎしっぽ http://www12.plala.or.jp/ksp/ —3

偏微分方程式を解く

この 2 式はそれぞれ変数が一つなのでもはや偏微分方程式でなく，常微分方程式です．変数分離定数 λ を導入することによって，偏微分方程式を二つの常微分方程式に分けることができたのですね．

分離した常微分方程式を解く
二つの常微分方程式に分けることができたので，今度はそれらを一つずつ解いて行きます．まずは x に関する微分方程 式 (5) を解きましょう． 式 (5) は少し変形すると

d2 f (x)= λf (x)dx2です．これに解として f (x) = eax を代入すると

d2 eax= λeaxdx2(6)になります．そして

d2 eaxdd deax=(aeax ) = a2 eax=2dx
dx dxdxより 式 (6) は∴ (a2 − λ)eax = 0

a2 eax = λeaxとなります．ここで eax = 0 とすればこの方程式を満たしますが，この場合は考えません．だってそうす
ると f (x) = 0 となって意味がないでしょう．ということは，

a2 − λ = 0(7)となりますね．この式から a を決定します．
ここで注意が必要です． λ はまだどんな値かさっぱり分からないので，λ が正，ゼロ，負の三つの場合で場合分けして考えないといけません．

λ が正の場合
式 (7) よりa2 = λですね．いまは λ が正の場合を考えているので，そのままルートを付けてかまいません．すると√a=± λです． eax を解として代入したのですから，解は√e λx ,√e− λxの二つです．一般解はこれらの線形結合で書けて√√

f (x) = Ae λx + Be− λx— 物理のかぎしっぽ http://www12.plala.or.jp/ksp/ —4

偏微分方程式を解くとなります（ A, B は任意定数）．ここで境界条件 u(0, t) = u(1, t) = 0 を使います．いまは u(x, t) =f (x)g(t) としてるから f (0) = f (1) = 0 です．したがって
f (0) = A + B = 0√√

f (1) = Ae λ + Be− λ = 0です．これから定数 A, B を求めてみます． e のなんとか乗は常に正ですから，これらを満たすのはA = B = 0 のときのみです．これらを f (x) に代入するとやっぱり f (x) = 0 となってしまいます．これを「恒等的にゼロである」といいます．せっかくさっき f (x) = 0 になる解を避けたのに，これでは意味がありません．
したがって λ が正のときは恒等的にゼロでない解，つまり非自明解は存在しません．他の場合を調べましょう．

λ がゼロの場合
λ = 0 なら 式 (6) はd2 f (x)=0dx2となります．つまり x で 2 階微分してゼロになるものが f (x) だといっています．2 階微分でゼロになる式はたくさんありますが，一般解は

f (x) = Ax + Bと書けますよね（実際に x で 2 階微分してみてください．ゼロになるでしょう？）．境界条件 f (0) =f (1) = 0 を適用すると
f (0) = B = 0
f (1) = A + B = 0です．これらを満たすのは A = B = 0 のときなので，またしても非自明解を持ちません．他の場合を調べます．といっても残りは λ が負の場合だけなので，今度こそ解があると期待できます．

λ が負の場合
a2 − λ = 0 より
a2 = λです．ここで λ は負なのでそのままルートを付けてはいけません！ a2 は λ に等しいから負の値なので，2 乗して負になるように虚数単位 i をルートの前につけておきます． λ はルートの中にあるので絶対値をつけます．したがって√
a = ±i |λ|です． eax を解として代入したんですから，解は√ei |λ| x ,√e−i |λ| x— 物理のかぎしっぽ http://www12.plala.or.jp/ksp/ —5偏微分方程式を解くの 2 つです．一般解はこれらの線形結合で，√√
f (x) = Aei |λ| x + Be−i |λ| xと表されます． e の肩に虚数単位 i が乗っていますね．オイラーの公式をつかいましょう．

√√√
ei |λ|x = cos( |λ| x) + i sin( |λ| x)√√√e−i |λ|x = cos( |λ| x) − i sin( |λ| x)
より}
{}{√√
√
√
f (x) = A cos( |λ| x) + i sin( |λ| x) + B cos( |λ| x) − i sin( |λ| x)√√= (A + B) cos( |λ| x) + i(A − B) sin( |λ| x)となって，ここで A + B = C, i(A − B) = D と新しい定数に置き換えると

f (x) = C cos(√√
|λ| x) + D sin( |λ| x)になります．境界条件 f (0) = 0 を適用すると

f (0) = C cos(0) + D sin(0) = C · 1 + D · 0 = C = 0です．したがって C = 0 ということが分かったので√
f (x) = D sin( |λ| x)と書けます．さらに上式に境界条件 f (1) = 0 を適用すると√√
f (1) = D sin( |λ| · 1) = D sin( |λ|) = 0ですね． D がゼロなら常に f (x) = 0 ですから，解をもつためには D ̸= 0 でないといけません．したがって，√sin( |λ|) = 0です． sin x は周期関数で， x = π, 2π, 3π, . . . でゼロになります．つまり√|λn | = nπ,(n = 1, 2, 3, . . . )がいえます．ここで λ を λn としているのは，それぞれの n に対応した定数 λ がいくつもあることを表しています．

√√|λn | が分かったので f (x) = D sin( |λ| x) はf (x) = D sin(nπx),(n = 1, 2, 3, . . . )となりますが， λ を λn にしたのと同様に，それぞれの n に対応する f (x) と D も多数存在するので，

fn (x) = Dn sin(nπx),(n = 1, 2, 3, . . . )(8)— 物理のかぎしっぽ http://www12.plala.or.jp/ksp/ —6

偏微分方程式を解くと書いておきます．やれやれ，ちょっと疲れましたがやっと境界条件を満たす f (x) の解が求まりました．

もう 1 つの常微分方程式を解く
最初に変数分離したとき，二つの常微分方程式に分かれました．そのうち一つはさっき考えたので，今度はもう一方の常微分方程 式 (4) を考えます． 式 (4) は少し変形すると

d2 g(t)= λc2 g(t)dt2(9)

です． 式 (5) が境界条件を満たす解を先ほど求めたので，それに対応する解を求めます．

√|λn | = nπ,(n = 1, 2, 3, . . . )|λn | = n2 π 2 ,(n = 1, 2, 3, . . . )λn = −n2 π 2 ,(n = 1, 2, 3, . . . )の両辺を 2 乗するとで， λ は負だったのでです．これを 式 (9) に代入すると

d2 g(t)= −(cnπ)2 g(t)dt2となります．解として ebt を代入すると

d2 ebt= −(cnπ)2 ebt
dt2∴ b2 = −(cnπ)2

ここで c, n, π はいずれも正ですから

b = ±icnπです．したがって解（特解）はeicnπt ,e−icnπtの 2 つで，一般解はこれらの線型結合でg(t) = Eeicnπt + F e−icnπtと書けます．例によってオイラーの公式で展開して整理すると

g(t) = (E + F ) cos(cnπt) + i(E − F ) sin(cnπt)
さらに E + F = G, i(E − F ) = H と置くとgn (t) = Gn cos(cnπt) + Hn sin(cnπt),(n = 1, 2, 3, . . . )(10)— 物理のかぎしっぽ http://www12.plala.or.jp/ksp/ —7

偏微分方程式を解くとなります．

2 つの常微分方程式を合わせる
さて，境界条件を考えることで変数分離した常微分方程式のそれぞれの解が分かりました．再掲すると

fn (x) = Dn sin(nπx),(n = 1, 2, 3, . . . )gn (t) = Gn cos(cnπt) + Hn sin(cnπt),(n = 1, 2, 3, . . . )です． f (x) と g(t) はもともと波動方程 式 (1) を変数分離するために u(x, t) = f (x)g(t) とおいたものでした．したがって

un (x, t) = fn (x)gn (t)= Dn sin(nπx) {Gn cos(cnπt) + Hn sin(cnπt)}= sin(nπx) {Dn Gn cos(cnπt) + Dn Hn sin(cnπt)} ,(n = 1, 2, 3, . . . )です．ここで Dn Gn = An , Dn Hn = Bn とおきます（ A と B は最初の方で使いましたが，ここでは新しい定数として使います）．さらに cnπ = ωn とおきます．すると

un (x, t) = sin(nπx) {An cos(ωn t) + Bn sin(ωn t)} ,(n = 1, 2, 3, . . . )(11)となります．これは波動方程 式 (1) の境界条件の 式 (2) を満たす解です．

初期条件を考える
境界条件を満たす解が分かりましたから，あとは初期条件を満たすようにするだけです．un (x, t) はい
ろいろな n に対する解がたくさんありますが，重ね合わせの原理よりそれらすべてを足したもの（線型結合）も解となります．すなわち

u(x, t) ==∞∑n=1∞∑un (x, t)sin(nπx) {An cos(ωn t) + Bn sin(ωn t)}(12)

n=1です．初期条件∂u(x,0)= 0 を当てはめるため，まず u(x, t) を t で偏微分しておきましょう．
∂t∞∂ ∑∂u(x, t)=sin(nπx) {An cos(ωn t) + Bn sin(ωn t)}∂t∂t n=1∞==∂ ∑{An sin(nπx) cos(ωn t) + Bn sin(nπx) sin(ωn t)}∂t n=1∞∑{An sin(nπx){−ωn sin(ωn t)} + Bn sin(nπx)ωn cos(ωn t)}

n=1— 物理のかぎしっぽ http://www12.plala.or.jp/ksp/ —8

偏微分方程式を解く

そして上式に∂u(x,0)= 0 を適用すると∂t∞∂u(x, 0) ∑={An sin(nπx){−ωn sin(ωn · 0)} + Bn sin(nπx)ωn cos(ωn · 0)}∂tn=1==∞∑{An sin(nπx) · 0 + Bn sin(nπx)ωn · 1}

n=1∞
∑Bn ωn sin(nπx) = 0

n=1となります．また，もう一つの初期条件 u(x, 0) = x(1 − x) よりu(x, 0) ===∞
∑n=1∞
∑n=1∞
∑sin(nπx){An cos(ωn · 0) + Bn sin(ωn · 0)}sin(nπx){An · 1 + 0}
An sin(nπx) = x(1 − x)n=1です．まとめると∞∑∂u(x, 0)=Bn ωn sin(nπx) = 0,∂tn=1u(x, 0) =∞∑(0 < x < 1)

An sin(nπx) = x(1 − x),(0 < x < 1)

n=1となります（ x の範囲 ( 0 < x < 1 ) はもともとの波動方程式のものです）．これはフーリエ正弦級数の形をしています．したがって 式 (12) が初期条件を満たすためには∂u(x,0)は 0 の半区間でのフーリエ正∂t弦級数， u(x, 0) は x(1 − x) の半区間でのフーリエ正弦級数でなければならないので，フーリエ級数の公式から

Bn =212
An =1∫ 10 · sin(nπx) dx0∫ 1x(1 − x) sin(nπx) dx0となり，

An =2{cos(nπ) + 1},nπ

Bn = 0が得られます．こうして係数 An と Bn が決まったので，これらを 式 (12) に代入すれば境界条件と初期条件を満たす 1 次元波動方程式の解

u(x, t) =∞∑2n=1nπsin(nπx){cos(nπ) + 1} cos(ωn t)が得られます．これでようやく解けました．お疲れさまです．

— 物理のかぎしっぽ http://www12.plala.or.jp/ksp/ —