matlab用lism求零输入响应,实验2 lti系统的时域分析.docx

实验2 lti系统的时域分析.docx

键入文字实验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表示输入信号数值向量及其时间向量。YLSIMSYS,X,T这种调用格式不绘出波形，而是返回响应的数值向量。3)连续时间系统的冲激响应与阶跃响应MATLAB提供了IMPULSE函数来求指定时间范围内由模型SYS描述的连续时间系统的单位冲激响应。IMPULSE函数的基本调用格式如下IMPULSESYS在默认的时间范围内绘出系统冲激响应的时域波形；IMPULSESYS,T绘出系统在0～T范围内冲激响应的时域波形；IMPULSESYS,TSTPTE绘出系统在TS～TE范围内，以TP为时间间隔取样的冲激响应波形。Y,TIMPULSE这种调用格式不绘出冲激响应波形，而是返回冲激响应的数值向量及其对应的时间向量。函数STEP用于求解单位阶跃响应，函数STEP同样也有如下几种调用格式STEPSYS；键入文字STEPSYS,T；STEPSYS,TSTPTE；Y,TSTEP。各种调用格式参数所代表的意思可参考上述IMPULSE函数。2离散时间系统时域分析的MATLAB实现1)离散时间系统的MATLAB表示LTI离散系统通常可以由系统差分方程描述，设描述系统的差分方程为𝑎0𝑦𝑛𝑎1𝑦𝑛1𝑎𝑁𝑦𝑛𝑁𝑏0𝑥𝑚𝑏1𝑥𝑚1𝑏𝑀𝑥𝑚𝑀则在MATLAB里，我们可以用如下两个向量来表示这个系统；𝑏𝑏0,𝑏1,,𝑏𝑀；𝑎𝑎0,𝑎1,,𝑎𝑁2)离散时间系统对任意输入的响应MATLAB提供了求LTI离散系统响应的专用函数FILTER，该函数用于求取由差分方程所描述的离散时间系统在执行时间范围内对输入序列所产生的响应，该函数基本调用格式为YFILTERB,A,X其中，X为输入序列，Y为输出序列，输出序列Y对应的时间区间与X对应的时间区间相同。3)离散时间系统的单位抽样响应MATLAB提供了函数IMPZ来求指定时间范围内，由向量B和A描述的离散时间系统的单位抽样响应，具体调用格式如下IMPZB,A在默认的时间范围内绘出系统单位抽样响应的时域波形；IMPZB,A,N绘出系统在0～N时间范围内单位抽样响应的时域波形；IMPZB,A,NSNE绘出系统在NS～NE范围内单位抽样响应的时域波形；Y,TIMPZ这种调用格式不绘出单位抽样响应波形，而是返回单位抽样响应的数值向量及其对应的时间向量。3．卷积和与卷积积分1)离散时间序列的卷积和卷积和是离散系统时域分析的基本方法之一，离散时间序列的卷积和𝑥1𝑛和𝑥2𝑛定义如下𝑥𝑛𝑥𝑛𝑥1𝑛𝑥2𝑛∞∑𝑘‒∞𝑥1𝑘𝑥2𝑛𝑘对于离散LTI系统，设其输入信号为，单位抽样响应为，则其零状态响应为𝑥𝑛H𝑛𝑦𝑛𝑦𝑛𝑥𝑛H𝑛即离散LTI系统的零状态响应可以表示出输入信号与单位抽样响应的卷积。𝑥𝑛H𝑛因此，离散时间序列的卷积和计算对于我们进行离散系统时域分析具有重要的意义。MATLAB的CONV函数可以用来求两个离散序列的卷积和，调用格式为XCONVX1,X2。2)连续时间信号的卷积积分卷积积分是连续系统时域分析的有效方法和工具之一，连续时间信号的卷𝑥1𝑡和𝑥2𝑡积积分定义如下𝑥𝑡𝑥𝑡𝑥1𝑡𝑥2𝑡∫∞∞𝑥1𝜏𝑥2𝑡𝜏𝑑𝜏对于连续LTI系统，设其输入信号为，单位抽样响应为，其零状态响应为，𝑥𝑡H𝑡𝑦𝑡则有键入文字𝑦𝑡𝑥𝑡H𝑡即连续LTI系统的零状态响应可以表示为输入信号与单位冲激响应的卷积。因此，连续时间信号卷积积分对连续系统的时域分析具有非常重要的意义。利用MATLAB可以采用数值计算的方法近似计算卷积积分。卷积积分可用求和运算来实现𝑥𝑡𝑥1𝑡𝑥2𝑡∫∞∞𝑥1𝜏𝑥2𝑡𝜏𝑑𝜏LIM∆→0∞∑𝑘‒∞𝑥1𝑘∆𝑥2𝑡𝑘∆∙∆现在考虑只求时，则有𝑡𝑛∆𝑥𝑡的值𝑥𝑛∆𝑥𝑛∆∆∙∞∑𝑘‒∞𝑥1𝑘∆𝑥2𝑛𝑘∆当足够小，的数值近似。我们可以利用计算离散卷积序列卷积和的∆𝑥𝑛∆就是𝑥𝑡CONV函数来计算卷积积分，具体步骤如下①将连续时间信号和以时间间隔Δ进行取样，得到离散序列和𝑥1𝑡𝑥2𝑡𝑥1𝑛∆；𝑥2𝑛∆②构造离散序列和对应的时间向量和；𝑥1𝑡𝑥2𝑡𝑡1𝑡2③调用函数CONV计算卷积积分在时的近似采样值；𝑡𝑛∆𝑥𝑛∆④构造离散序列对应的时间向量。𝑥𝑛∆𝑛三、实验内容(1)已知描述模拟低通、高通、带通和带阻滤波器的微分方程如下，试采用MATLAB绘出各系统的单位冲激响应和单位阶跃响应。①2TXYTY实验过程与实验结果MATLAB程序如下B1A1SQRT21SYSTFB,ASUBPLOT121IMPULSESYSSUBPLOT122STEPSYS上述程序的运行结果为键入文字0510010010203040506IMPULSERESPONSETIMESECONDSAMPLITUDE0246800204060811214STEPRESPONSETIMESECONDSAMPLITUDE②2TXYTY实验过程与实验结果MATLAB程序如下B100A1SQRT21SYSTFB,ASUBPLOT121IMPULSESYSSUBPLOT122STEPSYS上述程序的运行结果为02468161412108060402002IMPULSERESPONSETIMESECONDSAMPLITUDE024604020020406081STEPRESPONSETIMESECONDSAMPLITUDE③TXYTY实验过程与实验结果MATLAB程序如下B10A111SYSTFB,A键入文字SUBPLOT121IMPULSESYSSUBPLOT122STEPSYS上述程序的运行结果为051004020020406081IMPULSERESPONSETIMESECONDSAMPLITUDE0510010010203040506STEPRESPONS