实验二-离散时间信号与系统的Z变换分析

1、实验二离散时间信号与系统的Z变换分析一、 实验目的1、熟悉离散信号z变换的原理及性质2、熟悉常见信号的 Z变换3、了解正/反Z变换的MATLAB实现方法4、了解离散信号的Z变换与其对应的理想抽样信号的傅氏变换和拉氏变换之间的关系5、了解利用MATLAB实现离散系统的频率特性分析的方法二、实验原理1、正/反Z变换Z变换分析法是分析离散时间信号与系统的重要手段。如果以时间间隔Ts对连续时间信号f(t)进行理想抽样，那么，所得的理想抽样信号f (t)为：f(t) f(t)* Ts(t)f(t)* (t kTs)k理想抽样信号f (t)的双边拉普拉斯变换F为：stksTF (s) f (t)* (t

2、kTs) e dtf(kTs)e skk若令f(kTs)f(k) , z esTs , 那么f (t)的双边拉普拉斯变换 F (s)为：F (s) f(k)zkF(z)zesTsk则离散信号f(k)的Z变换定义为：kF(z) f(k)zk从上面关于Z变换的推导过程中可知，离散信号f(k)的Z变换F(z)与其对应的理想抽样信号f (t)的拉氏变换F (s)之间存在以下关系：F (s) F(z)zesTs同理，可以推出离散信号f(k)的Z变换F(z)和它对应的理想抽样信号f (t)的傅里叶变换之间的关系为F (j ) F(z)z如果已知信号的Z变换F(z),要求出所对应的原离散序列f(k),就需要

3、进行反Z变换,-1 一 k 1f(k)否？F(z)zk1dz反Z变换的定义为:下:其中，C为包围F(z)zk 1的所有极点的闭合积分路线。在MATLAB语言中有专门对信号进行正反Z变换的函数ztrans()和itrans()其调用格式分别如F=ztrans( f )F=ztrans(f,v)对f(n)进彳t Z变换，其结果为F(z) 对f(n)进彳t Z变换，其结果为F(v)F=ztrans(f,u,v)f=itrans ( F )对f(u)进彳T Z变换，其结果为F(v)对F(z)进彳T Z反变换，其结果为f(n)f=itrans(F,u)f=itrans(F,v,u)对F(z)进彳T Z反

4、变换，其结果为f(u)对F(v)进彳T Z反变换，其结果为f(u)注意：在调用函数ztran()及iztran()之前，要用syms命令对所有需要用到的变量行说明，即要将这些变量说明成符号变量。t,u,v,w )等进例.用MATLAB求出离散序列f(k) (0.5)k (k)的Z变换MATLAB程序如下:syms k zf=0.5Ak;%定义离散信号Fz=ztrans(f)运行结果如下：%对离散信号进行Z变换Fz =2*z/(2*z-1)例.已知一离散信号的Z变换式为F(z)2z2z 1,求出它所对应的离散信号f(k)MATLAB程序如下:syms k zFz=2* z/(2*z-1);%定义

