微分方程式:ステップごとの解答つき練習問題

図、ヒント、詳しい解答つきの無料の微分方程式練習問題。毎回新しい数値で、3 つのレベル:指数関数的な増加と減少, ニュートンの冷却法則, 1 階線形微分方程式。

無料の練習問題やさしい
ty1

ばね・質量・ダンパ系が y″ + 4.8y′ + 4y = 0 で表されます。y″ + 2ζωny′ + ωn²y = 0 と比べて、減衰比 ζ を求めなさい。

y″ + 4.8y′ + 4y = 0
次の問題 →
微分方程式を無料で練習今すぐ微分方程式を練習毎日 5 問まで無料。カード登録は不要です。

練習する内容

  1. やさしい指数関数的な増加と減少
  2. ふつうニュートンの冷却法則
  3. むずかしい1 階線形微分方程式

重要な公式

指数関数的な増加・減少
dy/dt = ky → y = y0ekt
半減期
k = ln 2 / t½
ニュートンの冷却法則
T(t) = Ta + (T0 − Ta)e−kt
積分因子
y′ + p(x)y = q(x): μ = e∫p dx, (μy)′ = μq
混合タンク
dm/dt = r·c − (r/V)·m

微分方程式をゼロから学ぶ

事前に学ぶこと: 微分, 積分

1微分方程式とは

微分方程式は、ある量とその変化率を結び付けます。それが何であるかではなく、どのように変化するかを記述します。

細菌は分裂するので、大きな集団ほど速く増えます。増加率は細胞数 N に比例します:

dN/dt = k · N

ここで dN/dt は N の時間 t についての導関数、k は速度定数です。k = 0.4 h⁻¹、N = 1000 個のとき、集団はその瞬間に 1 時間あたり 400 個増えます。

方程式を解くとは、関数 N(t) そのものを求めることです。N(0) = 200 のような初期値も必要です。これを初期条件といいます。

セルフチェック ある個体数が dN/dt = 0.2 · N に従う。この瞬間 N = 500 である。どれだけの速さで増えているか。
  1. 1 時間あたり 0.2
  2. 1 時間あたり 100
  3. 1 時間あたり 2500
  4. 1 時間あたり 500

答え: 1 時間あたり 100. 数を方程式に代入します:dN/dt = 0.2 · 500 = 100。500 を 0.2 で割ると 2500 になりますが、これは k を逆に使っています。

2指数関数的な増加と減衰

方程式 dy/dt = k·y の解は、常に同じ指数関数です:

y(t) = y0 · ekt

ここで y0 は t = 0 での値です。k が正なら増加を意味します。減衰の方程式は dy/dt = −k·y で、解は y = y0e−kt です。

放射性崩壊は、半減期 t½、つまり物質が半分になるまでの時間で表されることがよくあります。速度定数はそこから求まります:

k = ln 2 / t½

時間を求めるには、両辺の自然対数をとります:

t = ln(y / y0) / k

たとえば、k = 0.4 h⁻¹ で増える 200 個の細胞は、t = ln(5)/0.4 ≈ 4.02 h 後に 1000 個になります。

セルフチェック ある物質の半減期は 10 日である。80 mg から始める。20 日後には何 mg 残るか。
  1. 0 mg
  2. 20 mg
  3. 40 mg
  4. 60 mg

答え: 20 mg. 20 日は半減期 2 回分です:80 → 40 → 20 mg。0 mg は、半分が 2 つで全体になると考えた結果です。量は毎回半分になります。

3冷却と線形方程式

ニュートンの冷却の法則:熱い物体は、室温よりずっと熱いときのほうが速く冷えます。室温を Ta とすると:

dT/dt = −k · (T − Ta) ⇒ T(t) = Ta + (T0 − Ta) · e−kt

差 T − Ta は指数関数的に減衰するので、温度は室温に向かい、室温を下回ることはありません。

1 回の測定から k を求めるには、既知の値を代入して ln をとります:k = −ln[(T1 − Ta)/(T0 − Ta)] / t1。

1 階線形方程式は y′ + p(x)·y = q(x) の形をしています。これに積分因子 μ を掛けます:

μ = e∫p dx, (μ·y)′ = μ·q

次に両辺を積分し、μ で割ります。p = 2/x のとき、μ = e2 ln x = x² です。

セルフチェック 90 °C の紅茶を、20 °C の部屋に置いておく。十分に長い時間がたつと、温度は何度に近づくか。
  1. 0 °C
  2. 20 °C
  3. 55 °C
  4. 90 °C

答え: 20 °C. 項 (T0 − Ta)e−kt は 0 に近づくので、T は Ta = 20 °C に近づきます。室温を忘れると 0 °C になります。

