RL直列回路において、以下の小問に答える。 1) $LI' + RI = 0$ の一般解を求める。 2) $V(t) = V_0$ の場合の特殊解 $I_p$ を求め、方程式の一般解を求める。初期条件 $I(0) = 0$ を満たす解を求め、電流 $I(t)$ が時間に対してどのように変化するかを説明する。 3) $V(t) = V_0 \sin(\omega t)$ の場合の特殊解 $I_p$ を未定係数法で求める。$I_p(t) = A\sin(\omega t + \delta)$ とおいて、$A$ と $\delta$ の関係式を導く。 4) 前問の関係式において、特に $t=0$, $t=\frac{\pi}{2\omega}$ のときを考え、$\delta = -\tan^{-1}\left(\frac{\omega L}{R}\right)$、$A = \frac{V_0}{\sqrt{(\omega L)^2 + R^2}}$ となることを示す。

応用数学微分方程式RL直列回路電気回路初期条件一般解特殊解未定係数法過渡現象定常状態
2025/5/9

1. 問題の内容

RL直列回路において、以下の小問に答える。
1) LI+RI=0LI' + RI = 0 の一般解を求める。
2) V(t)=V0V(t) = V_0 の場合の特殊解 IpI_p を求め、方程式の一般解を求める。初期条件 I(0)=0I(0) = 0 を満たす解を求め、電流 I(t)I(t) が時間に対してどのように変化するかを説明する。
3) V(t)=V0sin(ωt)V(t) = V_0 \sin(\omega t) の場合の特殊解 IpI_p を未定係数法で求める。Ip(t)=Asin(ωt+δ)I_p(t) = A\sin(\omega t + \delta) とおいて、AAδ\delta の関係式を導く。
4) 前問の関係式において、特に t=0t=0, t=π2ωt=\frac{\pi}{2\omega} のときを考え、δ=tan1(ωLR)\delta = -\tan^{-1}\left(\frac{\omega L}{R}\right)A=V0(ωL)2+R2A = \frac{V_0}{\sqrt{(\omega L)^2 + R^2}} となることを示す。

2. 解き方の手順

