. . 实验六 离散信号与系统的Z变换分析 学院 班级 学号一、实验目的1.熟悉离散信号Z变换的原理及性质2.熟悉常见信号的Z变换3.了解正/反Z变换的MATLAB实现方法4.了解离散信号的Z变换与其对应的理想抽样信号的傅氏变换和拉氏变换之间的关系5.了解利用MATLAB实现离散系统的频率特性分析的方法二、实验原理1. 正/反Z变换Z变换分析法是分析离散时间信号与系统的重要手段。如果以时间间隔对连续时间信号f(t)进行理想抽样，那么，所得的理想抽样信号为：理想抽样信号的双边拉普拉斯变换Fd (s)为：若令 ，那么的双边拉普拉斯变换F (s)为：则离散信号f(k)的Z变换定义为： 从上面关于Z变换的推导过程中可知，离散信号f(k)的Z变换F(z)与其对应的理想抽样信号的拉氏变换Fd (s)之间存在以下关系： 同理，可以推出离散信号f(k)的Z变换F(z)和它对应的理想抽样信号的傅里叶变换之间的关系为 如果已知信号的Z变换F(z)，要求出所对应的原离散序列f(k)，就需要进行反Z变换，反Z变换的定义为： 其中，C为包围的所有极点的闭合积分路线。在MATLAB语言中有专门对信号进行正反Z变换的函数ztrans( ) 和itrans( )。其调用格式分别如下：l F=ztrans( f ) 对f(n)进行Z变换，其结果为F(z)l F=ztrans(f,v) 对f(n)进行Z变换，其结果为F(v)l F=ztrans(f,u,v) 对f(u)进行Z变换，其结果为F(v)l f=itrans ( F ) 对F(z)进行Z反变换，其结果为f(n)l f=itrans(F,u) 对F(z)进行Z反变换，其结果为f(u)l f=itrans(F,v,u ) 对F(v)进行Z反变换，其结果为f(u)注意： 在调用函数ztran( )及iztran( )之前，要用syms命令对所有需要用到的变量（如t,u,v,w）等进行说明，即要将这些变量说明成符号变量。例 用MATLAB求出离散序列 的Z变换。MATLAB程序如下：syms k zf=0.5k; %定义离散信号Fz=ztrans(f) %对离散信号进行Z变换运行结果如下：Fz = 2*z/(2*z-1)例 已知一离散信号的Z变换式为 ，求出它所对应的离散信号f(k)。MATLAB程序如下：syms k zFz=2* z/(2*z-1); %定义Z变换表达式fk=iztrans(Fz,k) %求反Z变换运行结果如下：fk = (1/2)k2.离散系统的频率特性同连续系统的系统函数H(s)类似，离散系统的系统函数H(z)也反映了系统本身固有的特性。对于离散系统来说，如果把其系统函数H(z)中的复变量z换成，那么所得的函数就是此离散系统的频率响应特性，即离散时间系统的频率响应为：其中， 称为离散系统的幅频特性，称为系统的相频特性。同连续系统一样，离散时间系统的幅频特性也是频率的偶函数，相频特性也是频率的齐函数。由于是频率的周期函数，所以离散系统的频率响应特性也是频率的周期函数，其周期为，或者频率周期为。实际上，这就是抽样系统的抽样频率，而其中的T则是系统的抽样周期。频率响应呈现周期性是离散系统特性区别于连续系统特性的重要特点。因此，只要分析在围的情况，便可分析出系统的整个频率特性。鉴于离散系统频率响应特性的特点，为了表示方便起见，我们通常将其中的用一个变量来代替，即令代入系统函数H(z)中，用函数来表示离散系统的频率响应特性。相应地，用表示幅频特性，而相频特性仍用来表示。应该特别注意的是，虽然这里的变量w仍然称为频率变量，但是它已经不是原来意义上的角频率概念，而实际上是表示角度的概念。我们称之为数字角频率。它与原来角频率的关系为：。也就是说，根据离散系统的系统函数H(z)，令其中的，并且代入0围不同的频率值（实际上是角度值），就可以逐个计算出不同频率时的响应，求出离散系统的频率响应特性。再利用离散系统频率特性的周期性特点（周期为），求出系统的整个频率特性。离散系统的幅频特性曲线和相频特性曲线能够直观地反映出系统对不同频率的输入序列的处理情况。在函数随的变换关系中，在=0附近，反映了系统对输入信号低频部分的处理情况，而在附近，则反映了系统对输入信号高频部分的处理情况。一般来说，分析离散系统频率响应特性就要绘制频率响应曲线，而这是相当麻烦的。虽然可以通过几何矢量法来定性画出频率响应特性曲线，但一般来说这也是很麻烦的。值得庆幸的是，MATLAB为我们提供了专门用于求解离散系统频率响应的函数freqz() ，其调用格式如下：lH,w=freqz(B,A,N) 其中，B和A分别是表示待分析的离散系统的系统函数的分子，分母多项式的向量，N为正整数，返回向量H则包含了离散系统频率响应函数在围的N个频率等分点的值。向量w则包含围的N个频率等分点。在默认情况下N=512。lH,w=freqqz(B,A,N,whole) 其中，B，A和N的意义同上，而返回向量H包含了频率响应函数在围N个频率等分点的值。由于调用freqz()函数只能求出离散系统频率响应的数值，不能直接绘制曲线图，因此，我们可以先用freqz()函数求出系统频率响应的值，然后再利用MATLAB的abs()和angle()函数以及plot()命令，即可绘制出系统在或围的幅频特性和相频特性曲线。例若离散系统的系统函数为，请用MATLAB计算频率围10个等分点的频率响应的样值。 MATLAB程序如下：A=1 0; %分母多项式系数向量B=1 -0.5; %分子多项式系数向量H,w=freqz(B,A,10) %求出对应围10个频率点的频率响应样值运行结果如下：H = 0.5000 0.5245 + 0.1545i 0.5955 + 0.2939i 0.7061 + 0.4045i 0.8455 + 0.4755i 1.0000 + 0.5000i 1.1545 + 0.4755i 1.2939 + 0.4045i 1.4045 + 0.2939i 1.4755 + 0.1545iw = 0 0.3142 0.6283 0.9425 1.2566 1.5708 1.8850 2.1991 2.5133 2.8274例 用MATLAB计算前面离散系统在频率围200个频率等分点的频率响应值，并绘出相应的幅频特性和相频特性曲线。MATLAB程序如下：A=1 0;B=1 -0.5;%H,w=freqz(B,A,200);H,w=freqz(B,A,200,whole); %求出对应围200个频率点的频率响%应样值HF=abs(H); %求出幅频特性值HX=angle(H); %求出相频特性值subplot(2,1,1);plot(w,HF) %画出幅频特性曲线subplot(2,1,2);plot(w,HX) %画出相频特性曲线运行结果如下：三、 实验容1 求出下列离散序列的Z变换（1） 解：syms k z;f = 0.5 k * cos(k * pi / 2);Fz = ztrans(f)结果：Fz =4*z2/(4*z2+1)（2）解：syms k zf = k*(k - 1) * (2/3)k;Fz = ztrans(f)结果：Fz =24*z/(3*z-2)3（3）解：syms k zf = Heaviside(k) - Heaviside(k - 5);Fz = ztrans(f)结果：Fz =1+1/z+1/z2+1/z3+1/z4（4）解：syms k zf = k*(k-1)*(Heaviside(k) - Heaviside(k-5);Fz = ztrans(f)结果：Fz =2/z4*(z2+3*z+6)2已知下列单边离散序列的z变换表达式，求其对应的原离散序列。（1） 解：syms k zFz = (z2 + z + 1)/(z2 + z - 2);f = iztrans(Fz,k)结果：f =-1/2*charfcn0(k)+1/2*(-2)k+1（2）解：syms z kFz = 1 + 1/z + 1/z2 + 1/z3 + 1/z4;f = iztrans(Fz,k)结果：f =charfcn0(k)+charfcn1(k)+charfcn2(k)+charfcn3(k)+charfcn4(k)（3） 解：syms z kFz = 2*(z2 + 3*z + 6)/z4;f = iztrans(Fz)结果：f =12*charfcn4(k)+6*charfcn3(k)+2*charfcn2(k) （4）解：syms z kFz = z*(z2 + z + 1)/(z + 1)*(z-2)*(z+3);f = iztrans(Fz,k)结果：f =7/15*2k-1/6*(-1)k+7/10*(-3)k3 已知离散系统的系统函数H (z)如下，请绘出系统的幅频和相频特性曲线，并说明系统的作用（1） 解：A = 1 7/6 1/3;B = 4 4;H w = freqz(B,A,200,whole);HF = abs(H);HX = angle(H);subplot(2,1,1);plot(w,HF)subplot(2,1,2);plot(w,HX) （2）解：A = 1 0 0.81;B = 1 0 -1;H w = freqz(B,A,200,whole);HF = abs(H);HX = angle(H);subplot(2,1,1);plot(w,HF)subplot(2,1,2);plot(w,HX)4已知描述离散系统的差分方程为：请绘出系统的幅频和相频特性曲线，并说明系统的作用。解：A = 1 0.25 0;B = 1 -1.2 0.35;H w = freqz(B,A,200,whole);HF = abs(H);HX = angle(H);subplot(2,1,1);plot(w,HF)subplot(2,1,2);plot(