うさぎでもわかる微分方程式 Part08 未定係数法を用いた定数係数線形微分方程式の特殊解の求め方

スポンサードリンク

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

前回は同次式の定数係数2階(n階)線形微分方程式について説明していきました。

今回は、非同次の定数係数線形微分方程式の4つの解き方

の中でも「未定係数法」を用いた方法について説明していきたいと思います。

未定係数法を含む残りの3つの方法の長所・短所も載せておくので、特殊解をどう求めようか迷った人はご覧ください。

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

www.momoyama-usagi.com

オイラー微分方程式に関する記事です。

スポンサードリンク

1.非同次式と特殊解

まずは、非同次式微分方程式の一般解がどのような形であったかを確認しましょう。

2階の非同次の線形微分方程式d2ydx2+adydx+b=R(x)の一般解は、

  • 同次方程式の一般解(R(x)=0 のときの一般解)
  • 非同次方程式を満たす1つの特殊解(なんでもいい)

2つの和で構成されています。

f:id:momoyama1192:20200411191930g:plainつまり、同次方程式d2ydx2+adydx+b=0の一般解を求め、さらにd2ydx2+adydx+b=R(x)を満たすような特殊解 y1つ求めることで非同次方程式の一般解を求めることができます。

なお、n階の非同次定数係数線形微分方程式andnydxn+an1dn1ydxn1++a1dydx+y=R(x)の場合も、2階のときと同じように同次方程式andnydxn+an1dn1ydxn1++a1dydx+y=0の一般解と「非同次方程式の特殊解1つ」の和で一般解を求めることができます。

スポンサードリンク

2.未定係数法とは

特殊解の求め方には、大まかにわけて

  • 未定係数法
  • 定数変化法
  • 微分演算子法

の3パターンあります。

今回は微分方程式の右辺 R(x) の部分から、特殊解の形を予測、仮定し、予想した特殊解の形を実際に代入することで、係数部分を求める未定係数法について説明していきたいと思います。

では、実際に未定係数法を用いて特殊解を求める流れを例題で説明していきたいと思います。

例題1

微分方程式d2ydx23dydx+2y=4xについて次の問いに答えなさい。

(1) 同次方程式d2ydx23dydx+2y=0の一般解を求めなさい。(復習)

(2) y=ax+b とおくことで特殊解を1つ見つけ、微分方程式の一般解を求めなさい。

解答1

(1)

これは復習ですね。

特性方程式k23k+2=0から、(k1)(k2)=0とし、k=1,2 となるので、任意定数 C1, C2 を用いてy=C1ex+C2e2xと表せます。

(2)

問題文の通り、y=ax+b とおいてみましょう。

すると、dydx=a,   d2ydx2=0となりますね。

