example3_12b.m

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

n=[0:1:99];x=cos(0.48*pi*n)+cos(0.52*pi*n);
subplot(2,1,1);stem(n,x);title('signal 0<=n<=99');
X=dft(x,100);
magX=abs(X(1:51));
k=0:1:50;
w=2*pi*k/100;
subplot(2,1,2);stem(w/pi,magX);
title('DFT');