离散系统变换域分析—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);