实验2 LTI系统的时域分析

[键入文字]实验2LTI系统的时域分析一、实验目的1.掌握利用MATLAB对系统进行时域分析的方法。2.掌握连续时间系统零状态相应、冲击响应和阶跃响应的求解方法。3.掌握求解离散时间系统响应、单位抽样响应的方法。4.加深对卷积积分和卷积和的理解。掌握利用计算机进行卷积积分和卷积和计算的方法。二、实验原理与方法1.连续时间系统时域分析的MATLAB实现1）连续时间系统的MATLAB表示LTI连续系统通常可以由系统微分方程描述，设描述系统的微分方程为：𝑎𝑁𝑦(𝑁)(𝑡)+𝑎𝑁−1𝑦𝑁−1(𝑡)+⋯+𝑎0𝑦(𝑡)=𝑏𝑀𝑥(𝑀)(𝑡)+𝑏𝑀−1𝑥(𝑀−1)(𝑡)+⋯+𝑏0𝑥(𝑡)则在MATLAB中，可以建立系统模型如下：𝑏=[𝑏𝑀,𝑏𝑀−1,⋯,𝑏0]；𝑎=[𝑎𝑁,𝑎𝑁−1,⋯,𝑎0]；𝑠𝑦𝑠=𝑡𝑓(𝑏,𝑎)；其中，tf是用于创建系统模型的函数，向量a和b是以微分方程求导的降幂次序来排列的，如果有缺项，应用0补齐。2）连续时间系统的零状态响应零状态响应指系统的初始状态为零，仅由输入信号所引起的响应。MATLAB提供了一个用于求解零状态响应的函数lism，其调用格式如下：lsim（sys,x,t）绘出输入信号及响应的波形，x和t表示输入信号数值向量及其时间向量。y=lsim(sys,x,t)这种调用格式不绘出波形，而是返回响应的数值向量。3）连续时间系统的冲激响应与阶跃响应MATLAB提供了impulse函数来求指定时间范围内由模型sys描述的连续时间系统的单位冲激响应。impulse函数的基本调用格式如下：impulse(sys)在默认的时间范围内绘出系统冲激响应的时域波形；impulse(sys,T)绘出系统在0～T范围内冲激响应的时域波形；impulse(sys,ts:tp:te)绘出系统在ts～te范围内，以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离散系统通常可以由系统差分方程描述，设描述系统的差分方程为：𝑎0𝑦(𝑛)+𝑎1𝑦(𝑛−1)+⋯+𝑎𝑁𝑦(𝑛−𝑁)=𝑏0𝑥(𝑚)+𝑏1𝑥(𝑚−1)+⋯+𝑏𝑀𝑥(𝑚−𝑀)则在MATLAB里，我们可以用如下两个向量来表示这个系统：𝑏=[𝑏0,𝑏1,…,𝑏𝑀]；𝑎=[𝑎0,𝑎1,…,𝑎𝑁]；2）离散时间系统对任意输入的响应MATLAB提供了求LTI离散系统响应的专用函数filter，该函数用于求取由差分方程所描述的离散时间系统在执行时间范围内对输入序列所产生的响应，该函数基本调用格式为y=filter(b,a,x)其中，x为输入序列，y为输出序列，输出序列y对应的时间区间与x对应的时间区间相同。3）离散时间系统的单位抽样响应MATLAB提供了函数impz来求指定时间范围内，由向量b和a描述的离散时间系统的单位抽样响应，具体调用格式如下：impz(b,a)在默认的时间范围内绘出系统单位抽样响应的时域波形；impz(b,a,N)绘出系统在0～N时间范围内单位抽样响应的时域波形；impz(b,a,ns:ne)绘出系统在ns～ne范围内单位抽样响应的时域波形；[y,t]=impz(…)这种调用格式不绘出单位抽样响应波形，而是返回单位抽样响应的数值向量及其对应的时间向量。3．卷积和与卷积积分1）离散时间序列的卷积和卷积和是离散系统时域分析的基本方法之一，离散时间序列𝑥1(𝑛)和𝑥2(𝑛)的卷积和𝑥(𝑛)定义如下：𝑥(𝑛)=𝑥1(𝑛)∗𝑥2(𝑛)=∑𝑥1(𝑘)𝑥2(𝑛−𝑘)+∞𝑘=−∞对于离散LTI系统，设其输入信号为𝑥(𝑛)，单位抽样响应为ℎ(𝑛)，则其零状态响应𝑦(𝑛)为𝑦(𝑛)=𝑥(𝑛)∗ℎ(𝑛)即离散LTI系统的零状态响应可以表示出输入信号𝑥(𝑛)与单位抽样响应ℎ(𝑛)的卷积。因此，离散时间序列的卷积和计算对于我们进行离散系统时域分析具有重要的意义。MATLAB的conv函数可以用来求两个离散序列的卷积和，调用格式为x=conv(x1,x2)。2）连续时间信号的卷积积分卷积积分是连续系统时域分析的有效方法和工具之一，连续时间信号𝑥1(𝑡)和𝑥2(𝑡)的卷积积分𝑥(𝑡)定义如下：𝑥(𝑡)=𝑥1(𝑡)∗𝑥2(𝑡)=∫𝑥1(𝜏)𝑥2(𝑡−𝜏)𝑑𝜏+∞−∞对于连续LTI系统，设其输入信号为𝑥(𝑡)，单位抽样响应为ℎ(𝑡)，其零状态响应为𝑦(𝑡)，则有𝑦(𝑡)=𝑥(𝑡)∗ℎ(𝑡)[键入文字]即连续LTI系统的零状态响应可以表示为输入信号与单位冲激响应的卷积。因此，连续时间信号卷积积分对连续系统的时域分析具有非常重要的意义。利用MATLAB可以采用数值计算的方法近似计算卷积积分。卷积积分可用求和运算来实现𝑥(𝑡)=𝑥1(𝑡)∗𝑥2(𝑡)=∫𝑥1(𝜏)𝑥2(𝑡−𝜏)𝑑𝜏+∞−∞=lim∆→0∑𝑥1(𝑘∆)𝑥2(𝑡−𝑘∆)∙∆+∞𝑘=−∞现在考虑只求𝑡=𝑛∆时𝑥(𝑡)的值𝑥(𝑛∆)，则有𝑥(𝑛∆)=∆∙∑𝑥1(𝑘∆)𝑥2[(𝑛−𝑘)∆]+∞𝑘=−∞当∆足够小，𝑥(𝑛∆)就是𝑥(𝑡)的数值近似。我们可以利用计算离散卷积序列卷积和的conv函数来计算卷积积分，具体步骤如下：①将连续时间信号𝑥1(𝑡)和𝑥2(𝑡)以时间间隔Δ进行取样，得到离散序列𝑥1(𝑛∆)和𝑥2(𝑛∆)；②构造离散序列𝑥1(𝑡)和𝑥2(𝑡)对应的时间向量t1和t2；③调用函数conv计算卷积积分在𝑡=𝑛∆时的近似采样值𝑥(𝑛∆)；④构造离散序列𝑥(𝑛∆)对应的时间向量𝑛。三、实验内容（1）已知描述模拟低通、高通、带通和带阻滤波器的微分方程如下，试采用MATLAB绘出各系统的单位冲激响应和单位阶跃响应。①)()()(2)(txtytyty实验过程与实验结果：matlab程序如下：b=[1];a=[1sqrt(2)1];sys=tf(b,a);subplot(121);impulse(sys);subplot(122);step(sys);上述程序的运行结果为：0510-0.100.10.20.30.40.50.6ImpulseResponseTime(seconds)Amplitude0246800.20.40.60.811.21.4StepResponseTime(seconds)Amplitude[键入文字]②)()()(2)(txtytyty实验过程与实验结果：matlab程序如下：b=[100];a=[1sqrt(2)1];sys=tf(b,a);subplot(121);impulse(sys);subplot(122);step(sys);上述程序的运行结果为：③)()()()(txtytyty实验过程与实验结果：matlab程序如下：b=[10];a=[111];sys=tf(b,a);subplot(121);impulse(sys);subplot(122);step(sys);上述程序的运行结果为：02468-1.6-1.4-1.2-1-0.8-0.6-0.4-0.200.2ImpulseResponseTime(seconds)Amplitude0246-0.4-0.200.20.40.60.81StepResponseTime(seconds)Amplitude[键入文字]④)()()()()(txtxtytyty实验过程与实验结果：matlab程序如下：b=[110];a=[111];sys=tf(b,a);subplot(121);impulse(sys);subplot(122);step(sys);上述程序的运行结果为：（2）已知某系统可以由如下微分方程描述)()(6)()(txtytyty①请利用MATLAB绘出该系统单位冲激响应和单位阶跃响应的时域波形；实验过程与实验结果：matlab程序如下：b=[1];a=[116];0510-0.4-0.200.20.40.60.81ImpulseResponseTime(seconds)Amplitude0510-0.100.10.20.30.40.50.6StepResponseTime(seconds)Amplitude0510-0.6-0.5-0.4-0.3-0.2-0.100.1ImpulseResponseTime(seconds)Amplitude0510-0.200.20.40.60.811.2StepResponseTime(seconds)Amplitude[键入文字]sys=tf(b,a);subplot(121);impulse(sys);subplot(122);step(sys);上述程序的运行结果为：②根据冲激响应的时域波形分析系统的稳定性；由波形图可知，t趋于无穷时，y(t)均趋于0，所以该系统是稳定的。③如果系统的输入为)()(tuetxt，求系统的零状态响应。实验过程与实验结果：matlab程序如下：b=[1];a=[116];sys=tf(b,a);t=0:0.01:10;x=exp(-t);lsim(sys,x,t);上述程序的运行结果为：0510-0.2-0.100.10.20.30.40.5ImpulseResponseTime(seconds)Amplitude051000.050.10.150.20.250.30.35StepResponseTime(seconds)Amplitude012345678910-0.200.20.40.60.811.2LinearSimulationResultsTime(seconds)Amplitude[键入文字]（3）已知描述离散系统的微分方程如下，试采用MATLAB绘出各系统的单位抽样响应，并根据单位抽样响应的时域波形分析系统的稳定性。①)()2(2)1(3)(nxnynyny实验过程与实验结果：matlab程序如下：b=[1];a=[132];impz(b,a,0:15);上述程序的运行结果为：分析：由图像可知，n趋于无穷时，y(n)趋于无穷，所以该系统是不稳定的。②)1(3)()2(8.0)1(5.0)(nxnxnynyny实验过程与实验结果：matlab程序如下：b=[1-3];a=[1-0.50.8];impz(b,a,0:10);上述程序的运行结果为：051015-8-6-4-2024x104n(samples)AmplitudeImpulseResponse012345678910-2.5-2-1.5-1-0.500.511.522.5n(samples)AmplitudeImpulseResponse[键入文字]分析：由图像可知，n趋于无穷时，y(n)有界，且趋于0，所以该系统是稳定的。（4）已知系统可以由如下差分方程描述)()2(25.0)1()(nxnynyny试采用MATLAB绘出该系统的单位抽样响应波形和单位阶跃响应波形。实验过程与实验结果：matlab程序如下：clear;b=[1];a=[110.25];subplot(121);impz(b,a,0:10);subplot(122);stepz(b,a);上述程序的运行结果为：（5