信号与系统实验二：连续时间LTI系统的时域分析

一、实验目的

1、学会使用符号法求解连续系统的零输入响应和零状态响应；

2、学会使用符号法求解连续系统的完全响应。

二、实验原理及实例分析

连续时间系统零输入响应和零状态响应的符号求解连续时间系统可以使用常系数微分方程来描述，其完全响应由零输入响应和零状态响应组成。MATLAB符号工具箱提供了 dsolve 函数，可以实现对常系数微分方程的符号求解，其调用格式为：

dsolve(‘eq1,eq2…’,’cond1,cond2,…’,’v’)其中参数eq表示各个微分方程，它与MATLAB符号表达式的输入基本相同，微分和导数的输入是使用 Dy，D2y，D3y 来表示y的一价导数，二阶导数，三阶导数；参数cond表示初始条件或者起始条件；参数v表示自变量，默认是变量t。通过使用dsolve函数可以求出系统微分方程的零输入响应和零状态响应，进而求出完全响应。

三、实验内容

1、求解[实例1]：

试用MATLAB命令求齐次微分方程的零输入响应，已知起始条件为

实验代码：

eq='D3y + 2*D2y + Dy = 0';

cond = 'y(0)=1, Dy(0)=1, D2y(0)=2';

yt = dsolve(eq, cond); simplify(yt);

实验结果：

>> ans

ans = 5 - 3*t*exp(-t) - 4*exp(-t)

2、求解[实例2]：

已知输入,试用MATLAB命令求解微分方程的零状态响应。

实验代码：

eq1='D3y + 4*D2y + 8*Dy = 3*Dx + 8*x';

eq2='x=heaviside(t)';

cond='y(-0.01)=0, Dy(-0.01)=0, D2y(-0.01)=0';

yt=dsolve(eq1, eq2, cond); simplify(yt.y);

实验结果：

 >> ans

ans = t/2 + abs(t)/2 + (3*heaviside(t))/8 - (3*cos(2*t)*exp(-2*t)*heaviside(t))/8 -    (3*sin(2*t)*exp(-2*t)*heaviside(t))/8 

3、求解零状态响应和零输入响应

已知系统的微分方程和激励信号如下，使用MATLAB画出系统的零状态响应和零输入响应。

实验代码：

%% 零输入响应：   
eq='D2r + 4*Dr + 4*r = 0';

cond='r(0)=0, Dr(0)=1';

rzi=dsolve(eq, cond);simplify(rzi);


%% 零状态响应：   
syms s t;

xt=exp(-1*t)*heaviside(t);

Hs=(s+3)/(s^2+4*s+4);

Xs=laplace(xt);  Ys=Hs*Xs;

yt=ilaplace(Ys);


%% 图像代码：       
subplot(2, 1 ,1);

eq='D2r + 4*Dr + 4*r = 0';

cond='r(0)=0, Dr(0)=1';

rzi=dsolve(eq, cond);simplify(rzi);

fplot(rzi, [-1, 10]);xlabel('Time');ylabel('rzi(t)');

legend('t*exp(-2*t)'); grid on; axis([-1, 8, -3, 0.5]);

subplot(2, 1 ,2)

sys=tf([1, 3],[1, 4, 4]);

t=-1:0.01:10; f=exp(-t).*heaviside(t); y=lsim(sys,f,t);

plot(t,y);xlabel('Time');ylabel('rzs(t)');

legend('(2*exp(t)-2-t)exp(-2*t)'); grid on; axis([-1, 8, 0, 0.35]);

suptitle('r(t)的零输入和零状态响应');

实验图像：

四、实验总结

最后一道练习题若采用实验指导所示方法求系统的零状态响应，计算结果为(exp(-2*t)*(t + 2*exp(t))*(sign(t) + 1))/2，即，这与用拉普拉斯变换的方式计算所得结果，笔算结果与后者相符，因此采用拉普拉斯变换的方式。

总的说来，此次使用MATLAB求解连续系统的零输入响应和零状态响应，我再次熟悉了求解各个相应地过程，只有MATLAB的使用者本身熟练掌握了零输入和零状态响应的相关概念，如初值、各个响应的含义，才能对程序显示的结果做出理性判断。同时，平时求这两个响应并未对图像有过深度思考，这一点也正是实验的目的。