打印

MSP430 FFT遇到的问题

[复制链接]
1012|3
手机看帖
扫描二维码
随时随地手机跟帖
跳转到指定楼层
楼主
ikemada|  楼主 | 2014-7-26 11:45 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式
先用 AD600HZ采样100HZ正弦波,得到128个值,没有问题;

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

//我的FFT ,已经倒位序完成。AD量给到Fft_Real【128】中
void FFT(void)
{
int NN=7; //这里因为128是2的7次方,如果是计算256点,则是2的8次方,N就是8,如果是512点则N=9,如此类推


unsigned int NUM_FFT=128; //这里要算多少点的fft就赋值多少,值只能是2的N次方
    for( i=1; i<=NN; i++)             /* for(1) */
    {
        b=1;
        b <<=(i-1);                 //蝶式运算,用于计算 隔多少行计算。例如第一级 1和2行计算,,,第二级
        for( j=0; j<=b-1; j++)       /* for (2) */
        {
            p=1;
            p <<= (NN-i);
            p = p*j;
            for( k2=j; k2<NUM_FFT; k2=k2+2*b)   /* for (3) 基二fft */
            {
                Temp_Real = Fft_Real[k2]; Temp_Imag = Fft_Image[k2]; temp = Fft_Real[k2+b];
                Fft_Real[k2] = Fft_Real[k2] + ((Fft_Real[k2+b]*COS_TAB[p])>>7) + ((Fft_Image[k2+b]*SIN_TAB[p])>>7);
                Fft_Image[k2] = Fft_Image[k2] - ((Fft_Real[k2+b]*SIN_TAB[p])>>7) + ((Fft_Image[k2+b]*COS_TAB[p])>>7);
                Fft_Real[k2+b] = Temp_Real - ((Fft_Real[k2+b]*COS_TAB[p])>>7) - ((Fft_Image[k2+b]*SIN_TAB[p])>>7);
                Fft_Image[k2+b] = Temp_Imag + ((temp*SIN_TAB[p])>>7) - ((Fft_Image[k2+b]*COS_TAB[p])>>7);
//                //移位,防止溢出。结果已经是本值的1/64
              Fft_Real[k2] >>= 1;
                Fft_Image[k2] >>= 1;
               Fft_Real[k2+b]  >>= 1;
                Fft_Image[k2+b]  >>= 1;
            }
        }
    }
   }
还有一个问题,计算模值:
for(n=1;n<128;n++)
        {
                //Fft_Real=conventionsignal()
                if(Fft_Real[n]<0) Fft_Real[n] = -Fft_Real[n];
                if(Fft_Image[n]<0) Fft_Image[n] = -Fft_Image[n];
                result[n]=sqrt((Fft_Real[n]*Fft_Real[n])+(Fft_Image[n]*Fft_Image[n]));
        }
我在FFT之后计算各频率模制,这里要不要都取正?

相关帖子

沙发
ikemada|  楼主 | 2014-7-26 11:46 | 只看该作者
本帖最后由 ikemada 于 2014-7-26 13:00 编辑

这是FFT转换好后实部和虚部的量,不知道对么?

222.jpg (61.93 KB )

实部

实部

333.jpg (60.39 KB )

虚部

虚部

111.jpg (63.75 KB )

adC采样数据

adC采样数据

444.jpg (109.97 KB )

实部^2+虚部^2

实部^2+虚部^2

使用特权

评论回复
板凳
firstblood| | 2014-7-26 17:03 | 只看该作者
最后计算模值的就可以的啊,然后找出最大分量的就可以的啊

使用特权

评论回复
地板
dirtwillfly| | 2014-7-26 20:05 | 只看该作者

使用特权

评论回复
发新帖 我要提问
您需要登录后才可以回帖 登录 | 注册

本版积分规则

2

主题

4

帖子

0

粉丝