MATLAB课程设计基于MATLAB离散系统Z域分析.docxVIP

课程设计任务书学生姓名：专业班级：电子科学与技术1001班指导教师：韩屏工作单位：信息工程学院题目:基于MATLAB的离散系统的Z域分析课题要求：利用MATLAB强大的图形处理功能，符号运算功能和数值计算功能，实现离散系统的Z域分析仿镇波形。课题内容：用MATLAB绘制离散系统极零图，根据极零图分布观察系统单位响应的时域特性并分析系统的稳定性。将极零图与h(k)对照起来画，看两者之间的关系。至少以六个例子说明。用MATLAB实现离散系统的频率特性分析以二个实例分别代表低通，高通滤波器，绘出极零图，幅频特性，相频特性。用MATLAB绘出梳状滤波器极零图与幅频特性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年 月 日1离散系统的Z域分析1.1 z变换z变换是离散信号与系统分析的重要方法和工具。z变换在离散信号与系统分析中的地位和作用，类似于连续信号与系统分析中的拉普拉斯变换，它将离散系统的数字模型——差分方程转化为简单的代数方程，使其求解过程得以简化。离散序列x(n)的z变换定义为：。在MATLAB中可以利用符号表达式计算一个因果序列的z变换。其命令格式为：syms n; 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为离散序列的符号表达式，输出参量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变换，即对于线性移不变系统，若n0时，h(n)=0，则系统为因果系统；若，则系统稳定。由于h(n)为因果序列，所以H(z)的收敛域为收敛圆外部区域，因此H(z)的收敛域为收敛圆外部区域时，系统为因果系统。因为，若z=1时H(z)收敛，即，则系统稳定，即H(z)的收敛域包括单位圆时，系统稳定。因此因果稳定系统应满足的条件为：，即系统函数H(z)的所有极点全部落在z平面的单位圆之内。2离散系统零极点图及零极点分析2.1离散系统零极点线性时不变离散系统可用线性常系数差分方程描述，即 （2-1）其中为系统的输出序列，为输入序列。将式（2-1）两边进行Z变换的 （2-2）将式（2-2）因式分解后有： （2-3）其中为常数，为的个零点，为的个极点。系统函数的零极点分布完全决定了系统的特性，若某系统函数的零极点已知，则系统函数便可确定下来。因此，系统函数的零极点分布对离散系统特性的分析具有非常重要意义。通过对系统函数零极点的分析，可以分析离散系统以下几个方面的特性：系统单位样值响应的时域特性；离散系统的稳定性；离散系统的频率特性；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);titl