FFT算法C语言实现

{// 倒叙
  unsigned short invert_temp;
Invert_Bdata invert_bdata_temp1,invert_bdata_temp2;
for(invert_temp = 0;invert_temp < 256;invert_temp ++)
{
invert_bdata_temp1.all = invert_temp;

invert_bdata_temp2.bit.bit0 = invert_bdata_temp1.bit.bit7;
invert_bdata_temp2.bit.bit1 = invert_bdata_temp1.bit.bit6;
invert_bdata_temp2.bit.bit2 = invert_bdata_temp1.bit.bit5;
invert_bdata_temp2.bit.bit3 = invert_bdata_temp1.bit.bit4;

invert_bdata_temp2.bit.bit4 = invert_bdata_temp1.bit.bit3;
invert_bdata_temp2.bit.bit5 = invert_bdata_temp1.bit.bit2;
invert_bdata_temp2.bit.bit6 = invert_bdata_temp1.bit.bit1;
invert_bdata_temp2.bit.bit7 = invert_bdata_temp1.bit.bit0;

ch1_run_process_varial.invert_data[invert_temp] = ch1_run_process_varial.fft_assy[invert_bdata_temp2.all];
}
}
{// 蝶形算法
unsigned int m,s1,s2,s3,i,j,k;
float wr,wi,qr,qi;

for(i = 0;i < 256;i ++)
{
ch1_run_process_varial.rek[i] = ch1_run_process_varial.invert_data[i]; 
ch1_run_process_varial.imk[i] = 0;
}