波動方程式
甲南大学理工学部機能分子化学科 山本雅博

This manuscript is modified on March 26, 2012 3 : 53 pm

この文章は基本的に有山の教科書 [1] を参考にして作成した。記号やわかりにくいところは一部改変してある。

1

波動とは波は、媒質のある物理量が周期的に変化しがそれが伝播（でんぱ）いくと考えられる。音波なら媒質は空気の密度が波の進行方向で粗密になって伝播していくし（縦波）、電磁波なら電場と磁場が波の進行方向と直角に変位して伝播していく (横波)。このように波動は，空間と時間で指定される関数によって表される。

Figure 1: 縦波 (longitudinal wave) と横波 (transverse wave). 波の進行方向に対して，変位が平行であるのが縦波で，変位が進行方向と垂直であるのが横波である。地震波では、縦波（P 波）は毎秒約 7km、横網（S 波・表面波）は毎秒約4km の速さで伝わる。この伝播速度差を利用して、震源に近い地点における P 波の観測に基づき、後から来る S 波の伝播を時系列的に予測し、震源からある程度以上（P 波と S 波の時間差が充分に開くほど）離れた地点に対しては、その到達前に予測を発表することができる。（緊急地震速報）
今、波の進行方向を x 方向とし、ある時間 t 場所 x での変位を u(x, t) とする。時間 t = t0 = 0 での変位を f (x) とすれば、f (x) = u(x, 0) である。今、この変位は一定の形を保ったまま一定の速度 v で +x 方向に伝播していくとする。
(Fig.2) (波が弱まったり形が変化することもあるが短い時間の範囲では近似的に成立するとする。)

u(x, δt) = f (x − vδt)

(1)u(x, −δt) = f (x + vδt)(2)

δt を t で置き換えると
u(x, t) = f (x − vt),

propagate + x direction−→(3)←−(4)となる。変位が −x 方向に進む場合は、v を −v に変えればよくu(x, t) = f (x + vt),

propagate − x directionとなる。(Fig.3)
Fig.2,3 では，f として波束（波がある部分に集中している）の例を挙げたが，どのような関数形でもよい。詳細はフーリエ級数やフーリエ変換の解説等を参考にしてほしいが，波束やその他の一般の関数は，周期的な波の重ね合わせで表現できる。例えば，Fig.4 にガウス波束の例を示した。ある波長を中心としてガウス分布の大きさで正弦波を重ね合わせるとガウス波束が再現できる。
従って，一般的な波の性質を調べるには，周期的な波についてその性質を調べておけば，一般の関数についても解析が可能になる。周期的な波は，一般にAei(kx−ωt) = A cos(kx − ωt) + iA sin(kx − ωt)1(5)と書ける。この Euler の式の証明は Appendix にある。cos, sin は周期 2π の関数である。
t = 0 とき，kx = 0, 2π で一周期となる。x = 0 から波長 (wavelength)λ で一周期となるので，kλ = 2π, k = 2π/λとなる。k は，波数 (wavenumber) と呼ばれる。
また，x = 0 のとき，ωt = 0, 2π で一周期となる。t = 0 から周期 (period)T で１周期となるので，T = 2π/ω となる。ω = 2πν で振動数 (frequency)ν を定義すれば，ν = 1/T となる。ω は角振動数 (angular frequency) とよばれる。
exp(iθ) の θ に相当する kx − ωt を位相 phase とよぶ。cos の波において一番振幅が正方向に大きいのは位相がゼロkx − ωt = 0 のところである。その位置は，x = (ω/k)t となり，+x 方向に進行する波となる。また，この時の波の移動速度を位相速度 phase velocity といい，v = ω/k となる。−x 方向に進行する波はしたがって，Aei(kx+ωt) となる。
(Fig. 5, 6, 7)

t0 + δtvδtt = t0vδtt0 − δtxFigure 2: 波束の伝播。波束は +x 方向に速度 v で形を変えずに移動している。

t0 + δtvδtt = t0vδtt0 − δtxFigure 3: 波束の伝播。波束は −x 方向に速度 v で形を変えずに移動している。

2波動方程式の導入u を，x および t で微分する。また，この時 X を以下の様に定義する。f は f (X) となる。
X∂u
∂x
∂2u
∂x2∂u
∂t
∂2u
∂t2
∂2u
∂t2

≡ x − vt, f (x − vt) = f (X)
∂X df
df
==
∂x dXdX∂ df∂X ∂ df
d2 f
===
∂x dX
∂x ∂X dX
dX 2
∂X df
df
== −v
∂t dX
dX
∂ df
∂X ∂ df
d2 f
= −v
= −v
= v2
∂t dX
∂t ∂X dX
dX 22∂
u= v2 2
∂x2

