物理学の中でも「振動」は非常に広い応用範囲を持つテーマであり、力学・電気・量子力学など多岐にわたる分野でその理解が求められます。
「振動数の方程式とは何か」「どのようにして導出されるのか」と疑問を持つ方も多いでしょう。
振動数を求める方程式は、微分方程式による運動方程式の解析を通じて導出されます。
単振動・調和振動子・減衰振動など、さまざまな振動モデルに応じた方程式の形を理解することが、物理学の深い理解への近道です。
本記事では、振動数の方程式の導出プロセスを段階的に丁寧に解説し、具体的な解法も合わせてお伝えします。
初学者の方でも理解できるよう、数学的な準備から丁寧に説明しますので、ぜひ最後までご覧ください。
目次
振動数の方程式は単振動の運動方程式から導出される
それではまず、振動数の方程式の本質とその導出方法について解説していきます。
結論として、振動数の方程式は単振動（調和振動）の運動方程式を解くことで導かれます。
その出発点となるのが、ニュートンの第二法則とフックの法則を組み合わせた微分方程式です。
単振動の運動方程式（最も基本的な振動数の方程式）
m(d²x/dt²) = -kx
整理すると：d²x/dt² = -(k/m)x = -ω²x
この方程式の解から振動数が求まります：ω = √(k/m)、f = ω/(2π) = (1/2π)√(k/m)
微分方程式としての振動方程式の意味
単振動の運動方程式 d²x/dt² + ω²x = 0 は、2階線形常微分方程式の代表的な形です。
この式の意味を言葉で説明すると、「位置xの加速度（位置の2階微分）が、位置xに比例した大きさで逆向きに働く」ということです。
つまり、物体が平衡点から遠ざかるほど、強い復元力が働いて平衡点へ引き戻そうとする運動を記述しています。
この「復元力が変位に比例する」という条件こそが、単振動の定義そのものです。
バネに取り付けられた質量ブロックがこの条件を満たす典型例であり、「調和振動子（ハーモニック・オシレーター）」とも呼ばれます。
微分方程式の解は x(t) = A cos(ωt + φ) という形を取り、Aは振幅、φは初期位相を表します。
一般解の導出過程を丁寧に追う
微分方程式 d²x/dt² + ω²x = 0 の一般解を導出する手順を確認しましょう。
この方程式の解として x = e^(λt) を仮定すると、代入により特性方程式 λ² + ω² = 0 が得られます。
これを解くと λ = ±iω（iは虚数単位）となります。
オイラーの公式 e^(iθ) = cosθ + isinθ を用いることで、一般解は以下のように表されます。
一般解の形
x(t) = C₁e^(iωt) + C₂e^(-iωt)
= A cos(ωt) + B sin(ωt) （A、Bは実数の任意定数）
= R cos(ωt + φ) （R：振幅、φ：初期位相）
ここで R = √(A²+B²)、tan(φ) = -B/A
この一般解が示す通り、単振動は角振動数ω、振幅R、初期位相φの3つのパラメータで完全に記述されます。
振動数fはf=ω/(2π)で求まり、周期TはT=2π/ωです。
初期条件を使った特殊解の求め方
物理問題では「t=0のときの位置と速度」が初期条件として与えられることが多くあります。
一般解 x(t) = A cos(ωt) + B sin(ωt) に初期条件を代入することで、A、Bの値が決まります。
【例題】t=0でx₀=0.1m、v₀=0 m/sの初期条件でばね振動子（ω=10rad/s）を解く
一般解：x(t) = A cos(10t) + B sin(10t)
初期条件①：x(0) = A = 0.1 → A = 0.1
速度：dx/dt = -10A sin(10t) + 10B cos(10t)
初期条件②：v(0) = 10B = 0 → B = 0
特殊解：x(t) = 0.1 cos(10t) （単位：m）
振動数：f = 10/(2π) ≈ 1.59Hz
初期条件を代入することで任意定数が決まり、一意の解（特殊解）が得られます。
物理的には「初期位置と初速度が運動を完全に決定する」ということを意味しているのです。
ばね-質量系・振り子・電気回路それぞれの振動数方程式
続いては、代表的な振動システムそれぞれの振動数方程式を確認していきます。
異なる物理系でも、同じ形の微分方程式に帰着することが「振動の普遍性」を示しています。
ばね-質量系（調和振動子）の振動数方程式
ばね定数k（N/m）のばねと質量m（kg）の物体からなる系の運動方程式は、フックの法則から導かれます。
ばねの復元力はF=-kxであり、ニュートンの第二法則ma=Fより以下の方程式が得られます。
ばね-質量系の振動数方程式
m(d²x/dt²) + kx = 0
→ d²x/dt² + (k/m)x = 0
角振動数：ω = √(k/m) rad/s
振動数：f = (1/2π)√(k/m) Hz
周期：T = 2π√(m/k) 秒
質量mが大きいほど振動が遅く（fが小さく）、ばね定数kが大きいほど振動が速い（fが大きい）ことが公式から直接読み取れます。
この直感は「硬いバネほど速く振動し、重い物体ほどゆっくり振動する」という日常的な感覚とも一致するでしょう。
単振り子の近似振動数方程式
長さLの振り子（小角近似：sinθ≈θ）の運動方程式は次のように表されます。
単振り子の振動数方程式（小角近似）
d²θ/dt² + (g/L)θ = 0
角振動数：ω = √(g/L) rad/s
振動数：f = (1/2π)√(g/L) Hz
周期：T = 2π√(L/g) 秒
注：この近似は振れ角が小さい（約5度以下）場合に成立します
単振り子の場合、振動数は振り子の長さLと重力加速度gのみに依存し、質量や振幅には依存しません。
これが「等時性」と呼ばれる性質であり、ガリレオが発見したとされる重要な物理現象です。
重力加速度gが場所によって異なるため、高精度な時計では振り子の長さを調整することで周期を正確に1秒に合わせる必要があります。
LC回路の電気振動方程式——力学との数学的類似
インダクタンスLのコイルとキャパシタンスCのコンデンサからなるLC回路の方程式は、力学的振動と数学的に全く同じ形をしています。
LC回路の振動数方程式
L(d²q/dt²) + (1/C)q = 0
→ d²q/dt² + (1/LC)q = 0
角振動数：ω = 1/√(LC) rad/s
振動数：f = 1/(2π√(LC)) Hz
周期：T = 2π√(LC) 秒
（qは電荷量、Lはインダクタンス（H）、Cはキャパシタンス（F））
力学系のmがLに、kがC⁻¹に対応しており、数学的構造が完全に一致しています。
このアナロジー（類比）は、異なる物理系を統一的に理解するうえで非常に強力な視点です。
ラジオの選局やフィルター回路の設計において、このLC回路の振動数方程式は今も現役で使われています。
減衰振動・強制振動における振動数方程式の拡張
続いては、より実際の物理現象に近い「減衰振動」と「強制振動」における振動数方程式について確認していきます。
実際の振動では摩擦や抵抗によってエネルギーが失われるため、単純な単振動よりも複雑な挙動を示します。
減衰振動の方程式と振動数への影響
空気抵抗や摩擦などの減衰力（速度に比例するとして-bv=-b(dx/dt)）を考慮した運動方程式は以下のようになります。
減衰振動の方程式
m(d²x/dt²) + b(dx/dt) + kx = 0
両辺をmで割ると：d²x/dt² + 2γ(dx/dt) + ω₀²x = 0
ここで γ = b/(2m)（減衰定数）、ω₀ = √(k/m)（固有角振動数）
減衰振動の角振動数：ω_d = √(ω₀² – γ²)
（γ
減衰が存在すると、振動数は固有振動数ω₀よりも小さくなります（ω_d
また、振幅は時間とともに指数関数的に減少していきます。
減衰が強すぎる場合（γ ≥ ω₀）は振動せずに指数関数的に平衡点へ近づく「過減衰」または「臨界減衰」となります。
強制振動と共振現象における振動数の重要性
外部から周期的な力が加えられる「強制振動」では、駆動振動数が系の固有振動数と一致したとき「共振（共鳴）」が起こります。
強制振動の方程式
m(d²x/dt²) + b(dx/dt) + kx = F₀cos(Ωt)
F₀：外力の振幅、Ω：外力の角振動数
共振条件：Ω ≈ ω₀ = √(k/m)
共振時の振幅は最大となり、減衰が弱いほど振幅が著しく大きくなります
共振現象は身の回りにも多く存在し、橋梁の崩壊事故（タコマナローズ橋）や楽器の音響設計、MRI装置の動作原理など、幅広い場面に関わっています。
エンジニアリングの現場では、構造物の固有振動数を正確に把握し、外部からの励振振動数と一致しないように設計することが安全設計の基本となります。
振動数方程式の解法まとめと活用場面
これまで見てきた各種振動数方程式の解法を整理しておきましょう。
| 振動の種類 | 方程式の形 | 振動数の式 | 主な応用例 |
|---|---|---|---|
| 単振動（非減衰） | d²x/dt² + ω²x = 0 | f = ω/(2π) | ばね、振り子、LC回路 |
| 減衰振動 | d²x/dt² + 2γẋ + ω₀²x = 0 | f_d = √(ω₀²-γ²)/(2π) | ショックアブソーバー |
| 強制振動 | d²x/dt² + 2γẋ + ω₀²x = F₀cos(Ωt)/m | 共振：Ω=ω₀ | 橋梁設計、音響 |
| 単振り子 | d²θ/dt² + (g/L)θ = 0 | f = (1/2π)√(g/L) | 時計、地震計 |
| LC回路 | d²q/dt² + q/(LC) = 0 | f = 1/(2π√(LC)) | ラジオ、フィルター |
各振動モデルに応じた方程式を正しく立て、解法（特性方程式・初期条件の代入）を実行することで振動数を求めることができます。
まず「どのタイプの振動か」を正確に把握することが、正しい方程式を選ぶための第一歩です。
振動数方程式の解法における数学的テクニック
続いては、振動数方程式を解くために必要な数学的テクニックについて確認していきます。
微分方程式の解法は複数ありますが、振動問題では特に有効な手法があります。
特性方程式を利用した解法
2階線形常微分方程式 d²x/dt² + px’/dt + qx = 0 の一般解を求めるために、「特性方程式」を利用する方法が最も標準的です。
x = e^(λt)を仮定して代入すると、λ²+pλ+q=0という代数方程式（特性方程式）が得られます。
判別式D=p²-4qの符号によって、解の形が異なります。
判別式による場合分け
D > 0：実数の異なる2解 → 非振動（過減衰）
D = 0：重解 → 臨界減衰
D
単振動（p=0）の場合：λ=±iω → D=-4ω²
振動が生じる条件はD
このとき固有値の虚部が角振動数ωに対応し、振動数f=ω/(2π)が決まります。
複素指数関数とオイラーの公式の活用
振動の数学的解析において、オイラーの公式 e^(iθ) = cosθ + isinθ は非常に重要な役割を担います。
複素指数関数を使うと、振動の計算が大幅に簡略化できます。
特に電気回路のインピーダンス計算やフーリエ解析では、複素数表現が標準的に用いられます。
複素表現の利点
位置 x(t) = Re[Ae^(iωt)]（Aは複素振幅）
微分：dx/dt = Re[iωAe^(iωt)]
これにより、微分演算が「iωをかける」という代数演算に変換されます
インピーダンス：Z = R + i(ωL – 1/(ωC))
複素表現を使いこなすことで、振動・波動の問題における計算効率が飛躍的に向上します。
大学物理・工学の演習では積極的に活用する価値があるテクニックです。
エネルギー的アプローチによる振動数の導出
微分方程式を解く以外にも、エネルギー保存の観点から振動数を導出することができます。
単振動では、運動エネルギーと弾性ポテンシャルエネルギーの和が一定に保たれます。
エネルギーによる振動数導出（ばね-質量系）
全エネルギー：E = (1/2)mv² + (1/2)kx² = 一定
x(t) = A cos(ωt) を代入すると
E = (1/2)mω²A²sin²(ωt) + (1/2)kA²cos²(ωt)
Eが時間によらず一定となるためには mω² = k が必要
→ ω = √(k/m)（微分方程式の解と一致）
エネルギー的アプローチは、複雑な系の固有振動数を見積もる際にも有効です。
「運動エネルギーの最大値=ポテンシャルエネルギーの最大値」という条件から振動数が決まるという直感は、様々な振動問題に応用できるでしょう。
まとめ
本記事では、振動数の方程式の導出と解法について、微分方程式の基礎から調和振動子、単振り子、LC回路、そして減衰振動・強制振動まで幅広く解説してきました。
振動数の方程式の出発点は、ニュートンの運動方程式にフックの法則（または類似の復元力）を代入した2階線形常微分方程式です。
特性方程式を解き、複素固有値の虚部から角振動数ω、さらに振動数f=ω/(2π)が求まります。
ばね-質量系、単振り子、LC回路など異なる物理系でも、同じ数学的構造の方程式に帰着することが振動現象の普遍性を示しています。
減衰や外部駆動力が加わると方程式は複雑になりますが、基本的な解法の流れは変わりません。
共振現象は工学設計において特に重要であり、固有振動数の正確な把握が安全な構造物や機器の設計に直結します。
ぜひ本記事を参考に、振動数方程式の導出と解法を自分の手で追いかけてみてください。
