matlab ljdt,MATLAB_离散系统z域分析

实验八 离散系统的Z 域分析

一、目的

(1)掌握利用MATLAB 绘制系统零极点图的方法 (2)掌握离散时间系统的零极点分析方法 (3)掌握逆Z 变换概念及MATLAB 实现方法

二、离散系统零极点

线性时不变离散系统可用线性常系数差分方程描述，即

()()N M

i

j

i j a y n i b x n j ==-=-∑∑ (8-1)

其中()y k 为系统的输出序列，()x k 为输入序列。

将式(8-1)两边进行Z 变换的

00

()

()

()()

()

M

j

j

j N

i

i i b z

Y z B z H z X z A z a z

-=-==

==

∑∑ (8-2) 将式(8-2)因式分解后有：

11

()

()()

M

j

j N

i

i z q H z C

z p ==-=-∏∏ (8-3)

其中C 为常数，(1,2,,)j q j M = 为()H z 的M 个零点，(1,2,,)i p i N = 为()H z 的N 个极点。

系统函数()H z 的零极点分布完全决定了系统的特性，若某系统函数的零极点已知，则系统函数便可确定下来。

因此，系统函数的零极点分布对离散系统特性的分析具有非常重要意义。通过对系统函数零极点的分析，可以分析离散系统以下几个方面的特性：

● 系统单位样值响应()h n 的时域特性； ● 离散系统的稳定性； ● 离散系统的频率特性；

三、离散系统零极点图及零极点分析 1．零极点图的绘制

设离散系统的系统函数为

()

()()

B z H z A z =

则系统的零极点可用MATLAB 的多项式求根函数roots()来实现，调用格式为：

p=roots(A)

其中A 为待求多项式的系数构成的行矩阵，返回向量p 则是包含多项式所有根的列向

量。如多项式为231

()48

B z z z =++，则求该多项式根的MATLAB 命令为：

A=[1 3/4 1/8]; P=roots(A)