数字信号处理matlab学习（四）

离散系统变换域分析—Z变换一

一.实验用函数

1、ztrans

功能：求无限长序列的z变换。

调用格式：

X=ztrans(x)；用于求解无限长序列的z变换。

2、iztrans

功能：求函数X(z)的z反变换。

调用格式：

x= iztrans (X)

3、syms

功能：定义符号对象。

调用格式：

sysm t,i,x：将变量t,i,x声明为符号变量。

4、residuez

功能：有理多项式的部分分式展开。

调用格式

[r,p,c]=residuez(b,a)；求解以系数向量b，a表示的系统函数的部分分式展开。

5、zplane

功能：绘制零极点分布图。

调用格式

zplane(z,p)：绘制由列向量z确定的零点、列向量p确定的极点构成的零极点分布图。

zplane(b,a)：绘制由行向量b,a构成的系统函数确定的零极点分布图。

[hz,hp,ht]=zplane(z,p)：获得三个句柄向量：hz为零点线句柄；hp为极点线句柄；ht为坐标轴、单位圆及文本对象的句柄。  

二.参考实例

1、求序列的z变换

syms n;
x = (n*(n-1))/2
X = ztrans(x)

2、求函数的z反变换

syms z;
X = z/(z-1)^2;
x = iztrans(X)

3.用部分分式法求z反变换（|z|>1）

a = [1 -1.5 0.5];
b = [1];
[r,p,c] = residuez(b,a)

4.、已知一个离散系统的函数和输入序列，求系统在变换域的响应Y(z)及时间域的响应y(n)。

syms z;
H = z^2/(z^2-1.5*z+0.5);
X = z/(z-1);
Y = H*X
y = iztrans(Y)

5.用部分分式法求系统函数的z反变换，并用图形与impz求得的结果相比较。

a = [1 0 0.34319 0 0.60439]; 
b = [0.1321 0 0.3963 0 0.3962];
[r,p,c] = residuez(b,a)
N = 50;
n = 0:N-1;
h = r(1)*p(1).^n+r(2)*p(2).^n+r(3)*p(3).^n+r(4)*p(4).^n+c(1).*[n==0];
subplot(1,2,1);  stem(n,real(h));
subplot(1,2,2);  impz(b,a,n);

6、系统极点的位置对系统响应的影响

1）研究z右半平面的实数极点对系统响应的影响。

已知系统的零-极点增益模型

a = [1 -0.85];
b = 1;
N = 40;
n = 1:N-1;
subplot(212);impz(b,a,n);
subplot(211);zplane(b,a);

2）研究z右半平面的复数极点对系统响应的影响

求系统的零极点分布图以及系统的冲激响应，并判断系统的稳定性。

b = [1 -0.3];
a = [1 -1 0.74];
N = 40;
n = 1:N-1;
subplot(212);impz(b,a,n);
subplot(211);zplane(b,a);