Excel VBA 数学教室ではアフィリエイトプログラムを利用して商品を紹介しています。

定数係数二階斉次方程式

二階線形微分方程式

一般に次のような形の微分方程式
 (A)d2ydx2+pdydx+qy=r(x)
定数係数二階線形微分方程式 とよびます (p,q は定数)。また上式で r(x)=0 とした方程式
 (B)d2ydx2+pdydx+qy=0
斉次方程式(homogeneous equation)といい、さらに p=0 とした方程式
 (C)d2ydx2+qy=0
を二階線形微分方程式の 標準形 (canonical form) とよびます。一般に斉次方程式 (B) は標準形 (C) に帰着できます。y=f(x)z(x) とおいて (A1) に代入すると
 fz+(2f+pf)z+(f+pf+qf)z=0
という方程式が得られます。ここで z の係数 2f+pf が 0 になるように f を定めます。すなわち
 2f+pf=0
を解いて
 f(x)=epx/2
が得られます。すると z の係数は
 f+pf+qf=(qp24)f
となるので、斉次方程式
 (B)d2ydx2+pdydx+qy=0
の一般解は
 y=epx/2z(x)
で表されます。ここに z(x) は標準形
 d2ydz2+(qp24)z=0
の解です。

標準形方程式の解

標準形方程式
 (C)d2ydx2+qy=0
の解はその形から指数関数であることが予測できます。実際、複素数 k を用いて y=ekx とおいて上式に代入してみると
 (k2+q)ekx=0
すなわち k2=q という関係式が得られます。
q の符号によって解は 2 種類の形に分類されます。m を正の実数とすると
 
q=m2<0、すなわち k=±m のとき、一般解は
 (C1)y=aemx+bemx
という実数の指数関数で与えられます。
 
q=m2>0、すなわち k=±im のとき、一般解は
 (C2)y=ceimx+deimx
というように複素指数関数で与えられます。

斉次方程式の解

斉次方程式
 (B)d2ydx2+pdydx+qy=0
の解の1つは y=eλx の形で与えられます。代入してみると
 (B1)λ2+pλ+q=0
という λ に関する 2 次方程式が得られますが、この式は 特性方程式 (characteristic equation) とよばれます。そこで (B) の一般解を
 y=eλ1xz(x)
とおいて代入してみると
 eλ1x[z+(2λ1+p)z+(λ12+pλ1+q)z)]=0
となりますが、特性方程式により
 λ12+pλ1+q=0
なので、
 z+(2λ1+p)z=0
という式が得られます。これに e(2λ1+p)x をかけると
 ddx[e(2λ1+p)xz=0]
という形に書き直せます。特性方程式
 (B1)λ2+pλ+q=0
の解と係数の関係 λ1+λ2=p を用いると
 ddx[e(λ1λ2)xz=0]
積分すると c を定数として
 e(λ1λ2)xz=c
これを解いて
 z={cλ2λ1e(λ2λ1)x(λ1λ2)cx+a(λ1=λ2)
となります。よって斉次方程式
 (B)d2ydx2+pdydx+qy=0
の一般解は y=eλ1xz(x) より
 (B2)y={aeλ1x+beλ2x(λ1λ2)(ax+b)eλx(λ1=λ2)
となります。ただし λ1,λ2 は特性方程式
 (B1)λ2+pλ+q=0
の判別式 D の値によって次のように分類されます。

D=p24q>0 のとき、r=p2/4q とおくと
 λ1=p2+r,λ2=p2r
D=p24q=0 のとき、
 λ1=λ2=p2
D=p24q<0 のとき、s=qp2/4 とおくと
 λ1=p2+is,λ2=p2is
斉次方程式は特性方程式を用いて簡単に解を求めることができます。たとえば
 y3y+2y=0
の特性方程式は
 λ23λ+2=0
です。これを解くと λ=1,2 なので一般解は
 y=aex+be2x
となります。また、
 y4y+4y=0
の特性方程式は
 λ24λ+4=0
です。これを解いて λ=2 、一般解は
 y=(a+bx)e2x
となります。

エクセルや数学に関するコメントをお寄せください