(实验四Z变换.docVIP

实验四 Z变换 【实验目的】 通过MATLAB仿真离散时间系统，研究其时频域特性，加深对离散系统的冲激响应，频率响应分析和零极点分布概念和理解。 【实验原理】 1.Z变换原理 （1）.Z变换 在数字信号处理的分析方法中，除时域分析方法外，还有变换域分析方法。后者通常指Z变换和傅里叶变换法。变换域分析的最大优点是将离散系统的差分方程转化为简单的代数方程，使其求解大大简化，也使得对系统的特性分析更为方便。 对于离散时间信号，设序列为x(n)，则其Z变换定义为： ，其中z为复变量，是一个以时部为横坐标，虚部为纵坐标构成的平面上的变量。Z变换记作，X(z)存在的z的集合称为收敛域（ROC），一般为 由于ROC是由定义的，因此一般为环形区域。根据ROC的特点，可以判定序列是右边序列、左边序列、双边序列等。 Z变换具有一些重要的特性，是傅里叶变换的推广，包括线性、时移特性、频移特性、尺度变换、共轭、翻褶、Z域微分、序列相乘、序列卷积等一系列性质。 （2）.系统函数 离散线性时不变（LTI）系统的系统函数H(z)定义为： H(z) = Z[h(n)] = (4.4) 若用差分方程表示系统，则有 如果系统起始状态为零，直接对上式的两边Z变换，并利用移位特性，有 因此，系统函数H(z)的分子和分母的系数正好等于差分方程的系数。归一化，即使得y(n)前的参数为1，此时可以对上式的分子、分母进行因式分解，可得 得到系统的增益函数K、零点、极点。它们分别由差分方程的系数决定，确定了系统的因果性和稳定性。 （3）.系统的频率响应 系统的频率响应定义在单位圆上() 上的系统函数，即 (4.8) 其中称为系统的幅度响应，称为系统的相位响应。 由于具有周期性，且周期为，通常只研究区间，相对于这一区间，“低频”就是靠近0处的频率，而“高频”就是靠近的区域。 （4）.结构 从频率角度来看，离散LTI系统对信号的响应过程实质上是对信号滤波的过程，这时离散LTI系统被称为数字滤波器。如果系统有极点，则为无限长单位冲激响应（IIR）数字滤波器，否则为有限长单位冲激响应（FIR）数字滤波器。前者分为直接I型，直接II型、级联型、并联型等。后者分为横截型、级联型等。 2.matlab实现 （1）. impz（）函数 如果把系数{},{}整合，构成矢量a、b，则这两个矢量构成了数字滤波器的所有参数，可以表征这个数字滤波器。因此，根据这两个变量，可以计算相应数字滤波器的冲激响应、频率响应等一系列内容。 在MATLAB中，有专门计算数字滤波器冲激响应的函数impz，并且能直接在当前图形窗口中画出离散数字系统的冲激响应图形,这在上次实验中也做过。 （2）f reqz（）函数 MATLAB中有专门计算数字滤波器频率响应的函数f reqz，能直接求出幅频特性和相频特性，并且直接在当前图形窗口中画出数字系统的幅频、相频特性曲线。 例4-2 对一离散LTI系统，有系统函数 求其幅频、相频特性曲线 程序如下： b = [0.2,0.3,1]; a = [1,0.4,1]; f reqz(b, a,128) 运行该程序后，直接得到下图，即该系统的幅频、相频特性曲线。 如果将程序最后一句语句改为[h,w] = f reqz(b,a, 128)，则不会产生幅频、相频特性曲线，而是向量h得到该数字滤波器的128点的幅频响应，这些点均匀地分布在上半单位圆(即)，而这些频率点将记录在向量w中。当然，点数也可以不取128，而选择其他任意值，只需要是大于0的整数，但最好能选2的幂次,方便FFT进行快速算法。如果缺省，默认取512个点。 （3）系统的零极点图和zplane（）函数 根据公式 ，可以得到零点{}和极点{}，零极点通常用二维坐标图来表示，其中极点用×表示，零点用○表示。同时还会绘出单位圆，以示参考。 在MATLAB中，有专门计算离散系统零极点的函数zplane，并且能直接在当前图形窗口中画出离散数字系统的零极点图形。 例4-3 同例4-1中的离散LTI系统，差分方程为 y(n)-1.1y(n - 1) + 1.5y(n - 2) – 0.7y(n - 3) + 0.3y(n - 4) = 0.2x(n) + 0.1x(n - 1) + 0.3x(n- 2) + 0.1x(n -3) + 0.2x(n -4) 求其零极点，并画出相应图形。 程序如下： b = [0.2, 0.1,0.3,0.1, 0.2]; a = [1,-1.1, 1.5,-0.7,0.3]; z