example3_4.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

n=[0:1:3];
x=[1,1,1,1];
N=4;
X=dft(x,N);
magX=abs(X),phaX=angle(X)*180/pi
subplot(2,1,1);stem(n,magX);
xlabel('k');ylabel('|X(k)|');title('amplitude of DFT');
subplot(2,1,2);stem(n,phaX);
xlabel('k');ylabel('degree');title('phase of DFT');