1) 斉次方程式 LI+RI=0LI' + RI = 0 を解く。
LI=RILI' = -RI
dII=RLdt\frac{dI}{I} = -\frac{R}{L} dt
dII=RLdt\int \frac{dI}{I} = \int -\frac{R}{L} dt
lnI=RLt+C\ln |I| = -\frac{R}{L} t + C
I(t)=eRLt+C=eCeRLt=KeRLtI(t) = e^{-\frac{R}{L}t + C} = e^C e^{-\frac{R}{L}t} = K e^{-\frac{R}{L}t}
ただし、KK は任意定数。
2) V(t)=V0V(t) = V_0 の場合を考える。
特殊解を Ip=AI_p = A (定数) と仮定する。
LIp+RIp=V0LI'_p + RI_p = V_0 に代入すると、L(0)+RA=V0L(0) + RA = V_0 より、A=V0RA = \frac{V_0}{R}
したがって、特殊解は Ip=V0RI_p = \frac{V_0}{R}
一般解は、I(t)=Ip+Ih=V0R+KeRLtI(t) = I_p + I_h = \frac{V_0}{R} + K e^{-\frac{R}{L}t}
初期条件 I(0)=0I(0) = 0 より、0=V0R+K0 = \frac{V_0}{R} + K なので、K=V0RK = -\frac{V_0}{R}
よって、I(t)=V0RV0ReRLt=V0R(1eRLt)I(t) = \frac{V_0}{R} - \frac{V_0}{R} e^{-\frac{R}{L}t} = \frac{V_0}{R}(1 - e^{-\frac{R}{L}t})
tt \rightarrow \infty で、I(t)V0RI(t) \rightarrow \frac{V_0}{R} に収束する。つまり、時間経過とともに電流は V0R\frac{V_0}{R} に近づく。
3) V(t)=V0sin(ωt)V(t) = V_0 \sin(\omega t) の場合を考える。
Ip(t)=Asin(ωt+δ)I_p(t) = A\sin(\omega t + \delta)LI+RI=V(t)LI' + RI = V(t) に代入する。
LIp=LAωcos(ωt+δ)LI'_p = LA\omega \cos(\omega t + \delta)
RIp=RAsin(ωt+δ)RI_p = RA\sin(\omega t + \delta)
LAωcos(ωt+δ)+RAsin(ωt+δ)=V0sin(ωt)LA\omega \cos(\omega t + \delta) + RA\sin(\omega t + \delta) = V_0 \sin(\omega t)
LAωcos(ωt+δ)=LAω[cos(ωt)cos(δ)sin(ωt)sin(δ)]LA\omega \cos(\omega t + \delta) = LA\omega [\cos(\omega t)\cos(\delta) - \sin(\omega t)\sin(\delta)]
RAsin(ωt+δ)=RA[sin(ωt)cos(δ)+cos(ωt)sin(δ)]RA\sin(\omega t + \delta) = RA[\sin(\omega t)\cos(\delta) + \cos(\omega t)\sin(\delta)]
A[Rcos(δ)sin(ωt)+Rsin(δ)cos(ωt)+Lωcos(δ)cos(ωt)Lωsin(δ)sin(ωt)]=V0sin(ωt)A[R\cos(\delta) \sin(\omega t) + R\sin(\delta) \cos(\omega t) + L\omega \cos(\delta) \cos(\omega t) - L\omega \sin(\delta) \sin(\omega t)] = V_0 \sin(\omega t)
A[(Rcos(δ)Lωsin(δ))sin(ωt)+(Rsin(δ)+Lωcos(δ))cos(ωt)]=V0sin(ωt)A[ (R\cos(\delta) - L\omega \sin(\delta))\sin(\omega t) + (R\sin(\delta) + L\omega \cos(\delta))\cos(\omega t) ] = V_0 \sin(\omega t)
sin(ωt)\sin(\omega t) の係数比較より、A(Rcos(δ)Lωsin(δ))=V0A(R\cos(\delta) - L\omega \sin(\delta)) = V_0
cos(ωt)\cos(\omega t) の係数比較より、A(Rsin(δ)+Lωcos(δ))=0A(R\sin(\delta) + L\omega \cos(\delta)) = 0
Rsin(δ)+Lωcos(δ)=0R\sin(\delta) + L\omega \cos(\delta) = 0
tan(δ)=LωR\tan(\delta) = -\frac{L\omega}{R}
δ=tan1(ωLR)\delta = -\tan^{-1}\left(\frac{\omega L}{R}\right)
A=V0Rcos(δ)Lωsin(δ)=V0Rcos(δ)+LωLωRcos(δ)=V0cos(δ)(R+(Lω)2R)=V0cos(δ)R2+(Lω)2RA = \frac{V_0}{R\cos(\delta) - L\omega \sin(\delta)} = \frac{V_0}{R\cos(\delta) + L\omega \frac{L\omega}{R}\cos(\delta)} = \frac{V_0}{\cos(\delta)(R + \frac{(L\omega)^2}{R})} = \frac{V_0}{\cos(\delta)\frac{R^2 + (L\omega)^2}{R}}
cos(δ)=11+tan2(δ)=11+(Lω)2R2=RR2+(Lω)2\cos(\delta) = \frac{1}{\sqrt{1+\tan^2(\delta)}} = \frac{1}{\sqrt{1 + \frac{(L\omega)^2}{R^2}}} = \frac{R}{\sqrt{R^2 + (L\omega)^2}}
A=V0RR2+(Lω)2R2+(Lω)2R=V0R2+(Lω)2A = \frac{V_0}{\frac{R}{\sqrt{R^2+(L\omega)^2}}\frac{R^2 + (L\omega)^2}{R}} = \frac{V_0}{\sqrt{R^2 + (L\omega)^2}}
4) t=0t=0 のとき、Ip(0)=Asin(δ)=Asin(tan1(ωLR))<0I_p(0) = A\sin(\delta) = A\sin(-\tan^{-1}(\frac{\omega L}{R})) < 0
t=π2ωt = \frac{\pi}{2\omega} のとき、Ip(π2ω)=Asin(ωπ2ω+δ)=Asin(π2tan1(ωLR))=Acos(tan1(ωLR))>0I_p(\frac{\pi}{2\omega}) = A\sin(\omega \frac{\pi}{2\omega} + \delta) = A\sin(\frac{\pi}{2} - \tan^{-1}(\frac{\omega L}{R})) = A\cos(-\tan^{-1}(\frac{\omega L}{R})) > 0
δ=tan1(ωLR)\delta = -\tan^{-1}\left(\frac{\omega L}{R}\right)
A=V0(ωL)2+R2A = \frac{V_0}{\sqrt{(\omega L)^2 + R^2}}

