搜索
bottom↓
回复: 7

MSP430 FFT遇到的问题

[复制链接]
(236133999)

出0入0汤圆

发表于 2014-7-26 11:37:22 | 显示全部楼层 |阅读模式
先用 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. unsigned int NUM_FFT=128; //这里要算多少点的fft就赋值多少,值只能是2的N次方
  5.     for( i=1; i<=NN; i++)             /* for(1) */
  6.     {
  7.         b=1;
  8.         b <<=(i-1);                 //蝶式运算,用于计算 隔多少行计算。例如第一级 1和2行计算,,,第二级
  9.         for( j=0; j<=b-1; j++)       /* for (2) */
  10.         {
  11.             p=1;
  12.             p <<= (NN-i);
  13.             p = p*j;
  14.             for( k2=j; k2<NUM_FFT; k2=k2+2*b)   /* for (3) 基二fft */
  15.             {
  16.                 Temp_Real = Fft_Real[k2]; Temp_Imag = Fft_Image[k2]; temp = Fft_Real[k2+b];
  17.                 Fft_Real[k2] = Fft_Real[k2] + ((Fft_Real[k2+b]*COS_TAB[p])>>7) + ((Fft_Image[k2+b]*SIN_TAB[p])>>7);
  18.                 Fft_Image[k2] = Fft_Image[k2] - ((Fft_Real[k2+b]*SIN_TAB[p])>>7) + ((Fft_Image[k2+b]*COS_TAB[p])>>7);
  19.                 Fft_Real[k2+b] = Temp_Real - ((Fft_Real[k2+b]*COS_TAB[p])>>7) - ((Fft_Image[k2+b]*SIN_TAB[p])>>7);
  20.                 Fft_Image[k2+b] = Temp_Imag + ((temp*SIN_TAB[p])>>7) - ((Fft_Image[k2+b]*COS_TAB[p])>>7);
  21. //                //移位,防止溢出。结果已经是本值的1/64
  22.               Fft_Real[k2] >>= 1;
  23.                 Fft_Image[k2] >>= 1;
  24.                Fft_Real[k2+b]  >>= 1;
  25.                 Fft_Image[k2+b]  >>= 1;
  26.             }
  27.         }
  28.     }
  29.    }
复制代码



还有一个问题,计算模值:
  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之后计算各频率模制,这里要不要都取正?

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册

x
(236133843)

出0入0汤圆

 楼主| 发表于 2014-7-26 11:39:58 | 显示全部楼层
附图给出FFT后实部和虚部的值,不知道哪里不对?
(236133283)

出0入0汤圆

发表于 2014-7-26 11:49:18 | 显示全部楼层
本帖最后由 charlesg 于 2014-7-26 11:52 编辑

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

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

出0入0汤圆

 楼主| 发表于 2014-7-26 12:52:29 | 显示全部楼层
本帖最后由 ikemada 于 2014-7-26 12:56 编辑
charlesg 发表于 2014-7-26 11:49
建议楼主传一个频域的图,我们可以看看你的fft结果和理想值的差异

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


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

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册

x
(236121653)

出0入0汤圆

发表于 2014-7-26 15:03:08 | 显示全部楼层
ikemada 发表于 2014-7-26 12:52
谢谢!麻烦啦~
RESULT没有开方,就是 实部^2+虚部^2;
表示完全看不到频率响应啊。。。 AD的数据倒很清楚 ...

还没求log,result怎么会出现负值?。。。
(236115697)

出0入0汤圆

 楼主| 发表于 2014-7-26 16:42:24 | 显示全部楼层
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]));
}
(235322127)

出0入0汤圆

发表于 2014-8-4 21:08:34 | 显示全部楼层
result...这么算,不应该是负的啊
(205591480)

出0入0汤圆

发表于 2015-7-14 23:39:21 | 显示全部楼层
楼主解决问题了吗
回帖提示: 反政府言论将被立即封锁ID 在按“提交”前,请自问一下:我这样表达会给举报吗,会给自己惹麻烦吗? 另外:尽量不要使用Mark、顶等没有意义的回复。不得大量使用大字体和彩色字。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

手机版|Archiver|amobbs.com 阿莫电子论坛 ( 公安交互式论坛备案:44190002001997 粤ICP备09047143号 )

GMT+8, 2022-1-18 12:24

© Since 2004 www.amobbs.com, 原www.ourdev.cn, 原www.ouravr.com

快速回复 返回顶部 返回列表