正弦電動勢 RLC 串聯電路的電流分析 | 線性二階微分方程應用

 
Current Analysis of A Series RLC Circuit Driven by A Sinusoidal Voltage Input


如圖一所示,為一正弦電動勢的 RLC 串聯電路。

正弦電動勢 RLC 串聯電路
圖一  正弦電動勢 RLC 串聯電路

假設迴路的電流為 $i$,則電阻 $R$ 的電壓降為 $Ri$,電感 $L$ 的電壓降為 $L\frac{di}{dt}$,電容 $C$ 的電壓降為 $\frac{1}{C}\int i\;dt$。根據 $KVL$,可知
$$Ri+L\frac{di}{dt}+\frac{1}{C}\int i\;dt=E\;\sin \omega t$$
上式等號兩邊對 $t$ 微分,可得
$$R\frac{di}{dt}+L\frac{d^2i}{dt^2}+\frac{i}{C}=E\omega\cos \omega t$$
整理得
$$\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\frac{d^2i}{dt^2}+\frac{R}{L}\frac{di}{dt}+\frac{i}{LC}=\frac {E}{L}\omega \cos \omega t\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(1)$$
為一線性二階非齊次微分方程 (Second-order Linear Non-homogeneous Differential Equation)。其通解 (General Solution) 的形式為 $i_h+i_P$,其中的 $i_h$ 為齊次解 (Homogeneous Solution),而 $i_p$ 則是特殊解 (Particular Solution)。


求齊次解 $i_h$

要找出式子 (1) 的齊次解,即是找出下面式子的的通解。
$$\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\frac{d^2i}{dt^2}+\frac{R}{L}\frac{di}{dt}+\frac{i}{LC}=0\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(2)$$
而 (2) 的特徵方程式 (Characteristic Equation) 為
$$r^2+\frac{R}{L}r+\frac{1}{LC}=0$$
可得
$$r_1=\frac{-\frac{R}{L}+\sqrt{(\frac{R}{L})^2-\frac{4}{LC}}}{2}\;\;\;\;$$
$$\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;=\frac{-R+\sqrt{R^2-\frac{4L}{C}}}{2L}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(3)$$
以及
$$\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;r_2=\frac{-R-\sqrt{R^2-\frac{4L}{C}}}{2L}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(4)$$
所以齊次解
$$\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;i_h=c_1e^{r_1t}+c_2e^{r_2t}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(5)$$

而由 (3) 與 (4) 可知,要決定 $r_1$ 與 $r_2$,有三種情況需要分析:

情況 1   過阻尼 (Overdamping) $(R^2-\frac{4L}{C} \gt 0)$

可知 $r_1$ 與 $r_2$ 皆是負值,所以 (5) 的 $c_1e^{r_1t}+c_2e^{r_2t}$ 值,會隨著時間 $t$ 增加而快速減少而最終趨近於 0 。因此,在此情況下 $i_h$ 稱為暫態解 (Transient Solution),而從後面分析可看到呈現週期性波形型式的特別解 $i_p$ 則稱為穩態解 (Steady-state Solution) 。

情況 2   臨界阻尼 (Critical Damping) $(R^2-\frac{4L}{C} = 0)$

可知 $r_1=r_2$,則
$$i_h=(c_1+tc_2)e^{\frac{-R}{2L}t}$$
同樣地,$(c_1+tc_2)e^{\frac{-R}{2L}t}$ 的值,會隨著 $t$ 增加而快速減少而最終趨近於 0 。

情況 3   欠阻尼 (Underdamping) $(R^2-\frac{4L}{C} \lt 0)$

令 $P=\frac {4L}{C}-R^2$,則
$$r_1=\frac {-R+jP}{2L}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;r_2=\frac {-R-jP}{2L}$$
代入 (5),可得
$$i_h=c_1e^{\frac {-R+jP}{2L}t}+c_2e^{\frac {-R-jP}{2L}t}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
$$=e^{-\frac{R}{2L}t}\left ( c_1e^{\frac {jP}{2L}t}+c_2e^{\frac {-jP}{2L}t}\right )\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
$$=e^{-\frac{R}{2L}t}\left \{c_1\left [cos\left(\frac{P}{2L}t\right)+j\;sin \left(\frac{P}{2L}t\right)\right ]+c_2\left [cos\left (\frac{P}{2L}t\right)-j\;sin\left (\frac{P}{2L}t\right)\right ] \right \}\;\;\;\;\;\;\;\;$$
$$=e^{-\frac{R}{2L}t}\left [ (c_1+c_2)cos \left (\frac{P}{2L}t\right )-j(c_2-c_1)sin \left (\frac {P}{2L}t\right ) \right ]\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
$$=Ke^{-\frac{R}{2L}t}cos\left(\frac{P}{2L}t+\theta \right )\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
其中的 $K$ 與 $\theta$ 均為常數。從上式可看出,$i_h$ 一樣會隨 $t$ 的增加而快速衰減至近乎 0。


