大学・微分方程式・微分・積分

2階線形微分方程式(定数係数)の練習問題

特性方程式から一般解を求め、振動・減衰の振る舞いが分かる。

考え方

ばねの振動や電気回路のように、加速度(2階導関数)を含む y′′+py′+qy=0y'' + py' + qy = 0 の形の方程式を、特性方程式から解けるようになります。揺れながら止まるのか、揺れずに戻るのかを、係数から判断できます。

特性方程式
y=eλty = e^{\lambda t} を代入して得られる λ\lambda の2次方程式 λ2+pλ+q=0\lambda^2 + p\lambda + q = 0。
一般解
2つの任意定数 C1C_1, C2C_2 を含む解。初期条件 y(0)y(0), y′(0)y'(0) で定数が決まる。
減衰振動・過減衰
揺れながら小さくなる動きと、揺れずにゆっくり戻る動き。

指数関数 eλte^{\lambda t} は、微分しても λ\lambda 倍されるだけで形が変わりません。そこで y=eλty = e^{\lambda t} とおいて y′′+py′+qy=0y'' + py' + qy = 0 に代入すると、eλt(λ2+pλ+q)=0e^{\lambda t}(\lambda^2 + p\lambda + q) = 0 となり、特性方程式を解く問題になります。

解が異なる2つの実数 λ1\lambda_1, λ2\lambda_2 なら、一般解は y=C1eλ1t+C2eλ2ty = C_1e^{\lambda_1 t} + C_2e^{\lambda_2 t} です。どちらも負なら、揺れずに 0 に戻ります(過減衰)。重解 λ\lambda なら y=(C1+C2t)eλty = (C_1 + C_2t)e^{\lambda t} です。

判別式が負で、解が α±βi\alpha \pm \beta i(複素数)のときは、y=eαt(C1cos⁡βt+C2sin⁡βt)y = e^{\alpha t}(C_1\cos\beta t + C_2\sin\beta t) になります。実部 α\alpha が減衰の速さ、虚部 β\beta が揺れの速さ(角振動数)を表し、α<0\alpha < 0 なら揺れながら小さくなります。

最後に初期条件で定数を決めます。異なる実数解なら、y(0)=C1+C2y(0) = C_1 + C_2、y′(0)=λ1C1+λ2C2y'(0) = \lambda_1C_1 + \lambda_2C_2 の連立方程式を解きます。

y = e^(-0.5t) cos 3t のグラフ。点線の包絡線 ±e^(-0.5t) の間で揺れながら、0 に近づいていく12345-11Oty包絡線 ±e^(-0.5t)
y=e−0.5tcos⁡3ty = e^{-0.5t}\cos 3t(特性方程式の解が −0.5±3i-0.5 \pm 3i)。点線の ±e−0.5t\pm e^{-0.5t} の間で揺れながら小さくなる減衰振動です。解が負の実数2つなら、揺れずに 0 に戻ります。

ポイント

特性方程式: y′′+py′+qy=0⇒λ2+pλ+q=0y'' + py' + qy = 0 \Rightarrow \lambda^2 + p\lambda + q = 0

異なる実数解 λ1,λ2\lambda_1, \lambda_2: y=C1eλ1t+C2eλ2ty = C_1e^{\lambda_1 t} + C_2e^{\lambda_2 t} 重解 λ\lambda: y=(C1+C2t)eλty = (C_1 + C_2t)e^{\lambda t}

複素数解 α±βi\alpha \pm \beta i: y=eαt(C1cos⁡βt+C2sin⁡βt)y = e^{\alpha t}(C_1\cos\beta t + C_2\sin\beta t)

解き方の手順

  1. y′′y'' を λ2\lambda^2、y′y' を λ\lambda、yy を 1 に置きかえて、特性方程式を作る。
  2. 因数分解か解の公式で λ\lambda を求める(判別式が負なら −m2=mi\sqrt{-m^2} = mi)。
  3. 解の種類(異なる実数・重解・複素数)に合った一般解の形を書く。
  4. 初期条件があれば y(0)y(0) と y′(0)y'(0) の式を作り、C1C_1, C2C_2 を求める。

