MathLabs

微分方程式と力学系

線形微分方程式

重ね合わせの原理を持つ高階方程式で、特性方程式によって解く。

直観直感:振動と重ね合わせの原理

ばねにつながれた質点、振幅の小さい振り子、RLC回路の電流はすべて同じ形の法則に従う:未知関数とその導関数が線形に組み合わさる。二つの別々の運動がそれぞれ法則を満たすなら、それらの和(または重み付き組み合わせ)も法則を満たす——これが重ね合わせの原理である。だからこそ、たった二つの基本解からあらゆる運動を組み立てられる。

特性多項式 r^2-3r+2 の放物線。r軸と r=1、r=2 で交わる。
減衰調和振動子 y′′+2ζy′+y=0y'' + 2\zeta y' + y = 0:解は指数関数的に減衰する包絡線 ±e−ζt\pm e^{-\zeta t} の内側で振動する。

中高正確な記述:定数係数2階線形方程式

定義: 斉次線形微分方程式

a,b,ca,b,c を定数、a≠0a\neq 0 として、ay′′+by′+cy=0ay'' + by' + cy = 0 の形の2階方程式を斉次線形方程式と呼ぶ。右辺がゼロで、すべての項で yy が1次でしか現れないからである。

ay′′+by′+cy=0ay'' + by' + cy = 0

erxe^{rx} を微分するとそれ自身に rr を掛けたものになるため、試行解 y=erxy=e^{rx} を用いる。方程式に代入し共通因子 erx≠0e^{rx}\neq 0 で割ると、微分方程式は rr に関する単純な2次方程式になる:

ar2+br+c=0ar^2+br+c=0
特性方程式ar2+br+c=0ar^2+br+c=0の3つの場合と対応する一般解
判別式根一般解
Δ>0\Delta>0二つの異なる実根 r1≠r2r_1\neq r_2y=c1er1x+c2er2xy = c_1 e^{r_1 x} + c_2 e^{r_2 x}
Δ=0\Delta=0一つの重根 rry=(c1+c2x)erxy = (c_1 + c_2 x) e^{r x}
Δ<0\Delta<0共役複素根 r=α±iβr = \alpha \pm i\betay=eαx(c1cos⁡(βx)+c2sin⁡(βx))y = e^{\alpha x}\left(c_1 \cos(\beta x) + c_2 \sin(\beta x)\right)

大学大学レベル:定理と証明

y1y_1 と y2y_2 がともにay′′+by′+cy=0ay'' + by' + cy = 0を満たすなら、任意の定数 c1,c2c_1,c_2 に対して組み合わせ y=c1y1+c2y2y=c_1y_1+c_2y_2 も方程式を満たす。さらにロンスキアンW(y1,y2)=y1y2′−y2y1′W(y_1,y_2) = y_1y_2' - y_2y_1'がある点で0でなければ、方程式のすべての解はこの形をとる。

なぜ正しいのか?

微分は線形演算である:和の導関数は導関数の和であり、定数は外に出せる。方程式の左辺は yy、y′y'、y′′y'' に定数を掛けて足し合わせたものだけでできているので、二つの解の組み合わせを代入すると、単に二つのゼロを足すことになる。

証明

