问题描述
要计算下边电路中电容两端电压随时间变化。(根据戴维南定理,可以用一个独立电压源与电阻串联代替其他一般线性电路)
设初始条件为:
- v c ( 0 ) = V 0 v_c(0) = V_0 vc(0)=V0,电容的初始电压。
- v i ( t ) = V I v_i(t) = V_I vi(t)=VI,直流电压源电压。
现在来看下 t > 0 t > 0 t>0 之后的电容电压 v c ( t ) v_c(t) vc(t) 随时间的变化。
数学推导
根据电容的伏安特性,有
i = C d v d t i = C \frac{dv}{dt} i=Cdtdv
根据节点电压法可以写出这个电路的节点方程
R C d v c ( t ) d t + v c ( t ) = v i ( t ) RC\frac{dv_c(t)}{dt} + v_c(t) = v_i(t) RCdtdvc(t)+vc(t)=vi(t)
因为设定的初始条件 $ v_i(t) = V_I $
带入后,方程为:
R C d v c ( t ) d t + v c ( t ) = V I RC\frac{dv_c(t)}{dt} + v_c(t) = V_I RCdtdvc(t)+vc(t)=VI
我们的目的是解出 v c ( t ) v_c(t) vc(t)
这里使用求特解加齐次解的方法解这个微分方程
v c ( t ) = v c p ( t ) + v c h ( t ) v_c(t) = v_{cp}(t) + v_{ch}(t) vc(t)=vcp(t)+vch(t)
v c p ( t ) v_{cp}(t) vcp(t) ,特解
v c h ( t ) v_{ch}(t) vch(t) ,齐次解
下面是求解的具体的三个步骤:
- 找到一个特殊解
- 找出对应的齐次方程的通解
- 根据初始条件,计算常数参数。
1.找特殊解
R C d v c p ( t ) d t + v c p ( t ) = V I RC\frac{dv_{cp}(t)}{dt} + v_{cp}(t) = V_I RCdtdvcp(t)+vcp(t)=VI
可以尝试 v c p = V I v_{cp} = V_I vcp=VI ,带入微分方程中进行验证
因为常数的微分为 0, 所以上式为:
$ 0 + V_I = V_I$
等式成立,可以确定 v c p = V I v_{cp} = V_I vcp=VI 是 $ v_c(t) $ 的一个特解,得到
v c p = V I v_{cp} = V_I vcp=VI 。
2.计算齐次方程的通解
R C d v c ( t ) d t + v c ( t ) = V I RC\frac{dv_c(t)}{dt} + v_c(t) = V_I RCdtdvc(t)+vc(t)=VI 对应的齐次方程为
R C d v c h ( t ) d t + v c h ( t ) = 0 RC\frac{dv_{ch}(t)}{dt} + v_{ch}(t) = 0 RCdtdvch(t)+vch(t)=0
因为 d e x d x = e x \frac{de^x}{dx} = e^x dxdex=ex
所以一般这种微分方程可以尝试使用 v c h ( t ) = A e s t v_{ch}(t) = Ae^{st} vch(t)=Aest 这种形式的通式来尝试,其中 A 、 s A 、 s A、s 都为待确定参数,
带入验证,
R C d A e s t d t + A e s t = 0 RC\frac{d{Ae^{st}}}{dt} + Ae^{st} = 0 RCdtdAest+Aest=0
可以化为
R C A s e s t + A e s t = 0 RCAse^{st} + Ae^{st} = 0 RCAsest+Aest=0
这里舍去 A = 0 A = 0 A=0 的解, 可以约掉 A e s t Ae^{st} Aest,化简为
R C s + 1 = 0 RCs + 1 = 0 RCs+1=0
即当 s = − 1 R C s = -\frac{1}{RC} s=−RC1 时, v c h ( t ) = A e s t v_{ch}(t) = Ae^{st} vch(t)=Aest 成立。
将 R C 记为 τ RC 记为 \tau RC记为τ , 读作 tào,单位是时间
v c h ( t ) = A e − t τ v_{ch}(t) = Ae^{-\frac{t}{\tau}} vch(t)=Ae−τt
3.带入初始条件,计算常数参数
根据前两步计算,我们还差一个未知常数参数 A 没有得到,这一步使用初始条件来计算 A ,目前我们已经得到了
-
v c p = V I v_{cp} = V_I vcp=VI
-
v c h ( t ) = A e − t τ v_{ch}(t) = Ae^{-\frac{t}{\tau}} vch(t)=Ae−τt
带入 v c ( t ) = v c p ( t ) + v c h ( t ) v_c(t) = v_{cp}(t) + v_{ch}(t) vc(t)=vcp(t)+vch(t) 可得
v c ( t ) = V I + A e − t τ , τ = R C v_c(t) = V_I + Ae^{-\frac{t}{\tau}},\tau = RC vc(t)=VI+Ae−τt,τ=RC
已知 t = 0 t = 0 t=0 时, v c ( t ) = V 0 v_c(t) = V_0 vc(t)=V0
可得
V 0 = V I + A e − 0 τ V_0 = V_I + Ae^{-\frac{0}{\tau}} V0=VI+Ae−τ0
V 0 = V I + A V_0 = V_I + A V0=VI+A
A = V 0 − V I A = V_0 - V_I A=V0−VI
最终结果
通过以上三步,可以得到最终的结果为:
v c ( t ) = V I + ( V 0 − V I ) e − t τ , τ = R C v_c(t) = V_I + (V_0 - V_I)e^{-\frac{t}{\tau}},\tau = RC vc(t)=VI+(V0−VI)e−τt,τ=RC
举例验证
通过设定一些具体的参数进行验证,设
- V 0 = 5 v V_0 = 5v V0=5v ,电容上携带的初始电压
- V I = 3.4 v V_I = 3.4v VI=3.4v ,电路中独立电压源的直流电压
- R = 100 k Ω R = 100k\Omega R=100kΩ ,电阻阻值
- C = 10 u F C = 10uF C=10uF ,电容容值
将这些具体数值带入函数中,可得
v c ( t ) = 3.4 + ( 5 − 3.4 ) e − t 1 v_c(t) = 3.4 + (5 - 3.4)e^{-\frac{t}{1}} vc(t)=3.4+(5−3.4)e−1t
绘制曲线
使用绘图软件,将下边这个函数绘制出来
v c ( t ) = 3.4 + ( 5 − 3.4 ) e − t 1 v_c(t) = 3.4 + (5 - 3.4)e^{-\frac{t}{1}} vc(t)=3.4+(5−3.4)e−1t
可以得到如下曲线,只考虑 t > 0 的部分
LTspice仿真
下边是仿真结果
元件实测
实验测试通过设置下降沿单次触发得到的波形曲线
对电阻电容的测量,误差不大,仿真就按标称值进行的
电脑电源改装的电压源,手动插拔切换的电压
用面包板和直插的元件搭建的测试电路
参考
麻省理工学院公开课——电路与电子学 p13: https://www.bilibili.com/video/BV1ts411v7Ep?p=13&vd_source=55a350b8650fa259fcd33c6529c36cf5