例

y′′+8y′+12y=0y'' + 8y' + 12y = 0、y(0)=3y(0) = 3、y′(0)=−10y'(0) = -10 を解いてください。

  1. 特性方程式: λ2+8λ+12=(λ+2)(λ+6)=0\lambda^2 + 8\lambda + 12 = (\lambda + 2)(\lambda + 6) = 0 より λ=−2,−6\lambda = -2, -6
  2. 一般解: y=C1e−2t+C2e−6ty = C_1e^{-2t} + C_2e^{-6t}
  3. y(0)=C1+C2=3y(0) = C_1 + C_2 = 3、y′(0)=−2C1−6C2=−10y'(0) = -2C_1 - 6C_2 = -10 を解いて C1=2C_1 = 2, C2=1C_2 = 1
  4. y=2e−2t+e−6ty = 2e^{-2t} + e^{-6t}(揺れずに 0 に戻る過減衰)

見る

計算例

ばねにつながれたおもりを、粘り気の強い液体の中で動かします。つり合いの位置からの変位 yy(cm)は、時刻 tt(秒)について y′′+8y′+15y=0y'' + 8y' + 15y = 0 をみたします(抵抗が強く、振動せずに戻る「過減衰」)。

問い最初の変位が y(0)=−6y(0) = -6、速度が y′(0)=26y'(0) = 26 のとき、yy を tt の式で求めてください。

途中の式

  1. y=eλty = e^{\lambda t} を代入すると、特性方程式 λ2+8λ+15=0\lambda ^{2} + 8\lambda + 15 = 0 が得られます。
  2. (λ+3)(λ+5)=0(\lambda+3)(\lambda+5) = 0 より λ=−3,−5\lambda = -3, -5
  3. 一般解: y=C1e−3t+C2e−5ty = C_1e^{-3t} + C_2e^{-5t}
  4. y(0)=C1+C2=−6y(0) = C_1 + C_2 = -6、y′(0)=−3C1−5C2=26y'(0) = -3C_1 - 5C_2 = 26 を解いて C1=−2C_1 = -2, C2=−4C_2 = -4
  5. y=−2e−3t−4e−5ty = -2e^{-3t} - 4e^{-5t}。どちらの項も時間とともに 0 に近づき、おもりはつり合いの位置に戻ります。

答えy(t の式。例: y = 7e^(-2t) - e^(-6t)。指数はかっこで囲む)は −2e−3t−4e−5t-2e^{-3t} - 4e^{-5t}

一緒にやる

小さく分けて、一緒に解く

1つずつ答えを確かめながら進みます。間違えても大丈夫です。

ばねにつながれたおもりを、粘り気の強い液体の中で動かします。つり合いの位置からの変位 yy(cm)は、時刻 tt(秒)について y′′+6y′+8y=0y'' + 6y' + 8y = 0 をみたします(抵抗が強く、振動せずに戻る「過減衰」)。

最初の変位が y(0)=2y(0) = 2、速度が y′(0)=−12y'(0) = -12 のとき、yy を tt の式で求めてください。

小さな問いに分けて、1つずつ一緒に進めます。

ステップ 1 / 3

特性方程式 λ2+6λ+8=0\lambda ^{2} + 6\lambda + 8 = 0 の解 λ\lambda をすべて求めてください。

解が複数あるときはカンマで区切ります(λ=2, -3)。「λ=」は省略できます。解がないときは「解なし」。

読み取り:(入力すると、ここに数式で表示されます)

自分でやる

練習問題

数字は毎回変わります。2〜3問続けて解けたら、次の単元や元のレッスンに進みましょう。

問題を準備しています…

学び直しマップ

前後の単元

分からないところがあったら前提まで戻り、分かったら使う先へ進みましょう。

学び直しマップで、この単元のまわりを見る