2階定数係数線形微分方程式の解の公式|特性方程式の3つの場合の証明・未定係数法の例題【大学院入試・微分方程式】

2階定数係数線形微分方程式の解の公式|特性方程式の3つの場合の証明・未定係数法の例題【大学院入試・微分方程式】

合格大作戦 算数・数学公式辞典 › 関数

大学院入試

y′′+ay′+by=0,λ2+aλ+b=0 ⟹ y=C1eλ1x+C2eλ2xy''+ay'+by=0,\quad\lambda^2+a\lambda+b=0\ \Longrightarrow\ y=C_1e^{\lambda_1x}+C_2e^{\lambda_2x}

重解なら (C1+C2x)eλx(C_1+C_2x)e^{\lambda x}、虚数解 α±βi\alpha\pm\beta i なら eαx(C1cos⁡βx+C2sin⁡βx)e^{\alpha x}(C_1\cos\beta x+C_2\sin\beta x)

この公式のポイント

  • y=eλxy=e^{\lambda x} を入れると特性方程式 λ2+aλ+b=0\lambda^2+a\lambda+b=0
  • 実数解・重解・虚数解の3つで一般解の形が決まる
  • 重解は xeλxxe^{\lambda x}、虚数解はオイラーの公式で三角関数に
  • 非同次は同次の一般解+特殊解
うかるくん

うかるくん

特性方程式が重解のとき、C1eλx+C2eλxC_1e^{\lambda x}+C_2e^{\lambda x} じゃだめなの?

合格先生

合格先生

それは (C1+C2)eλx(C_1+C_2)e^{\lambda x} で、定数が実質1つだ。1階の方程式2つに分けて解くと、(C1+C2x)eλx(C_1+C_2x)e^{\lambda x} が自然に出てくるぞ。

2階定数係数線形微分方程式の解の公式

y′′+ay′+by=0y''+ay'+by=0(a, ba,\ b は定数)の形の微分方程式は、y=eλxy=e^{\lambda x} を入れてできる2次方程式(特性方程式)の解で、一般解がすべて書けます。ばねの振動や電気回路(RLC回路)を表す、もっとも大事な微分方程式です。

同次方程式 y″+ay′+by=0 の一般解

特性方程式 λ2+aλ+b=0\lambda^2+a\lambda+b=0 の解によって

① 異なる2つの実数解 λ1, λ2\lambda_1,\ \lambda_2:y=C1eλ1x+C2eλ2xy=C_1e^{\lambda_1x}+C_2e^{\lambda_2x}

② 重解 λ\lambda:y=(C1+C2x)eλxy=(C_1+C_2x)e^{\lambda x}

③ 虚数解 α±βi\alpha\pm\beta i(β≠0\beta\ne0):y=eαx(C1cos⁡βx+C2sin⁡βx)y=e^{\alpha x}(C_1\cos\beta x+C_2\sin\beta x)

非同次方程式 y′′+ay′+by=f(x)y''+ay'+by=f(x) の一般解は

y=(同次方程式の一般解)+(非同次方程式の特殊解1つ)y=(\text{同次方程式の一般解})+(\text{非同次方程式の特殊解1つ})
xy0123451y″+3y′+2y=0:λ=−1, −2(実数解2つ)y″+2y′+y=0:λ=−1(重解)y″+2y′+5y=0:λ=−1±2i(虚数解)どれも y(0)=1, y′(0)=0 から出発
同じ初期条件でも、特性方程式の解の型で動きがちがう(③だけ振動する)

証明

① 2つの1階の方程式に分ける

特性方程式の解を λ1, λ2\lambda_1,\ \lambda_2 とすると、解と係数の関係より a=−(λ1+λ2)a=-(\lambda_1+\lambda_2)、b=λ1λ2b=\lambda_1\lambda_2。そこで z=y′−λ2yz=y'-\lambda_2y とおくと

z′−λ1z=y′′−(λ1+λ2)y′+λ1λ2y=y′′+ay′+by=0z'-\lambda_1z=y''-(\lambda_1+\lambda_2)y'+\lambda_1\lambda_2y=y''+ay'+by=0

よって z′=λ1zz'=\lambda_1z で、z=Aeλ1xz=Ae^{\lambda_1x}(1階線形微分方程式)。

② もう一度1階の方程式を解く

y′−λ2y=Aeλ1xy'-\lambda_2y=Ae^{\lambda_1x} に積分因子 e−λ2xe^{-\lambda_2x} をかけると (e−λ2xy)′=Ae(λ1−λ2)x(e^{-\lambda_2x}y)'=Ae^{(\lambda_1-\lambda_2)x}。

  • λ1≠λ2\lambda_1\ne\lambda_2 なら、積分して e−λ2xy=Aλ1−λ2e(λ1−λ2)x+C2e^{-\lambda_2x}y=\dfrac{A}{\lambda_1-\lambda_2}e^{(\lambda_1-\lambda_2)x}+C_2、つまり y=C1eλ1x+C2eλ2xy=C_1e^{\lambda_1x}+C_2e^{\lambda_2x}。
  • λ1=λ2=λ\lambda_1=\lambda_2=\lambda なら、右辺は定数 AA なので e−λxy=Ax+C1e^{-\lambda x}y=Ax+C_1、つまり y=(C1+C2x)eλxy=(C_1+C_2x)e^{\lambda x}。

