void IFFT()
{ int16 temp;
fft.ipcbptr = ipcb; // FFT computation buffer
fft.magptr = mag; // Store mag. square in separate buff
fft.winptr = (long *)win; // Window coefficient array
fft.init(&fft); // Twiddle factor pointer initialization
for(temp=0;temp<FFT_N;temp++)
{
ipcb[temp*2+1]=-ipcb[temp*2+1]; //求共轭复根
}
/***************再求FFT变换*************************/
CFFT32_brev2(ipcb,ipcb,FFT_N);
fft.izero(&fft);
fft.calc(&fft);
/****************************************************/
for(temp=0;temp<FFT_N;temp++)
{
ipcb[temp*2]=((ipcb[temp*2]));
ipcb[temp*2+1]=(-ipcb[temp*2+1]);
}
}
或者谁帮我看看这个程序错在哪里？？？
|
楼主
