离散系统的频率响应和零极点分布

实验目的

本实验结合理论教学内容。学习和掌握利用变换分析系统特性的方法，加深对系统函数、零极点分布对系统特性的影响、以及离散系统的频率响应分析的理解。

本实验结合理论教材教学内容。学习利用系统的频率响应特性，分析线性时不变（LTI）系统对信号的响应，加深对滤波概念的理解。

主要函数

y= ztrans(x):Z变换函数。输入参数为进行Z变换的表达式。输出参数为Z变换的结果。

y= iztrans(x):Z逆变换函数。输入参数为进行Z逆变换的表达式。输出参数为Z逆变换的结果。

[z,p,k]= residuez(b,a) :有理多项式的部分分式展开函数。输人参数b、a分别表示有理多项式的分子、分母系数。输出参数z、p、k分别表示展开后的零点向量、极点向量及增益向量。

zplane(z,p,k):绘制零极点图。输人参数x、p、k分别表示系统函数的零点向量、极点向量和增益向量。

zplane(b,a):绘制零极点图。输人参数b、a分别表示系统函数的分子、分母系数。

Y= dlsim(b,a,x):计算离散时间系统H(z)的系统响应。输人参数b、a分别表示H(z)的分子、分母系数,x表示输人信号。输出参数y表示系统的响应。当此函数不带输出参数时，可绘制出系统的输出响应曲线。

h=dimpulse(b,a):计算离散时间系统H(z)的单位冲激响应。输人参数b、a分别是系统函数分子、分母系数。输出参数h是离散时间系统的单位冲激响应。此函数当不带输出变量时，可绘制出单位冲激响应曲线。

[h, T]= dimpulse(b,a):计算离散时间系统H(z)的单位冲激响应函数。输人参数b、a分别是系统函数分子、分母系数。输出参数h是单位冲激响应的列向量;T是时间变量。

[h,T]= impz(b,a,N):计算离散时间系统H(z)的单位冲激响应。输人参数b、a分别表是系统函数分子、分母系数。当N为整数向量时，仅计算N指定的整数点上的h(n)。输出参数h是单位冲激响应的列向量;T是时间变量。

[H ,W]= freqz(b,a,N):求离散线性时不变系统的频率响应函数。输人参数b、a分别表示离散线性时不变系统的分子、分母系数,N为在0~π之间的频率抽样点数，默认N=512输出参数H、W分别表示频率响应向量和相应的频率。

[y,x,n]=dinitial(S,B,C,D)计算高散时间LTI系统由初始状态工。所引起的零输人响应，和状态响应x,取样点数由函数自动选取。n为仿真所用的点数。A.B.C.D表示系统的状态空间矢量。

[y.x]= dstep(A,B,C,D) .返回离散LTI系统的单位阶跃响应y向量和时间状态历史记录上向量。A.B.C、D表示系统的状态空间矢量。

[H,W]- freqz(b,a,N):计算出N个频率点上的频率响应存放在H向量中,N个频率存放在向量w. freqz 函数自动将这N个频点均匀设置在频率范围[0,]上。默认w和N时,freqz函数自动选取512个频率点计算。b和a分别为数字滤波器系统函数H(z)的分子和分母多项式系数向量。

y= fiter(b,a,x,xic);对向量x中的数据进行滤波处理，产生输出序列向量y. b和a分别为数字滤波器系统函数H(z)的分F和分母多项式系数向量。xic为系统的等效初始状态输人数组。

xic= filtic(b,a,Y,X):由于初始条件Y和X计算系统的等效于初始状态输入数组

实验内容

 

 

 实验代码

1.clear,clc
2.a=[1,-0.9];
3.b=[1];
4.figure(1);
5.subplot(2,1,1);
6.zplane(b,a)
7.[z,p,k]=residuez(b,a);
8.% h=impz(b,a)
9.% stem(h)
10.
11.figure;
12.[H, W] = freqz(b, a);
13.subplot(2, 1, 1);
14.plot(W / pi, abs(H));
15.title('幅度曲线');
16.xlabel('x/\pi');
17.ylabel('H(ejw)');
18.subplot(2, 1, 2);
19.plot(W / pi, angle(H));
20.title('相位');
21.xlabel('x/\pi');
22.ylabel('相角');
23.
24.
25.figure;
26.impz(b, a);

 

 

 

 

1.clear,clc
2.a = [1, -0.9];
3.b = [1, 0.9];
4.[r,p,k] = residuez(b,a);
5.figure;
6.h = impz(b, a);
7.stem(h);
8.title("单位冲激响应");
9.
10.figure;
11.[H, W] = freqz(b, a);
12.subplot(2, 1, 1);
13.plot(W / pi, abs(H));
14.title('幅度曲线');
15.xlabel('x/\pi');
16.ylabel('H(ejw)');
17.subplot(2, 1, 2);
18.plot(W / pi, angle(H));
19.title('相位');
20.xlabel('x/\pi');
21.ylabel('相角');