5、Z变换表达式fk=iztrans(Fz,k)运行结果如下：%求反Z变换fk =(1/2件例：求序列f (k)(k1) (t 4)的Z变换.clc;clear all syms n hn=sym( kroneckerDelta(n, 1) + kroneckerDelta(n, 2)+ kroneckerDelta(n, 3)Hz=ztrans(hn)Hz=simplify(Hz)2、离散系统的频率特性同连续系统的系统函数 H(s)类似，离散系统的系统函数 H(z)也反映了系统本身固有的特性。对于离 散系统来说，如果把其系统函数 H(z)中的复变量z换成ej ej Ts (其中Ts),那么所得的

6、函数H (ej )就是此离散系统的频率响应特性，即离散时间系统的频率响应为：H(ej) |H(ej )gej ( ) H(矶 ej其中，H(ej )称为离散系统的幅频特性，()称为系统的相频特性。同连续系统一样，离散时间系统的幅频特性也是频率的偶函数，相频特性也是频率的齐函数。由于ej是频率的周期函数，所以离散系统的频率响应特性也是频率的周期函数，其周期为2 ,2或者角频率周期为t 。实际上，这就是抽样系统的抽样频率，而其中的T则是系统的抽样周期。Ts频率响应呈现周期性是离散系统特性区别于连续系统特性的重要特点。因此，只要分析H (ej )在 2范围内的情况，便可分析出系统的整个频率特性。H

7、(ej )函数来表示离散系统的频率响应特性，H (ej )表示幅频特性，而相频特性仍用()来表示。应该特别 注立的是，虽然这里的变量仍然称为频率变量， 但是它已经不是原来意义上的角频率概念,而实际上是表示 鱼直的概念。我们称之为 数字频率。它与原来角频率的关系为：Ts。也就是说，根据离散系统的系统函数H(z),令其中的z ej ,并且代入。2范围内不同的频率值(实际上是角度值)，就可以逐个计算出不同频率时的响应，求出离散系统的频率响应特性。再利用离散系统频率特性的周期性 特点(周期为2 ),求出系统的整个频率特性。离散系统的幅频特性曲线和相频特性曲线能够直观地反映出系统对不同频率的输入序列的处

8、理情况。 在函数H(ej )随 的变换关系中，在=0附近，反映了系统对输入信号低频部分的处理情况，而在 =附近，则反映了系统对输入信号高频部分的处理情况。一般来说，分析离散系统频率响应特性就要绘制频率响应曲线，而这是相当麻烦的。虽然可以通过 几何矢量法来定性画出频率响应特性曲线，但一般来说这也是很麻烦的。值得庆幸的是，MATLAB为我们提供了专门用于求解离散系统频率响应的函数freqz(),其调用格式如下：H,w=freqz(B,A,N) 其中，B和A分别是表示待分析的离散系统的系统函数的分子，分母多 项式的向量，N为正整数，返回向量 H则包含了离散系统频率响应函数H(ej )在0 范围内的N

9、个频率等分点的值。向量 则包含0 范围内的N个频率等分点。在默认情况下N=512 。H,w=freqz(B,A,N,whole)其中，B, A和N的意义同上，而返回向量 H包含了频率响应函数H(ej )在02范围内N个频率等分点的值。由于调用freqz()函数只能求出离散系统频率响应的数值，不能直接绘制曲线图，因此，我们可以先 用freqz()函数求出系统频率响应的值，然后再利用MATLAB的abs()和angle()函数以及plot()命令，即可绘制出系统在0 或02 范围内的幅频特性和相频特性曲线。例.若离散系统的系统函数为H(z)z 0.5，请用MATLAB计算0 频率范围内10个等分点

10、的频率响应H(ej )的样值。%分母多项式系数向量%分子多项式系数向量%求出对应0 范围内10个频率点的频率响应样值MATLAB程序如下：A=1 0;B=1 -0.5;H,w=freqz(B,A,10)运行结果如下：H =0.50000.5245 + 0.1545i0.5955 + 0.2939i0.7061 + 0.4045i0.8455 + 0.4755i1.0000 + 0.5000i1.1545 + 0.4755i1.2939 + 0.4045i1.4045 + 0.2939i1.4755 + 0.1545i w =00.31420.62830.94251.25661.57081.88

11、502.19912.51332.8274例.用MATLAB计算前面离散系统在 02频率范围内200个频率等分点的频率响应值，并绘出相应的幅频特性和相频特性曲线。MATLAB程序如下：A=1 0;B=1 -0.5;H,w=freqz(B,A,200);H,w=freqz(B,A,200,whole);%求出对应02 范围内200个频率点白频率响 应样值HF=abs(H);%求出幅频特性值HX=angle(H);%求出相频特性值subplot(2,1,1);plot(w,HF)subplot(2,1,2);plot(w,HX)%画出幅频特性曲线%画出相频特性曲线运行结果如下：运行结果分析：从该系统

12、的幅频特性曲线可以看出，该系统呈高通特性，是一阶高通滤波器。实验内容1 .求出下列离散序列的Z变换fi(k) (2)kcos(7) (k)f2(k) k(k 1)(|)k (k)f3(k)(k) (k 5)f4(k) k(k 1) (k) (k 5)已知下列单边离散序列的 z变换表达式，求其对应的原离散序列。 Fi(z)2z-2 z F2(z) F3(z)2(z23z 6) F4(z)z(z2z 1)(z 1)(z 2)(z 3)已知离散系统的系统函数H (z)如下，请绘出系统的幅频和相频特性曲线，并说明系统的作用 H(z) H(z)4z 4(z 2)(z i)z2 1z2 0.81已知描述离

13、散系统的差分方程为：y(k) 1.2y(k 1) 0.35y(k 2) e(k) 0.25e(k 1)请绘出系统的幅频和相频特性曲线，并说明系统的作用。预习要求、 熟悉正反 z 变换的意义及用 MATLAB 软件实现的方法、 熟悉离散系统的频率响应特性及用 MATLAB 软件实现的方法、 编写 MATLAB 程序实验报告要求、简述实验目的及实验原理、计算相应 z 变换或反 z 变换的理论值，并与实验结果进行比较、记录离散系统的频率响应特性曲线，分析系统作用、写出程序清单、收获与建议% 参考程序三1.clc;clear all% 定义离散信号%对离散信号进行Z 变换% 定义离散信号%对离散信号进行Z 变换syms k zf1=0.5Ak*cos(k*pi/2);Fz1=ztrans(f1)实验二 1. f2=k*(k-1)*(2/3)Ak;Fz2=ztrans(f2)实验二 1. f3=sym(kroneckerDel