资源描述:
《matlab课程设计---基于matlab的离散系统的z域分析》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、Matlab应用实践课程设计课程设计任务书学生姓名:专业班级:电子科学与技术1001班指导教师:韩屏工作单位:信息工程学院题目:基于MATLAB的离散系统的Z域分析课题要求:利用MATLAB强大的图形处理功能,符号运算功能和数值计算功能,实现离散系统的Z域分析仿镇波形。课题内容:一.用MATLAB绘制离散系统极零图,根据极零图分布观察系统单位响应的时域特性并分析系统的稳定性。将极零图与h(k)对照起来画,看两者之间的关系。至少以六个例子说明。二.用MATLAB实现离散系统的频率特性分析1.以二个实例分别代表低通,高通滤波器,绘出极零图,幅频特性,相频特性。2.用MATLAB绘出梳状滤
2、波器极零图与幅频特性FIR型IIR型设N=8,a=0.8,0.9,0.98三.用MATLAB实现巴特沃兹滤波器分析1.用MATLAB绘制巴特沃兹滤波器频率特性曲线(w,n作为参数变化)2.用MATLAB绘制巴特沃兹滤波器的极零点分布图(w,n作为参数变化)将两种图对照起来看极点分布与频率特性之间的关系。时间安排:学习MATLAB语言的概况第1天学习MATLAB语言的基本知识第2、3天学习MATLAB语言的应用环境,调试命令,绘图能力第4、5天课程设计第6-9天答辩第10天指导教师签名:2013年月日系主任(或责任教师)签名:2013年月20Matlab应用实践课程设计日目录1离散系统
3、的Z域分析31.1z变换31.2利用MATLAB的符号运算实现z变换31.3离散系统的系统函数及因果稳定的系统应满足的条件32离散系统零极点图及零极点分析42.1离散系统零极点42.2零极点的绘制53MATLAB实现离散系统的频率特性分析113.1低通滤波器113.2高通滤波器123.3梳状滤波器的特性分析134MATLAB实现巴特沃兹滤波器分析175总结体会196参考文献1920Matlab应用实践课程设计1离散系统的Z域分析1.1z变换z变换是离散信号与系统分析的重要方法和工具。z变换在离散信号与系统分析中的地位和作用,类似于连续信号与系统分析中的拉普拉斯变换,它将离散系统的数字
4、模型——差分方程转化为简单的代数方程,使其求解过程得以简化。离散序列x(n)的z变换定义为:。在MATLAB中可以利用符号表达式计算一个因果序列的z变换。其命令格式为:symsn;f=(1/2)^n+(1/3)^n;ztrans(f)1.2利用MATLAB的符号运算实现z变换如果离散序列()xn可以用符号表达式,则可以直接用MATLAB的ztrans函数来求离散序列的单边变换。调用ztrans函数的命令格式为①Z=ztrans(X)②Z=ztrans(X,w)格式①中输如参量X为离散序列的符号表达式,输出参量Z为返回默认符号自变量为n的关于X的变换的符号表达式。格式②中输如参量X为离
5、散序列的符号表达式,输出参量Z为返回符号自变量为w的关于X的变换的符号表达式。1.3离散系统的系统函数及因果稳定的系统应满足的条件一个线性移不变离散系统可以用它的单位抽样响应h(n)来表示其输入与输出关系,即y(n)=x(n)*h(n)对该式两边取z变换,得:Y(z)=X(z)·H(z)则:将H(z)定义为系统函数,它是单位抽样响应h(n)的z变换,即对于线性移不变系统,若n<0时,h(n)=0,则系统为因果系统;若,则系统稳定。由于h(n)为因果序列,所以H(z)的收敛域为收敛圆外部区域,因20Matlab应用实践课程设计此H(z)的收敛域为收敛圆外部区域时,系统为因果系统。因为,
6、若z=1时H(z)收敛,即,则系统稳定,即H(z)的收敛域包括单位圆时,系统稳定。因此因果稳定系统应满足的条件为:,即系统函数H(z)的所有极点全部落在z平面的单位圆之内。2离散系统零极点图及零极点分析2.1离散系统零极点线性时不变离散系统可用线性常系数差分方程描述,即(2-1)其中为系统的输出序列,为输入序列。将式(2-1)两边进行Z变换的(2-2)将式(2-2)因式分解后有:(2-3)其中为常数,为的个零点,为的个极点。系统函数的零极点分布完全决定了系统的特性,若某系统函数的零极点已知,则系统函数便可确定下来。因此,系统函数的零极点分布对离散系统特性的分析具有非常重要意义。通过对
7、系统函数零极点的分析,可以分析离散系统以下几个方面的特性:l系统单位样值响应的时域特性;l离散系统的稳定性;l离散系统的频率特性;20Matlab应用实践课程设计2.2零极点的绘制(1)MATLAB程序实现:b=[0,1,-1];a=[1,-2.5,1];rz=roots(b)rp=roots(a)subplot(1,2,1),zplane(b,a);title('系统的零极点分布图');subplot(1,2,2),impz(b,a,20);title(