의 스펙트럼 정보를 구하여 보자.

여기서는 유한개의 샘플을 이용하여 이산시간 퓨리에 변환을 실행한다. 이산시간 퓨리에 변환을 하는 fft 함수를 사용한다. 명령어를 아래와 같이 편집한다.

n=[0:99:1];  
x=cos(0.48*pi*n)+cos(0.52*pi*n); 
n3=[0:99:1]; 
y3=[x(1:10:1) zeros(1,90)]; 
subplot(2,1,1); 
stem(n3,y3); 
title("signal x(n), 0 <= n <= 9 + 90 zeros"); 
xtitle("n")  
Y3=fft(y3); 
magY3=abs(Y3(1:51:1)); 
k3=0:50:1; 
w3=2*pi/100*k3; 
subplot(2,1,2); 
plot(w3/pi,magY3); 
title("DTFT Magnitude"); 
xtitle("frequency in pi units") 

위와 같은 명령을 실행하면 의 스펙트럼 정보를 다음 그림과 같이 확인할 수 있다.