实验二 z 变换和系统频域特性
一、画出因果系统 y(n)=0.9y(n-1)+x(n)的零极点图和频率响应曲线,并
求单位抽样响应.
N=50;
b=[1 0];
a=[1 -0.9];
subplot(2,2,1);zplane(b,a);xlabel('Real part');ylabel('Imaginary part');grid
on;%根据系统函数返回零极点图
x=[1 zeros(1,N-1)];
n=0:N-1;
subplot(2,2,2);stem(n,filter(b,a,x));xlabel('n');ylabel('h(n)');grid on;%根据
系统函数返回单位冲激响应(离散的)
[H,w]=freqz(b,a,N);%根据系统函数返回 0~2pi 间等间隔的 N 个频率
w 相应的频率响应 H
subplot(2,2,3);plot(w,abs(H));xlabel('w');ylabel('|H(exp(jw))|');grid on;
subplot(2,2,4);plot(w,angle(H));xlabel('w');ylabel('arg(H(exp(jw)))');grid
on;
二、求
���
��
� z
zz
z
zX 1,
143
)(
2
的 z 反变换.
b=[0 1 0];%X(z)的分子多项式系数
a=[3 -4 1];%X(z)的分母多项式系数
[h n]=impz(b,a,30)%由系统函数的分子多项式 b 和分母多项式 a 生成
评论0