Skip to content
Published on

工業数学シリーズ 第2回:1次微分方程式の概念

Authors

工業数学シリーズ 第2回:1次微分方程式の概念

1次微分方程式は工業数学の出発点です。名前が馴染みないように見えますが、実は「ある量の変化速度がわかっているとき、元の量の動きを復元する問題」と考えればずっと身近になります。

問題意識

現実のシステムでは、ほとんどの場合、静的な値よりも変化の過程が重要です。現在の温度、現在の電圧、現在の速度よりも「どれだけ速く変わっているか」がシステムの性質を決めることが多いです。微分方程式はまさにこの変化の規則を数式で書き留めたものです。

1次微分方程式とは、最高次の微分が1回だけ現れる式を意味します。一般的に次のように書きます。

F(x,y,y)=0F(x, y, y') = 0

またはより直接的に

dydx=f(x,y)\frac{dy}{dx} = f(x, y)

ここでxxは独立変数、yyは未知関数、yy'はその変化率です。

解を求めるとは何か

1次微分方程式の解を求めるとは、式を満たす関数y(x)y(x)を見つけるという意味です。

例えば

dydx=2x\frac{dy}{dx} = 2x

を満たす関数は

y=x2+Cy = x^2 + C

です。微分してみると確かに2x2xになるからです。つまり微分方程式の解は一つの数ではなく、通常は関数全体です。

なぜ積分定数が生まれるのか

微分は情報の一部を消します。x2x^2を微分しても2x2xx2+5x^2 + 5を微分しても2x2xです。したがって微分方程式を逆に解く過程では、失われた情報を定数CCで復元する必要があります。

このため一般解は通常

y=何らかの式+Cy = \text{何らかの式} + C

の形になります。そして初期条件や境界条件が与えられると、その時CCが決まります。

初期値問題の意味

現実の問題では、現在の状態が一緒に与えられることが多いです。例えば

dydx=2x,y(0)=3\frac{dy}{dx} = 2x, \quad y(0) = 3

であれば、一般解y=x2+Cy = x^2 + Cy(0)=3y(0) = 3を代入してC=3C = 3を得て、

y=x2+3y = x^2 + 3

という特殊解を得ます。

このように微分方程式と初期条件が一緒に与えられる問題を初期値問題と呼びます。工学では、回路を入れたときの初期電圧がいくらか、システム開始時の位置と速度がいくらかなどの情報がこの役割を果たします。

傾き場による直観

dydx=f(x,y)\frac{dy}{dx} = f(x, y)

は座標平面の各点での傾きを教えてくれます。この情報を矢印や小さな線分で描くと傾き場になり、解はその傾き場に沿って進む曲線になります。

例えば

dydx=y\frac{dy}{dx} = y

では、yyが大きいところほど傾きが急で、y=0y=0付近ではほぼ平らです。そのため解は指数関数の形になります。

この直観は、後に正確な解を求めるのが難しい問題でも非常に重要です。手で完全に解けなくても、解が増加しているのか減少しているのか、どの値に引き込まれているのかくらいは読み取れるようにならなければなりません。

代表例題

次の問題を見てみましょう。

dydt=0.5y,y(0)=10\frac{dy}{dt} = -0.5y, \quad y(0) = 10

この式は現在の量の半分の割合で減少する現象を表しています。一般解は

y(t)=Ce0.5ty(t) = Ce^{-0.5t}

であり、初期条件を代入すると

10=Ce0=C10 = Ce^0 = C

なので

y(t)=10e0.5ty(t) = 10e^{-0.5t}

です。

この解は、時間が経つにつれて値が0に近づく減衰現象を示しています。

工学応用

冷却問題

ある物体の温度が周囲温度に近づく過程は、しばしば

dTdt=k(TTs)\frac{dT}{dt} = -k(T - T_s)

のようにモデリングします。現在の温度と周囲温度の差が大きいほど速く冷えるという意味です。

簡単なRC回路

コンデンサの電圧変化も1次微分方程式から始まります。回路解析で2次、3次システムに進む前に最初に見る典型的な例です。

サービストラフィック減衰モデル

開発者の観点では、キャッシュTTL後に減少するトラフィック、徐々に安定するキュー長、簡単なフィードバックシステムも似た構造で見ることができます。

よくある間違い

式と解を区別できない

微分方程式は関数が満たすべき規則であり、解はその規則を満たす関数です。両者は同じものではありません。

初期条件を入れるのが遅すぎる

一般解を求めた後に特殊解を得るために初期条件を入れるという流れを忘れると、計算がよく混乱します。

変化率と関数値を混同する

yyyy'は全く異なる情報です。値が大きいからといって変化率が必ず大きいわけではなく、逆に値が小さくても変化率は大きいことがあります。

一行まとめ

1次微分方程式は「変化率の規則」から「関数の時間に沿った動き」を復元する問題です。

次回予告

次の記事では、1次微分方程式を実際に解く代表的な方法である変数分離法、積分因子法、完全微分方程式を順番に整理します。

参考資料

  • Erwin Kreyszig, Advanced Engineering Mathematics, 10th Edition
  • Dennis G. Zill, A First Course in Differential Equations
  • Paul Dawkins, Differential Equations Notes