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

声明：华南理工大学 线性系统作业 苏老师 本人特别菜，不研究相关的方向，差点挂科，这个作业的内容仅供交流。

第二次作业

题目

针对一个两阶线性系统，分析其零输入响应的特点，给出一组解答的基底，用仿真说明其特点。

解答思路

对一个具体的二阶的电路系统，求解零输入响应，给定一组基底，分析零输入响应的特点，利用MATLAB仿真验证零输入响应的特点。

系统的描述

设各组成元件的参数值为已知，取电源电压$e(t)$为输入变量，将电容电压$u_c$和电容电流$i_L$作为输出。

图1 电路原理图

计算

回路方程的计算：将电容电压$u_c$和电感电流$i_L$作为状态变量，根据回路基尔霍夫定律，得到
$$u_c+R_{2}C\frac{d{u}_c}{dt}=L\frac{d{i}_L}{dt}$$
$$(C\frac{d{u}_c}{dt}+i_L)R_1+L\frac{d{i}_L}{dt}=e(t)$$

转化为状态方程和输出方程的形式为：
$${\dot{u}}_c=(-\frac{1}{(R_2+R_1)C})u_c+(-\frac{R_1}{(R_2+R_1)C})i_L+(\frac{1}{(R_1+R_2)C})e(t)$$
$${\dot{i}}_L=(\frac{R_1}{(R_1+R_2)L})u_c+(-\frac{R_1{R}_2}{(R_1+R_2)L})i_L+(\frac{R_2}{(R_1+R_2)L})e(t)$$
引入参数矩阵的表示：

状态矩阵和输出方程可以表达为：

根据状态方程求解该系统的零输入响应，为了方便求解我们设 $R_1=200k\Omega$， $R_2=800k\Omega$，$C=10^{-6}F$，$L=1\times 10^{-6}H$
所以

两个特征值${\lambda }_1=-1.25,{\lambda }_2=-1.6\times 10^{11}$对应的特征向量分别为

将初始状态分解

我们将

作为一组基底去表示初始状态，初始状态可以用$v_1$ , $v_2$线性表示，系数分别为${\alpha }_1=(u_c(t_0)-1.25\times 10^{-6}{i}_L(t_0))$ ，${\alpha }_2=(-1.25\times 10^{-6}{u}_c(t_0)+i_L(t_0))$
将零输入响应分解：

从上述分解的结果中可以得到，零输入响应的结果是:
$$x(t)={\alpha }_1{e}^{{\lambda }_{1}t}{v}_1+{\alpha }_2{e}^{{\lambda }_{2}t}{v}_2$$

分析可以得到两阶线性系统的零输入响应特点
当初始状态可以用$x(t_0)={\alpha }_1{v}_1+{\alpha }_2{v}_2$ ，$A$对应的特征向量$v_1$ 和 $v_2$，对应的特征根${\lambda }_1,{\lambda }_2$ ，零输入响应可以表示为$x(t)={\alpha }_1{e}^{{\lambda }_{1}t}{v}_1+{\alpha }_2{e}^{{\lambda }_{2}t}{v}_2$。
设初始电压$u_c=100000V$ ，$i_L=0.01A$ ，$t_0=0$ ,图中直接用实线表示
设初始电压$u_c=300000V$，$i_L=0.03A$ ，$t_0=0$ ,图中用x表示
设初始电压$u_c=500000V$ ，$i_L=0.05A$ ，$t_0=0$ ，图中用。表示
带入计算响应，画图可以得到：

图2 电感电流计算结果画图

图3 电容电压计算结果画图

电路仿真

设初始电压$u_c=100000V$ ，$i_L=0.01A$ ，$t_0=0$ ， 进行MATLAB仿真，选择仿真时间7s。设初始电压$u_c=300000V$，$i_L=0.03A$ ，$t_0=0$ ， 进行MATLAB仿真，选择仿真时间7s。设初始电压 $u_c=500000V$ ，$i_L=0.05A$ ，$t_0=0$ ， 进行MATLAB仿真，选择仿真时间7s。

图4 MTALB电路仿真内部设计图

图5 MATLAB仿真总图

仿真的结果如下：

图6 电容两端电压仿真结果

图7 电感电流仿真结果

从图中可以看出零输入相应输出和初始状态成正比和计算的结果相符，从示波器的仿真结果也可以看出仿真的结果和计算画图的结果基本一致。
验证：
取仿真$t=1s$时，从仿真示波器结果中读取
$i_{L\_1}(1)=0.036759$，$i_{L\_1}(1)=0.10645$ ，$i_{L\_1}(1)=0.17786$
$u_{c\_1}(1)=28419.1483$ , $u_{c\_2}(1)=85367.6259$，$u_{c\_3}(1)=142885.5883$
将数值代入验证这些数值都可以用基底线性表示。
例如：

因此电路仿真也基本可以验证二阶线性系统零输入响应的特点，当初始状态可以用$x(t_0)={\alpha }_1{v}_1+{\alpha }_2{v}_2$， 对应的特征向量$v_1$和$v_2$ ，对应的特征根${\lambda }_1,{\lambda }_2$ ，零输入响应可以表示为 $x(t)={\alpha }_1{e}^{{\lambda }_{1}t}{v}_1+{\alpha }_2{e}^{{\lambda }_{2}t}{v}_2$

附录

%%%计算响应画图的代码%%%%%
figure
x=linspace(0,7,100);
y1=(exp(-1.25x)(110^5-1.25*10-60.01)+1.2510^{-6*exp(-1.6*10}11x)(-1.2510^-6*1*105+0.01));
y2=(exp(-1.25x)(510^5-1.25*10-60.05)+1.2510^{-6*exp(-1.6*10}11x)(-1.2510^-6*5*105+0.05));
y3=(exp(-1.25x)(310^5-1.25*10-60.03)+1.2510^{-6*exp(-1.6*10}11x)(-1.2510^-6*3*105+0.03));
plot(x,y1,‘b:’,x,y2,’-o’,x,y3,’-rx’);
xlabel(‘t/s’)
ylabel(‘U_c/V’)
figure
x=linspace(0,7,100);
y1=1.2510^{-6*(exp(-1.25*x)*(1*10}5-1.2510^{-6*0.01)+exp(-1.6*10}11x)(-1.2510^-6*1*105+0.01));
y2=1.2510^{-6*(exp(-1.25*x)*(5*10}5-1.2510^{-6*0.05)+exp(-1.6*10}11x)(-1.2510^-6*5*105+0.05));
y3=1.2510^{-6*(exp(-1.25*x)*(3*10}5-1.2510^{-6*0.03)+exp(-1.6*10}11x)(-1.25*10^-6*3*105+0.03));
plot(x,y1,‘b:’,x,y2,’-o’,x,y3,’-rx’);
xlabel(‘t/s’)
ylabel(‘i_L/A’)