x4x52.m

来自「dsp学习中matlab试验源代码」· M 代码 · 共 17 行

n=0:16; 
m=0:15;
a=cos((pi/4)*n);    
b=sin((pi/8)*n);   
c=a+b*j;
g=fft(c,16);
for i=0:15
  z(i+1)=g(16-i);
end
p=conj(z);
d=0.5*(g+p); 
e=-0.5*j*(g-p); 
f=ifft(d,16);
subplot(2,1,1);stem(m,f);
h=ifft(e,16);
subplot(2,1,2);stem(m,h);