离散时间系统分析与MATLAB实现.docVIP

离散时间系统分析及MATLAB实现 摘自：张登奇,周婷,梁莺. HYPERLINK /kns55/detail.aspx?dbName=CJFD2009FileName=YYSF200903012DBCode=CJFQ \t _blank 离散时间系统分析及MATLAB实现[J]. HYPERLINK /kns55/oldNavi/Bridge.aspx?LinkType=BaseLinkDBCode=cjfqTableName=cjfdbaseinfoField=BaseIDValue=YYSF \t _blank 湖南理工学院学报(自然科学版), HYPERLINK /kns55/oldNavi/Bridge.aspx?LinkType=IssueLinkDBCode=cjfqTableName=cjfdyearinfoShowField=cnameField=BaseID*year*issueValue=YYSF*2009*03 \t _blank 2009,(03) 摘 要：线性时不变离散时间系统是最基本的数字系统，差分方程和系统函数是描述系统的常用数学模型，单位脉冲响应和频率响应是描述系统特性的主要特征参数，零状态响应和因果稳定性是系统分析的重要内容。文章从系统的分析流程、系统模型的创建、时域分析、频域分析和因果稳定性分析等方面，介绍了线性时不变离散时间系统的基本分析方法，并以实例形式列举了MATLAB实现程序。 关键词：MATLAB；离散时间系统；系统分析；传输函数 离散时间系统是将一个序列变换成另一序列的系统，它有多种类型，其中线性时不变离散时间系统是最基本、最重要的系统。差分方程反映了系统输入与输出的运动状态，是在时域描述系统的通用数学模型；系统函数是零状态下系统输出与输入的Z变换之比，在时域与频域之间起桥梁作用。分析系统就是在已知系统结构或系统模型条件下，从时域和频域两方面分析系统输入与输出的关系，前者重点研究系统的时间特性，后者主要研究系统的频率特性。下面从系统分析流程、系统模型创建、系统时域分析、系统频域分析和因果稳定性分析等方面，介绍线性时不变离散时间系统的基本分析方法，并以实例形式列举MATLAB在系统分析过程中的具体应用。 3 系统的时域分析 系统的时域分析主要包括求表征系统时间特性的单位脉冲响应特征量和在时域求解系统输出两方面内容。 3.1 单位脉冲响应的计算 根据差分方程求解单位脉冲激励下系统的零状态响应，或将系统函数进行Z反变换都可算出系统的单位脉冲响应，具体算法可参见参考文献[3]。在MATLAB中描述系统的差分方程或系统函数都是用系数向量表示，调用impz函数就可直接算出系统的单位脉冲响应。如实例1描述的系统，其单位脉冲响应的计算及显示程序如下： b=[0.3,0.06,0,0]; %系数向量不齐后面补0 a=[1,-1.1,0.55,-0.125]; %系数向量不齐后面补0 [hn,n]=impz(b,a,16), %列向求出16点单位脉冲响应 stem(n,hn,.); grid; %绘制点状图并加网格 xlabel(n);ylabel(hn);title(单位脉冲响应); 若要写出闭环形式，可调用residuez函数将系统函数展开成部分分式形式，再通过查表求Z反变换即可。 3.2 系统输出的时域计算 在时域上计算离散时间系统的输出，实际上就是直接求解差分方程或作卷积运算。参考文献[3]列举了迭代法、时域经典法、卷积法等常用方法及应用实例。考虑到分析系统的目的在于综合，系统设计时不存在初始问题，因此，分析系统响应重点分析零状态响应。只要掌握了分析系统的概念、原理和方法，繁杂的计算可由MATLAB完成。 实例2：试计算实例1中，当输入序列分别为单位脉冲、单位阶跃和一般序列时，系统的输出响应。 方法1：调用filter函数实现 b=[0.3,0.06,0,0]; a=[1,-1.1,0.55,-0.125]; x1=[1,zeros(1,15)]; %产生16点单位脉冲序列 x2=ones(1,16); %产生16点单位阶跃序列 x3=exp(-0.2*[0:15]); %用指数序列代表一般序列 y1=filter(b,a,x1), %计算单位脉冲响应 y2=filter(b,a,x2), %计算单位阶跃响应 y3=filter(b,a,x3), %计算一般序列响应 方法2：调用conv函数实现 b=[0.3,0.06,0,0]; a=[1,-1.1,0.55,-0.125]; [hn,n]=impz(b,a,16); %求出16点单位脉冲响应 x=exp(-0.2*(0:15)); %输入或产生一般序列 y1=conv(hn,x), %用线性卷积求系统响应 y2=filter(