求特殊解 $i_p$

觀察 (1) 要找到特殊解,可令 
$$i_p=A\cos \omega t+B\sin \omega t\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(6)$$
則
$$i_p'=-\omega A\sin \omega t+B\omega \cos \omega t$$
且
$$i_p''=-\omega^2A\cos \omega t-B\omega^2\sin \omega t$$
將上面三式代入 (1),可得
$$-A\omega^2\cos \omega t-B\omega^2\sin \omega t+\frac{R}{L}(-\omega A\sin \omega t+B\omega\cos \omega t)+\frac{A}{LC}\cos \omega t+\frac{B}{LC}\sin \omega t=\frac{E\omega}{L}\cos \omega t$$
整理得
$$\cos \omega t \left ( -A\omega^2+\frac{BR\omega}{L}+\frac{A}{LC}\right )+\sin \omega t \left (-B\omega^2-\frac{AR\omega }{L}+\frac{B}{LC}\right ) = \frac{E\omega}{L}\cos \omega t$$
即
$$-A \left (\omega^2-\frac{1}{LC}\right )+B \left (\frac{R \omega}{L}\right )=\frac{E \omega}{L}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(7)$$
$$-A\left (\frac{R \omega}{L}\right )+B \left (\frac{1}{LC}-\omega^2 \right )=0\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(8)$$
為簡化上面的聯立方程式,令
$$m=\omega^2-\frac{1}{LC}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(9)$$
且
$$n=\frac{R\omega}{L}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(10)$$
則 (7) 與 (8) 可改寫如下:
$$-Am+Bn=\frac{En}{R}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(11)$$
$$-An-Bm=0\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(12)$$
解 (11) 與 (12) 聯立方程式,可得
$$A=-\frac{Enm}{R(n^2+m^2)}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(13)$$
以及
$$B=\frac{En^2}{R(n^2+m^2)}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(14)$$
將 (9)、(10) 代入 (13) 及 (14),可知
$$A=\frac{-E\frac{R\omega}{L}(\omega^2-\frac{1}{LC})}{R \left [(\frac{R\omega}{L})^2+(\omega^2-\frac{1}{LC})^2 \right ]}$$
$$\;=\frac{-E\frac{\omega^2}{L^2}(\omega L-\frac{1}{\omega C})}{\frac{\omega^2}{L^2}R^2+\frac{\omega^2}{L^2}(\omega L-\frac{1}{\omega C})^2}$$
$$\;=\frac{-E(\omega L-\frac{1}{\omega C})}{R^2+(\omega L-\frac{1}{\omega C})^2}\;\;\;\;\;\;\;\;\;$$
令
$$X=\omega L-\frac{1}{\omega C}$$
則
$$A=\frac{-EX}{R^2+X^2}$$
又
$$B=\frac{E(\frac{R\omega}{L})^2}{R \left [(\frac{R\omega}{L})^2+(\omega^2-\frac{1}{LC})^2 \right ]}$$
$$\;\;\;\;=\frac{(\frac{\omega}{L})^2ER}{(\frac{\omega}{L})^2\left [R^2+(\omega L-\frac{1}{\omega C})^2\right ]}$$
$$=\frac{ER}{R^2+X^2}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
接著再將 $A$ 與 $B$ 代入 (6),即
$$i_p=\frac{-EX}{R^2+X^2}\cos \omega t+\frac{ER}{R^2+X^2}\sin \omega t\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
$$\;\;\;\;\;\;=\frac{E}{\sqrt{R^2+X^2}}\left (\frac{-X}{\sqrt {R^2+X^2}}\cos \omega t+\frac{R}{\sqrt {R^2+X^2}}\sin \omega t\right )$$
令
$$\sin \theta=\frac{X}{\sqrt {R^2+X^2}},\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;cos\;\theta=\frac{R}{\sqrt {R^2+X^2}}$$
則
$$i_p=\frac{E}{\sqrt {R^2+X^2}}\sin(\omega t-\tan^{-1}\frac{X}{R})$$
因為當時間 $t$ 增加,很快地 $i_h$ 就會衰減至趨近於 $0$,整個 $i$ 等式,則會趨近於 $i_p$,即穩態解。而此穩態解還可以採用工程上常用複數方法與相量式的方式快速地得到,說明如下:


以複數方法 (Complex Method) 求取 $i_p$

