振動する弦のモデル化
前回の復習
前回の変数分離法の計算で、
u(x,t)=12{f∗(x−vt)+f∗(x+vt)}
という形にu(x,t)が書けることを見ました。イメージとしては下図のようになるのでしたね。
図4 式(35)の物理的解釈（再掲）
図のように、f∗(x−vt)のグラフはf∗(x)のグラフを右へvtだけ移動させて得られます。従って、f∗(x−vt) (v>0)は時間tの増加と共に右へ進行する波をあらわしています。同様に、f∗(x+vt)は左へ進行する波をあらわしています。つまり、u(x,t)はこれら2つの進行波の重ね合わせです。
これらの内容を基にして、今回はd’ Alembert の解法による解の求積を見ていきます。
波動論と音楽
d’ Alembert の解法による解の求積を見ていく前に、少し波動論と音楽の関係について見てみたいと思います。
これまでの議論から、波動方程式の解と高校物理における波動論の話の間に密接な関係があることは少なからず実感でしょう。波動論は音楽との関連も非常に深いものです。ここでは休憩も兼ねて、音楽への応用について簡単に述べることにしましょう。これまでの議論より、振動数は
fn=λn2π=vn2L=n2L√Tρ
で与えられます。振動数が高ければ高いほど、高い音を出す楽器となるので、張力Tを変えれば調律ができることが分かるでしょう。故に、ギターやヴァイオリン、コントラバスなどの弦楽器は調律の際に張っている弦の張力を調整する仕組みになっているはずです。
しかし、張力は無限に大きくすることができないので別の工夫も必要となりますね。振動数を大きくするためには、張力Tを大きくするほかに、線密度ρを小さくする（弦を細くする）、あるいは弦を張る距離Lを小さくするという2つの方法も考えられます。
ヴァイオリンとコントラバスを思い出してみましょう。ヴァイオリンは高い音を出せるように弦を張る距離（楽器の大きさ）は小さくなっていて、弦も細いものが用いられています。一方で、コントラバスは低い音を出せるように楽器の大きさは大きく、弦も太いものが用いられています。もしこれらの楽器を見る機会があれば、今回の内容を思い出しながら大きさに注目してみてください。
d’ Alembert の解法による解の求積
波動方程式(8)の解が式(35)で与えられました。この節では、式(8)を次の方法で変形すると式(35)が直接に得られることを示しましょう。新しい独立変数
y=x+vt、 z=x−vt
を導入すると、uがyとzの関数になり、式(8)の偏導関数がyとzについての偏導関数であらわされることになります。連鎖法則から、
∂u∂x=∂u∂y∂y∂x+∂u∂z∂z∂x
となります。この式の右辺に再び連鎖法則を適用しましょう。ここで現れる全ての偏微分が連続とすると∂2u/(∂z∂y)=∂2u/(∂y∂z)です。∂y/∂x=1と∂z/∂x=1を用いると、
∂2u∂x2=∂∂x(∂u∂y+∂u∂z)=∂∂y(∂u∂y+∂u∂z)∂y∂x+∂∂z(∂u∂y+∂u∂z)∂z∂x=∂2u∂y2+2∂2u∂y∂z+∂2u∂z2
となります。同様のやり方でtについての2階偏導関数も計算すると
∂2u∂t2=v2(∂2u∂y2−2∂2u∂y∂z+∂2u∂z2)
となります。これら2つの結果を波動方程式(8)に代入すると、
∂2u∂y∂z=0
が得られる。d’ Alembert の解法の利点は、式(38)が積分を2回繰り返すことで簡単に解けるという点です。まず、zについて積分すると、
∂u∂y=h(y)
となります。但し、h(y)はyの任意関数です。これをyについて積分すると、
u(y、z)=∫h(y)dy+ψ(z)
となります。但し、ψ(z)はzの任意関数である。積分はyのみの関数なので、この積分結果をϕ(y)とすると、uの解はu=ϕ(y)+ψ(z)の形になります。式(37)を用いると、
u(x,t)=ϕ(y)+ψ(z)=ϕ(x+vt)+ψ(x−vt)
となります。これが波動方程式(8)のd’ Alembert の解です。
初期条件を満たすd’ Alembert の解
ここでは以下の初期条件を課したd’ Alembert の解について議論します。
u(x,0)=f(x)
∂u∂t|t=0=g(x)
これは前に考えた初期条件と同じものです。式(39)を微分すると、
∂u∂t=vϕ′(x+vt)−vψ′(x−vt)
が得られます。但し、プライム′は引数x+vtとx−vtについての微分です。式(39-42)から、
u(x,0)=ϕ(x)+ψ(x)=f(x)
∂u∂t|t=0=vϕ′(x)−vϕ′(x)=g(x)
となります。式(44)をvで割ってxで積分すると、
ϕ(x)−ψ(x)=k(x0)+1v∫xx0g(s)ds where k(x0)=ϕ(x0)−ψ(x0)
となります。上式と式(43)を加えるとψが消去されます。これを2で割って以下の式を得ることができます。
ϕ(x)=12f(x)+12v∫xx0g(s)ds+12k(x0)
同様に式(43)から式(45)を引いて2で割ると
ψ(x)=12f(x)−12v∫xx0g(s)ds−12k(x0)
となります。式(46)でxをx+vtでおきかえると、積分はx0からx+vtまでになります。同様に式(47)でxをx−vtでおきかえると、積分はx0からx−vtまでになりますが、積分の符号を変えると積分範囲はx−vtからx0までになります。このようにしておいてϕ(x+vt)とψ(x−vt)を加えると、式(39)からu(x,t)が以下で与えられることが分かります。
u(x,t)=12{f(x+vt)+f(x−ct)}+12v∫x+vtx−vtg(s)ds
特に初速度が0ならば、上式は、
u(x,t)=12{f(x+vt)+f(x−ct)}
となります。これは式(35)と同じ結論です。境界条件(9)から、関数fは奇関数で周期2Lを持つことも明らかでしょう。このように、2つの初期条件と境界条件から解が一意的に決定できます。
このd’ Alembert の解の導出は変数分離の方法よりも洗練されているだろう。しかし、d’ Alembert の解法は限定的で、Fourier 級数による解法は様々な偏微分方程式に利用することができます。