3. 最終的な答え

1) I(t)=KeRLtI(t) = K e^{-\frac{R}{L}t}
2) I(t)=V0R(1eRLt)I(t) = \frac{V_0}{R}(1 - e^{-\frac{R}{L}t})
3) δ=tan1(ωLR)\delta = -\tan^{-1}\left(\frac{\omega L}{R}\right), A=V0(ωL)2+R2A = \frac{V_0}{\sqrt{(\omega L)^2 + R^2}}
4) δ=tan1(ωLR)\delta = -\tan^{-1}\left(\frac{\omega L}{R}\right), A=V0(ωL)2+R2A = \frac{V_0}{\sqrt{(\omega L)^2 + R^2}} (上記参照)

「応用数学」の関連問題

$\alpha$, $\beta$, $\gamma$ は鋭角で、$\tan \alpha = 2$, $\tan \beta = 5$, $\tan \gamma = 8$ であるとき、以下の値を求...

三角関数加法定理角度
2025/6/5

質量 $m$ の質点が、ばね定数 $k$ のばねに繋がれており、水平面上を運動する。ばねの自然長からの質点の変位を $x$ としたとき、以下の問いに答える。 (1) 質点の運動方程式を求める。 (2)...

力学微分方程式単振動運動方程式
2025/6/5

質量 $m$ の質点が $x$ 軸上を運動しており、速度 $v$ に比例する力 $-bv$ が作用している。 (1) 運動方程式を求める。 (2) 運動方程式の解となり得るものを選択する。 (3) 初...

運動方程式微分方程式力学積分
2025/6/5

スカラー場 $\phi = 2xy - y^2z^2$ が与えられたとき、$\Delta \phi$ を求めよ。ここで $\Delta$ はラプラシアンを表す。

ベクトル解析ラプラシアン偏微分
2025/6/5

質量 $m$ の質点がx軸上を運動しており、x軸と逆向きに一定の大きさ $F$ の力が作用している。以下の3つの問いに答える。 (1) 運動方程式を求める。 (2) 運動方程式の解となり得るものをすべ...

力学運動方程式積分微分
2025/6/5

2次元ベクトル場 $\vec{F} = 2y \hat{i} + 2x \hat{j}$ の回転を計算し、その結果の意味を説明する。

ベクトル解析回転ベクトル場偏微分
2025/6/5

ベクトル $\vec{A} = 4\vec{i} + \vec{j} + 2\vec{k}$, $\vec{B} = \vec{i} - 2\vec{j} + 3\vec{k}$, $\vec{C} ...

ベクトルベクトル演算内積外積単位ベクトル三角形の面積
2025/6/5

質量$m$のおもりが長さ$l$の糸でつるされ、水平面内で等速円運動をしている。糸が鉛直方向との角度$\theta$を保っているとき、この円運動の周期$T$とおもりの速さ$v$を求める。

力学円運動物理
2025/6/5

質量 $m$ のおもりが長さ $l$ の糸でつるされ、水平面内で等速円運動をしている。糸が鉛直方向との角度 $\theta$ を保っているとき、この円運動の周期 $T$ とおもりの速さ $v$ を求め...

力学円運動物理三角関数周期速さ
2025/6/5

(1) 水平面となす角$\theta$で初速度$v_1$で投げ上げられた物体Aについて、最高点に達するまでの時間、最高点の座標、地上に落ちるまでの時間、最大水平到達距離を求める。 (2) 座標$(s,...

力学放物運動運動方程式衝突
2025/6/5