由 (1) 可知,等號右邊為 $\frac{E\omega }{L}\cos \omega t$,而由尤拉公式 (Euler's Formula) 
$$e^{j\omega t}=\cos \omega t+j\sin \omega t$$
可知
$$\cos \omega t=Re \left \{e^{j\omega t}\right\}$$
其中的 $Re$ 表示實部之意。因此,為求 (1) 的穩態解,我們先考慮如下的式子:
$$\frac{d^2i}{dt^2}+\frac{R}{L}\frac{di}{dt}+\frac{i}{LC}=\frac{E\omega}{L}e^{j\omega t}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(15)$$
令複數穩態解
$$i_p^*(t)=Ke^{j\omega t}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(16)$$ 
則
$$i_p^{*'}(t)=j\omega Ke^{j\omega t}$$
且
$$i_p^{*''}(t)=-\omega^2K e^{j\omega t}$$
代入 (15) 可得
$$\left (-\omega^2K+\frac{R}{L}j\omega K+\frac{K}{LC}\right )e^{j\omega t}=\frac{E\omega}{L}e^{j\omega t}$$
$$\left ( \frac{1}{LC}-\omega^2+j\frac{R\omega}{L}\right )K=\frac{E\omega}{L}$$
可知
$$K=\frac{E\omega}{L\left (\frac{1}{LC}-\omega^2+j\frac{R\omega}{L}\right )}$$
$$=\frac{E\omega}{\frac{1}{C}-L\omega^2+jR\omega}\;\;$$
$$=\frac{E}{\frac{1}{\omega C}-\omega L+jR}\;\;\;\;$$
令 
$$X=\omega L-\frac {1}{\omega C}$$
則
$$K=\frac{E}{-X+jR}$$
$$\;\;\;\;\;\;\;\;\;\;=\frac{E(-X-jR)}{S^2+R^2}$$
代回 (16),可得
$$i_p^{*}=\frac{E\left (-X-jR\right )}{X^2+R^2}e^{j\omega t}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
$$=\frac{E(-X-jR)}{X^2+R^2}(\cos \omega t+j\sin \omega t)\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
$$=\frac{E}{X^2+R^2} [(-X\cos \omega t+R\sin \omega t)+j(-X\sin \omega t-R\cos \omega t) ]\;\;$$
因此
$$i_p=\frac{E}{X^2+R^2} (-X\cos \omega t+R\sin \omega t)\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
$$=\frac{E\sqrt{(-X)^2+R^2}}{X^2+R^2}\left (\frac{-X}{\sqrt {(-X)^2+R^2}}\cos \omega t+\frac{R}{\sqrt {(-X)^2+R^2}}\sin \omega t \right )$$
$$=\frac{E}{\sqrt {R^2+X^2}}\sin \left (\omega t-\tan^{-1}\frac{X}{R}\right )\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$


以相量 (Phasor) 方法求取 $i_p$

已知電源電壓 $E\sin\;\omega t$ 的相量式
$$\overline E=E\;\;\;\;\;\;\;\;\;\;(直角坐標表示法)$$
$$\;\;\;\;\;\;=E\angle  0^\circ\;\;\;\;\;\;\;\;\;\;(極坐標表示法)$$
電阻 $R$ 的相量式
$$\overline R=R\;\;\;\;\;\;\;\;\;\;(直角坐標表示法)$$
$$\;\;\;\;\;\;\;=R\angle  0^\circ\;\;\;\;\;\;\;\;\;\;(極坐標表示法)$$
電容抗 $X_C$ 的相量式
$$\overline X_C=-j\;X_C=-j\frac{1}{\omega C}\;\;\;\;\;\;\;\;\;\;(直角坐標表示法)$$
$$\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;=X_C\angle  -90^\circ=\frac {1}{\omega C}\angle  -90^\circ\;\;\;\;\;\;\;\;\;\;(極坐標表示法)$$
電感抗 $X_L$ 的相量式
$$\overline X_L=j\;X_L=j\;\omega L\;\;\;\;\;\;\;\;\;\;(直角坐標表示法)$$
$$\;\;\;\;\;\;\;\;\;\;\;\;\;\;=X_L\angle  90^\circ=\omega L\angle  90^\circ\;\;\;\;\;\;\;\;\;\;(極坐標表示法)$$
因此,根據歐姆定律 (Ohm's Law)
$$\overline I_p=\frac{\overline E}{\overline R+\overline X_C+\overline X_L}\;\;\;\;\;\;\;$$
$$=\frac{E}{R-j\frac{1}{\omega C}+j\;\omega L}$$
$$=\frac{E}{R+j(\omega L-\frac{1}{\omega C})}$$
令 $X=\omega L-\frac{1}{\omega C}$,$X$ 在此稱為電抗 (Reactance),則
$$\overline I_p=\frac{E}{R+jX}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
$$\;\;\;\;\;\;=\frac{E\angle 0^\circ}{\sqrt{R^2+X^2}\angle \tan^{-1}\frac{X}{R}}$$
$$\;\;\;\;\;\;\;\;\;\;\;=\frac{E}{\sqrt {R^2+X^2}}\angle -\tan^{-1}\frac{X}{R}$$
故
$$i_p=\frac{E}{\sqrt {R^2+X^2}}\sin \left(\omega t-\tan^{-1}\frac{X}{R}\right )$$




最後通解 $i(t)$

$$i=i_h+i_p=c_1e^{r_1t}+c_2e^{r_2t}+\frac{E}{\sqrt {R^2+X^2}}\sin \left(\omega t-tan^{-1}\frac{X}{R}\right )$$
其中 $r_1=\frac{-R+\sqrt{R^2-\frac{4L}{C}}}{2L}$,$r_2=\frac{-R-\sqrt{R^2-\frac{4L}{C}}}{2L}$,且 $X=\omega L-\frac{1}{\omega C}$。

而上式的唯二未知常數 $c_1$ 與 $c_2$ 的求取,則須借助於電路的兩個初始條件。如果初始電流 $i(0)=0$,且電容的初始電荷 $Q(0)=0$,則迴路應用 $KVL$ 所得到的式子
$$Ri+L\frac{di}{dt}+\frac{1}{C}\int i\;dt=E\sin \omega t$$
可改寫成
$$Ri+L\frac{di}{dt}+\frac{Q}{C}=E\sin \omega t$$
當 $t=0$ 時,
$$Ri(0)+Li'(0)+\frac{Q(0)}{C}=0$$
因此,可知
$$i'(0)=0$$
而有了 $i(0)=0$ 以及 $i'(0)=0$ 這兩個初始條件,就可以找出 $c_1$ 與 $c_2$ 的值。


範例:如圖一所示 RLC 串聯電路,若電源電壓 $v(t)=100\sin (1000t)\;V$,R=4 $\Omega$,C=250$\mu$ F,L=1m H,試求穩態時的 $i(t)$、$v_L(t)$ 以及 $v_C(t)$。

解:
$$\overline V=100\angle 0^\circ\;V$$
$$X_L=\omega L=1000\times 1m=1\;\Omega$$
$$X_C=\frac{1}{\omega C}=\frac{1}{1000\times 250\mu}=\frac{1}{1000\times 250\times 10^{-6}}=4\;\Omega$$
又
$$\overline I=\frac{\overline V}{\overline R+\overline X_L+\overline X_C}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
$$=\frac{\overline V}{R+jX_L-jX_C}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
$$\;=\frac{\overline V}{R+j(X_L-X_C)}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
$$\;\;\;\;\;\;\;\;\;\;\;\;\;\;=\frac{V\angle 0^\circ}{\sqrt {R^2+(X_L-X_C)^2}\;\angle \tan^{-1}\frac{X_L-X_C}{R}}$$
$$\;\;=\frac{100\angle 0^\circ}{\sqrt {4^2+(1-4)^2}\;\angle \tan^{-1}\frac{1-4}{4}}$$
$$=\frac {100\angle 0^\circ}{5\angle -36.9^\circ}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
$$=20\angle 36.9^\circ\;A\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
因此
$$i(t)=20\sin (1000t+36.9^\circ)\;A$$
又
$$\overline V_L=\overline I\;\overline X_L\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
$$\;\;=20\angle 36.9^\circ\times1\angle 90^\circ$$
$$=20\angle 126.9^\circ\;V\;\;\;\;\;\;$$
所以
$$v_L(t)=20\sin (1000t+126.9^\circ)\;V$$
而
$$\overline V_C=\overline I\;\overline X_C\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;$$
$$\;\;=20\angle 36.9^\circ\times 4\angle -90^\circ$$
$$=80\angle -53.1^\circ\;V\;\;\;\;\;\;\;\;$$
故
$$v_C(t)=80\sin (1000t-53.1^\circ)\;V$$





留言

這個網誌中的熱門文章

三段式電子開關電路

陣列 (C++)

分壓偏壓 BJT 放大電路的直流分析及其近似解的條件

首數、尾數與位數

MOSFET 共汲極放大電路 (源極隨耦器) 小訊號分析

為什麼理想的 OPA 電壓放大器有虛短路與虛斷路現象 | OPA 電壓放大器增益推導

從 0 積到 1,sqrt (x^2+1) 的積分 | 除法定則推導

OP-Amp 方波產生電路 | OP-Amp 方波產生電路的數學分析