线性系统粗浅认识——第五次作业

声明：本人特别菜，不研究相关的方向，差点挂科，这个作业的内容仅供交流。

题目

仿真说明LC电路不是BIBO稳定。设计一个输入信号使得电容电压趋向无穷，说明该信号与LC的关系。
我进行了三个方面的分析

解答

一、串联RLC系统

串联RLC系统如下：

图1 RLC电路电路图

输出时电容电压$U_c$，输入是电源电压$U_s$ ，电容大小为$C$，电感为$L$，电阻大小为$R$。
$$\begin{array}{l}{U}_s(t)=Ri(t)+L\frac{di}{dt}+U_c(t)\\ i(t)=C\frac{d{U}_c}{dt}\end{array}$$
可以得到
$$U_s(t)=RC\frac{d{U}_c}{dt}+LC\frac{d{U}_c^2}{d{t}^2}+U_c$$
两边进行拉普拉斯变换可以得到
$$U_s(s)=RCs{U}_c(s)+LC{s}^2{U}_c(s)+U_c(s)$$
所以
$$\frac{U_c(s)}{U_s(s)}=\frac{1}{LC{s}^2+RCs+1}=\frac{1}{LC(s-{\lambda }_1)(s-{\lambda }_2)}$$

$\mathrm{Re}({\lambda }_1)<0,$且$\mathrm{Re}({\lambda }_2)<0$ ，因此所有的根都位于$s$的开的左半平面，因此系统一定是BIBO稳定。
当$(RC{)}^2-4LC=0$时候，${\lambda }_1={\lambda }_2$ ，该系统具有重根
下面令$C=100F,R=100\Omega ,L=4H$，满足 $(RC{)}^2-4LC=0$
利用MATLAB的sinmulink模块进行仿真
仿真电路图如下所示：

图2 RLC电路仿真图

图3 RLC电路电压仿真图

从图中看出仿真结果输出是有界的，与上述结论该系统BIBO稳定的相符合。

二、串联LC系统

从RLC系统中去掉$R$，输出时电容电压$U_c$，输入是电源电压$U_s$，电容大小为$C$，电感为$L$。
$$\begin{array}{l}{U}_s(t)=L\frac{di}{dt}+U_c(t)\\ i(t)=C\frac{d{U}_c}{dt}\end{array}$$
可以得到
$$U_s(t)=LC\frac{d{U}_c^2}{d{t}^2}+U_c$$
两边进行拉普拉斯变换可以得到
$$G(s)=\frac{U_c(s)}{U_s(s)}=\frac{1}{LC{s}^2+1}=\frac{\frac{1}{LC}}{s^2+\frac{1}{LC}}=\frac{1}{\sqrt{LC}}\frac{\frac{1}{\sqrt{LC}}}{s^2+{(\frac{1}{\sqrt{LC}})}^2}$$
所以
$$g(t)=\frac{1}{\sqrt{LC}}\sin (\frac{1}{\sqrt{LC}}t)$$
进行反拉式变换
$$g(t)=\frac{1}{\sqrt{LC}}\sin (\frac{1}{\sqrt{LC}}t)$$

输入为直流电压源

$$U_s(t)=U_s<M_1$$
$$y(t)={\int }_0^t\frac{1}{\sqrt{LC}}\sin (\frac{1}{\sqrt{LC}}(t-\tau ))U_{s}d\tau =U_s(1-\cos (\frac{1}{\sqrt{LC}}t))$$
$$\underset{t\to \infty }{\lim }\left|y(t)\right|=\underset{t\to \infty }{\lim }\left|U_s(1-\cos (\frac{1}{\sqrt{LC}}t))\right|<2U_s=M$$

此时满足$\left|y(t)\right|<M$ ，$M$为任意的正实数，因此该LC系统满足BIBO稳定。$L=10^{-2}H,C=10^{-5}C,U_s=10V$，仿真的电路图为：

图4 直流电源 LC电路仿真的电路图

仿真时间是6s,示波器输出的结果为：

图５　直流电压电源输入LC电路电压输出

从图中观察，$y(t)$有界，和上面的公式分析一致，此时有界输入能产生有界的输出。

输入为正弦波交流电压源

