常微分2.2

3.如左图所示的PL电路，试求：

（1） 当开关$S_1$合上10s后，电感$L$上的电流

（2）$S_1$合上10s后再将$S_2$合上20s后，电感$L$上的电流

4.试求右图所示的RL电路电感上电流$I(t)$的变化规律，并解释其物理意义，设$t=0$时，$I=0$

3.解：
（1） 设图中电路上的电流为$I(t)$，其中$R_1=10\Omega ,L=2H,E=50V$，且有前提条件为$I(0)=0$，当开关合上时电路方程：
$$\frac{dI}{dt}+\frac{R_1}{L}I=\frac{E}{L},即\frac{dI}{dt}=-5I+25,I(0)=0$$

下面展示 Matlab求解代码

syms I(t);
dsolve(diff(I)==-5*I+25,I(0)==0)

求解结果为：

绘制图像：

I=5 - 5*exp(-5*t)  %dsolve函数求解结果

t=0:0.01:12;
I=eval(I);
plot(t,I)

故而最终求得电感上的电流方程为：$I\left(t\right)=5-5e^{-5t}$
当$t=10s$时，$I\left(10\right)=5-5e^{-50}\approx 5A$

（2） $S_1$合上10s后再合上$S_2$,所以有$\frac{1}{R}=\frac{1}{R_1}+\frac{1}{R_2}=\frac{3}{20}$，此时电路方程为$$\frac{dI}{dt}=-\frac{10}{3}I+25,I(0)=5$$
下面展示 Matlab求解代码

syms I(t);
dsolve(diff(I)==-10/3*I+25,I(0)==5)

求解结果为：

绘制图像：

I=15/2 - (5*exp(-(10*t)/3))/2  %dsolve函数求解结果

t=0:0.01:20;
I=eval(I);
plot(t,I)

故而最终求得电感上的电流方程为：$I\left(t\right)=7.5-2.5e^{-\frac{10}{3}t}$
当$t=20s$时，$I\left(20\right)=7.5-2.5e^{\frac{200}{3}}\approx 7.5A$

4.解：
根据$RI+L\frac{dI}{dt}=E=U_m\sin wt$可得$$\frac{dI}{dt}=-\frac{R}{L}I+\frac{U_m}{L}\sin wt,且I\left(0\right)=0$$

下面展示 Matlab求解代码

syms I(t);
sym R;
syms L;
syms U;
syms w;
dsolve(diff(I)== -R/L*I+U+L*sin(w*t),I(0)==0)

求解结果为：

故而最终求得电感上的电流方程为：$$I\left(t\right)=\frac{L\left(R^{2}U+L{R}^2\sin wt+L^{2}U{W}^2-L^{2}Rw\cos wt\right)}{R\left(L^2{w}^2+R^2\right)}-e^{-\frac{Rt}{L}}\frac{-L^{3}Rw+U{L}^3{w}^2+UL{R}^2}{L^{2}R{w}^2+R^3}$$
物理意义：电路上电流分两部分，一部分随时间逐渐衰减最后为0；一部分稳定，与电动势周期相同。