Python时域仿真RC并联电路异常：高RC值下电压过冲问题排查

方程推导 #

RC并联电路的KCL方程推导正确：

由电容电流 ( C\frac{dU}{dt} ) 与电阻电流 ( \frac{U}{R} ) 之和等于输入电流 ( I )，整理得：
dU/dt + U/(R*C) = I/C

齐次解与特解推导无误，最终通解为：
u(t) = (U_ref - R*I) * exp(-(t - t_offset)/(R*C)) + R*I
其中 ( U_ref ) 为初始电压，( t_offset ) 为初始时刻，( I ) 为当前输入电流。

Python实现逻辑 #

• 定位电流突变的时间索引（如电流在第9到10个时间步间突变，记录第10个索引）
• 遍历时间序列，逐点计算电压：U_total[i] = (U_ref - R*I) * np.exp(-(time[i] - offset_time)/(R*C)) + R*I
• 电流突变时，更新 ( U_ref ) 为前一时刻电压 ( U_total[i-1] )，更新 ( offset_time ) 为当前时间 ( time[i] )，并切换为当前电流 ( current[i] )

问题现象 #

• 电流变化缓慢时，无论R、C取值，仿真结果正常
• R、C均小于1且电流快速变化时，仿真结果正常
• 相同电流波形下，R、C大于1时，( U_total[i] ) 出现严重过冲，仅跟随 ( R*I ) 变化

疑问解答 #

1. 方程求解无问题
   推导过程完全符合RC并联电路的KCL定律，通解公式在直流输入、阶跃电流输入场景下的理论正确性已验证，部分场景仿真正常也佐证了这一点。

2. 电流突变时的初始条件更新方式在高RC值、快速变流场景下失效
   当RC值较大（时间常数大，电压响应慢）且电流快速突变时，若直接用前一时刻的电压 ( U_total[i-1] ) 作为新的初始电压 ( U_ref )，并将当前时间 ( time[i] ) 作为新的起始时刻 ( offset_time )，会忽略电流突变时刻到当前时间步之间的电压响应过程：

   • 电流突变发生在两个时间步之间，而非当前时间步的起始点
   • 高RC值下，电压本应缓慢趋近新的稳态 ( R*I )，但这种近似方式相当于强制让电压在当前时间步直接跳变到稳态趋势，从而产生过冲
   • 当RC值较小时，电压响应快，这种时间步内的近似误差可以忽略；但RC值大时，误差被放大，导致仿真结果失真

正确的处理方式应是：在检测到电流突变时，先计算突变时刻的电压值（以突变前的电流和初始条件计算到突变时刻的电压），再以此电压作为新的初始条件，从突变时刻开始计算后续的电压响应。

内容的提问来源于stack exchange，提问作者Kerdiorp