基于Z变换的离散系统分析.doc

皖西学院数字信号处理课程设计报告数字信号处理课程设计题目： 基于Z变换的离散系统分析 学院： 信息工程学院 专业： 通信工程 班级： 学号： 姓名： 指导教师： 目 录摘要1第一章 背景31.1背景知识31.2软件介绍31.3MATLAB软件功能简介4第二章 设计目的及要求62.1设计目的62.2设计要求6第三章 设计原理73.1Z变换的定义73.2Z变换的定义83.3几种不同特性的序列的介绍93.4频率响应函数与系统函数的定义113.5系统函数的零极点分析123.6离散时间LTI系统的频率特性分析133.7系统函数的极点分布与系统因果性、稳定性的关系13第四章 设计过程与结果分析15收获与体会21参考文献22皖西学院数字信号处理课程设计报告摘要离散时间系统是将一个序列变换成另一序列的系统,它有多种类型,其中线性时不变离散时间系统是最基本、最重要的系统. Z变换是离散信号与系统分析的重要方法和工具，是对离散序列进行的一种数学变换。它在离散时间系统中的地位，如同拉普拉斯变换在连续时间系统中的地位。离散时间信号的变换已成为分析线性时不变离散时间系统问题的重要工具。在数字信号处理、计算机控制系统等领域有广泛的应用。Matlab语言是一种广泛应用于工程计算及数值分析领域的新型高级语言，Matlab功能强大、简单易学、编程效率高。将Matlab引入到线性时不变离散时间系统的Z域分析，通过运用Matlab的LST对象和符号运算功能分析离散系统的有关问题，给出相应的计算程序和运算结果，在运用图形分析功能，将结果用图形表示出来，从而是离散系统的分析得到直观可视化效果。关键词:离散时间系统；Z变换；MatlabAbstractDiscrete time system is a sequence of transform into another sequence of system, it has a variety of types, in which the linear time-invariant discrete-time system is the most basic, most important system. Z transform is a discrete signal and system analysis is an important method and tool, is a discrete sequence of a mathematical transformation. It is in discrete time system in the position, as Laplasse transform in continuous time system status. Discrete time signal transformation has become the analysis of linear time invariant discrete time systems are important tools. In digital signal processing, computer control systems are widely used in the field of.The Matlab language is widely used in engineering calculation and the numerical analysis in the field of advanced language, Matlab powerful, easy programming, high efficiency. The Matlab is introduced to linear time-invariant discrete time system Z domain analysis, through the use of Matlab LST objects and symbolic computation of discrete system function analysis of related issues, the corresponding calculation program and results of operations, in the use of graphical analysis function, the results using the graph expresses, from but discrete system analysis is intuitive visual effects.Keywords: Discrete time system; Z transform; Matlab21第一章 背景1.1 背景知识数字信号处理（Digital Signal Processing,简称DSP）是一门设计许多学科而又广泛应用于许多领域的新兴学科。DSP有两种含义：Digital Signal Processing（数字信号处理）、Digital Signal Processor（数字信号处理器）。我们常说的DSP指的是数字信号处理器。数字信号处理器是一种适合完成数字信号处理运算的处理器。20世纪60年代以来，随着计算机和信息技术的飞速发展，数字信号处理技术应运而生并得到迅速的发展。在过去的二十多年的时间里，数字信号处理已经在通信等领域得到极为广泛的应用。数字信号处理是利用计算机专用处理设备，以数字形式对信号进行采集、变换、滤波、估值增强、压缩、识别等处理，以得到符合人们需要的信号形式。它是以众多学科为理论基础的，它所涉及的范围极其广泛。例如，在数学领域，微积分、概率统计、随即过程、数值分析等都是数字信号处理的基本工具，与网络理论、信号与系统、控制论、通信理论、故障诊断等也密切相关，近来新兴的一些学科，如人工智能、模式识别、神经网络等，都与数字信号处理密不可分。可以说，数字信号处理是把许多经典的理论体系作为自己的理论基础，同时又使自己成为一系列新兴学科的理论基础。1.2 软件介绍 MATLAB是矩阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品的商业数字数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、经融建模设计与分析等领域。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解答问题要比用C,FORTRAN等语言完成相同的事情简捷的多，并且mathworks也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C,FORTRAN，C+，JAVA的支持。可以直接调用，用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。1.3 MATLAB软件功能简介MATLAB的名称源自Matrix Laboratory,1984年由美国Mathworks公司推向市场。它是一种科学计算软件，专门以矩阵的形式处理数据。MATLAB将高性能的数值计算和可视化集成在一起，并提供了大量的内置系数，从而被广泛的应用于科学计算、控制系统和信息处理等领域的分析、仿真和设计过程。MATLAB软件包括五大通用功能（Nemeric）、符号运算功能（Symbolic）、数据可视化（Graphic）、数字图形文字统一处理功能（Notebook）和建模仿真可视化功能（Simulink）。其中，符号运算功能的实现是通过请求MAPLE内核计算并将结果返回到MATLAB命令窗口。该软件有三大特点，一是功能强大；二是界面友善、语言自然；三是开放性强。目前，Mathworks公司已推出30多个应用工具箱。MATLAB在线性代数、矩阵分析、数值及优化、数值统计和随机信号分析、电路与系统、系统动力学、图像处理、控制理论分析和系统设计、过程控制、建模和仿真、通信系统以及财政金融等众多领域的理论研究和工程设计中得到了广泛应用。MATLAB在信号与系统中的应用主要包括符号运算和数值计算仿真分析。由于信号与系统课程的许多内容都是基于公式演算，而MATLAB借助符号数学工具箱提供的符号运算功能，能基本满足信号与系统课程的需求。例如解微分方程、傅立叶正反变换、拉普拉斯正反变换和z正反变换等。MATLAB在信号与系统中的另一主要应用是数值计算与仿真分析，主要包括函数波形绘制、函数运算、冲级响应与阶跃响应仿真分析、信号的时域分析、信号的频谱分析、系统的S域分析和零极点图绘制等内容。数值计算仿真分析可以帮助学生更深入地理解理论知识，并为将来使用MATLAB进行信号处理领域的各种分析和实际应用打下基础。第二章 设计目的及要求2.1 设计目的1. 理解并掌握序列的Z变换以及序列特性与收敛域之间的关系，2. 学会利用Z变换求解离散时间系统的系统函数的零、极点，3. 学会分析离散时间系统的系统函数的零、极点分布，4. 学会利用系统的零极点分布分析系统的频率响应特性，5. 学会运用MATLAB分析离散时间系统的系统函数的零、极点分布，6. 学会运用MATLAB进行离散时间系统的频率特性分析，7. 学会运用MATLAB分析离散系统的系统函数零极点分布与系统稳定性，因果性的关系。2.2 设计要求1. 巩固所学的专业技术知识；2. 提高综合运用所学理论知识独立分析和解决问题的能力；3. 进一步提高程序设计及调试能力；4. 更好地将理论与实践相结合；5. 学习和掌握科学研究资料检索的方法，学习对已有资料进行消化总结的方法6. 学习撰写科学报告的基本方法；7. 本设计要求分组合作完成；8. 熟悉使用课程设计实验平台，掌握其基本的操作方法；9. 了解课程设计实验平台的源代码，掌握其程序结构及在此平台上添加处理程序的方法；10. 设计过程中详细记录产生的图形、参数、数据等，用于编写课程设计报告。第三章 设计原理3.1 Z变换的定义 一个离散序列x(n)的Z变换定义为 （3.1-1）式子中Z是一个复变量，它所在的平面称为Z平面。注意在定义中，对n求和是在负无穷，正无穷之间求和。可以称为双边Z变换。还有一种称为单边Z变换的定义，如下式：（3.1-2）这种单边Z变换的求和限是从零到无限大，因此对于因果序列，用两种Z变换定义计算的结果是一样的。注意：Z变换存在的条件是等号右边级数收敛，要求级数绝对可和，即（3.1-3）Z的变量取值的域称为收敛域，收敛域一般取环状域，即常用的Z变换是一个有理函数：X(z)的零点：P(z)的根， X(z)的极点：Q(z)的根。收敛域中没有极点， 收敛域总是用极点限定其边界。3.2 几种不同特性的序列的介绍Z平面上收敛域的位置和序列有着密切的关系：(1) 有限长序列 即序列从到的序列值不全为零，此范围之外序列值为零，这样的序列称为有限长序列。其Z变换为其收敛情况：(2) 右边序列右边序列是指只在时有值。右序列的Z变换表示为右边序列Z变换的收敛域为u 因果序列其Z变换收敛域包括是因果序列的特征。(3) 左边序列 左边序列是指在时有值左边序列Z变换的收敛域为 如果，收敛域应包括，即(4) 双边序列一个双边序列可看作一个右边序列和一个左边序列之和，其Z变换表示为第一部分收敛域为；第二部分收敛域为收敛域为如果，则存在公共收敛区域：；如果，两个收敛域没有交集，则X(z)没有收敛域，那么X(z)不存在。3.3 频率响应函数与系统函数的定义在时域中，一个线性时不变系统完全由它的单位脉冲响应h(n)来表示。 其Z变换定义为线性时不变系统的系统函数。系统函数：频响函数：在单位圆上()的系统函数就是系统的频率响应。3.4 系统函数的零极点分析离散时间系统的系统函数定义为系统零状态响应的Z变换与激励的Z变换之比，即如果系统函数的有理函数表示式为那么，在MATLAB中系统函数的零极点就可通过函数roots得到，也可借助函数tf2zp得到，tf2zp的语句格式为Z,P,K=tf2zp(B,A)其中，B与A分别表示的分子与分母多项式的系数向量。它的作用是将的有理分式表示式转换为零极点增益形式，即若要获得系统函数的零极点分布图，可直接应用zplane函数，其语句格式为zplane(B,A)其中，B与A分别表示的分子和分母多项式的系数向量。它的作用是在Z平面上画出单位圆、零点与极点。3.5 离散时间LTI系统的频率特性分析对于因果稳定的离散时间系统，如果激励序列为正弦序列，则系统的稳态响应为。其中，通常是复数。离散时间系统的频率响应定义为其中，称为离散时间系统的幅频特性；称为离散时间系统的相频特性；是以（，若零，）为周期的周期函数。因此，只要分析在范围内的情况，便可分析出系统的整个频率特性。MATLAB提供了求离散时间系统频响特性的函数freqz，调用freqz的格式主要有两种。一种形式为H,w=freqz(B,A,N)其中，B与A分别表示的分子和分母多项式的系数向量；N为正整数，默认值为512；返回值w包含范围内的N个频率等分点；返回值H则是离散时间系统频率响应在范围内N个频率处的值。另一种形式为H,w=freqz(B,A,N,whole)与第一种方式不同之处在于角频率的范围由扩展到。3.6 系统函数的极点分布与系统因果性、稳定性的关系单位脉冲响应h(n)为因果序列的系统称为因果系统， 因果系统的系统函数H(z)具有包括z=点的收敛域，即一个线性时不变系统稳定的充分必要条件而Z变换的收敛域由满足稳定系统的系统函数H(z)必须在单位圆上收敛，即收敛域包括单位圆|z|=1，H(ej)存在。u 因果稳定系统它的系统函数H(z)必须在从单位圆到的整个Z域内收敛，也就是说，系统函数的全部极点必须在单位圆内。具体系统的因果性和稳定性可由系统函数H(z)的极点分布和收敛域确定。第四章 设计过程与结果分析1. 自行设计以下几种序列：有限长序列、右序列、左序列和双边序列，分别求其Z变换和收敛域，分析Z变换收敛域的规律。(1) 对序列x(n)=RN(n)的Z变换及其收敛域:解 因为所以这是一个因果的有限长序列，因此收敛域为(2) 对序列, 求其Z变换及收敛域。解 这是一个因果序列，其Z变换为 收敛域： 由于X(z)只在z=a处有一极点，整个收敛域应该在极点所在的圆外。 结果分析：右边序列的Z变换如果有N个有限极点z1,z2,zN， 那么收敛域一定在模最大的极点所在的圆外(3) 对序列, 求其Z变换及收敛域。解 这是一个左边序列。其Z变换为此等比级数在，即收敛。由于函数只在z=a处有一极点，整个收敛域应该在极点所在的圆内。 结果分析：对于左边序列，如果序列Z变换有N个有限极点z1, z2, , zN，那么收敛域一定在模最小的极点所在的圆内结论：一个左边序列与一个右边序列的变换表达式是完全一样的。所以，只给出Z变换的闭合表达式不能正确得到原序列，需要已知收敛域。(4) 对序列x(n)=a|n|, a为实数，求其Z变换及收敛域。解 这是一个双边序列，其Z变换为 若，则存在公共收敛域 若，则无公共收敛域,序列两端都发散2. 已知某离散系统的系统函数，要求：（1）参数、由运行时输入；（2）画出系统的零、极点分布图；（3）判断系统的稳定性；（4）判断系统的因果性；（5）绘出系统的幅频响应曲线；（6）分析零、极点分布对系统频率响应特性的影响。解 将系统函数H(z)变换得，运用MATLAB计算零极点，绘制频率响应曲线以及判断系统稳定性，程序如下：c1=1;c2=3;d1=4;d2=2;%定义系统初始变量x1=c1+c2;x2=d1+d2;y1=c1*c2;y2=d1*d2; %定义系统初始变量B=1 x1 y1;A=1 x2 y2;%设置系统函数系数向量B和Asubplot(2,2,1);zplane(B,A);%绘制零极点图title(零极点分布图)H,w=freqz(B,A);%计算频率响应subplot(2,2,2);plot(w/pi,abs(H);%绘制幅频响应曲线xlabel(omega/pi);ylabel(|H(ejomega)|);title(系统函数幅频特性)subplot(2,2,3);plot(w/pi,angle(H); %绘制相频响应曲线xlabel(omega/pi);ylabel(phi(omega);title(系统函数相频特性)disp(系统零点为:)Z=roots(B) %求H(z)的零点，并显示disp(系统极点为:)P=roots(A) %求H(z)的极点，并显示disp(系统极点模的最大值为:)M=max(abs(P) %求所有极点模的最大值，并显示if M1 disp(系统稳定),else disp(系统不稳定),end %通过判断极点模的最大值M,如果M1,则显示“系统稳定”，否则显示“系统不稳定”程序运行结果如下：图1 系统函数的零极点分布及幅频、相频特性系统零点为:Z = -3 -1系统极点为:P = -4 -2系统极点模的最大值为:M = 4系统不稳定由于H(z)