$U_s(t)=U_s\sin (w_{1}t+\phi )<\left|U_s\right|\le M_1$
利用积化和差三角公式
$$\begin{gathered} y(t)={\int }_0^t\frac{1}{\sqrt{LC}}\sin (\frac{1}{\sqrt{LC}}(t-\tau ))U_s\sin (w_1\tau +\phi )d\tau ={\int }_0^t\frac{U_s}{\sqrt{LC}}\sin (\frac{1}{\sqrt{LC}}(t-\tau ))\sin (w_1\tau +\phi )d\tau \\ ={\int }_0^t\frac{U_s}{2\sqrt{LC}}\left[\cos (\frac{1}{\sqrt{LC}}(t-\tau )-(w_1\tau +\phi ))-\cos (\frac{1}{\sqrt{LC}}(t-\tau )+(w_1\tau +\phi ))\right]d\tau \\ ={\int }_0^t\frac{U_s}{2\sqrt{LC}}\left[\cos (-(\frac{1}{\sqrt{LC}}+w_1)\tau +\frac{1}{\sqrt{LC}}t-\phi ))-\cos ((w_1-\frac{1}{\sqrt{LC}})\tau +\frac{1}{\sqrt{LC}}t+\phi ))\right]d\tau \\ ={\frac{U_s}{2\sqrt{LC}}\left[\frac{1}{-(\frac{1}{\sqrt{LC}}+w_1)}\sin (-(\frac{1}{\sqrt{LC}}+w_1)\tau +\frac{1}{\sqrt{LC}}t-\phi ))-\frac{1}{(w_1-\frac{1}{\sqrt{LC}})}\sin ((w_1-\frac{1}{\sqrt{LC}})\tau +\frac{1}{\sqrt{LC}}t+\phi ))\right]|}_0^t \\ =\frac{U_s}{2\sqrt{LC}}\frac{1}{-(\frac{1}{\sqrt{LC}}+w_1)}\left[\sin (-w_{1}t-\phi )-\sin (\frac{1}{\sqrt{LC}}t-\phi )\right]+\frac{U_s}{2\sqrt{LC}}\frac{1}{\frac{1}{\sqrt{LC}}-w_1}\left[\sin (w_{1}t+\phi )-\sin (\frac{1}{\sqrt{LC}}t+\phi )\right] \end{gathered}$$
根据和差化积的三角公式可以得到
$$y(t)=\frac{U_s}{2\sqrt{LC}}\frac{1}{-(\frac{1}{\sqrt{LC}}+w_1)}\left[\sin (-w_{1}t-\phi )-\sin (\frac{1}{\sqrt{LC}}t-\phi )\right]+\frac{U_s}{2\sqrt{LC}}\frac{1}{\frac{1}{\sqrt{LC}}-w_1}\left[\sin (w_{1}t+\phi )-\sin (\frac{1}{\sqrt{LC}}t+\phi )\right]$$
$$\begin{gathered} \underset{t\to \infty }{\lim }\left|y(t)\right|=\underset{t\to \infty }{\lim }\left|\frac{U_s}{(1+w_1\sqrt{LC})}+\frac{U_s}{(1-w_1\sqrt{LC})}\right| \\ <\left|\frac{U_s}{(1+w_1\sqrt{LC})}\right|+\left|\frac{U_s}{(1-w_1\sqrt{LC})}\right|=M \end{gathered}$$

因此，输入有界的正弦交流电时，该系统的输出时有界的$L=10^{-2}H,C=10^{-2}C,U_s=20\sin (10^{2}t)$，仿真电路图为：

图６ 交流电源 LC电路仿真的电路图

示波器输出结果

图７　交流电源输入LC电路电压输出

仿真的输出结果有界，和分析结果一致。

输入为sgn函数控制电压源

此时没有考虑到相同频率的情况，相同频率类似于并联分析
$U_s(t)=U_s\mathrm{sgn}(\sin (\frac{1}{\sqrt{LC}}(t-\tau )))<\left|U_s\right|\le M_1$
利用积化和差三角公式
$$y(t)={\int }_0^t\frac{1}{\sqrt{LC}}\sin (\frac{1}{\sqrt{LC}}(t-\tau ))U_s\mathrm{sgn}(\sin (\frac{1}{\sqrt{LC}}(t-\tau )))d\tau ={\int }_0^t\left|\frac{U_s}{\sqrt{LC}}\sin (\frac{1}{\sqrt{LC}}(t-\tau ))\right|d\tau$$
$\left[\frac{t}{2\pi \sqrt{LC}}\right]$表示$\frac{t}{2\pi \sqrt{LC}}$ 向下取整
$\underset{t\to \infty }{\lim }\left|y(t)\right|$ 不是有界，有界的输入不能产生有界的输出，该系统非BIBO稳定$L=10^{-2}H,C=10^{-2}C,$用方波进行仿真，仿真电路图为：

图8 方波控制电压源LC电路仿真的电路图

示波器输出结果

图9　方波控制电压源输入LC电路电压输出

方波仿真的输出结果有界，方波并不是sgn函数信号，在仿真中我不知如何构造sgn函数，或者说这个输入信号已经改变系统的结构，可能不存在。

三、对并联的LC电路进行分析

输出时电容电压$U_c$，输入是电源电流$I_s$，电容大小为$C$，电感为$L$。
$$\begin{array}{l}{U}_c(t)=L\frac{d{i}_L}{dt}\\ i_c(t)=C\frac{d{U}_c}{dt}\\ i_c(t)+i_c(t)=I_s\end{array}$$
可以得到
$$C\frac{d{U}_c}{dt}+{\int }_0^t\frac{U_c}{L}dt=I_s$$
两边进行拉普拉斯变换可以得到
$$sC{U}_c(s)+\frac{U_c(s)}{sL}=I(s)$$

所以
$$G(s)=\frac{U_c(s)}{I_s(s)}=\frac{sL}{s^{2}CL+1}=\frac{1}{c}\frac{s}{s^2+\frac{1}{CL}}$$
进行反拉式变换
$$g(t)=\frac{1}{C}\cos (\frac{1}{\sqrt{LC}}t)$$
输入为交流电流源
$$I_s(t)=I_s\sin (wt)<M_1$$
$$I_s(s)=I_s\frac{w}{s^2+w^2}$$
$$y(s)=I_s\frac{w}{s^2+w^2}\frac{1}{c}\frac{s}{s^2+\frac{1}{CL}}$$

当$w^2=\frac{1}{CL}$时，反拉氏变换
$$y(t)=-\frac{I_s}{2}t\sin \frac{1}{\sqrt{LC}}t$$

当$t\to \infty$的时候， $y(t)=\frac{I_s}{2}t\sin \frac{1}{\sqrt{LC}}t$为振幅趋向于无穷的正弦震荡，因此输出无界，满足题目要求，系统不再BIBO稳定.
$L=1H,C=1C,I_s=\sin t$，仿真的电路图为：

图10 交流电流源 LC并联电路仿真的电路图

仿真时间是600s,示波器输出的结果为：

图11　交流电流源 LC并联电路输入电压输出

从图中观察，$y(t)$输出无界，和上面的公式分析一致，此时有界输入能产生无界的输出，输入信号是正弦信号，频率满足$w=1/\sqrt{LC}$