1、理论基础
L
C
d
2
U
C
d
t
2
+
R
C
d
U
C
d
t
+
U
c
=
U
o
p
L C \frac{d^2{U_C }}{dt^2 }+RC\frac{ dU_C }{dt}+U_c=U_{op}
LCdt2d2UC+RCdtdUC+Uc=Uop
2、算法
运用龙格库塔法解上述微分方程
f1 = f(UC,I,t)
g1 = g(UC,I,t)
f2 = f(UC+f1*h/2,I+h/2*g1,t+h/2)
g2 = g(UC+f1*h/2,I+h/2*g1,t+h/2)
f3 = f(UC+f2*h/2,I+h/2*g2,t+h/2)
g3 = g(UC+f2*h/2,I+h/2*g2,t+h/2)
f4 = f(UC+f3*h,I+g3*h,t+h)
g4 = g(UC+f3*h,I+g3*h,t+h)
UC = UC +1/6*(f1+f2*2+f3*2+f4)*h
I = I +1/6*(g1+g2*2+g3*2+g4)*h
3、结果