4ステップごとの方法

  1. 型を見分けます:dy/dt = ky(増加または減衰)、ニュートンの冷却、または y′ + p(x)y = q(x)。
  2. その型の標準的な解を書きます。
  3. 与えられたデータで定数を求めます:y0 は初期値から、k は半減期または 2 回目の測定から。
  4. 問われたことに答えます:時間を代入して値を求める、または ln を使って時間を求めます。
  5. k の単位(h⁻¹、min⁻¹)を t の単位と照らして確認します。

線形方程式では:μ を求め、(μy)′ = μq を積分し、定数 C を加え、初期条件で C を求めます。

5例題

すべてのステップをたどりましょう。開く前に次のステップを予想してみましょう。

例題

放射性トレーサーの半減期は 8 日である。試料は最初 60 mg を含む。崩壊は dm/dt = −k·m に従う。12 日後に残る量 m はいくらか。

1

型を見分ける

これは指数関数的な減衰なので、解は次のとおりです:

m(t) = m0 · e−kt, m0 = 60 mg
2

半減期から k を求める

k = ln 2 / 8 = 0.6931 / 8 = 0.08664 日⁻¹
3

時間を代入する

k · t = 0.08664 · 12 = 1.0397
m = 60 · e−1.0397 = 60 · 0.3536 = 21.21 mg
4

答えを確認する

12 日は半減期 1.5 回分です。8 日後には 30 mg、16 日後には 15 mg が残るので、21.2 mg はその間にあり、予想どおりです。

60 · 0.51.5 = 21.21 mg
tm608301615247.51221.2
答え

12 日後には m ≈ 21.2 mg が残ります。

6よくある間違い

  • ニュートンの冷却の法則で、室温 Ta を引き、また足し戻すのを忘れる。
  • 減衰に正の指数を使い、量が減らずに増えてしまう。
  • 自然対数 ln の代わりに log10 を使う。
  • 時間の単位を混同する。たとえば、min⁻¹ の k に、時間で表した t を使う。
  • 積分因子を使うとき、左辺にだけ μ を掛ける。
セルフチェック 80 °C の部品が、20 °C の部屋で冷える。ある時刻 t に、e−kt = 0.5 になる。そのときの温度はいくらか。
  1. 30 °C
  2. 40 °C
  3. 50 °C
  4. 60 °C

答え: 50 °C. T = 20 + (80 − 20) · 0.5 = 20 + 30 = 50 °C です。80 をそのまま半分にすると 40 °C になりますが、室温を忘れています。

7試験問題の例

  • やさしい:指数関数的な増加と減衰。ある時間後の量、または、ある量に達するまでの時間を、半減期から求めることが多くあります。
  • ふつう:ニュートンの冷却の法則。1 回の測定で k を求め、その後で温度または時間を求めます。
  • むずかしい:y(1) = 1 の y′ + 2y/x = 3 のような 1 階線形方程式。積分因子を使い、別の x での y を求めます。

簡単な確認:減衰と冷却の答えは、初期値と最終値(0、または室温)の間にあるはずです。

アプリでステップごとに学ぶアプリでステップごとに学ぶ

解答例

例題 1やさしい
rr² + br + c?

方程式 y″ + 4y′ + 3y = 0 について、特性方程式の根のうち 0 に最も近いもの r を求めなさい。それが解の減衰の遅さを決めます。

y″ + 4y′ + 3y = 0

求めるもの: r

ヒントを表示
  1. y = ert を代入します。
  2. r² + 4r + 3 = 0.
  3. 因数分解するか、解の公式を使います。
答えとステップごとの解答を表示

2 階線形方程式:y″ + b·y′ + c·y = 0 では y = ert とおきます。これから特性方程式 r² + b·r + c = 0 が得られます。異なる 2 つの実根 r₁ ≠ r₂ のときは y = A·er₁t + B·er₂t(過減衰)、重根 r のときは y = (A + Bt)·ert(臨界減衰)、複素根 −α ± jω のときは y = e−αt(A cos ωt + B sin ωt)(不足減衰、振動)となります。初期値 y(0) と y′(0) から A と B が決まります。y″ + 2ζωny′ + ωn²y = 0 の形に書くと、ωn は固有角周波数、ζ は減衰比で、ζ < 1 なら振動します。入力が一定の y″ + by′ + cy = c·K では、解は定常値 K に自由応答を加えたものです。

固有方程式

r² + 4r + 3 = (r + 1)(r + 3) = 0

根

r = −1, r = −3
0 に最も近い: r = −1

答え

r = −1

一般解は y = A·e−1t + B·e−3t です。

例題 2ふつう
t (min)T (°C)5Ta​73 °C50 °C

81 °C のコーヒーを 21 °C の部屋に置いておきます。 5 分後の温度は 73 °C です。 温度はニュートンの冷却法則 dT/dt = −k·(T − Tₐ) に従います。 開始から 50 °C に達するまでの時間 t を求めなさい。

