応用数学 第4回 (4) 身近なところにある2階線形微分方程式
RLC 回路 ( 教科書 p.70, 例 4.4 )
2階線形微分方程式は身近なところにある、というお話で、まずは電気回路から。
設定 次の4つを直列につないだ電気回路を考えます:
- 抵抗値 R の抵抗
- インダクタンス L のコイル
- 容量 C のコンデンサ
- 起電力 f(t) の電源 ( t は時刻 )
ただし、
R,
L,
C は定数で、
f(t) は時間変化するとします。
問 この回路に流れる電流量 i(t) の満たす微分方程式を求めよ。
解 ポイントは、各素子がもたらす電圧への影響です。
- 抵抗値 R の抵抗は、電流量 i(t) に比例する R×i(t) の電圧降下をもたらします。
- インダクタンス L のコイルは、電流量の変化に比例する L×i′(t) の電圧降下をもたらします。
- 容量 C のコンデンサは、コンデンサに蓄えられる電荷量 q(t) に比例する C×q(t) の電圧降下をもたらします。
以上から
Cq(t)+Ri(t)+Li′(t)=f(t) ⋯⋯ (1)
という式が得られます。
更に、電荷量
q(t) の変化が電流量になりますので、
i(t)=q′(t).
従って
(1) を
t で微分することにより
Li″+Ri′+Ci=f′
これが
i=i(t) の満たす微分方程式で、
定数係数2階線形微分方程式という形です。
鉛直ばね振り子 ( 教科書 p.26, 例 8.4 )
次はばねの運動です。
設定 天井から吊り下げられたばねに、空気抵抗の大きな物体が吊り下げられているとします。
物体に掛かる力は次の3つです:
- 重力
- ばねの力
- 空気抵抗
ばねに何も吊るさないときの先端の位置を
x=0 とし、この物体の時刻
t での位置を
x(t) とします。
問 x(t) の満たす微分方程式を求めよ。
解
- 重力は、物体の質量を m, 重力加速度を g として −mg です。
- ばねの力は、何も吊るさないときの位置からの変位に比例します。
比例定数を k として −kx(t) と表されます。
- 空気抵抗は速度 x′(t) に比例します。比例定数を R として −Rx′(t) と表されます。
以上を、ニュートンの運動方程式
F=ma=mx″(t)
へ入れると
−mg−kx−Rx′=mx″
すなわち
x″+Rmx′+kmx=−g
やはり定数係数2階線形微分方程式になりました。
なおこれは、
X=x+mgk とおけば
X″+RmX′+kmX=0
となり、定数係数2階線形同次微分方程式に変換できます。