L[y]=ay′′+by′+cyL[y] = ay''+by'+cy とする。微分は定数 kk に対して (u+v)′=u′+v′(u+v)'=u'+v' および (ku)′=ku′(ku)'=ku' を満たすので、演算子 LL は線形である:L[c1y1+c2y2]=a(c1y1+c2y2)′′+b(c1y1+c2y2)′+c(c1y1+c2y2)=c1(ay1′′+by1′+cy1)+c2(ay2′′+by2′+cy2)=c1L[y1]+c2L[y2]L[c_1y_1+c_2y_2] = a(c_1y_1+c_2y_2)'' + b(c_1y_1+c_2y_2)' + c(c_1y_1+c_2y_2) = c_1(ay_1''+by_1'+cy_1) + c_2(ay_2''+by_2'+cy_2) = c_1L[y_1]+c_2L[y_2]。

y1,y2y_1,y_2 は方程式の解なので L[y1]=0L[y_1]=0 かつ L[y2]=0L[y_2]=0 であり、したがって L[c1y1+c2y2]=c1⋅0+c2⋅0=0L[c_1y_1+c_2y_2]=c_1\cdot 0+c_2\cdot 0=0:組み合わせも解である。これで前半が示された。

後半については、線形常微分方程式の理論により2階方程式の解空間はちょうど2次元であることが保証される(初期条件 y(x0)=y0y(x_0)=y_0、y′(x0)=y0′y'(x_0)=y_0' が唯一の解を定め、二つの自由パラメータを与える)。ある点 x0x_0 でロンスキアン W(y1,y2)=y1y2′−y2y1′W(y_1,y_2)=y_1y_2'-y_2y_1' が0でないことは、連立方程式 c1y1(x0)+c2y2(x0)=y0c_1y_1(x_0)+c_2y_2(x_0)=y_0、c1y1′(x0)+c2y2′(x0)=y0′c_1y_1'(x_0)+c_2y_2'(x_0)=y_0' がどんな初期値に対しても一意な解 (c1,c2)(c_1,c_2) を持つことを意味する。よってすべての解はある c1y1+c2y2c_1y_1+c_2y_2 に一致し、{y1,y2}\{y_1,y_2\} が2次元の解空間全体を張ることになる。

ay′′+by′+cy=0ay'' + by' + cy = 0 の一般解は特性方程式 ar2+br+c=0ar^2+br+c=0 の根によって完全に決まる:異なる実根なら y=c1er1x+c2er2xy = c_1 e^{r_1 x} + c_2 e^{r_2 x}、重根なら y=(c1+c2x)erxy = (c_1 + c_2 x) e^{r x}、共役複素根 r=α±iβr = \alpha \pm i\beta なら y=eαx(c1cos⁡(βx)+c2sin⁡(βx))y = e^{\alpha x}\left(c_1 \cos(\beta x) + c_2 \sin(\beta x)\right) となる。

なぜ正しいのか?

指数関数が自然な構成要素となるのは、erxe^{rx} を微分すると単に rr 倍されるだけだからである。この試行解を定数係数の線形方程式に入れると、微積分が代数——多項式の根を求める問題——に帰着する。

証明

試行解 y=erxy=e^{rx} を用いる。このとき y′=rerxy'=re^{rx}、y′′=r2erxy''=r^2e^{rx} であり、ay′′+by′+cy=0ay'' + by' + cy = 0 に代入すると (ar2+br+c)erx=0(ar^2+br+c)e^{rx}=0 となる。erxe^{rx} は決して0にならないので、特性方程式 ar2+br+c=0ar^2+br+c=0 が導かれる。2次方程式の解の公式より r=−b±b2−4ac2ar = \dfrac{-b \pm \sqrt{b^2-4ac}}{2a}。

Δ>0\Delta>0 の場合:二つの根 r1≠r2r_1\neq r_2 から解 er1xe^{r_1x} と er2xe^{r_2x} が得られる。そのロンスキアンは (r2−r1)e(r1+r2)x≠0(r_2-r_1)e^{(r_1+r_2)x}\neq 0 であり、重ね合わせの定理により一般解は y=c1er1x+c2er2xy = c_1 e^{r_1 x} + c_2 e^{r_2 x} となる。

Δ=0\Delta=0 の場合:重根 r=−b/(2a)r=-b/(2a) が一つだけあり、指数解 erxe^{rx} は一つしか得られない。もう一つの独立解を見つけるため y2=v(x)erxy_2=v(x)e^{rx}(階数低下法)を試して方程式に代入すると、rr が ar2+br+car^2+br+c の重根であるため vv と v′v' の項が打ち消し合い、av′′erx=0av''e^{rx}=0、すなわち v′′=0v''=0 が残るので v=c1+c2xv=c_1+c_2x。これにより第二の解 xerxxe^{rx} が得られ、一般解は y=(c1+c2x)erxy = (c_1 + c_2 x) e^{r x} となる。

Δ<0\Delta<0 の場合:根は共役複素数 r=α±iβr = \alpha \pm i\beta であり α=−b/(2a)\alpha=-b/(2a)、β=−Δ/(2a)\beta=\sqrt{-\Delta}/(2a)。オイラーの公式 eiβx=cos⁡(βx)+isin⁡(βx)e^{i\beta x}=\cos(\beta x)+i\sin(\beta x) により、複素解 e(α±iβ)xe^{(\alpha\pm i\beta)x} は二つの独立な実解 eαxcos⁡(βx)e^{\alpha x}\cos(\beta x) と eαxsin⁡(βx)e^{\alpha x}\sin(\beta x) になる(LL の係数が実数であるため、実部と虚部もそれぞれ解になることによる)。よって一般解は y=eαx(c1cos⁡(βx)+c2sin⁡(βx))y = e^{\alpha x}\left(c_1 \cos(\beta x) + c_2 \sin(\beta x)\right) となる。

大学実世界での応用と具体例

特性方程式の3つの場合は単なる代数ではない——減衰ばね質量系 my′′+cy′+ky=0my'' + cy' + ky = 0(質量 mm、減衰係数 cc、剛性 kk)の3つの物理的に異なる挙動を表している:振動しながら減衰する不足減衰(Δ<0\Delta<0)、まったく振動せずに戻る過減衰(Δ>0\Delta>0)、そして振動せずに戻る最速の臨界減衰(Δ=0\Delta=0)——これはまさに自動車のショックアブソーバーやドアクローザーの設計目標である。

例: 一般解を求める:異なる実根

y′′−5y′+6y=0y'' - 5y' + 6y = 0 を解け。

解答

特性方程式は r2−5r+6=0r^2-5r+6=0 である。

因数分解すると (r−2)(r−3)=0(r-2)(r-3)=0 となり、根は r1=2r_1=2 と r2=3r_2=3:異なる二つの実根である。

上の定理(Δ>0\Delta>0 の場合)により、一般解は y=c1e2x+c2e3xy = c_1e^{2x}+c_2e^{3x} であり、c1,c2c_1,c_2 は初期条件によって定まる任意定数である。

例: 非同次方程に対する未定係数法

y′′−y=e2xy'' - y = e^{2x} の特殊解を求めよ。

解答

斉次方程 y′′−y=0y''-y=0 の特性方程式は r2−1=0r^2-1=0、根は r=±1r=\pm 1 である。強制項 e2xe^{2x} の指数は 22 でこれらの根とは異なるため、共鳴は起きない。

試行解 yp=Ae2xy_p = Ae^{2x} を用いる。すると yp′′=4Ae2xy_p''=4Ae^{2x} となり、代入すると 4Ae2x−Ae2x=3Ae2x4Ae^{2x}-Ae^{2x}=3Ae^{2x} となり、これは e2xe^{2x} に等しくなければならない。

係数を比較すると 3A=13A=1 なので A=13A=\tfrac{1}{3} となり、特殊解 yp=13e2xy_p=\tfrac{1}{3}e^{2x} が得られる。完全な一般解は y=c1ex+c2e−x+13e2xy=c_1e^{x}+c_2e^{-x}+\tfrac{1}{3}e^{2x} である。

y′′−7y′+12y=0y''-7y'+12y=0 の一般解は何か。

y′′+2y′+5y=0y''+2y'+5y=0 について、特性方程式の根はどのような種類で、一般解は何か。

自動車のショックアブソーバーは my′′+cy′+ky=0my'' + cy' + ky = 0 でモデル化され、まったく跳ねることなく車体をできるだけ速く静止状態に戻すように設計される。技術者は特性方程式のどの場合を目指すべきか。

y′′−y=e2xy''-y=e^{2x} の特殊解を未定係数法で求めるには、どの試行解を使うべきか。

参考文献

  1. William E. Boyce, Richard C. DiPrima (2017). Elementary Differential Equations and Boundary Value Problems
  2. Lawrence Perko (2001). Differential Equations and Dynamical Systems