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

1階常微分方程式の練習問題

変数分離形・1階線形微分方程式を解き、成長・減衰のモデルに使える。

考え方

「変化の速さが今の量に比例する」のような関係(微分方程式)から、量を時間の式で表せるようになります。在庫の減り方、放射性物質の減衰、人口や預金の成長、薬の血中濃度のモデルです。

微分方程式
未知の関数 y(t)y(t) と、その導関数 dydt\dfrac{dy}{dt} の関係を表す式。
一般解・特殊解
任意の定数を含むすべての解(一般解)と、初期条件 y(0)=y0y(0) = y_0 などで定数を決めた解(特殊解)。
変数分離形
dydt=f(t)g(y)\dfrac{dy}{dt} = f(t)g(y) のように、yy の部分と tt の部分に分けられる形。

在庫が「残っている量に比例する速さで」売れていくなら、dydt=−ky\dfrac{dy}{dt} = -ky(k>0k > 0)と書けます。量が多いうちは速く減り、少なくなるとゆっくり減ります。このような「変化の速さ」の式から yy を求めるのが、微分方程式を解くことです。

変数分離形は、yy を含む部分を左辺、tt を含む部分を右辺に集めてから、両辺を積分します。1y dy=−k dt\dfrac{1}{y}\,dy = -k\,dt を積分すると、log⁡∣y∣=−kt+C\log|y| = -kt + C です。

指数に戻すと ∣y∣=e−kt+C=eCe−kt|y| = e^{-kt + C} = e^C e^{-kt} なので、y=Ae−kty = Ae^{-kt}(AA は定数)が一般解です。e−kt+Ce^{-kt + C} を e−kt+Ce^{-kt} + C としないよう注意します。最後に t=0t = 0 を代入すると y(0)=Ay(0) = A なので、AA は最初の量です。

1階線形微分方程式 dydt+P(t)y=Q(t)\dfrac{dy}{dt} + P(t)y = Q(t) は、両辺に e∫P(t) dte^{\int P(t)\,dt}(積分因子)をかけると左辺が1つの積の微分になり、積分して解けます。

y = 4e^(-0.5t) の減っていく曲線と、t = 0 と t = 3 での接線(点線)。量が多い t = 0 では接線が急で、量が少ない t = 3 では緩やか1234561234Oty量が多い → 速く減る量が少ない → ゆっくり減る
y=4e−0.5ty = 4e^{-0.5t}(dydt=−0.5y\dfrac{dy}{dt} = -0.5y の解)。点線は接線で、その傾きが減る速さです。量が多い t=0t = 0 では速く、量が少ない t=3t = 3 ではゆっくり減ります。

ポイント

dydt=ky\dfrac{dy}{dt} = ky の一般解: y=Aekty = Ae^{kt}(A=y(0)A = y(0))

変数分離: dydt=f(t)g(y)⇒∫dyg(y)=∫f(t) dt\dfrac{dy}{dt} = f(t)g(y) \Rightarrow \displaystyle\int\frac{dy}{g(y)} = \int f(t)\,dt

1階線形: dydt+P(t)y=Q(t)\dfrac{dy}{dt} + P(t)y = Q(t) は、両辺に e∫P(t) dte^{\int P(t)\,dt} をかけて積分する

解き方の手順

  1. yy の部分を左辺、tt の部分を右辺に分ける(1y dy=−k dt\dfrac{1}{y}\,dy = -k\,dt)。
  2. 両辺を積分し、積分定数 CC をつける(log⁡∣y∣=−kt+C\log|y| = -kt + C)。
  3. 指数の形に直して y=Ae−kty = Ae^{-kt} とする(ea+C=eCeae^{a + C} = e^Ce^a)。
  4. t=0t = 0 を代入して AA を決め、yy を tt の式で答える。

例

細菌の数 yy が dydt=0.2y\dfrac{dy}{dt} = 0.2y(tt は時間)に従い、最初は 500 個です。yy を tt の式で表してください。

  1. 変数分離: 1y dy=0.2 dt\dfrac{1}{y}\,dy = 0.2\,dt
  2. 両辺を積分: log⁡∣y∣=0.2t+C\log|y| = 0.2t + C
  3. y=Ae0.2ty = Ae^{0.2t}(A=±eCA = \pm e^C)
  4. y(0)=A=500y(0) = A = 500 より y=500e0.2ty = 500e^{0.2t}

見る

計算例

季節商品の在庫 yy(個)は、残っている在庫に比例する速さで売れていき、dydt=−0.1y\dfrac{dy}{dt} = -0.1y(tt は週)を満たすとします。最初(t=0t = 0)の在庫は 1100 個です。

最初の在庫
1,100個

問い微分方程式 dydt=−0.1y\dfrac{dy}{dt} = -0.1y、y(0)=1100y(0) = 1100 を解き、yy を tt の式で答えてください。

途中の式

  1. 変数分離: 1y dy=−0.1 dt\dfrac{1}{y}\,dy = -0.1\,dt
  2. 両辺を積分: log⁡∣y∣=−0.1t+C\log|y| = -0.1t + C
  3. ∣y∣=e−0.1t+C=eCe−0.1t|y| = e^{-0.1t + C} = e^{C}e^{-0.1t} なので y=Ae−0.1ty = Ae^{-0.1t}(A=±eCA = \pm e^{C})
  4. y(0)=A=1100y(0) = A = 1100 より y=1100e−0.1ty = 1100e^{-0.1t}
  5. 注意: e−0.1t+Ce^{-0.1t + C} は e−0.1t+Ce^{-0.1t} + C ではなく eCe−0.1te^{C}e^{-0.1t} です(指数法則 ep+q=epeqe^{p + q} = e^p e^q)。

答えy(t の式。例: y = 3e^(2t))は 1100e−0.1t1100e^{-0.1t}

一緒にやる

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

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

季節商品の在庫 yy(個)は、残っている在庫に比例する速さで売れていき、dydt=−0.4y\dfrac{dy}{dt} = -0.4y(tt は週)を満たすとします。最初(t=0t = 0)の在庫は 1600 個です。

微分方程式 dydt=−0.4y\dfrac{dy}{dt} = -0.4y、y(0)=1600y(0) = 1600 を解き、yy を tt の式で答えてください。

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

ステップ 1 / 3

変数を分離すると 1y dy=−0.4 dt\dfrac{1}{y}\,dy = -0.4\,dt。右辺 −0.4-0.4 を tt で積分すると?(積分定数 CC も付けて)

累乗は t^3、分数は 1/3 のように入力し、積分定数 C を忘れずに付けます。

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

自分でやる

練習問題

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

問題を準備しています…

学び直しマップ

前後の単元

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

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