example3_12b.m

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

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');