ex31b.m

来自「数字信号分析处理方面的好教材,有DSP的MATLAB源程序.」· M 代码 · 共 21 行

%
% b) 4-point DFT
x = [1,1,1,1];
N = 4; w1 = 2*pi/N; k = 0:N-1;
X = dft(x,N);
magX = abs(X), phaX = angle(X)*180/pi

w = [0:1:500]*2*pi/500;  %取Δt=1/500
[H] = freqz(x,1,w);
magH = abs(H); phaH = angle(H); phaH(126)=-47.5841*pi/180;

subplot(2,1,1);plot(w*N/(2*pi),magH,'--'); 
axis([-0.1,4.1,-1,5]); hold on
stem(k,magX);
ylabel('|X(k)|'); title('DFT的幅度: N=4');text(4.3,-1,'k')
hold off
subplot(2,1,2);plot(w*N/(2*pi),phaH*180/pi,'--');
axis([-0.1,4.1,-200,200]); hold on
stem(k,phaX);
ylabel('度'); title('DFT相角: N=4');text(4.3,-200,'k')