T0 = 81 °CTa = 21 °CT(5) = 73 °CT = 50 °C

求めるもの: t (min)

ヒントを表示
  1. まず t = 5 min での測定から k を求めます。
  2. k = ln(60/52)/5 = 0.02862 min⁻¹.
  3. t = ln((T₀ − Tₐ)/(T − Tₐ))/k.
答えとステップごとの解答を表示

1 階微分方程式:dy/dt = k·y の解は y = y₀·ekt です(k > 0 なら増加、k < 0 なら減少)。半減期 T からは k = ln 2/T が得られます。ニュートンの冷却法則 dT/dt = −k·(T − Tₐ) からは T(t) = Tₐ + (T₀ − Tₐ)·e−kt が得られます。線形方程式 y′ + p(x)·y = q(x) は積分因子 μ = e∫p dx を使って解きます:(μy)′ = μq なので、μy = ∫μq dx + C です。

冷却定数

k = ln(60/52)/5 = 0.02862 min⁻¹

目標温度に達するまでの時間

t = ln(60/29)/0.02862 = 25.4 min

答え

t = 25.4 min
例題 3むずかしい
ty1y(1) = ?

y′ + y = 4e−2t、y(0) = 4 を解きなさい。y(1) を求めなさい。

y′ + y = 4e−2ty(0) = 4

求めるもの: y(1)

ヒントを表示
  1. 積分因子は μ = et です。
  2. (ety)′ = 4e−t → ety = 4/(−1)·e−t + C.
  3. y = −4·e−2t + C·e−t. y(0) を使って C を求めます。
答えとステップごとの解答を表示

1 階微分方程式:dy/dt = k·y の解は y = y₀·ekt です(k > 0 なら増加、k < 0 なら減少)。半減期 T からは k = ln 2/T が得られます。ニュートンの冷却法則 dT/dt = −k·(T − Tₐ) からは T(t) = Tₐ + (T₀ − Tₐ)·e−kt が得られます。線形方程式 y′ + p(x)·y = q(x) は積分因子 μ = e∫p dx を使って解きます:(μy)′ = μq なので、μy = ∫μq dx + C です。

積分因子

μ = et
(ety)′ = 4e−t

一般解

y = 4/(−1)·e−2t + C·e−t = −4·e−2t + C·e−t

初期条件

4 = −4 + C → C = 8

計算する

y(1) = −4·0.1353 + 8·0.3679 = 2.402

答え

y(1) = 2.402

質問と回答

微分方程式とは?

細菌は分裂するので、大きな集団ほど速く増えます。増加率は細胞数 N に比例します:

ここで dN/dt は N の時間 t についての導関数、k は速度定数です。k = 0.4 h⁻¹、N = 1000 個のとき、集団はその瞬間に 1 時間あたり 400 個増えます。

微分方程式の問題の解き方は?(手順)

2 階線形方程式:y″ + b·y′ + c·y = 0 では y = ert とおきます。これから特性方程式 r² + b·r + c = 0 が得られます。異なる 2 つの実根 r₁ ≠ r₂ のときは y = A·er₁t + B·er₂t(過減衰)、重根 r のときは y = (A + Bt)·ert(臨界減衰)、複素根 −α ± jω のときは y = e−αt(A cos ωt + B sin ωt)(不足減衰、振動)となります。初期値 y(0) と y′(0) から A と B が決まります。y″ + 2ζωny′ + ωn²y = 0 の形に書くと、ωn は固有角周波数、ζ は減衰比で、ζ < 1 なら振動します。入力が一定の y″ + by′ + cy = c·K では、解は定常値 K に自由応答を加えたものです。

微分方程式でよくある間違いは?
  • それはもう一方の根で、0 からより遠く、その項はより速く減衰します。
  • b と c が正の r² + br + c = 0 の根は負です。
  • −b/2 は根の中点であり、根ではありません。
微分方程式の練習では何を学ぶ?
  • やさしい: 指数関数的な増加と減少
  • ふつう: ニュートンの冷却法則
  • むずかしい: 1 階線形微分方程式

問題ごとに数値が変わるので、答えではなく解き方が身につきます。 行き詰まったら、まず小さなヒントを開き、必要なときだけ詳しい解説を見ましょう。

練習のしかた

  • 問題ごとに数値が変わるので、答えではなく解き方が身につきます。
  • 行き詰まったら、まず小さなヒントを開き、必要なときだけ詳しい解説を見ましょう。
  • 進捗、間違えた問題、保存した問題は記録されるので、試験前に復習できます。

このトピックの前後のトピック

先に学ぶ

次に学ぶ

計算ツール

このトピックを含む科目

数学・統計学のほかのトピック

印刷用ワークシート

微分方程式を無料で練習今すぐ微分方程式を練習毎日 5 問まで無料。カード登録は不要です。