一般RL串联电路的微分方程计算

电感电压 Vl=di/dt，然后利用KVL等式
U=Ri+Ldi/dt
然后变换一下

你可以用wxMaxima，创建一个总和那些正弦信号的信号
s(t) := sum((1/(2*n-1))*sin(2*%pi*(2*n-1)*f*t),n,1,3);
我做了一个50Hz的你可以看下

可以假设是通过一个串联的开关来施加电压，电流的等式就是i(t)=V/R(1-exp[(-R/L)t]，其中V就是电流里面的U