文档介绍：实验三 连续时间LTI系统的时域分析
一、实验目的
1．学会用MATLAB求解连续系统的零状态响应；
2. 学会用MATLAB求解冲激响应及阶跃响应；
3．学会用MATLAB实现连续信号卷积的方法；
二、实验原理
1．连续时间系统零状态响应的数值计算
我们知道，LTI连续系统可用如下所示的线性常系数微分方程来描述，
在MATLAB中，控制系统工具箱提供了一个用于求解零初始条件微分方程数值解的函数lsim。其调用格式
y=lsim(sys,f,t)
式中，t表示计算系统响应的抽样点向量，f是系统输入信号向量，sys是LTI系统模型，用来表示微分方程，差分方程或状态方程。其调用格式
sys=tf(b,a)
式中，b和a分别是微分方程的右端和左端系数向量。例如，对于以下方程:
可用 获得其LTI模型。
注意，如果微分方程的左端或右端表达式中有缺项，则其向量a或b中的对应元素应为零，不能省略不写，否则出错。
例3-1 已知某LTI系统的微分方程为
y’’(t)+ 2y’(t)+100y(t)=f(t)
其中，，求系统的输出y(t).
解：显然，这是一个求系统零状态响应的问题。其MATLAB计算程序如下：
ts=0;te=5;dt=;
sys=tf([1],[1,2,100]);
t=ts:dt:te;
f=10*sin(2*pi*t);
y=lsim(sys,f,t);
plot(t,y);
xlabel('Time(sec)');
ylabel('y(t)');
2．连续时间系统冲激响应和阶跃响应的求解
在MATLAB中，对于连续LTI系统的冲激响应和阶跃响应，可分别用控制系统工具箱提供的函数impluse和step来求解。其调用格式为
y=impluse(sys,t)
y=step(sys,t)
式中，t表示计算系统响应的抽样点向量，sys是LTI系统模型。
例3-2已知某LTI系统的微分方程为
y’’(t)+ 2y’(t)+100y(t)=10f(t)
求系统的冲激响应和阶跃响应的波形.
解：ts=0;te=5;dt=;
sys=tf([10],[1,2,100]);
t=ts:dt:te;
h=impulse(sys,t);
figure;
plot(t,h);
xlabel('Time(sec)');
ylabel('h(t)');
g=step(sys,t);
figure;
plot(t,g);
xlabel('Time(sec)');
ylabel('g(t)');
3. 用MATLAB实现连续时间信号的卷积
信号的卷积运算有符号算法和数值算法，此处采用数值计算法，需调用MATLAB的conv( )函数近似计算信号的卷积积分。连续信号的卷积积分定义是
如果对连续信号和进行等时间间隔均匀抽样，则和分别变为离散时间信号和。其中，为整数。当足够小时，和既为连续时间信号和。因此连续时间信号卷积积分可表示为
采用数值计算时，只求当时卷积积分的值，其中，n为整数，既
其中，实际就是离散序列和的卷积和。当足够小时，序列就是连续信号的数值近似，既
上式表明，连续信号和的卷积，可用各自抽样后的离散时间序列的卷积再乘以抽样间隔。抽样间隔越小，误差越小。
例3-3用数值计算法求与的卷积积分。
解：因为是一个持续时间无限长的信号，而计算机数值计算不可能计算真正的无限长信号，所以在进行的抽样离散化时，所取的时间围让衰减到足够小就可以了，本例取。程序是
dt=; t=-1:dt:;
f1=(t>=0);
f2=(t>=2);
f3=f1-f2
f4=exp(-3*t).*f1;
f=conv(f3,f4)*dt
n=length(f); tt=(0:n-1)*dt-2;
subplot(221), plot(t,f3), grid on;
axis([-1,,-,]); title('f1(t)'); xlabel('t')
subplot(222), plot(t,f4), grid on;
axis([-1,,-,]); title('f2(t)'); xlabel('t')
subplot(223), plot(tt,f), grid on;
title('f(t)=f1(t)*f2(t)'); xl