MSP430 FFT遇到的问题

ikemada

先用 AD600HZ采样100HZ正弦波，得到128个值，没有问题；

查了很多FFT的算法，也看懂了，但是算出来的值还是不对！
求高手们帮忙，看看哪里出了问题：~

//我的FFT ，已经倒位序完成。AD量给到Fft_Real【128】中

1. void FFT(void)
   
2. {
   
3. int NN=7; //这里因为128是2的7次方，如果是计算256点，则是2的8次方，N就是8，如果是512点则N＝9，如此类推
   
4. 
   
5. 
   
6. unsigned int NUM_FFT=128; //这里要算多少点的fft就赋值多少，值只能是2的N次方
   
7.     for( i=1; i<=NN; i++)             /* for(1) */
   
8.     {
   
9.         b=1;
   
10.         b <<=(i-1);                 //蝶式运算，用于计算 隔多少行计算。例如第一级 1和2行计算，，，第二级
    
11.         for( j=0; j<=b-1; j++)       /* for (2) */
    
12.         {
    
13.             p=1;
    
14.             p <<= (NN-i);
    
15.             p = p*j;
    
16.             for( k2=j; k2<NUM_FFT; k2=k2+2*b)   /* for (3) 基二fft */
    
17.             {
    
18.                 Temp_Real = Fft_Real[k2]; Temp_Imag = Fft_Image[k2]; temp = Fft_Real[k2+b];
    
19.                 Fft_Real[k2] = Fft_Real[k2] + ((Fft_Real[k2+b]*COS_TAB[p])>>7) + ((Fft_Image[k2+b]*SIN_TAB[p])>>7);
    
20.                 Fft_Image[k2] = Fft_Image[k2] - ((Fft_Real[k2+b]*SIN_TAB[p])>>7) + ((Fft_Image[k2+b]*COS_TAB[p])>>7);
    
21.                 Fft_Real[k2+b] = Temp_Real - ((Fft_Real[k2+b]*COS_TAB[p])>>7) - ((Fft_Image[k2+b]*SIN_TAB[p])>>7);
    
22.                 Fft_Image[k2+b] = Temp_Imag + ((temp*SIN_TAB[p])>>7) - ((Fft_Image[k2+b]*COS_TAB[p])>>7);
    
23. //                //移位，防止溢出。结果已经是本值的1/64
    
24.               Fft_Real[k2] >>= 1;
    
25.                 Fft_Image[k2] >>= 1;
    
26.                Fft_Real[k2+b]  >>= 1;
    
27.                 Fft_Image[k2+b]  >>= 1;
    
28.             }
    
29.         }
    
30.     }
    
31.    }

复制代码

还有一个问题，计算模值：

1. for(n=1;n<128;n++)
   
2.         {
   
3.                 //Fft_Real=conventionsignal()
   
4.                 if(Fft_Real[n]<0) Fft_Real[n] = -Fft_Real[n];
   
5.                 if(Fft_Image[n]<0) Fft_Image[n] = -Fft_Image[n];
   
6.                 result[n]=sqrt((Fft_Real[n]*Fft_Real[n])+(Fft_Image[n]*Fft_Image[n]));
   
7.         }

复制代码

我在FFT之后计算各频率模制，这里要不要都取正？

ikemada

附图给出FFT后实部和虚部的值，不知道哪里不对？

charlesg

建议楼主传一个频域的图，我们可以看看你的fft结果和理想值的差异

计算模值不要取正，因为如果是负值，你取正再和原来乘，又变成负的了。要的是一个正的，就是平方效果就好

ikemada

charlesg 发表于 2014-7-26 11:49
建议楼主传一个频域的图，我们可以看看你的fft结果和理想值的差异

计算模值不要取正，因为如果是负值，你 ...

谢谢！麻烦啦~
RESULT没有开方，就是 实部^2+虚部^2；
表示完全看不到频率响应啊。。。 AD的数据倒很清楚 0 0

charlesg

ikemada 发表于 2014-7-26 12:52
谢谢！麻烦啦~
RESULT没有开方，就是 实部^2+虚部^2；
表示完全看不到频率响应啊。。。 AD的数据倒很清楚 ...

还没求log，result怎么会出现负值？。。。

ikemada

charlesg 发表于 2014-7-26 15:03
还没求log，result怎么会出现负值？。。。

我是这么写的
for(n=1;n<128;n++)
{
result[n]=((Fft_Real[n]*Fft_Real[n])+(Fft_Image[n]*Fft_Image[n]))；
}

yulee

result...这么算，不应该是负的啊

谦虚的行者

楼主解决问题了吗