(6)
(7)
(8)
(9)
(10)
(11)

)./01$/2/-%-&-&34561*/3))7)') """""))7)') """"))7)') """!""!#"!"$"!#&'#()*' ()+*,' -&- """"$%"

$ """%"

xFigure 4: 波束はいろいろな波の重ね合わせである。f (x) = A exp(−x2 /2/∆2 ) exp(ikx) で表される Gaussian∫∞2A∆dk 0 exp[− ∆2 (k − k 0 )2 + ik 0 x] と書ける。すなわち実数部だけ書く
wavepacket（ガウス波束）は，f (x) = √2π −∞][∫ +∞ 0A∆∆20 2と，Ref (x) = A exp(−x2 /2/∆2 ) cos(kx) = √(k−k)
cos(k 0 x) となる。積分ではなくと和のdkexp−2−∞[ 2π2]∑∆A∆2形に書くと，Ref (x) = √cos[(k − k”)x], k” = k − k 0 = ndk”,(n = −N, −N +k” (δk”) exp − 2 k”
2π1, ...., −2,[ −1, 0, 1,]2, ...., N − 1, N ) となる。すなわち，ガウス波束は，いろいろな波長の波 cos[(k − k”)x] を重み2√∆ exp − ∆ k”2 をかけて和をとったもので表される。
22π−2π−πkx

π2πFigure 5: 波の伝播。赤線は cos(kx) である。 δ > 0 とすると，青線は cos(kx − δ), 緑線は cos(kx + δ) である。２つの山の間隔から kλ = 2π がわかる。ある時間後に +x 方向に伝播する場合は cos(kx − δ),−x 方向に伝播する場合は，cos(kx + δ)となる。

3Figure 6: 波の位相ei(kx+ωt)
Figure 7: +x 方向への進行波の頂点に乗ってサーフィンする。頂点は，波の位相がゼロである。この条件から位相速度v = ω/k が得られる。波の最高到達地点の間隔は，波長 λ である。

cos(ωt)ωt = 0, ±2π, ...
T = 2π/ωT = 2π/ω = 1/νv = ω/k = 2πν/(2π/λ) = λνc = ω/k = λνFigure 8: ある場所における波の振動，周期，振動数，角振動数4が得られる。f (x + vt) あるいは f1 (x − vt) + f2 (x + vt) についても同様な式が得られる。(宿題：f (x + vt) あるいはf1 (x − vt) + f2 (x + vt) についても同様な式が得られることを証明せよ。) 式 (11) は，波動方程式 (wave equation) とよばれ，時間および空間に関して 2 階の偏微分方程式 となっている。また、それぞれの解の足し算もまた波動方程式の解になっていることは、線形性あるいは「重ね合わせの原理」が成り立つことを意味する。
逆に，式 (11) の方程式をみたす関数は，f (x − vt), f (x + vt), f1 (x − vt) + f2 (x + vt) となることも示すことができる。以下のように定義しよう。

X≡ x − vt,

Y ≡ x + vt

(12)

∂ 2 u/∂t2 と ∂ 2 u/∂x2 を u の X, Y についての微分で表してみる。
x =
∂u∂t
∂ ∂u
∂t ∂t

===

∂u
∂x
∂ ∂u
∂x ∂x===X +Y−X + Y, t=22v∂X ∂u∂Y ∂u∂u∂u+= −v+v
∂t ∂X
∂t ∂Y
∂X∂Y∂ ∂u∂ ∂u∂X ∂∂Y ∂ ∂u∂X ∂∂Y ∂ ∂u−v+v= −v(+)+ v(+)
∂t ∂X
∂t ∂Y
∂t ∂X
∂t ∂Y ∂X∂t ∂X∂t ∂Y ∂Y
22
∂2u
2 ∂ u
2 ∂ uv2−2v+v∂X 2∂X∂Y
∂Y 2
∂X ∂u
∂Y ∂u
∂u
∂u+=+∂x ∂X
∂x ∂Y
∂X
∂Y
∂ ∂u
∂ ∂u
∂X ∂∂Y ∂ ∂u∂X ∂∂Y ∂ ∂u+=(+)+(+)
∂x ∂X
∂x ∂Y
∂x ∂X
∂x ∂Y ∂X∂x ∂X
∂x ∂Y ∂Y∂2u∂2u∂2u+2+2∂X∂X∂Y∂Y 2(13)(14)(15)
(16)
(17)
(18)
(19)

波動方程式が成立するには，∂2u∂X∂Y= 0(20)となる。
u = f1 (X) + f2 (Y ) ならこの式は成立する。従って，u(x, t) = f1 (x − vt) + f2 (x + vt)
この式 Y について積分すると，∂u∂X= f10 (X)f10 (X) は、Y に依存しない任意の X の関数の微分である。この式をさらに X で積分すると、∫u =dXf10 (X) + f2 (Y ) = f1 (X) + f2 (Y )

(21)

(22)

ここで、f2 (Y ) は、X に依存しない Y の任意の関数である。以上の関係により、波動方程式を満たす関数 u(x, t) は、u(x, t) = f1 (x − vt) + f2 (x + vt) の関数によりあらわされることが示された。
以上では、波の一般的な性質から波動方程式を導いてきたが、実際には u がもつ物理的な性質から導かれるものである。ロープの張力やばねで連結された質点の例がよく使われる。

3３次元の波3.1

平面波 (plane wave)

波の進む方向に向きをもち大きさ 2π/λ のベクトルを波数ベクトル k としよう。１次元の波 cos(kx − ωt) で波の最大振幅は位相がゼロである線が ±x 方向に進行するとみなせたように、３次元では位相がゼロである面が ±k 方向に進行すると見なせる。この波を平面波という。波数ベクトルに垂直な面内のベクトル r を表す方程式は、

k·r=0(23)kx x + ky y + kz z = 0(24)である。したがって、３次元の平面波はu(r, t)= A exp[i(k · r ∓ ωt)]とあらわされる。

5(25)Figure 9: 平面波

位相 X と波束 u = f (X) より波動方程式を導こう。

X∂u∂t
∂ ∂u
∂t ∂t
∂u
∂x
∂ ∂u
∂x ∂x
∂2u
∂y 2
∇2 u
kx2 + ky2 + kz22

∂ f
∂X 2
∂2u
∂t2

= kx x + ky y + kz z ∓ ωt
∂X ∂f
∂f
== ∓ω
∂t ∂X
∂X
∂X ∂ ∂f∂2f= ∓ω= ω2∂t ∂X ∂X
∂X 2
∂X ∂f
∂f
== kx
∂x ∂X
∂X
∂X ∂ ∂f∂2f= kx= kx2∂x ∂X ∂X
∂X 2222∂ f∂ u
2 ∂= ky2,=kz∂X 2∂z 2∂X 22∂ f= k2, ∇2 = (∂ 2 /∂x2 ) + (∂ 2 /∂y 2 ) + (∂ 2 /∂z 2 )∂X 2= k2

(26)
(27)
(28)
(29)
(30)
(31)
(32)
(33)

21 ∂ u1= 2 ∇2 uω 2 ∂t2k( ω )2=∇2 = v 2 ∇2 u
k=(34)(35)

この式が 3 次元での波動方程式となる。ここで、v = ω/k である。

3.2球面波

3 次元の波動方程式は，平面波だけではなく，点を波源としてあらゆる方向に伝播する球面波（例：X 線の原子により散乱）も解としてもつ。球座標系で，Laplacian は∇2q u=1 ∂2ucos θ ∂u12 ∂u ∂ 2 u∂2u+ 2 + 2 2 + 2+ 2 2r ∂r∂rr ∂θr sin θ ∂θr sin θ ∂φ2(36)u に θ と φ の依存性がないとし，ru(r, t) = U (r, t) とおけば，2 ∂u ∂ 2 u+ 2r ∂r
∂r2∂U
∂2 U
1 ∂2U= v2 (+ 2 ) = v2,r ∂r r
∂r rr ∂r2

∇2q u =
1 ∂2Ur ∂t2

6(37)2∂2U2∂ U=v∂t2
∂r2(38)

となり，１次元の波動方程式と同じ形になる。１次元の時と同じ考え方で，U (r, t) = f1 (r − vt) + f2 (r + vt) とすることができるので，u(r, t) =U (r, t)
f1 (r − vt) + f2 (r + vt)=rr(39)と書くことができる。図に，球面波 (spherical wave) の典型例を示す。
OUTIN

Figure 10: 球面波。波は IN 方向に球の中心に向かって進行し f2 (r + vt)，OUT 方向に向かって球の中心から外側に進行する f1 (r − vt)。

4Appendixeiθ

= exp(iθ) = cos θ + i sin θ(40)はそれぞれを Taylor 展開すればよい。
If the function f (x) can be approximated by the polynominals around x = af (x)' a0 + a1 (x − a) + a2 (x − a)2 + a3 (x − a)3 + a4 (x − a)4 + ... + an (x − a)n(41)

If |x − a| << 1, the above equations can be converged. We can get the 1st to n-th derivative such as
f 0 (x) =f 00 (x) =a1 + 2a2 (x − a) + 3a3 (x − a)2 + 4a4 (x − a)3 + ... + nan (x − a)n−1
2a2 + 3 × 2a3 (x − a) + 4 × 3a4 (x − a)2 + ... + n(n − 1)an (x − a)n−2(42)(43)

f 000 (x) = 3 × 2a3 + 4 × 3 × 2a4 (x − a) + ... + n(n − 1)(n − 1)an (x − a)n−3f 0000 (x) = 4 × 3 × 2a4 + ... + n(n − 1)(n − 2)(n − 3)an (x − a)n−4
f (n) (x) = n!an

(44)
(45)
(46)

1 000f (a),
3!

(47)

Then we have
1 (n)
f (a)
n!
1
1
1f (x) ' f (a) + f 0 (a)(x − a) + f 00 (a)(x − a)2 + f 000 (a)(x − a)3 + ... + f (n) (a)(x − a)n2!
3!
n!
a1

= f 0 (a),

a2 =

1 00
f (a),
2!

a3 =This is the Taylor expansion around x = a.

7

a4 =

1 0000f (a),
4!

an =(48)

exeixcos(x)1 2
11x + x3 ... + xn + ...
2!
3!
n!
111
1
= 1 + ix − x2 − ix3 + x4 + ix5 ...
2!
3!
4!
5!
11
1
1
= 1 − x2 + x4 ..., sin(x) = x − x3 + x5 ...,
2!
4!
3!
5!
= 1+x+(49)(50)
then we have eix = cos x + i sin x(51)(52)

上の eix の展開では，x << 1 のところで展開した。より一般に x = x0 の回りでの展開ではどうなるのであろうか？
x0 = x1 + dx, dx << 1 とする。

eix0= ei(x1 +dx) = eix1 eidx = eix1 [1 + idx −1(dx)2 − ...] = eix1 [cos(dx) + i sin(dx)]2!となる。同様に，x1 = x2 + dx, x2 = x3 + dx, ...., xn−1 = xn + dx, 0 = xn = x0 − ndx

eix1= eix2 [cos(dx) + i sin(dx)],

eix0 = eix2 [cos(dx) + i sin(dx)]2

eixn−1...
= eixn [cos(dx) + i sin(dx)],

eix0 = eixn [cos(dx) + i sin(dx)]n = [cos(dx) + i sin(dx)]n[cos(dx) + i sin(dx)]2[cos(dx) + i sin(dx)]3[cos(dx) + i sin(dx)]n
eix0

=
=
=

cos2 (dx) − sin2 (dx) + 2i cos(dx) sin(dx) = cos(2dx) + i sin(2dx)[cos(2dx) + i sin(2dx)][cos(dx) + i sin(dx)]cos(2dx) cos(dx) − sin(2dx) sin(dx) + i[sin(2dx) cos(dx) + sin(dx) cos(2dx)]=

cos(3dx) + i sin(3dx)
...

=
=

cos(ndx) + i sin(ndx) = cos(x0 ) + i sin(x0 )
cos(x0 ) + i sin(x0 )(53)(54)となり，一般の x に対して，eix = cos x + i sin x が示された。

複素平面

虚数軸i = √-1z = r eiθrθ

実数軸→：長さrで実軸となす角θz = r e = r (cosθ ＋ i sinθ）iθ

z = Re(z) + i Im(z),Re(z) = r cosθ,Im(z) = r sinθFigure 11: 複素平面における複素数の表示。横軸が実数軸で縦軸が虚数軸である。複素数 z を z = reiθ (ここで，r は原点からの距離，θ は実軸からの角度) で表すと, 実数成分は Rez = r cos θ，虚数成分は Imz = r sin θ となる。

8ei(kx−ωt)t=0t = +∆tx = −∆xx=0x = +∆x−2π−πkxπ2πFigure 12: +x 方向への進行波 ei(kx−ωt) . x = 0(中央の図) および t = 0 で, e0 = 1 すなわち赤の矢印の位置にあったものが，t = +∆t 後には，e−iω∆t = cos(ω∆t) − i sin(ω∆t) すなわち青の矢印の位置にある。この波は，複素平面上を半径1 で時計回りに角速度 ω で回転しているものと考えることができる。矢印の実軸上への射影を右図の赤線と青線で示した。青線は，赤線から見て +x 方向への進行波となっている。

ei(kx+ωt)t = +∆tt=0x = −∆xx=0x = +∆x−2π−πkxπ2πFigure 13: −x 方向への進行波 ei(kx+ωt) x = 0(中央の図) および t = 0 で, e0 = 1 すなわち赤の矢印の位置にあったもの
が，t = +∆t 後には，eiω∆t = cos(ω∆t) + i sin(ω∆t) すなわち緑の矢印の位置にある。この波は，複素平面上を半径 1で反時計回りに角速度 ω で回転しているものと考えることができる。矢印の実軸上への射影を右図の赤線と緑線で示した。緑線は，赤線から見て −x 方向への進行波となっている。

9References[1] 有山 正孝, 「振動・波動」 (基礎物理学選書 (8)) 1986,10

裳華房