← 戻る ライブラリ › 微分方程式と力学系 › 常微分方程式 微分方程式と力学系
線形微分方程式 重ね合わせの原理を持つ高階方程式で、特性方程式によって解く。
直観 直感:振動と重ね合わせの原理 ばねにつながれた質点、振幅の小さい振り子、RLC回路の電流はすべて同じ形の法則に従う:未知関数とその導関数が線形に 組み合わさる。二つの別々の運動がそれぞれ法則を満たすなら、それらの和(または重み付き組み合わせ)も法則を満たす——これが重ね合わせの原理 である。だからこそ、たった二つの基本解からあらゆる運動を組み立てられる。
減衰調和振動子 y ′ ′ + 2 ζ y ′ + y = 0 y'' + 2\zeta y' + y = 0 y ′′ + 2 ζ y ′ + y = 0 :解は指数関数的に減衰する包絡線 ± e − ζ t \pm e^{-\zeta t} ± e − ζ t の内側で振動する。 中高 正確な記述:定数係数2階線形方程式 定義: 斉次線形微分方程式
a , b , c a,b,c a , b , c を定数、a ≠ 0 a\neq 0 a = 0 として、a y ′ ′ + b y ′ + c y = 0 ay'' + by' + cy = 0 a y ′′ + b y ′ + cy = 0 の形の2階方程式を斉次線形 方程式と呼ぶ。右辺がゼロで、すべての項で y y y が1次でしか現れないからである。
a y ′ ′ + b y ′ + c y = 0 ay'' + by' + cy = 0 a y ′′ + b y ′ + cy = 0 e r x e^{rx} e r x を微分するとそれ自身に r r r を掛けたものになるため、試行解 y = e r x y=e^{rx} y = e r x を用いる。方程式に代入し共通因子 e r x ≠ 0 e^{rx}\neq 0 e r x = 0 で割ると、微分方程式は r r r に関する単純な2次方程式になる:
a r 2 + b r + c = 0 ar^2+br+c=0 a r 2 + b r + c = 0 特性方程式a r 2 + b r + c = 0 ar^2+br+c=0 a r 2 + b r + c = 0 の3つの場合と対応する一般解 判別式 根 一般解 Δ > 0 \Delta>0 Δ > 0 二つの異なる実根 r 1 ≠ r 2 r_1\neq r_2 r 1 = r 2 y = c 1 e r 1 x + c 2 e r 2 x y = c_1 e^{r_1 x} + c_2 e^{r_2 x} y = c 1 e r 1 x + c 2 e r 2 x Δ = 0 \Delta=0 Δ = 0 一つの重根 r r r y = ( c 1 + c 2 x ) e r x y = (c_1 + c_2 x) e^{r x} y = ( c 1 + c 2 x ) e r x Δ < 0 \Delta<0 Δ < 0 共役複素根 r = α ± i β r = \alpha \pm i\beta r = α ± i β y = e α x ( c 1 cos ( β x ) + c 2 sin ( β x ) ) y = e^{\alpha x}\left(c_1 \cos(\beta x) + c_2 \sin(\beta x)\right) y = e α x ( c 1 cos ( β x ) + c 2 sin ( β x ) )
大学 大学レベル:定理と証明 y 1 y_1 y 1 と y 2 y_2 y 2 がともにa y ′ ′ + b y ′ + c y = 0 ay'' + by' + cy = 0 a y ′′ + b y ′ + cy = 0 を満たすなら、任意の定数 c 1 , c 2 c_1,c_2 c 1 , c 2 に対して組み合わせ y = c 1 y 1 + c 2 y 2 y=c_1y_1+c_2y_2 y = c 1 y 1 + c 2 y 2 も方程式を満たす。さらにロンスキアンW ( y 1 , y 2 ) = y 1 y 2 ′ − y 2 y 1 ′ W(y_1,y_2) = y_1y_2' - y_2y_1' W ( y 1 , y 2 ) = y 1 y 2 ′ − y 2 y 1 ′ がある点で0でなければ、方程式のすべての解はこの形をとる。
なぜ正しいのか? 微分は線形演算である:和の導関数は導関数の和であり、定数は外に出せる。方程式の左辺は y y y 、y ′ y' y ′ 、y ′ ′ y'' y ′′ に定数を掛けて足し合わせたものだけでできているので、二つの解の組み合わせを代入すると、単に二つのゼロを足すことになる。
証明 L [ y ] = a y ′ ′ + b y ′ + c y L[y] = ay''+by'+cy L [ y ] = a y ′′ + b y ′ + cy とする。微分は定数 k k k に対して ( u + v ) ′ = u ′ + v ′ (u+v)'=u'+v' ( u + v ) ′ = u ′ + v ′ および ( k u ) ′ = k u ′ (ku)'=ku' ( k u ) ′ = k u ′ を満たすので、演算子 L L L は線形である:L [ c 1 y 1 + c 2 y 2 ] = a ( c 1 y 1 + c 2 y 2 ) ′ ′ + b ( c 1 y 1 + c 2 y 2 ) ′ + c ( c 1 y 1 + c 2 y 2 ) = c 1 ( a y 1 ′ ′ + b y 1 ′ + c y 1 ) + c 2 ( a y 2 ′ ′ + b y 2 ′ + c y 2 ) = c 1 L [ y 1 ] + c 2 L [ y 2 ] 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] L [ c 1 y 1 + c 2 y 2 ] = a ( c 1 y 1 + c 2 y 2 ) ′′ + b ( c 1 y 1 + c 2 y 2 ) ′ + c ( c 1 y 1 + c 2 y 2 ) = c 1 ( a y 1 ′′ + b y 1 ′ + c y 1 ) + c 2 ( a y 2 ′′ + b y 2 ′ + c y 2 ) = c 1 L [ y 1 ] + c 2 L [ y 2 ] 。
y 1 , y 2 y_1,y_2 y 1 , y 2 は方程式の解なので L [ y 1 ] = 0 L[y_1]=0 L [ y 1 ] = 0 かつ L [ y 2 ] = 0 L[y_2]=0 L [ y 2 ] = 0 であり、したがって L [ c 1 y 1 + c 2 y 2 ] = c 1 ⋅ 0 + c 2 ⋅ 0 = 0 L[c_1y_1+c_2y_2]=c_1\cdot 0+c_2\cdot 0=0 L [ c 1 y 1 + c 2 y 2 ] = c 1 ⋅ 0 + c 2 ⋅ 0 = 0 :組み合わせも解である。これで前半が示された。
後半については、線形常微分方程式の理論により2階方程式の解空間はちょうど2次元であることが保証される(初期条件 y ( x 0 ) = y 0 y(x_0)=y_0 y ( x 0 ) = y 0 、y ′ ( x 0 ) = y 0 ′ y'(x_0)=y_0' y ′ ( x 0 ) = y 0 ′ が唯一の解を定め、二つの自由パラメータを与える)。ある点 x 0 x_0 x 0 でロンスキアン W ( y 1 , y 2 ) = y 1 y 2 ′ − y 2 y 1 ′ W(y_1,y_2)=y_1y_2'-y_2y_1' W ( y 1 , y 2 ) = y 1 y 2 ′ − y 2 y 1 ′ が0でないことは、連立方程式 c 1 y 1 ( x 0 ) + c 2 y 2 ( x 0 ) = y 0 c_1y_1(x_0)+c_2y_2(x_0)=y_0 c 1 y 1 ( x 0 ) + c 2 y 2 ( x 0 ) = y 0 、c 1 y 1 ′ ( x 0 ) + c 2 y 2 ′ ( x 0 ) = y 0 ′ c_1y_1'(x_0)+c_2y_2'(x_0)=y_0' c 1 y 1 ′ ( x 0 ) + c 2 y 2 ′ ( x 0 ) = y 0 ′ がどんな初期値に対しても一意な解 ( c 1 , c 2 ) (c_1,c_2) ( c 1 , c 2 ) を持つことを意味する。よってすべての解はある c 1 y 1 + c 2 y 2 c_1y_1+c_2y_2 c 1 y 1 + c 2 y 2 に一致し、{ y 1 , y 2 } \{y_1,y_2\} { y 1 , y 2 } が2次元の解空間全体を張ることになる。
a y ′ ′ + b y ′ + c y = 0 ay'' + by' + cy = 0 a y ′′ + b y ′ + cy = 0 の一般解は特性方程式 a r 2 + b r + c = 0 ar^2+br+c=0 a r 2 + b r + c = 0 の根によって完全に決まる:異なる実根なら y = c 1 e r 1 x + c 2 e r 2 x y = c_1 e^{r_1 x} + c_2 e^{r_2 x} y = c 1 e r 1 x + c 2 e r 2 x 、重根なら y = ( c 1 + c 2 x ) e r x y = (c_1 + c_2 x) e^{r x} y = ( c 1 + c 2 x ) e r x 、共役複素根 r = α ± i β r = \alpha \pm i\beta r = α ± i β なら y = e α x ( c 1 cos ( β x ) + c 2 sin ( β x ) ) y = e^{\alpha x}\left(c_1 \cos(\beta x) + c_2 \sin(\beta x)\right) y = e α x ( c 1 cos ( β x ) + c 2 sin ( β x ) ) となる。
なぜ正しいのか? 指数関数が自然な構成要素となるのは、e r x e^{rx} e r x を微分すると単に r r r 倍されるだけだからである。この試行解を定数係数の線形方程式に入れると、微積分が代数——多項式の根を求める問題——に帰着する。
証明 試行解 y = e r x y=e^{rx} y = e r x を用いる。このとき y ′ = r e r x y'=re^{rx} y ′ = r e r x 、y ′ ′ = r 2 e r x y''=r^2e^{rx} y ′′ = r 2 e r x であり、a y ′ ′ + b y ′ + c y = 0 ay'' + by' + cy = 0 a y ′′ + b y ′ + cy = 0 に代入すると ( a r 2 + b r + c ) e r x = 0 (ar^2+br+c)e^{rx}=0 ( a r 2 + b r + c ) e r x = 0 となる。e r x e^{rx} e r x は決して0にならないので、特性方程式 a r 2 + b r + c = 0 ar^2+br+c=0 a r 2 + b r + c = 0 が導かれる。2次方程式の解の公式より r = − b ± b 2 − 4 a c 2 a r = \dfrac{-b \pm \sqrt{b^2-4ac}}{2a} r = 2 a − b ± b 2 − 4 a c 。
Δ > 0 \Delta>0 Δ > 0 の場合:二つの根 r 1 ≠ r 2 r_1\neq r_2 r 1 = r 2 から解 e r 1 x e^{r_1x} e r 1 x と e r 2 x e^{r_2x} e r 2 x が得られる。そのロンスキアンは ( r 2 − r 1 ) e ( r 1 + r 2 ) x ≠ 0 (r_2-r_1)e^{(r_1+r_2)x}\neq 0 ( r 2 − r 1 ) e ( r 1 + r 2 ) x = 0 であり、重ね合わせの定理により一般解は y = c 1 e r 1 x + c 2 e r 2 x y = c_1 e^{r_1 x} + c_2 e^{r_2 x} y = c 1 e r 1 x + c 2 e r 2 x となる。
Δ = 0 \Delta=0 Δ = 0 の場合:重根 r = − b / ( 2 a ) r=-b/(2a) r = − b / ( 2 a ) が一つだけあり、指数解 e r x e^{rx} e r x は一つしか得られない。もう一つの独立解を見つけるため y 2 = v ( x ) e r x y_2=v(x)e^{rx} y 2 = v ( x ) e r x (階数低下法)を試して方程式に代入すると、r r r が a r 2 + b r + c ar^2+br+c a r 2 + b r + c の重根であるため v v v と v ′ v' v ′ の項が打ち消し合い、a v ′ ′ e r x = 0 av''e^{rx}=0 a v ′′ e r x = 0 、すなわち v ′ ′ = 0 v''=0 v ′′ = 0 が残るので v = c 1 + c 2 x v=c_1+c_2x v = c 1 + c 2 x 。これにより第二の解 x e r x xe^{rx} x e r x が得られ、一般解は y = ( c 1 + c 2 x ) e r x y = (c_1 + c_2 x) e^{r x} y = ( c 1 + c 2 x ) e r x となる。
Δ < 0 \Delta<0 Δ < 0 の場合:根は共役複素数 r = α ± i β r = \alpha \pm i\beta r = α ± i β であり α = − b / ( 2 a ) \alpha=-b/(2a) α = − b / ( 2 a ) 、β = − Δ / ( 2 a ) \beta=\sqrt{-\Delta}/(2a) β = − Δ / ( 2 a ) 。オイラーの公式 e i β x = cos ( β x ) + i sin ( β x ) e^{i\beta x}=\cos(\beta x)+i\sin(\beta x) e i β x = cos ( β x ) + i sin ( β x ) により、複素解 e ( α ± i β ) x e^{(\alpha\pm i\beta)x} e ( α ± i β ) x は二つの独立な実解 e α x cos ( β x ) e^{\alpha x}\cos(\beta x) e α x cos ( β x ) と e α x sin ( β x ) e^{\alpha x}\sin(\beta x) e α x sin ( β x ) になる(L L L の係数が実数であるため、実部と虚部もそれぞれ解になることによる)。よって一般解は y = e α x ( c 1 cos ( β x ) + c 2 sin ( β x ) ) y = e^{\alpha x}\left(c_1 \cos(\beta x) + c_2 \sin(\beta x)\right) y = e α x ( c 1 cos ( β x ) + c 2 sin ( β x ) ) となる。
大学 実世界での応用と具体例 特性方程式の3つの場合は単なる代数ではない——減衰ばね質量系 m y ′ ′ + c y ′ + k y = 0 my'' + cy' + ky = 0 m y ′′ + c y ′ + k y = 0 (質量 m m m 、減衰係数 c c c 、剛性 k k k )の3つの物理的に異なる挙動を表している:振動しながら減衰する不足減衰 (Δ < 0 \Delta<0 Δ < 0 )、まったく振動せずに戻る過減衰 (Δ > 0 \Delta>0 Δ > 0 )、そして振動せずに戻る最速の臨界減衰 (Δ = 0 \Delta=0 Δ = 0 )——これはまさに自動車のショックアブソーバーやドアクローザーの設計目標である。
例: 一般解を求める:異なる実根
y ′ ′ − 5 y ′ + 6 y = 0 y'' - 5y' + 6y = 0 y ′′ − 5 y ′ + 6 y = 0 を解け。
解答 特性方程式は r 2 − 5 r + 6 = 0 r^2-5r+6=0 r 2 − 5 r + 6 = 0 である。
因数分解すると ( r − 2 ) ( r − 3 ) = 0 (r-2)(r-3)=0 ( r − 2 ) ( r − 3 ) = 0 となり、根は r 1 = 2 r_1=2 r 1 = 2 と r 2 = 3 r_2=3 r 2 = 3 :異なる二つの実根である。
上の定理(Δ > 0 \Delta>0 Δ > 0 の場合)により、一般解は y = c 1 e 2 x + c 2 e 3 x y = c_1e^{2x}+c_2e^{3x} y = c 1 e 2 x + c 2 e 3 x であり、c 1 , c 2 c_1,c_2 c 1 , c 2 は初期条件によって定まる任意定数である。
例: 非同次方程に対する未定係数法
y ′ ′ − y = e 2 x y'' - y = e^{2x} y ′′ − y = e 2 x の特殊解を求めよ。
解答 斉次方程 y ′ ′ − y = 0 y''-y=0 y ′′ − y = 0 の特性方程式は r 2 − 1 = 0 r^2-1=0 r 2 − 1 = 0 、根は r = ± 1 r=\pm 1 r = ± 1 である。強制項 e 2 x e^{2x} e 2 x の指数は 2 2 2 でこれらの根とは異なるため、共鳴は起きない。
試行解 y p = A e 2 x y_p = Ae^{2x} y p = A e 2 x を用いる。すると y p ′ ′ = 4 A e 2 x y_p''=4Ae^{2x} y p ′′ = 4 A e 2 x となり、代入すると 4 A e 2 x − A e 2 x = 3 A e 2 x 4Ae^{2x}-Ae^{2x}=3Ae^{2x} 4 A e 2 x − A e 2 x = 3 A e 2 x となり、これは e 2 x e^{2x} e 2 x に等しくなければならない。
係数を比較すると 3 A = 1 3A=1 3 A = 1 なので A = 1 3 A=\tfrac{1}{3} A = 3 1 となり、特殊解 y p = 1 3 e 2 x y_p=\tfrac{1}{3}e^{2x} y p = 3 1 e 2 x が得られる。完全な一般解は y = c 1 e x + c 2 e − x + 1 3 e 2 x y=c_1e^{x}+c_2e^{-x}+\tfrac{1}{3}e^{2x} y = c 1 e x + c 2 e − x + 3 1 e 2 x である。
よくある誤り. 強制項がそれ自身、斉次方程式の解になっている場合——例えば y ′ ′ − 4 y ′ + 4 y = e 2 x y'' - 4y' + 4y = e^{2x} y ′′ − 4 y ′ + 4 y = e 2 x では、斉次の特性方程式 r 2 − 4 r + 4 = 0 r^2-4r+4=0 r 2 − 4 r + 4 = 0 が重根 r = 2 r=2 r = 2 を持ち、強制項が e 2 x e^{2x} e 2 x である——単純な試行解 A e 2 x Ae^{2x} A e 2 x は失敗する(それは斉次方程式の解であるため、代入すると 0 = e 2 x 0=e^{2x} 0 = e 2 x という矛盾が生じる)。この共鳴 のケースでは、一致する根の重複度を k k k として試行解に x k x^k x k を掛ける必要がある:ここでは y p = A x 2 e 2 x y_p = Ax^2e^{2x} y p = A x 2 e 2 x 。 歴史的ノート
1743年、レオンハルト・オイラーは定数係数の線形方程式に y = e r x y=e^{rx} y = e r x を代入すると、r r r に関する多項式方程式——ここで使われている特性方程式——に帰着することを示し、微分方程式を解く問題を代数に変えた。その過程で振動する場合を説明する等式 e i β x = cos ( β x ) + i sin ( β x ) e^{i\beta x}=\cos(\beta x)+i\sin(\beta x) e i β x = cos ( β x ) + i sin ( β x ) も発見した。
レオンハルト・オイラー
y ′ ′ − 7 y ′ + 12 y = 0 y''-7y'+12y=0 y ′′ − 7 y ′ + 12 y = 0 の一般解は何か。
y = c 1 e 3 x + c 2 e 4 x y=c_1e^{3x}+c_2e^{4x} y = c 1 e 3 x + c 2 e 4 x y = c 1 e 7 x + c 2 e 12 x y=c_1e^{7x}+c_2e^{12x} y = c 1 e 7 x + c 2 e 12 x y = ( c 1 + c 2 x ) e 3 x y=(c_1+c_2x)e^{3x} y = ( c 1 + c 2 x ) e 3 x y = e 3 x ( c 1 cos ( 4 x ) + c 2 sin ( 4 x ) ) y=e^{3x}(c_1\cos(4x)+c_2\sin(4x)) y = e 3 x ( c 1 cos ( 4 x ) + c 2 sin ( 4 x )) y ′ ′ + 2 y ′ + 5 y = 0 y''+2y'+5y=0 y ′′ + 2 y ′ + 5 y = 0 について、特性方程式の根はどのような種類で、一般解は何か。
共役複素根;y = e − x ( c 1 cos ( 2 x ) + c 2 sin ( 2 x ) ) y=e^{-x}(c_1\cos(2x)+c_2\sin(2x)) y = e − x ( c 1 cos ( 2 x ) + c 2 sin ( 2 x )) 異なる二つの実根;y = c 1 e − x + c 2 e 5 x y=c_1e^{-x}+c_2e^{5x} y = c 1 e − x + c 2 e 5 x 重根;y = ( c 1 + c 2 x ) e − x y=(c_1+c_2x)e^{-x} y = ( c 1 + c 2 x ) e − x 共役複素根;y = e 2 x ( c 1 cos ( − x ) + c 2 sin ( − x ) ) y=e^{2x}(c_1\cos(-x)+c_2\sin(-x)) y = e 2 x ( c 1 cos ( − x ) + c 2 sin ( − x )) 自動車のショックアブソーバーは m y ′ ′ + c y ′ + k y = 0 my'' + cy' + ky = 0 m y ′′ + c y ′ + k y = 0 でモデル化され、まったく跳ねることなく車体をできるだけ速く静止状態に戻すように設計される。技術者は特性方程式のどの場合を目指すべきか。
Δ = 0 \Delta=0 Δ = 0 、臨界減衰:c 2 = 4 m k c^2=4mk c 2 = 4 mk Δ < 0 \Delta<0 Δ < 0 、不足減衰:c 2 < 4 m k c^2<4mk c 2 < 4 mk Δ > 0 \Delta>0 Δ > 0 、c c c をできるだけ大きくした過減衰減衰なし、c = 0 c=0 c = 0 y ′ ′ − y = e 2 x y''-y=e^{2x} y ′′ − y = e 2 x の特殊解を未定係数法で求めるには、どの試行解を使うべきか。
y p = A e 2 x y_p=Ae^{2x} y p = A e 2 x y p = A x e 2 x y_p=Axe^{2x} y p = A x e 2 x y p = A x 2 e 2 x y_p=Ax^2e^{2x} y p = A x 2 e 2 x y p = A cos ( 2 x ) + B sin ( 2 x ) y_p=A\cos(2x)+B\sin(2x) y p = A cos ( 2 x ) + B sin ( 2 x )