homew3.m

来自「西交大刘树棠写的《数字信号处理》那本书的所有源代码」· M 代码 · 共 11 行

n=[0:11];
x=cos(n*pi/6);
k=-200:200;w=(pi/100)*k;
X=x*exp(-j*pi/100.*(n'*k));%用矩阵的变换完成对序列x的DTFT，取两百个点观察
magX=abs(X);angX=angle(X);
subplot(211);plot(w/pi,magX);grid;
title('幅频响应');
xlabel('Frequency');ylabel('amplitude');
subplot(212);plot(w/pi,angX);grid;
title('相频响应');
xlabel('Frequency');ylabel('angle');