この計算はすべての解を求めているので、解はこれで全部です(任意の定数は2つ)。

③ 虚数解のときは実数の形に直す

λ=α±βi\lambda=\alpha\pm\beta i のとき、①の式は複素数のまま成り立ちます。オイラーの公式 e(α+βi)x=eαx(cos⁡βx+isin⁡βx)e^{(\alpha+\beta i)x}=e^{\alpha x}(\cos\beta x+i\sin\beta x) を使い、C1, C2C_1,\ C_2 をまとめなおすと y=eαx(C1cos⁡βx+C2sin⁡βx)y=e^{\alpha x}(C_1\cos\beta x+C_2\sin\beta x)。実際に入れて確かめることもできます。

④ 非同次方程式

ypy_p を特殊解とすると、どんな解 yy も差 y−ypy-y_p は同次方程式をみたします(左辺が1次式だから)。よって一般解=同次の一般解+ypy_p です。

y″+ay′+by=0特性方程式 λ²+aλ+b=0D>0実数解 λ₁≠λ₂C₁e^(λ₁x)+C₂e^(λ₂x)D=0重解 λ(C₁+C₂x)×e^(λx)D<0虚数解 α±βie^(αx)(C₁cosβx+C₂sinβx)D=a²−4b。C₁, C₂ は2つの初期条件で決まる
特性方程式の判別式 D=a2−4bD=a^2-4b で、一般解の形が3つに分かれる

大学院入試では

院試では、特性方程式から一般解を書き、初期条件で C1, C2C_1,\ C_2 を決める問題と、非同次方程式の特殊解(未定係数法)を求める問題がよく出ます。

大学院入試

例題1(異なる実数解)

y′′−5y′+6y=0, y(0)=1, y′(0)=0y^{\prime\prime}-5y^{\prime}+6y=0,\ y(0)=1,\ y^{\prime}(0)=0 を解きなさい。

解き方

λ2−5λ+6=0\lambda^2-5\lambda+6=0 より λ=2, 3\lambda=2,\ 3。y=C1e2x+C2e3xy=C_1e^{2x}+C_2e^{3x}。

y(0)=C1+C2=1y(0)=C_1+C_2=1、y′(0)=2C1+3C2=0y^{\prime}(0)=2C_1+3C_2=0 より C1=3, C2=−2C_1=3,\ C_2=-2。

答え y=3e2x−2e3xy=3e^{2x}-2e^{3x}

大学院入試

例題2(虚数解:単振動)

y′′+4y=0, y(0)=0, y′(0)=2y^{\prime\prime}+4y=0,\ y(0)=0,\ y^{\prime}(0)=2 を解きなさい。

解き方

λ2+4=0\lambda^2+4=0 より λ=±2i\lambda=\pm2i(α=0, β=2\alpha=0,\ \beta=2)。y=C1cos⁡2x+C2sin⁡2xy=C_1\cos2x+C_2\sin2x。

y(0)=C1=0y(0)=C_1=0、y′(0)=2C2=2y^{\prime}(0)=2C_2=2 より C2=1C_2=1。

答え y=sin⁡2xy=\sin2x

大学院入試

例題3(非同次:未定係数法)

y′′−3y′+2y=2xy^{\prime\prime}-3y^{\prime}+2y=2x の一般解と、y′′−3y′+2y=exy^{\prime\prime}-3y^{\prime}+2y=e^x の特殊解を求めなさい。

解き方

同次の特性方程式 λ2−3λ+2=0\lambda^2-3\lambda+2=0 より λ=1, 2\lambda=1,\ 2。

右辺 2x2x:yp=Ax+By_p=Ax+B とおくと −3A+2(Ax+B)=2x-3A+2(Ax+B)=2x。2A=2, −3A+2B=02A=2,\ -3A+2B=0 より A=1, B=32A=1,\ B=\dfrac32。一般解は C1ex+C2e2x+x+32C_1e^x+C_2e^{2x}+x+\dfrac32。

右辺 exe^x:λ=1\lambda=1 が特性方程式の解なので AexAe^x では左辺が 00 になり失敗。yp=Axexy_p=Axe^x とおくと、yp′=A(1+x)exy_p^{\prime}=A(1+x)e^x、yp′′=A(2+x)exy_p^{\prime\prime}=A(2+x)e^x で、左辺 =A{(2+x)−3(1+x)+2x}ex=−Aex=A\{(2+x)-3(1+x)+2x\}e^x=-Ae^x。A=−1A=-1。

