n=0:16; 
m=0:15;
a=cos((pi/4)*n);    
b=sin((pi/8)*n);   
c=a+b;
g=fft(c,16);
subplot(4,1,2); 
stem(m,abs(g))  ; title('X4(t)+X5(t)的序列16點(diǎn)的FFT');

d=real(g); 
e=imag(g); 
f=ifft(d);
subplot(2,1,1);stem(m,f);
title('用X(n)=X4(n)+X5(n)經(jīng)FFT后還原的X4(n)');
i=e*j;
h=ifft(i);
subplot(2,1,2);stem(m,h);
title('用X(n)=X4(n)+X5(n)經(jīng)FFT后還原的X5(n)');