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

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

一、 实验目的

1、 熟悉离散信号Z 变换的原理及性质

2、 熟悉常见信号的 Z 变换

3、 了解正/反Z 变换的MATLAB 实现方法

4、 了解离散信号的Z 变换与其对应的理想抽样信号的傅氏变换和拉氏变换之间的关系

5、 了解利用MATLAB 实现离散系统的频率特性分析的方法

二、 实验原理

1、正/反Z 变换

Z 变换分析法是分析离散时间信号与系统的重要手段。

如果以时间间隔T s 对连续时间信号f (t)进行理

想抽样，那么，所得的理想抽样信号

f (t)为：

则离散信号f (k )的Z 变换定义为:

k

F(z) f(k)z

k

从上面关于Z 变换的推导过程中可知，离散信号 f (k )的Z 变换 拉氏变换F (s)之间存在以下关系:

F (s) F(z)

同理，可以推出离散信号 f (k )的Z 变换F(z)和它对应的理想抽样信号

F(j ) F(z)

z e j Ts

f (t)

f(t)* Ts (t) f (t)

(t kT s )

理想抽样信号 f (t)的双边拉普拉斯变换

F (s)为:

F (s)

f(t)*

k

(t kT s )

e st

dt f (kT s )e ksT s

k

若令f (kT s )

f(k) ， z e sTi ， 那么f (t)的双边拉普拉斯变换

F (s)为:

F (s)

f(k)z k

FO zesI

F(z)与其对应的理想抽样信号

f (t)的

f (t)的傅里叶变换之间的关系为

如果已知信号的Z变换F(z),要求出所对应的原离散序列f(k),就需要进行反Z变换,

其中，C为包围F(z)z k1的所有极点的闭合积分路线。

在MATLAB语言中有专门对信号进行正反Z变换的函数ztrans()和itrans() 下:

F=ztrans( f ) 对f(n)进行Z变换，其结果为 F(z)

F=ztrans(f,v) 对f(n)进行Z变换，其结果为F(v)

F=ztrans(f,u,v) 对f(u)进行Z变换，其结果为F(v)

f=itrans ( F ) 对F(z)进行Z反变换，其结果为f(n)

f=itrans(F,u) 对F(z)进行Z反变换，其结果为 f(u) f=itrans(F,v,u ) 对F(v)进

行Z反变换，其结果为 f(u)

注意：在调用函数ztran()及iztran()之前，要用syms命令对所有需要用到的变量

行说明，即要将这些变量说明成符号变量。

k

例①.用MATLAB求出离散序列f(k) (0.5) (k)的Z变换

MATLAB程序如下:

syms k z

f=0.5A k; %定义离散信号

Fz =

2*z/(2*z-1)

clc;clear all syms n hn=sym( 'kroneckerDelta(n, 1) + kroneckerDelta(n, 2)+ kroneckerDelta(n, 3)'

Hz=ztra ns(h n)

Hz=simplify(Hz)反Z变换的定义为:

f(k) 21

j?F(z)z k1dz

其调用格式分别如

t,u,v,w )等进

Fz=ztra

ns(f)

运行结果如下：

%对离散信号进行Z变换

例②.已知一离散信号的Z变换式为F(z)

2z

2z 1

，求出它所对应的离散信号f(k)

MATLAB程序如下：

syms k z

Fz=2* z/(2*z-1);

fk=iztra ns(F z,k)

运行结果如下：

fk = %定义Z变换表达式%求反Z变换

例③：求序列f (k)(k 1) (t 4)的Z 变换.

2、离散系统的频率特性

同连续系统的系统函数 H (s)类似，离散系统的系统函数

散系统来说，如果把其系统函数

H (z )中的复变量z 换成e j

e

j T s

(其中

H (e j )就是此离散系统的频率响应特性，即离散时间系统的频率响应为:

其中， H(e j

)称为离散系统的幅频特性，

()称为系统的相频特性。同连续系统一样，离散时间系统

的幅频特性也是频率的偶函数，相频特性也是频率的齐函数。

由于e j

是频率的周期函数，所以离散系统的频率响应特性也是频率

的周期函数，其周期为2 ， 2

或者角频率周期为

T

——。实际上，这就是抽样系统的抽样频率，而其中的

T 则是系统的抽样周期。

T

s

频率响应呈现周期性是离散系统特性区别于连续系统特性的重要特点。 因此，只要分析H(e j

)在| | 2

范围内的情况，便可分析出系统的整个频率特性。

H (e j )函数来表示离散系统的频率响应特性，

|H(e j )表示幅频特性，而相频特性仍用

()来表

示。应该特别 注意的是,虽然这里的变量

仍然称为频率变量， 但是它已经不是原来意义上的

角频率概念,

而实际上是表示 角度的概念。我们称之为 数字频率。它与原来角频率的关系为:

T s 。也就是说，根

据离散系统的系统函数 H (z ),令其中的z e j

,并且代入0〜2范围内不同的频率值(实际上是角度值)， 就可以逐个计算出不同频率时的响应，求出离散系统的频率响应特性。再利用离散系统频率特性的周期性 特点(周期为2

)，求出系统的整个频率特性。

离散系统的幅频特性曲线和相频特性曲线能够直观地反映出系统对不同频率的输入序列的处理情况。 在函数H(e j

)随 的变换关系中，在=0附近，反映了系统对输入信号低频部分的处理情况， 而在 =

附近，则反映了系统对输入信号高频部分的处理情况。

一般来说，分析离散系统频率响应特性就要绘制频率响应曲线，而这是相当麻烦的。虽然可以通过 几何矢量法来定性画出频率响应特性曲线， 但一般来说这也是很麻烦的。值得庆幸的是，MATLAB 为我们

提供了专门用于求解离散系统频率响应的函数

freqz()，其调用格式如下：

[H ,w]=freqz(B,A,N)

其中，B 和A 分别是表示待分析的离散系统的系统函数的分子，分母多

项式的向量，N 为正整数，返回向量 H 则包含了离散系统频率响应函数 H(e j )在0~ 范围内的N 个频

率等分点的值。向量

则包含0~ 范围内的N 个频率等分点。在默认情况下

N=512。

[H ,w]=freqz(B,A,N,'whole')

其中，B, A 和N 的意义同上，而返回向量 H 包含了频率响应

函数H(e j

)在0~2范围内N 个频率等分点的值。

由于调用freqz()函数只能求出离散系统频率响应的数值，不能直接绘制曲线图，因此，我们可以先 用freqz()函数求出系统频率响应的值，然后再利用

MATLAB 的abs()和angle()函数以及plot()命令，即

可绘制出系统在0 ~ 或0 ~ 2 范围内的幅频特性和相频特性曲线。

H(e j ) |H(e j ) ge j ()

H(z )z

H (z )也反映了系统本身固有的特性。对于离

T s )，那么所得的函数