うさぎでもわかる微分方程式 Part06 定数係数線形微分方程式(同次式編)

スポンサードリンク

こんにちは、ももやまです。

今回から2階以上の線形微分方程式(基本は2階)の解き方について説明していきたいと思います。

まずは、2階以上の線形微分方程式の中でも、

  • 定数係数
  • 同次式

の解き方について説明していきたいと思います。

前回の微分方程式の記事はこちら!

www.momoyama-usagi.com

2階線形微分方程式ってどんなのだっけ、基本解ってなんだっけとなった人は必ず上にある前回の記事で復習しておきましょう。

スポンサードリンク

1.定数係数2階線形微分方程式とは

定数係数2階線形微分方程式とは、2階線形微分方程式d2ydx2+P(x)dydx+Q(x)y=R(x)P(x), Q(x) の部分が定数になっている微分方程式、つまりd2ydx2+adydx+by=R(x)となっている微分方程式のことを表します。(a, b は定数)

さらに、R(x)=0 のとき、つまりd2ydx2+adydx+by=0の形の場合、同次式の2階線形微分方程式となります。

今回はこの同次式の定数係数2階線形微分方程式について解いていきたいと思います。

スポンサードリンク

2.同次式定数係数2階線形微分方程式の一般解

(1) まずは同次式定数係数1階微分方程式の一般解から

まずは、復習として、同次式の定数係数1階微分方程式dydx+ay=0の一般解を出してみましょう。

この微分方程式は、両辺を y で割り、1ydydx=aと変数分離形にできるので、任意定数 C1 を用いて1y dy=a dxlog|y|=ax+C1log|y|=logeax+logeC1|y|=eaxeC1y=±eC1eaxと一般解を求めることができます。さらに、C=eC1 とすることで、一般解をy=Ceaxの形にすることができます。

(2) 定数係数2階線形微分方程式の基本解

1階の同次式の定数係数1階微分方程式dydx+ay=0の一般解は y=Ceax と、ekx の形になりましたね。

なので、2階の同次の定数係数微分方程式も y=ekx の形になりそうですよね!

実際に y=ekx と仮定して、同次式の2階線形微分方程式d2ydx2+adydx+by=0に代入してみましょう。dydx=kekx,   d2ydx2=k2ekxなので、k2ekx+akekx+bekx=0(k2+ak+b)ekx=0となりますね。

なので、赤色部分 k2+ak+b が0となるような k、つまり(1)k2+ak+b=0を満たすような k のとき、y=ekx は解になるといえますね。

この方程式 k2+ak+b=0 のことを特性方程式と呼びます。

特性方程式 k2+ak+b=02次の方程式なので、方程式を満たす k は重解のときを除いて必ず2つありますね。なので、解の数も2つ存在します。

このように、特性方程式から導き出せる1次独立な解のことを基本解と呼びます。

(基本解がもしわからない人 or 忘れた人はこちらの記事で復習しましょう。)

1つ例として、微分方程式d2ydx25dydx+6y=0を満たすような基本解を特性方程式から求めてみましょう。

特性方程式は、k25k+6=0となるので、(k2)(k3)=0より、特性方程式の解は k=2,3 となるので、y=e2x, y=e3x の2つが基本解となります。

スポンサードリンク

3.3つの解パターン

同次式定数係数2階線形微分方程式d2ydx2+adydx+by=0の一般解の形は、特性方程式k2+ak+b=0解の形(判別式)によって、3パターンに分けられます。具体的には、

  1. 特性方程式が異なる2つの実数解をもつとき
  2. 特性方程式が重解をもつとき
  3. 特性方程式が異なる2つの虚数解をもつとき

の3つに分けられます。

先に、解のパターンがどうなるかを書きましょう。

同次式の定数係数2階線形微分方程式の3つの解パターン同次式の定数係数2階線形微分方程式d2ydx2+adydx+by=0の解のパターンは、特性方程式k2+ak+b=0の2つの解 k1, k2 の形、つまり判別式D=a24bによって決まる。

(i) k1, k2相異なる2つの実数解のとき
D>0

  • 基本解:ek1x, ek2x
  • 一般解:C1ek1x+C2ek2x

(ii) k1, k2重解k1=k2)のとき
D=0

  • 基本解:ek1x, xek1x
  • 一般解:C1ek1x+xC2ek1x

(iii) k1, k2相異なる2つの(共役な)虚数解のとき
D>0

