离散系统的变换域分析【数字信号处理四】

一、求系统H(z)的零、极点、幅频响应和相位响应

num=[0.0528,0.0797,0.1295,0.1295,0.797,0.0528];
den=[1,-1.8107,2.4947,-1.8801,0.9537,-0.2336];
k=2048;
w=0:pi/k:pi;
h=freqz(num,den,w);
figure(1);
zplane(num,den); 
title("零极点图")
figure(2);
subplot(2,1,1);
plot(w/pi,abs(h));
title('幅度谱'); 
xlabel('\omega/\pi');ylabel('幅值');
subplot(2,1,2);
plot(w/pi,angle(h));
title('相位谱'); 
xlabel('\omega/\pi');ylabel('弧度');

二、用Matlab验证DFT运算的对称性质

N=32;
n=0:N-1;
nn=mod(-n,N)+1;
xn=cos(2*pi*6*n/N)+1i*sin(2*pi*6.8*n/N);
% x[n]  -->  X[k]
Xk=fft(xn);
subplot(2,7,1)
stem