問題文の微分方程式に代入すると、03a+2(ax+b)=4x2ax+(2b3a)=4xとなります。よって、{2a=43a+2b=0{a=2b=3となるので、特殊解の1つがy=2x+3とわかります。

よって、非同次方程式の一般解は、「同次方程式の一般解」と「特殊解の1つの和」で表されるので、y=C1ex+C2e2x+2x+3となる。

このように右辺 R(x) が多項式となっている場合、多項式の次数に合わせて特殊解を仮定します。

例1:R(x)=4 → 定数なので y=a と仮定
例2:R(x)=6x+4 → 1次の式なので y=ax+b と仮定
例3:R(x)=2x2+1 → 2次の式なので y=ax2+bx+c と仮定

多項式を微分しても次数は増えることがない(というか必ず1つ減る)ため、R(x) の次数がそのまま特殊解の次数となります*1

例題2 右辺が e^kx の場合

微分方程式d2ydx2dydx2y=e3xについて次の問いに答えなさい。

(1) 同次方程式d2ydx2dydx2y=0の一般解を求めなさい。(復習)

(2) 特殊解を1つ見つけ、微分方程式の一般解を求めなさい。

解説2

(1)

これは復習ですね。

特性方程式k2k2=0から、(k2)(k+1)=0とし、k=2,1 となるので、任意定数 C1, C2 を用いてy=C1e2x+C2exと表せます。

(2)

ekx を微分しても、ekx の形は崩れませんよね。

(微分していくと kekx, k2ekx, … となる)

そこで、右辺が R(x)=ekx の場合は、y=aekx とおきます

今回は、R(x)=e3x なので、y=ae3x とおきましょう。

すると、dydx=3ae3x,   d2ydx2=9ae3xとなりますね。

問題文の微分方程式に代入すると、9ae3x3ae3x2ae3x=e3x4ae3x=e3xとなります。よって、a=14となるので、特殊解の1つがy=14e3xとわかります。

よって、非同次方程式の一般解は、「同次方程式の一般解」と「特殊解の1つ」の和で表されるので、y=C1e2x+C2ex+14e3xとなります。

例題3 右辺が e^kx かつ同次方程式に基本解が含まれている場合

微分方程式d2ydx2dydx6y=e3xについて次の問いに答えなさい。

(1) 同次方程式d2ydx2dydx6y=0の一般解を求めなさい。(復習)

(2) 特殊解を1つ見つけ、微分方程式の一般解を求めなさい。

解説3

(1)

さすがに慣れてきたでしょう。

特性方程式k2k6=0から、(k3)(k+2)=0とし、k=3,2 となるので、任意定数 C1, C2 を用いてy=C1e3x+C2e2xと表せます。

(2)

ダメな例

ekx 型なので、y=ae3x とおきます。

すると、dydx=3ae3x,   d2ydx2=9ae3xとなりますね。

問題文の微分方程式に代入すると、9ae3x3ae3x6ae3x=e3x0=e3x

( ̄ー ̄?).....??ありゃ??

意味の分からない数式 0=e3x が出てしまいました。

実は、y=ae3x とおいてしまうと、一般解はy=C1e3x+C2e2x+ae3xのように、同次形の基本解の e3x と相殺されてしまい、y=C1e3x+C2e2xとなってしまいます。

正しい例

そこで、ekx 以外で、微分しても ekx の形が崩れない関数として、xekx を考えます。

今回の微分方程式の場合、R(x)=e3x なので、y=axe3x とおきます。

すると、dydx=ae3x+3axe3x d2ydx2=3ae3x+3ae3x+9axe3x=6ae3x+9axe3xとなりますね。

問題文の微分方程式に代入すると、6ae3x+9axe3x(ae3x+3axe3x)6axe3x=e3x(9a3a6a)xe3x+(6aa)e3x=e3x5ae3x=e3xとなり、うまく左辺の xe3x の項が消えてくれます。5a=1となるので、特殊解の1つがy=15xe3xとわかります。

よって、非同次方程式の一般解は、同次方程式の一般解と特殊解の1つの和で表されるので、y=C1e3x+C2e2x+15xe3xとなります。

例題4 右辺が sin kx / cos kx の形の場合

微分方程式d2ydx25dydx+6y=sinxについて次の問いに答えなさい。

(1) 同次方程式d2ydx25dydx+6y=0の一般解を求めなさい。(復習)

(2) 特殊解を1つ見つけ、微分方程式の一般解を求めなさい。

解説4

(1)

これは復習ですね。

特性方程式k25k+6=0から、(k2)(k3)=0とし、k=2,3 となるので、任意定数 C1, C2 を用いてy=C1e2x+C2e3xと表せます。

(2)

sinkx を微分していくと、kcoskx, k2sinkx, k3coskx, k4sinkx, … と「sinkxcoskx の形」になりますね。

そこで、右辺が R(x)=sinkx もしくは R(x)=coskx の場合は、y=asinkx+bcoskx とおきます

今回は、R(x)=sinx なので、y=asinx+bcosx とおきましょうか。

すると、dydx=acosxbsinxd2ydx2=asinxbcosxとなりますね。

問題文の微分方程式に代入すると、asinxbcosx5(acosxbsinx)+6(asinx+bcosx)=sinx(a+5b+6a)sinx+(b5a+6b)cosx=sinx(5a+5b)sinx+(5a+5b)cosx=sinxとなります。

よって、{5a+5b=15a+5b=0{a=110b=110となるので、特殊解の1つがy=110sinx+110cosxとわかります。

よって、非同次方程式の一般解は、「同次方程式の一般解」と「特殊解の1つの和」で表されるので、y=C1e2x+C2e3x+110sinx+110cosxとなります。

なお、右辺がR(x)=sinkx+coskxのように sinkxcoskx が両方含まれていた場合でも、R(x)=asinkx+bcoskxとおくことで特殊解を求められます。

例題5 右辺が sin kx / cos kx かつ同次方程式に基本解が含まれている場合

あまりないパターンですが、念のためやっておきましょう。

微分方程式d2ydx2+y=sinxについて次の問いに答えなさい。

(1) 同次方程式d2ydx2+y=0の一般解を求めなさい。(復習)

(2) 特殊解を1つ見つけ、微分方程式の一般解を求めなさい。

解説5

(1)

特性方程式k2+1=0から、(k+i)(ki)=0とし、k=±i となるので、任意定数 C1, C2 を用いてy=C1sinx+C2cosxと表せます。

(2)

R(x)=sinx なので、y=asinx+bcosx とおきたくなりますが、同次方程式の基本解が sinx, cosx なので、そのままおいてしまうと、一般解y=C1sinx+C2cosx+asinx+bcosxC1sinxC2cosx に相殺されてしまいます。

そこで、微分しても sinkx / coskx の形が崩れないように、特殊解をy=axsinx+bxcosxとおきます。すると、dydx=asinx+axcosx+bcosxbxsinx=(abx)sinx+(ax+b)cosxdydx=asinx+axcosx+bcosxbxsinx=(abx)sinx+(ax+b)cosxd2ydx2=bsinx+(abx)cosx+acosx(ax+b)sinx=(baxb)sinx+(abx+a)cosx=(ax2b)sinx+(2abx)cosxとなるので、代入すると、(ax2b)sinx+(2abx)cosx+axsinx+bxcosx=sinx2bsinx+2acosx=sinxとなります。

よって、{2b=12a=0{a=0b=12となるので、特殊解の1つがy=12xcosxとわかります。

よって、非同次方程式の一般解はy=C1sinx+C2cosx12xcosxとなります。

ここで、特殊解のおき方をまとめておきましょう。

非同次の定数係数線形微分方程式andnydxn+an1dn1ydxn1++a1dydx+y=R(x)の特殊解のおきかたのコツ

(1) R(x) が多項式の場合

→ 特殊解を R(x) の次数に合わせて特殊解を仮定

例1:R(x) が1次なら ax+b
例2:R(x) が3次なら ax3+bx2+cx+d

(2) R(x)=ekx の形の場合

→ 特殊解を aekx とおく

→ ただし、ekx が同次方程式の基本解に含まれていれば x 倍して axekx とおく

xekx も同次方程式の基本解に含まれていれば、さらに x 倍して ax2ekx とおく。

→ (以下略)

(3) R(x)=sinkx+coskx の形の場合
R(x)sinkx, coskx 単品でもOK)

→ 特殊解を asinkx+bcoskx とおく

→ ただし、sinkx もしくは coskx が同次方程式の基本解に含まれていれば x 倍して axsinkx+bxcoskx とおく

→ さらに、xsinkx もしくは xcoskx も同次方程式の基本解に含まれていれば、さらに x 倍して ax2sinkx+bx2coskx とおく

(以下略)

未定係数法による特殊解のおき方

未定係数法は、

  • 特殊解の形の予想が容易
  • 定数係数線形微分方程式である

場合にしか使えませんが、次回紹介する定数変化法、微分演算子法に比べて単純かつ簡単に特殊解を求めることができます。

スポンサードリンク

3.オイラー微分方程式への適用(応用)

前回同次式のオイラー微分方程式x2d2ydx2+axdydx+b=0を定数係数2階線形微分方程式d2ydt2+dydt+b=0にして解く方法を説明しましたね。

(オイラー微分方程式について興味がある人はこちらをご覧ください!)

オイラー微分方程式の場合でも、定数係数2階線形微分方程式に変形してあげることで未定係数法を適用させることができます。

1問だけ例題で説明しましょう。

例題6

微分方程式x2d2ydx23dydx+4y=x3について、次の問いに答えなさい。

(1) x=logt とおき、定数係数2階線形微分方程式の形にしなさい。

(2) (1)で求めた微分方程式の一般解を求めなさい。

(3) 与えられた微分方程式の一般解を求めなさい。

解説6

(1)

x=et、つまり t=logx とおくことで、d2ydt24dydt+4y=(x3)t=e3tと定数係数2階微分方程式にすることができる。

(2)

まず同次方程式d2ydt24dydt+4y=0の一般解を求める。特性方程式は、k24k+4=0となるので、(k2)2=0とすることで k=2 の2重解となり、同次方程式の一般解は任意定数 C1, C2 を用いてy=C1e2t+C2te2tとなる。

つぎに、特殊解を1つ探す。

右辺が R(x)=e3tekt の形になっているため、特殊解をy=ae3tとおく。

すると、dydx=3ae3t,   d2ydx2=9ae3tとなる。

問題文の微分方程式に代入すると、9ae3t12ae3t+4ae3tae3t=e3tとなります。よって、a=1となるので、特殊解の1つがy=e3tと求められる。

よって、非同次方程式の一般解はy=C1e2t+C2te2t+e3tとなる。

(3)

t=logx を代入することで、y=y(x) の形に戻す。

よって、x2d2ydx23dydx+4y=x3の一般解は、y=C1e2t+C2te2t+e3t=C1e2logx+C2logxe2logx+e3logx=C1elogx2+C2logxelogx2+elogx3=C1x2+C2x2logx+x3と求められる。

このように、定数係数線形微分方程式ではなかった場合でも、(2階の)オイラーの微分方程式であれば、定数係数におきかえることで未定係数法を適用することができます。

4.練習問題

では、実際に5問ほど練習してみましょう。

なお、オイラーの微分方程式は練習問題には入れていません。

練習1

微分方程式dydx3y=exの一般解を求めなさい。

(定数変化法を使わずに解いてみましょう。)

練習2

微分方程式d2ydx2+dydx2y=x2の一般解を求めなさい。

練習3

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

練習4

微分方程式d2ydx26dydx+9y=e3xの一般解を求めなさい。

練習5

微分方程式d2ydx2+4dydx+4y=cos2xの一般解を求めなさい。

練習6

微分方程式d3ydx32d2ydx25dydx+6y=3x25x+6の一般解を求めなさい。

5.練習問題の答え

解答1

まず同次方程式dydx3y=0の一般解を求める。特性方程式は、k3=0となるので、k=3 となり、同次方程式の一般解は任意定数 C を用いてy=Ce3xとなる。

つぎに、特殊解を1つ探す。

右辺が R(x)=exekx の形になっているため、特殊解をy=aexとおく。

(このとき ex が基本解に含まれていないかどうかを必ず確認すること!!)

すると、dydx=aexとなる。

問題文の微分方程式に代入すると、aex3aex=e3x2aex=e3xとなります。よって、a=12となるので、特殊解の1つがy=12exと求められる。

よって、一般解は、「同次形の一般解」+「特殊解の1つ」となるので、y=Ce3x12exが答え。

解答2

まず同次方程式d2ydx2+dydx2y=0の一般解を求める。特性方程式は、k2+k2=0となるので、(k+2)(k1)=0とすることで k=1,2 となり、同次方程式の一般解は任意定数 C1, C2 を用いてy=C1ex+C2e2xとなる。

つぎに、特殊解を1つ探す。

右辺が R(x)=x2 となっているため、特殊解をy=ax2+bx+cとおく。

すると、dydx=2ax+b,   d2ydx2=2aとなる。

問題文の微分方程式に代入すると、2a+(2ax+b)2(ax2+bx+c)=x22ax2+(2a2b)x+(2a+b2c)=x2となります。よって、{2a=12a2b=02a+b2c=0{a=12b=12c=34となるので、特殊解の1つがy=12x212x34と求められる。

よって、一般解は、「同次形の一般解」+「特殊解の1つ」となるので、y=C1ex+C2e2x12x212x34が答え。

解答3

まず同次方程式d2ydx2+2dydx+5y=0の一般解を求める。特性方程式は、k2+2k+5=0となるので、k=2±152=1±2iとすることで k=1±2i となり、同次方程式の一般解は任意定数 C1, C2 を用いてy=C1exsin2x+C2excos2xとなる。

つぎに、特殊解を1つ探す。

右辺が R(x)=e5xekx の形になっているため、特殊解をy=ae5xとおく。

(このとき e5x が基本解に含まれていないかどうかを必ず確認すること!!)

すると、dydx=5ae5x,   d2ydx2=25ae5xとなる。

問題文の微分方程式に代入すると、25ae5x+25ae5x+5ae5x=e5x40ae5x=e5xとなります。よって、a=140となるので、特殊解の1つがy=140e5xと求められる。

よって、一般解は、「同次形の一般解」+「特殊解の1つ」となるので、y=C1exsin2x+C2excos2x+140e5xが答え。

解答4

まず同次方程式d2ydx26dydx+9y=e3xの一般解を求める。特性方程式は、k26k+9=0となるので、(k3)2=0とすることでk=3(2重解)となり、同次方程式の一般解は任意定数 C1, C2 を用いてy=C1e3x+C2xe3xとなる。

つぎに、特殊解を1つ探す。

右辺が R(x)=e3x となっているため、特殊解をy=ae3xとおきたいが、e3x, xe3x も同次方程式の基本解となっているので、特殊解をy=ax2e3xとおく。

すると、dydx=2axe3x+3ax2e3x

d2ydx2=2ae3x+6axe3x+6axe3x+9ax2e3x=2ae3x+12axe3x+9ax2e3xとなる。

問題文の微分方程式に代入すると、2ae3x+12axe3x+9ax2e3x6(2axe3x+3ax2e3x)+9ax2e3x=e3x(9a18a+9a)x2e3x+(12a12a)xe3x+2ae3x=e3x2ae3x=e3xとなります。よって、a=12となるので、特殊解の1つがy=12x2e3xと求められる。

よって、一般解は、「同次形の一般解」+「特殊解の1つ」となるので、y=C1e3x+C2xe3x+12x2e3xが答え。

解答5

まず同次方程式d2ydx2+4dydx+4y=0の一般解を求める。特性方程式は、k2+4k+4=0となるので、(k+2)2=0とすることでk=2(2重解)となり、同次方程式の一般解は任意定数 C1, C2 を用いてy=C1e2x+C2xe2xとなる。

つぎに、特殊解を1つ探す。

右辺が R(x)=cos2xcoskx の形になっているため、特殊解をy=asin2x+bcos2xとおく。

すると、dydx=2acos2x2bsin2xd2ydx2=4asin2x4bcos2xとなる。

問題文の微分方程式に代入すると、4asin2x4bcos2x+4(2acos2x2bsin2x)+4(asin2x+bcos2x)=cos2x(4a8b+4a)sin2x+(4b+8a+4b)cos2x=cos2x8bsin2x+8acos2x=cos2xとなります。よって、{8a=18b=0{a=18b=0となるので、特殊解の1つがy=18sin2xと求められる。

よって、一般解は、「同次形の一般解」+「特殊解の1つ」となるので、y=C1e2x+C2xe2x+18sin2xが答え。

解答6

まず同次方程式d3ydx32d2ydx25dydx+6y=0の一般解を求める。特性方程式は、k32k25k+6=0となるので、(k1)(k2k+6)=0(k1)(k+2)(k3)=0とすることでk=1,2,3 となり、同次方程式の一般解は任意定数 C1, C2, C3 を用いてy=C1ex+C2e2x+C3e3xとなる。

つぎに、特殊解を1つ探す。

右辺が R(x)=3x25x+6 と2次の多項式になっているため、特殊解をy=ax2+bx+cとおく。

すると、dydx=2ax+b  d2ydx2=2ad3ydx3=0となる。

問題文の微分方程式に代入すると、022a5(2ax+b)+6(ax2+bx+c)=3x25x+66ax2+(10a+6b)x+(4a5b+6c)=3x25x+6となります。よって、{6a=310a+6b=54a5b+6c=6{a=12b=0c=43となるので、特殊解の1つがy=12x2+43と求められる。

よって、一般解は、「同次形の一般解」+「特殊解の1つ」となるので、y=C1ex+C2e2x+C3e3x+12x2+43が答え。

6.さいごに

今回は、非同次の定数係数線形微分方程式を未定係数法を用いて解く方法について説明しました。

特殊解の形が簡単に推測できる場合は、今回のように未定係数法を使って特殊解を求めましょう。

次回は、非同次方程式の解き方の2つ目として、定数変化法について説明していきたいと思います。

*1:例えば、R(x) が3次式なら、特殊解 y=y(x) は3次式、dydx は2次式、d2ydx2 は1次式となりますね。そのため、y=ax3+bx2+cx+d と3次の形でおいて微分方程式に代入して恒等式を解くことで特殊解が求められます。

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

おすすめの記事