2つの解をそれぞれ k1=α+βi, k2=αβi と実部虚部に分けることで

  • 基本解:eαxcosβx, eαxsinβx
  • 一般解:eαx(C1cosβx+C2sinβx)

の3パターンにわけられる。

1つずつ見ていきましょう。

(1) 特性方程式が異なる2つの実数解をもつとき

最も基本的なパターンです。

特性方程式k2+ak+b=0が異なる2つの実数解 k1, k2 をもつとき、つまり判別式 DD=a24b>0となるときの基本解は ek1x, ek2x となります。

また、一般解は2つの任意定数 C1, C2 を用いてy=C1ek1x+C2ek2xと表すことができます。

例題1

微分方程式d2ydx23dydx10y=0の基本解と一般解を求めなさい。

解説1

特性方程式は、k23k10=0となるので、(k5)(k+2)=0とすることで特性方程式の解は k=5,2 となる。

よって、基本解は y=e5x, y=e2x となる。

また、それぞれの基本解を定数倍したものの和も解となるので、一般解は任意定数 C1, C2 を用いてy=C1e5x+C2e2xとなる。

(2) 特性方程式が重解となるとき

特性方程式m2+am+b=0は、2次方程式なので、基本的には2つの解を持ちます。しかし、判別式D=a24b=0となるときは重解となるので、解は1つしか持ちませんね。

(この重解を k としましょう。)

すると、基本解は ekx となります。

(1)と同じようにすると、一般解は、y=C1ekx+C2ekxとなりますが、この形は1つにまとめられてy=Cekxとなってしまいます。

あれ、2階の微分方程式なので、任意定数は2つ出てくるはずなのに1つになっちゃいました。

ということで ekx になるべく近い形(微分しても ekx の形を保つような関数)かつ ekx の定数倍ではないもう1つの基本解を探してみましょう。

ここで、y=ekxx 倍した y=xekx を考えてみましょう。

すると、dydx=ekx+k1xekx=(kx+1)ekx d2ydx2=kek+kekx+k2xekx=2kek+k12xek=(k2x+2k)ekとなるので、d2ydx2+adydx+by=0に代入すると、(k2x+2k)ek+a(kx+1)ekx+bxekx=0((k2+ak+b)x+(2k+a))ekx=0と変形できます。

ここで、ekx>0 なので、y=xekx が解を持つためには(k2+ak+b)x+(2k+a)=0を満たす必要があります。

ekx が基本解のとき、k2+ak+b=0は必ず満たすため(特性方程式)、ekx が基本解のとき、xekx が基本解になるかどうかは2k+a=0、つまり a=2k を満たすかどうか確認すればOKです。

特性方程式m2+am+b=0が重解をもつとき、判別式D=a24b=0の関係 b=14a2 が成立するので、a=2k を代入すると、b=14(2k)2=144k2=k2となるので、m2+am+b=m22km+k2=(mk)2=0となり、特性方程式が重解をもつときだけ、ekx が解であれば xekx も解となることが示せましたね。

例題2

微分方程式d2ydx24dydx+4y=0の基本解と一般解を求めなさい。

解説2

微分方程式の特性方程式は、k24k+4=0となるので、(k2)2=0とすることで特性方程式の解は k=2(2重解)となる。

よって、基本解は y=e2x, y=xe2x となる。

また、それぞれの基本解を定数倍したものの和も解となるので、一般解は任意定数 C1, C2 を用いてy=C1e2x+C2xe2xとなる。

(3) 特性方程式が虚数解となるとき

特性方程式m2+am+b=0の判別式が負になるとき、つまりD=a24b<0となるときは虚数解となりますよね。

このときの特性方程式の解は、a±D2=a2+D2iと実部と虚部で表すことができますね。

ここで、特性方程式の解を k1=α+βi, k2=αβi と実部と虚部で分けて表現します。

すると、微分方程式の基本解は y=e(α+βi)xy=e(α+βi)xとなり、一般解はy=A1e(α+βi)x+A2e(α+βi)xとなります。

でも、微分方程式の一般解に虚数解があるのはちょっと…と思いますよね。

なんとか実数だけで表現することはできないでしょうか。

ということで、虚数で表された一般解をなんとか実数範囲だけで表すことを考えてみましょう。

ここで、オイラーの公式eiθ=cosθ+isinθを思い出してみましょう。

