Question

单片机FFT算法，256点，看看哪里不对

for( i=1; i<=7; i++) /* for(1) */
{
b=1;
b <<=(i-1); //碟式运算，用于计算隔多少行计算例如 第一极 1和2行计算
for( j=0; j<=b-1; j++) /* for (2) */
{
p=1;
p <<= (7-i);
p = p*j;
for( k=j; k<128; k=k+2*b) /* for (3) */
{
Temp_Real=Fft_Real[k];
Temp_Imag=Fft_Image[k];
temp=Fft_Real[k+b];
Fft_Real[k]=Fft_Real[k]+((Fft_Real[k+b]*COS_TAB[p])>>7)+((Fft_Image[k+b]*SIN_TAB[p])>>7);
Fft_Image[k]=Fft_Image[k]-((Fft_Real[k+b]*SIN_TAB[p])>>7)+((Fft_Image[k+b]*COS_TAB[p])>>7);
Fft_Real[k+b]=Temp_Real-((Fft_Real[k+b]*COS_TAB[p])>>7)-((Fft_Image[k+b]*SIN_TAB[p])>>7);
Fft_Image[k+b]=Temp_Imag+((temp*SIN_TAB[p])>>7)-((Fft_Image[k+b]*COS_TAB[p])>>7); // 移位.防止溢出. 结果已经是本值的 1/64
Fft_Real[k] >>= 1;
Fft_Image[k] >>= 1;
Fft_Real[k+b] >>= 1;
Fft_Image[k+b] >>= 1;
}
}
}

Answer 1

结果不能直接左移吧？如果都是浮点数的话，那该用除法。

追问

不是浮点数，采样的时候强制INT了。这个程序，在运算256个点的时候，不是要将程序中的7改为8,128改为256，这样理解对不？

追答

整形的计算出的结果根本没法用啊……我以前试验过，结果放弃了。

追问

好的 我试试