《实验2LTI系统的时域分析》由会员分享，可在线阅读，更多相关《实验2LTI系统的时域分析（11页珍藏版）》请在金锄头文库上搜索。
1、实验2 LTI系统的时域分析一、实验目的1. 掌握利用MATLAB对系统进行时域分析的方法。2. 掌握连续时间系统零状态相应、冲击响应和阶跃响应的求解方法。3. 掌握求解离散时间系统响应、单位抽样响应的方法。4.加深对卷积积分和卷积和的理解。掌握利用计算机进行卷积积分和卷积和计算的方法。二、实验原理与方法1. 连续时间系统时域分析的MATLAB实现1）连续时间系统的MATLAB表示LTI连续系统通常可以由系统微分方程描述，设描述系统的微分方程为：则在MATLAB中，可以建立系统模型如下：；其中，tf是用于创建系统模型的函数，向量a和b是以微分方程求导的降幂次序来排列的，如果有缺项，应用0补齐。
2、2）连续时间系统的零状态响应零状态响应指系统的初始状态为零，仅由输入信号所引起的响应。MATLAB提供了一个用于求解零状态响应的函数lism，其调用格式如下：lsim（sys,x,t）绘出输入信号及响应的波形，x和t表示输入信号数值向量及其时间向量。y=lsim(sys,x,t)这种调用格式不绘出波形，而是返回响应的数值向量。3）连续时间系统的冲激响应与阶跃响应MATLAB提供了impulse函数来求指定时间范围内由模型sys描述的连续时间系统的单位冲激响应。impulse函数的基本调用格式如下：impulse(sys)在默认的时间范围内绘出系统冲激响应的时域波形；impulse(sys,T)
3、绘出系统在0T范围内冲激响应的时域波形；impulse(sys,ts:tp:te)绘出系统在tste范围内，以tp为时间间隔取样的冲激响应波形。y,t=impulse()这种调用格式不绘出冲激响应波形，而是返回冲激响应的数值向量及其对应的时间向量。函数step用于求解单位阶跃响应，函数step同样也有如下几种调用格式：step(sys)；step(sys,T)；step(sys,ts:tp:te)；y,t=step()。各种调用格式参数所代表的意思可参考上述impulse函数。2.离散时间系统时域分析的MATLAB实现1）离散时间系统的MATLAB表示LTI离散系统通常可以由系统差分方程描述，
4、设描述系统的差分方程为：则在MATLAB里，我们可以用如下两个向量来表示这个系统：；2）离散时间系统对任意输入的响应MATLAB提供了求LTI离散系统响应的专用函数filter，该函数用于求取由差分方程所描述的离散时间系统在执行时间范围内对输入序列所产生的响应，该函数基本调用格式为y=filter(b,a,x)其中，x为输入序列，y为输出序列，输出序列y对应的时间区间与x对应的时间区间相同。3）离散时间系统的单位抽样响应MATLAB提供了函数impz来求指定时间范围内，由向量b和a描述的离散时间系统的单位抽样响应，具体调用格式如下：impz(b,a)在默认的时间范围内绘出系统单位抽样响应的时域
5、波形；impz(b,a,N)绘出系统在0N时间范围内单位抽样响应的时域波形；impz(b,a,ns:ne)绘出系统在nsne范围内单位抽样响应的时域波形；y,t=impz()这种调用格式不绘出单位抽样响应波形，而是返回单位抽样响应的数值向量及其对应的时间向量。3卷积和与卷积积分1）离散时间序列的卷积和卷积和是离散系统时域分析的基本方法之一，离散时间序列的卷积和定义如下：对于离散LTI系统，设其输入信号为，单位抽样响应为，则其零状态响应为即离散LTI系统的零状态响应可以表示出输入信号与单位抽样响应的卷积。因此，离散时间序列的卷积和计算对于我们进行离散系统时域分析具有重要的意义。MATLAB的co
6、nv函数可以用来求两个离散序列的卷积和，调用格式为x=conv(x1,x2)。2）连续时间信号的卷积积分卷积积分是连续系统时域分析的有效方法和工具之一，连续时间信号的卷积积分定义如下：对于连续LTI系统，设其输入信号为，单位抽样响应为，其零状态响应为，则有即连续LTI系统的零状态响应可以表示为输入信号与单位冲激响应的卷积。因此，连续时间信号卷积积分对连续系统的时域分析具有非常重要的意义。利用MATLAB可以采用数值计算的方法近似计算卷积积分。卷积积分可用求和运算来实现现在考虑只求时，则有当足够小，的数值近似。我们可以利用计算离散卷积序列卷积和的conv函数来计算卷积积分，具体步骤如下：将连续时
7、间信号和以时间间隔进行取样，得到离散序列和；构造离散序列和对应的时间向量和；调用函数conv计算卷积积分在时的近似采样值；构造离散序列对应的时间向量。三、实验内容（1）已知描述模拟低通、高通、带通和带阻滤波器的微分方程如下，试采用MATLAB绘出各系统的单位冲激响应和单位阶跃响应。实验过程与实验结果：matlab程序如下：b=1;a=1 sqrt(2) 1;sys=tf(b,a);subplot(121);impulse(sys);subplot(122);step(sys);上述程序的运行结果为：实验过程与实验结果：matlab程序如下：b=1 0 0;a=1 sqrt(2) 1;sys=t
8、f(b,a);subplot(121);impulse(sys);subplot(122);step(sys);上述程序的运行结果为：实验过程与实验结果：matlab程序如下：b=1 0;a=1 1 1;sys=tf(b,a);subplot(121);impulse(sys);subplot(122);step(sys);上述程序的运行结果为：实验过程与实验结果：matlab程序如下：b=1 1 0;a=1 1 1;sys=tf(b,a);subplot(121);impulse(sys);subplot(122);step(sys);上述程序的运行结果为：（2）已知某系统可以由如下微分方程
9、描述请利用MATLAB绘出该系统单位冲激响应和单位阶跃响应的时域波形；实验过程与实验结果：matlab程序如下：b=1;a=1 1 6;sys=tf(b,a);subplot(121);impulse(sys);subplot(122);step(sys);上述程序的运行结果为：根据冲激响应的时域波形分析系统的稳定性；由波形图可知，t趋于无穷时，y(t)均趋于0，所以该系统是稳定的。如果系统的输入为，求系统的零状态响应。实验过程与实验结果：matlab程序如下：b=1;a=1 1 6;sys=tf(b,a);t=0:0.01:10;x=exp(-t);lsim(sys,x,t);上述程序的运行
10、结果为：（3）已知描述离散系统的微分方程如下，试采用MATLAB绘出各系统的单位抽样响应，并根据单位抽样响应的时域波形分析系统的稳定性。实验过程与实验结果：matlab程序如下：b=1;a=1 3 2;impz(b,a,0:15);上述程序的运行结果为：分析：由图像可知，n趋于无穷时，y(n)趋于无穷，所以该系统是不稳定的。实验过程与实验结果：matlab程序如下：b=1 -3;a=1 -0.5 0.8;impz(b,a,0:10);上述程序的运行结果为：分析：由图像可知，n趋于无穷时，y(n)有界，且趋于0，所以该系统是稳定的。（4）已知系统可以由如下差分方程描述试采用MATLAB绘出该系统
11、的单位抽样响应波形和单位阶跃响应波形。实验过程与实验结果：matlab程序如下：clear;b=1;a=1 1 0.25;subplot(121);impz(b,a,0:10);subplot(122);stepz(b,a);上述程序的运行结果为：（5）采用MATLAB计算如下两个序列的卷积，并绘出图形 实验过程与实验结果：matlab程序如下：x1=1 2 1 1;x2=1 1 1 1 1;x=conv(x1,x2)n1=-3:4;stem(n1,x,filled);上述程序的运行结果为：（6）已知某LTI离散时间系统，其单位抽样响应，系统的输入为，计算当时系统的零状态响应，绘出，和时域波形
12、。实验过程与实验结果：matlab程序如下：n=0:40;n1=0:80;x=sin(0.2*n).*heaviside(n);h=sin(0.5*n).*heaviside(n);y=conv(x,h)subplot(221);stem(n,x,filled);xlabel(n);title(x(n);subplot(222);stem(n,h,filled);xlabel(n);title(h(n);subplot(212);stem(n1,y,filled);xlabel(n);title(y(n);上述程序的运行结果为：（7）已知两个连续时间信号,其中为幅值为2，门宽为2的门函数， 为
13、幅值为1，门宽为4的门函数，试采用MATLAB求这两个信号的卷积。实验过程与实验结果：matlab程序如下：t=-10:10;dt=0.01;t1=-1:0.01:1;t2=-2:0.01:2;x1=2.*rectpuls(t1,2);x2=rectpuls(t2,4);x,t=sconv(x1,x2,t1,t2,dt);plot(t,x);xlabel(t);title(x1*x2);注：sconv为自行给出的计算卷积积分的通用函数，与书上169页8-10前的函数内容相同，在此不再列出。上述程序的运行结果为：四、体会和建议这次的实验内容比较简单，书上的代码可以看懂，我也模糊知道MATLAB里面的数是以矩阵的方式存储的，所以在对方程分析时十分快捷。这次实验也让我掌握利用MATLAB对系统进行时域分析的方法，学会如何利用MATLAB求解连续时间系统零状态相应、冲击响应和阶跃响应、离散时间系统响应、单位抽样响应。同时对于卷积这一重要概念和方法进一步熟悉和理解，学会如何利用计算机进行卷积积分和卷积和计算，为以后的学习和分析提供了有力的工具保障。
