物理・工学・数学を学んでいると必ず登場する「線形微分方程式」。
「1階と2階で何が違うの？」「同次と非同次の解法の違いは？」「線形代数とどう関係するの？」と疑問を感じる方も多いでしょう。
本記事では、線形微分方程式の意味・1階・2階の解法・同次・非同次の違いから、線形代数との関係まで、わかりやすく解説していきます。
具体的な解法ステップと例を交えながら丁寧に説明しますので、ぜひ最後までご覧ください。
線形微分方程式とは？基本的な意味と分類
それではまず、線形微分方程式の基本的な意味と分類について解説していきます。
線形微分方程式とは、未知関数とその導関数が1次式（線形）の形で表された微分方程式のことです。
一般的な n 階線形微分方程式は次の形をしています。
【n階線形微分方程式の一般形】
aₙ(x)yⁿ + aₙ₋₁(x)y⁽ⁿ⁻¹⁾ + … + a₁(x)y’ + a₀(x)y = f(x)
左辺の係数 aᵢ(x) はxの関数（定数でも可）
f(x) = 0 → 同次方程式
f(x) ≠ 0 → 非同次方程式
線形微分方程式の最大の特徴は、「重ね合わせの原理」が成立することです。
y₁ と y₂ がともに同次方程式の解であれば、c₁y₁ + c₂y₂ も解になります。
線形微分方程式の分類
| 分類 | 条件 | 例 |
|---|---|---|
| 1階同次 | n=1, f(x)=0 | y’ + 2y = 0 |
| 1階非同次 | n=1, f(x)≠0 | y’ + 2y = x |
| 2階同次（定係数） | n=2, f=0, 定係数 | y” + 3y’ + 2y = 0 |
| 2階非同次（定係数） | n=2, f≠0, 定係数 | y” + 3y’ + 2y = eˣ |
| 変数係数 | 係数がxの関数 | xy” + y = 0 |
実用上最も重要なのは「定係数2階線形微分方程式」であり、バネ・マス・ダンパー系・電気回路などの物理問題で頻繁に登場します。
1階線形同次微分方程式の解き方
1階線形同次微分方程式 y’ + p(x)y = 0 は変数分離法で解けます。
【1階同次の解法】
dy/y = -p(x)dx と変数分離
両辺を積分して ln|y| = -∫p(x)dx + C
→ y = Ae^{-∫p(x)dx}（A は任意定数）
たとえば y’ + 2y = 0 の解は y = Ae^{-2x} となります。
1階非同次線形微分方程式の解き方
1階線形非同次方程式 y’ + p(x)y = q(x) は「積分因子法」で解けます。
【積分因子法の手順】
①積分因子 μ(x) = e^{∫p(x)dx} を計算
②両辺に μ を掛けると左辺が (μy)’ の形になる
③両辺を積分して y を求める
積分因子法は1階非同次方程式の標準的な解法として広く使われます。
2階定係数線形微分方程式の解法
続いては、最も実用的な2階定係数線形微分方程式の解法を確認していきます。
2階定係数同次方程式 ay” + by’ + cy = 0 の解は「特性方程式」を使って求めます。
特性方程式と解の形
y = e^{λx} とおいて代入すると、特性方程式 aλ² + bλ + c = 0 が得られます。
特性方程式の解によって、一般解の形が次のように変わります。
| 特性根の種類 | 一般解の形 |
|---|---|
| 異なる2実根 λ₁, λ₂ | y = C₁e^{λ₁x} + C₂e^{λ₂x} |
| 重根 λ（重複） | y = (C₁ + C₂x)e^{λx} |
| 複素根 α ± βi | y = e^{αx}(C₁cosβx + C₂sinβx) |
複素根が出る場合は振動する解（正弦・余弦）が現れるため、バネの振動や電気回路のLC振動と対応します。
非同次方程式の解法：未定係数法
非同次方程式 ay” + by’ + cy = f(x) の解は、同次方程式の一般解に「特解」を加えたものです。
特解を求める代表的な方法が「未定係数法」です。
f(x) の形に応じて特解の形を仮定し、方程式に代入してパラメータを決定します。
【未定係数法の特解の仮定形】
f(x) = eˣ → y_p = Aeˣ（λ=1が特性根でない場合）
f(x) = sin(x) → y_p = Acos(x) + Bsin(x)
f(x) = x² → y_p = Ax² + Bx + C
線形代数との関係
線形微分方程式と線形代数は深い関係を持ちます。
微分演算子 D を使うと、微分方程式を「行列の演算」のように扱えます。
線形微分方程式の解空間は線形空間を形成し、その次元は方程式の階数に等しくなります。
行列の固有値・固有ベクトルは線形微分方程式系の解の構造を理解するために直接使われます。
線形代数の「固有値・固有ベクトル」の理論が線形微分方程式系を完全に解く鍵となっています。
線形微分方程式の応用
続いては、線形微分方程式の実際の応用例を確認していきます。
バネ・マス・ダンパー系（力学）
物理の古典的な例として、バネ・マス・ダンパー系の運動方程式があります。
【バネ・マス系の運動方程式】
m × x” + c × x’ + k × x = F(t)
m：質量、c：減衰係数、k：バネ定数
F(t) = 0 → 自由振動（同次方程式）
F(t) ≠ 0 → 強制振動（非同次方程式）
この方程式は2階定係数非同次線形微分方程式であり、特性方程式を使って解を求めます。
橋・建物・機械部品の振動解析でこの方程式が活躍しています。
電気回路（LRC回路）への応用
電気回路のLRC回路（インダクタ・抵抗・コンデンサ）の電流・電圧の変化も2階線形微分方程式で記述できます。
L × q” + R × q’ + (1/C) × q = V(t) という形をとり、バネ・マス系と数学的に同じ構造を持ちます。
このような「アナロジー（類推）」は異なる物理現象が同じ数学的構造で記述できるという線形微分方程式の普遍性を示しています。
人口動態・感染症モデルへの応用
生物学・疫学でも線形微分方程式が使われます。
単純な人口増加モデル dN/dt = rN（指数成長）は1階線形同次方程式です。
SIR感染症モデルは非線形ですが、線形化近似によって線形微分方程式系として解析することもあります。
線形微分方程式は自然科学・社会科学の広い分野に応用される普遍的なツールと言えるでしょう。
まとめ
本記事では、線形微分方程式の意味・分類・1階同次・非同次の解法・2階定係数の解法・線形代数との関係・実際の応用まで詳しく解説しました。
線形微分方程式とは未知関数とその導関数が線形の形で表された方程式で、重ね合わせの原理が成立する美しい構造を持ちます。
特性方程式・積分因子法・未定係数法などの解法を習得することで、多くの物理・工学問題を解く力が身につきます。
線形微分方程式の理解は数学・物理・工学・データサイエンスの深い学習に不可欠な基礎となるでしょう。
ぜひ具体的な問題を繰り返し解いて、解法の感覚を身につけてみてください。