答え 一般解 y=C1ex+C2e2x+x+32y=C_1e^x+C_2e^{2x}+x+\dfrac32、特殊解 yp=−xexy_p=-xe^x

学部の授業とのつながり・院試での注意

  • どの授業で習うか:大学1〜2年の「微分方程式」、物理の力学(ばねの振動・減衰振動・強制振動)、電気回路で習います。図1の3つの型は、物理では過減衰・臨界減衰・減衰振動とよばれます。
  • 答案で示すこと:①特性方程式とその解、②一般解の形(どの型か)、③非同次なら特殊解の求め方(おき方の理由)、④初期条件から C1, C2C_1,\ C_2 を決める連立方程式。
  • よく出る形:右辺が多項式・ekxe^{kx}・sin⁡kx\sin kx・cos⁡kx\cos kx の未定係数法、kk が特性方程式の解と重なる「共鳴」の場合(xx をかける)、連立1階方程式を行列の固有値で解く問題、ラプラス変換で初期値問題を解く問題。

よくある間違い

腕でバツを作る合格先生

合格先生

右辺の ekxe^{kx} の kk が特性方程式の解なら、AekxAe^{kx} ではうまくいかない。AxekxAxe^{kx} とおくんだぞ。

  • 重解で C1eλx+C2eλxC_1e^{\lambda x}+C_2e^{\lambda x} と書く:これは (C1+C2)eλx(C_1+C_2)e^{\lambda x} で、任意の定数が実質1つしかありません。重解では xx をかけた xeλxxe^{\lambda x} が必要です。
  • 共鳴を見落とす:右辺の ekxe^{kx} の kk が特性方程式の解なら、AekxAe^{kx} ではなく AxekxAxe^{kx}(重解なら Ax2ekxAx^2e^{kx})とおきます。
  • 初期条件を先に入れる:特殊解を足す前の同次の一般解に初期条件を入れるのは誤り。一般解=同次+特殊解をそろえてから C1, C2C_1,\ C_2 を決めます。

練習問題

体で丸を作る合格先生

合格先生

特性方程式の解の型を見きわめれば、もう迷わないぞ。

  1. 大学院入試 y′′+2y′+y=0, y(0)=1, y′(0)=1y^{\prime\prime}+2y^{\prime}+y=0,\ y(0)=1,\ y^{\prime}(0)=1 を解きなさい。
    答えと解説を見るλ=−1\lambda=-1(重解)で y=(C1+C2x)e−xy=(C_1+C_2x)e^{-x}。y(0)=C1=1y(0)=C_1=1、y′(0)=C2−C1=1y^{\prime}(0)=C_2-C_1=1 より C2=2C_2=2。答え y=(1+2x)e−xy=(1+2x)e^{-x}
  2. 大学院入試 y′′+2y′+5y=0y^{\prime\prime}+2y^{\prime}+5y=0 の一般解を求めなさい。
    答えと解説を見るλ2+2λ+5=0\lambda^2+2\lambda+5=0 より λ=−1±2i\lambda=-1\pm2i。答え y=e−x(C1cos⁡2x+C2sin⁡2x)y=e^{-x}(C_1\cos2x+C_2\sin2x)
  3. 大学院入試 y′′+y=x2y^{\prime\prime}+y=x^2 の一般解を求めなさい。
    答えと解説を見る同次は λ=±i\lambda=\pm i で C1cos⁡x+C2sin⁡xC_1\cos x+C_2\sin x。yp=Ax2+Bx+Cy_p=Ax^2+Bx+C とおくと 2A+Ax2+Bx+C=x22A+Ax^2+Bx+C=x^2 より A=1, B=0, C=−2A=1,\ B=0,\ C=-2。答え y=C1cos⁡x+C2sin⁡x+x2−2y=C_1\cos x+C_2\sin x+x^2-2

関連する公式

よくある質問

特性方程式とは?

y″+ay′+by=0 に y=e^(λx) を入れると e^(λx)(λ²+aλ+b)=0 となることから得られる2次方程式 λ²+aλ+b=0 のことです。

重解のとき、なぜ x をかけるのですか?

e^(λx) だけでは任意の定数が1つしかなく、解が足りないからです。1階の方程式2つに分けて解くと、自然に (C₁+C₂x)e^(λx) が出てきます。

虚数解のとき三角関数になるのはなぜ?

オイラーの公式 e^((α+βi)x)=e^(αx)(cosβx+i sinβx) を使って、実数の解に組み直すと e^(αx)(C₁cosβx+C₂sinβx) になるからです。

非同次方程式はどう解きますか?

同次方程式の一般解に、非同次方程式の特殊解を1つ足します。特殊解は右辺の形に合わせて未定係数法で求めます。

あわせて読みたい

出典・参考

  • 大学初年級の微分積分学(常微分方程式)

最終更新:2026年10月11日/作成:大学受験合格大作戦