(オイラーの公式についてはこちらの記事をご覧ください。)

すると、eiβx=cosβx+isinβxとなるので、一般解はy=A1e(α+βi)x+A2e(αβi)x=A1eαxeiβx+A2eαxeiβx=A1eαx(cosβx+isinβx)+A2eαx{cos(βx)+isin(βx)}=A1eαx(cosβx+isinβx)+A2eαx(cosβxisinβx)=eαx{(A1+A2)cosβx+i(A1A2)sinβx}と表すことができる。

さらに、任意定数を改めて C1=A1+A2, C2=i(A1A2) とすると、一般解をeαx(C1cosβx+C2sinβx)と虚数 i を使わずに表現することができます!

(基本解は eaxcosβxeaxsinβx となる。)

例題3

微分方程式d2ydx2+y=0の基本解と一般解を求めなさい。

解説3

微分方程式の特性方程式は、k2+1=0となるので、k2=1とすることで特性方程式の解は k=±i(2重解)となる。

よって、基本解は y=e0cosx=cosx, y=e0sinx=sinx となる。

また、それぞれの基本解を定数倍したものの和も解となるので、一般解は任意定数 C1, C2 を用いてy=C1cosx+C2sinxとなる。

4.ばね振り子の単振動(応用)

実際に同次式の定数係数2階線形微分方程式の応用して、ばね振り子の単振動について少しだけですが見ていきましょう。

下の図に示すような軽い(質量を無視できる)ばねに重さ m のおもりを付け合うと、左側のようにある地点でつりあいます。

f:id:momoyama1192:20200411150302g:plain

このつり合った点を0として、x だけ引っ張った状態を考えましょう。

すると、つりあった点に戻ろうとする力(復元力)が働きますね。この力 f は、ばね定数 k、引っ張った距離 x を用いてf=kxと表すことができます。(変位とは逆向きなのでマイナスが必要なことに注意!)

すると、運動方程式は加速度 d2xdt2 を用いてmd2xdt2=kxとなりますね。早速この時刻 t に対する微分方程式を解いていくことで、ある時間 t に対するおもりがいる場所(変位) x を求めてみましょう。

まず、両辺を m で割り、さらに移項させてd2xdt2+kmx=0とし、さらにω2=kmとすると、d2xdt2+ω2x=0となります。

すると、特性方程式は、k2+ω2=0となるので、k=±iω となりますね。

ここで、特性方程式の解が k=α±βi となるときの基本解は eatcosβt, eatsinβt 、一般解は任意定数 C1, C2 を用いてx=eαt(C1cosβt+C2sinβt)となるのでしたね。

なので、微分方程式の基本解は cosωt, sinωt となり、一般解はx=e0C1cosωt+e0C2sinωt=C1cosωt+C2sinωtと求めることができます。

(角速度が ω となるので、単振動の周期 TT=2πω=2πmkとなります。)

さらに、C=C12+C22とし、C1=Asinα,   C2=Acosαと任意定数を C, α におきかえることで、x=Asinαcosωt+Acosαsinωt=A(sinαcosωt+cosαsinωt)=Asin(α+ωt)=Asin(ωt+α)と表現することもできます。

