实验2lti系统的时域分析

1、-. z.实验2 LTI系统的时域分析一、实验目的1. 掌握利用MATLAB对系统进展时域分析的方法。2. 掌握连续时间系统零状态相应、冲击响应和阶跃响应的求解方法。3. 掌握求解离散时间系统响应、单位抽样响应的方法。4.加深对卷积积分和卷积和的理解。掌握利用计算机进展卷积积分和卷积和计算的方法。二、实验原理与方法1. 连续时间系统时域分析的MATLAB实现1连续时间系统的MATLAB表示LTI连续系统通常可以由系统微分方程描述，设描述系统的微分方程为：则在MATLAB中，可以建立系统模型如下： QUOTE ； QUOTE ； QUOTE ；其中，tf是用于创立系统模型的函数，向量a和b是以微

2、分方程求导的降幂次序来排列的，如果有缺项，应用0补齐。2连续时间系统的零状态响应零状态响应指系统的初始状态为零，仅由输入信号所引起的响应。MATLAB提供了一个用于求解零状态响应的函数lism，其调用格式如下：lsimsys,*,t绘出输入信号及响应的波形，*和t表示输入信号数值向量及其时间向量。y=lsim(sys,*,t)这种调用格式不绘出波形，而是返回响应的数值向量。3连续时间系统的冲激响应与阶跃响应MATLAB提供了impulse函数来求指定时间围由模型sys描述的连续时间系统的单位冲激响应。impulse函数的根本调用格式如下：impulse(sys)在默认的时间围绘出系统冲激响应的

3、时域波形；impulse(sys,T)绘出系统在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里，我们可以用如下两个向量来表示这个系统：QUOTE ； QUOTE ；2离散时间系统对任意输入的响应MATLAB提供了求LTI离散系统响应的专用函数filter，该函数用于求取由差分方程所描述的离散时间系统在执行时间围对输入序列所产生的响应，该函数根本调用格式为y=filter(b,a,*)其中，*为输入序列，y为输出序列，输出序列y对应的时间区间与*对应的时间区间一样。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离散时间序列的卷积和卷积和是离散系统时域分析的根本方法之一，离散时间序列的卷积和 QUOTE 定义如下：对于离散LTI系统，设其输入信号为 QUOTE ，单位抽样响应为，则其零状态响应 QUOTE 为即离散LTI系统的零状态响应可以表示出输入信号 QUOTE 与单位抽样响应的卷积。因此，

6、离散时间序列的卷积和计算对于我们进展离散系统时域分析具有重要的意义。MATLAB的conv函数可以用来求两个离散序列的卷积和，调用格式为*=conv(*1,*2)。2连续时间信号的卷积积分卷积积分是连续系统时域分析的有效方法和工具之一，连续时间信号的卷积积分 QUOTE 定义如下：对于连续LTI系统，设其输入信号为 QUOTE ，单位抽样响应为，其零状态响应为 QUOTE ，则有即连续LTI系统的零状态响应可以表示为输入信号与单位冲激响应的卷积。因此，连续时间信号卷积积分对连续系统的时域分析具有非常重要的意义。利用MATLAB可以采用数值计算的方法近似计算卷积积分。卷积积分可用求和运算来实现现

7、在考虑只求 QUOTE 时 QUOTE ，则有当 QUOTE 足够小， QUOTE 的数值近似。我们可以利用计算离散卷积序列卷积和的conv函数来计算卷积积分，具体步骤如下：将连续时间信号 QUOTE 和 QUOTE 以时间间隔进展取样，得到离散序列 QUOTE 和 QUOTE ；构造离散序列 QUOTE 和 QUOTE 对应的时间向量 QUOTE 和 QUOTE ；调用函数conv计算卷积积分在 QUOTE 时的近似采样值 QUOTE ；构造离散序列 QUOTE 对应的时间向量 QUOTE 。三、实验容1描述模拟低通、高通、带通和带阻滤波器的微分方程如下，试采用MATLAB绘出各系统的单位冲

8、激响应和单位阶跃响应。实验过程与实验结果：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=tf(b,a);subplot(121);impulse(sys);subplot(122);step(sys);上述程序的运行结果为：实验过程与实验结果：matlab程序如下：b=1 0;a=1 1 1;sys=tf(b,a);subplot(121);imp

9、ulse(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*系统可以由如下微分方程描述请利用MATLAB绘出该系统单位冲激响应和单位阶跃响应的时域波形；实验过程与实验结果：matlab程序如下：b=1;a=1 1 6;sys=tf(b,a);subplot(121);impulse(sys);subplot(122);step(sys);上述程序的

10、运行结果为：根据冲激响应的时域波形分析系统的稳定性；由波形图可知，t趋于无穷时，y(t)均趋于0，所以该系统是稳定的。如果系统的输入为，求系统的零状态响应。实验过程与实验结果：matlab程序如下：b=1;a=1 1 6;sys=tf(b,a);t=0:0.01:10;*=e*p(-t);lsim(sys,*,t);上述程序的运行结果为：3描述离散系统的微分方程如下，试采用MATLAB绘出各系统的单位抽样响应，并根据单位抽样响应的时域波形分析系统的稳定性。实验过程与实验结果：matlab程序如下：b=1;a=1 3 2;impz(b,a,0:15);上述程序的运行结果为：分析：由图像可知，n趋

11、于无穷时，y(n)趋于无穷，所以该系统是不稳定的。实验过程与实验结果：matlab程序如下：b=1 -3;a=1 -0.5 0.8;impz(b,a,0:10);上述程序的运行结果为：分析：由图像可知，n趋于无穷时，y(n)有界，且趋于0，所以该系统是稳定的。4系统可以由如下差分方程描述试采用MATLAB绘出该系统的单位抽样响应波形和单位阶跃响应波形。实验过程与实验结果：matlab程序如下：clear;b=1;a=1 1 0.25;subplot(121);impz(b,a,0:10);subplot(122);stepz(b,a);上述程序的运行结果为：5采用MATLAB计算如下两个序列的

12、卷积，并绘出图形实验过程与实验结果：matlab程序如下：*1=1 2 1 1;*2=1 1 1 1 1;*=conv(*1,*2)n1=-3:4;stem(n1,*,filled);上述程序的运行结果为：6*LTI离散时间系统，其单位抽样响应，系统的输入为，计算当时系统的零状态响应，绘出，和时域波形。实验过程与实验结果：matlab程序如下：n=0:40;n1=0:80;*=sin(0.2*n).*heaviside(n);h=sin(0.5*n).*heaviside(n);y=conv(*,h)subplot(221);stem(n,*,filled);*label(n);title(*

13、(n);subplot(222);stem(n,h,filled);*label(n);title(h(n);subplot(212);stem(n1,y,filled);*label(n);title(y(n);上述程序的运行结果为：7两个连续时间信号,其中为幅值为2，门宽为2的门函数，为幅值为1，门宽为4的门函数，试采用MATLAB求这两个信号的卷积。实验过程与实验结果：matlab程序如下：t=-10:10;dt=0.01;t1=-1:0.01:1;t2=-2:0.01:2;*1=2.*rectpuls(t1,2);*2=rectpuls(t2,4);*,t=sconv(*1,*2,t1,t2,dt);plot(t,*);*label(t);title(*1*2);注：sconv为自行给出的计算卷积积分的通用函数，与书上169页8-10前的函数容一样，在此不再列出。上述程序的运行结果为：四、体会和建议