(ちなみに任意定数のつけかたを変えて、一般解をx=Acos(ωt+α)としてもOKです*1

実際の力学の問題では、一般解に加え、2つの初期条件が加わります。

初期条件の例としては、

  • 距離4だけ引っ張った状態でスタート(t=0 のとき x=4
  • 静かにおもりを離した(初速度は0、dxdt=0

のようなものがあります。

実際に ω=2 としたときの上の初期条件を満たすような特解を求めてみましょう。

まず、x(0)=4 より、x=C1cos0+C2sin0=C1=4が成立し、さらにdxdt=2C1sin0+2C2cos0=0となるので、C2=0 となるので、特解をx=4cos2tと求めることができます。

5.高階定数係数微分方程式(応用)

同次式の定数係数2階線形微分方程式の解の形は、特性方程式の解の形によって求められると説明しました。

3階以上の同次式高階定数係数微分方程式の解の形も、特性方程式を解くことで一般解を求めることができます。

同次式の定数係数n階線形微分方程式の3つの解パターン同次式の定数係数 n 次線形微分方程式andnydxn+an1dn1ydxn1++a1dydx+a0y=0は、以下のステップによって求めることができる。

  1. 特性方程式ankn+an1kn1++a1k+a0=0を解いて n 個の基本解 y1, y2, … yn を求める。
  2. それぞれの基本解ごとに別々の任意定数を掛け、C1y1, C2y2, …, Cnyn としてからすべてを足し、y=C1y1+C2y2++Cnynとしたものが一般解となる。

基本解の組み合わせは、下の4パターンに分けることができます。

n 次線形微分方程式の基本解のパターンは、特性方程式の解の形で以下の4パターンに分けられる。

(1) 特性方程式の解 k重解ではない実数解のとき

→ 基本解は ekx となる。

(2) 特性方程式の解 km 重解(ただし mn)のとき

→ 基本解は、ekx, xekx, x2ekx, …, xm1ekx となる。(m 個の基本解が存在する)

(3) 特定方程式の解のうちの2つ k1, k2相異なる共役な虚数解のとき

k1=a+bi, k2=abi とすることで基本解は eaxcosbx, eaxsinbx となる。

(4) 特定方程式の解のうちの2つ k1, k2m 重解の共役な虚数解になるとき
[4次以上の微分方程式で可能性あり]

k1=a+bi, k2=abi とすることで基本解は

  • eaxcosbx, xeaxcosbx, …, xm1eaxcosbx
  • eaxsinbx, xeaxsinbx, …, xm1eaxsinbx

となる。(2m 個の基本解が存在する)

同次式の定数係数n階線形微分方程式の基本解の4パターン

3問ほど例題で練習してみましょう。

例題4(重解ではない実数解になる場合)

微分方程式d3ydx36d2ydx2+11dydx6y=0の基本解と一般解を求めなさい。

解説4

特性方程式は、k36k+11k6=0となるので、(k1)(k25k+6)=0(k1)(k2)(k3)=0とすることで特性方程式の解は k=1,2,3 となる。

よって、基本解は ex, e2x, e3x となる。

また、一般解は任意定数 C1, C2, C3 を用いてy=C1ex+C2e2x+C3e3xとなる。

例題5(重解の実数解が含まれる場合)

微分方程式d3ydx3d2ydx28dydx+12y=0の基本解と一般解を求めなさい。

解説5

特性方程式は、k3k28k+12=0となるので、(k2)(k2+k6)=0(k2)2(k+3)=0とすることで特性方程式の解は k=3 と、k=2 の2重解となる。

よって、基本解は e3x と e2x, xe2x となる。

また、一般解は任意定数 C1, C2, C3 を用いてy=C1e3x+C2e2x+C3xe2x=C1e3x+(C2+C3x)e2xとなる。

例題6(虚数解が出てくる場合)

微分方程式d3ydx3+d2ydx2+4dydx+4y=0の基本解と一般解を求めなさい。

解説6

特性方程式は、k3+k2+4k+4=0となるので、(k+1)(k2+4)=0(k+1)(k+2i)(k2i)=0とすることで特性方程式の解は k=1 と、k=2i, k=2i の2重解となる。

よって、基本解は ex, cos2x, sin2x となる。

また、一般解は任意定数 C1, C2, C3 を用いてy=C1ex+C2cos2x+C3sin2xとなる。

6.練習問題

では、定数係数線形微分方程式の基本解と一般解を実際に求める練習をしましょう。

全部で10問用意しています。後ろのほうほど階数が増えていきます。

練習1

微分方程式dydx10y=0の基本解と一般解を求めなさい。

練習2

微分方程式d2ydx23dydx=0の基本解と一般解を求めなさい。

練習3

微分方程式d2ydx22dydx8y=0の基本解と一般解を求めなさい。

練習4

微分方程式d2ydx2+8dydx+16y=0の基本解と一般解を求めなさい。

練習5

微分方程式d2ydx24dydx+13y=0の基本解と一般解を求めなさい。

練習6

微分方程式d3ydx32d2ydx2d3ydx3+2y=0の基本解と一般解を求めなさい。

練習7

微分方程式d3ydx3+2d2ydx215dydx36y=0の基本解と一般解を求めなさい。

練習8

微分方程式d3ydx39d2ydx2+27d3ydx327y=0の基本解と一般解を求めなさい。

練習9

微分方程式d4ydx42d3ydx3+5d2ydx2=0の基本解と一般解を求めなさい。

練習10

微分方程式d4ydx4+2d2ydx2+y=0の一般解を求めなさい。

7.練習問題の答え

解答1

特性方程式は、k10=0となるので、特性方程式の解は k=10 となる。

よって、基本解は e10xなる。

さらに、一般解は任意定数 C を用いてy=C1e10xとなる。

同次式の定数係数の1階線形微分方程式は、特性方程式で一気に求めましょう!

解答2

特性方程式は、k23k=0となるので、k(k3)=0とすることで特性方程式の解は k=0, k=3 となる。

よって、基本解は e0=1, e3xなる。

さらに、一般解は任意定数 C1, C2 を用いてy=C1+C2e3xとなる。

解答3

特性方程式は、k22k8=0となるので、(k4)(k+2)=0とすることで特性方程式の解は k=4, k=2 となる。

よって、基本解は e4x, e2xなる。

さらに、一般解は任意定数 C1, C2 を用いてy=C1e4x+C2e2xとなる。

解答4

特性方程式は、k2+8k+16=0となるので、(k+4)2=0とすることで特性方程式の解は k=4 の2重解となる。

よって、基本解は e4x, xe4xなる。

さらに、一般解は任意定数 C1, C2 を用いてy=C1e4x+C2xe4x=(C1+C2x)e4xとなる。

解答5

特性方程式は、k24k+13=0となるので、k=4±16522=4±362=4±6i2=2±3iとすることで特性方程式の解は k=2±3i となる。

よって、基本解は e2xcos3x, e2xsin3xなる。

さらに、一般解は任意定数 C1, C2 を用いてy=C1e2xcos3x+C2e2xsin3xとなる。

解答6

特性方程式は、k32k2k+2=0となるので、(k1)(k2k2)=0(k1)(k+1)(k2)=0とすることで特性方程式の解は k=1, k=1, k=2 となる。

よって、基本解は ex, ex, e2xなる。

さらに、一般解は任意定数 C1, C2, C3 を用いてy=C1ex+C2ex+C3e2xとなる。

解答7

特性方程式は、k2+2k15k36=0となるので、(t+3)(t2t12)=0(t4)(t+3)2=0とすることで特性方程式の解は k=3k=4 の2重解となる。

よって、基本解は e3x, e4x, xe4xなる。

さらに、一般解は任意定数 C1, C2, C3 を用いてy=C1e3x+C2e4x+C3xe4x=C1e3x+(C2+C3x)e4xとなる。

解答8

特性方程式は、k39k2+27k27=0となるので、(k3)(k26k+9)=0(k3)3=0とすることで特性方程式の解は k=3 の3重解となる。

よって、基本解は e3x, xe3x, x2e3xなる。

さらに、一般解は任意定数 C1, C2, C3 を用いてy=C1e3x+C2xe3x+C3x2e3x=(C1+C2x+C3x2)e3xとなる。

解答9

特性方程式は、k42k3+5k2=0となるので、k2(k22k+5)=0とすることで特性方程式の解は k=0 (2重解)と、k=1±2i となる。

よって、基本解は e0=1, xe0=x, excos2x, exsin2x となる。

さらに、一般解は任意定数 C1, C2, C3, C4 を用いてy=C1+xC2+C3excos2x+C4exsin2x=C1+xC2+ex(C3cos2x+C4sin2x)となる。

解答10

特性方程式は、k4+2k2+1=0となるので、t=k2 とすると、t2+2t+1=0となるので、(t+1)2=0となり、t=1 の2重解となる。

さらに、k2=1 を解くと、k=±i となるので、特性方程式の解は k=±i の2重解となる。

よって、基本解は cosx, sinx, xcosx, xsinxとなる。

さらに、一般解は任意定数 C1, C2, C3, C4 を用いてy=C1cosx+C2xcosx+C3sinx+C4xsinx=(C1+C2x)cosx+(C3+C4x)sinxとなる。

8.さいごに

今回は、

  • 2階を中心とした同次式の定数係数線形微分方程式の解き方
  • 基本解とロンスキアン

について解説していきました。

次回は、少し特殊な形をした同次式の「2階オイラー微分方程式」を定数係数微分方程式の形に持ち込んで解く方法について説明したいと思います。

*1:Asin(ωt+α)の任意定数を改めて β=α+π/2 とおけば、Asin(ωt+β+π2)=Acos(ωt+β)となり、どちらも同じ一般解になることが示せる。

関連広告・スポンサードリンク

おすすめの記事