• 孙杨、宁泽涛扛起中国游泳大旗 2019-05-11
  • 那你就写个帖子驳呀!告诉你,我是普通老百姓就是错了也没有啥,你可是论坛“真理化身”呀! 2019-05-11
  • 守住青山不放松 护好绿水不辞难——在渝全国人大代表聚焦“共抓大保护、不搞大开发”专题调研记略 2019-05-10
  • 巴厘岛被爆众多海域垃圾成灾 潜水常被塑料瓶塑料袋包围 2019-05-04
  • 两个不能否定,怎会是停留在30年前? 2019-04-27
  • 国际在线:向世界报道中国,向中国报道世界 2019-04-27
  • 【理上网来喜迎十九大】外媒记者:稳步推进生态文明建设的中国更加美丽 2019-04-26
  • 董卿海霞李思思康辉撒贝宁 央视名嘴童年照谁更萌 2019-04-26
  • “龙江二号”微卫星传回地月合影 2019-04-18
  • 中国驻泰国大使考察合艾国光中学孔子课堂 2019-04-18
  • 【清园20年】感恩大回馈!半价!半价! 2019-04-18
  • 中央纪委紧盯节点 4年来通报曝光问题近4000起 2019-04-13
  • 2013环球企业领袖宁夏圆桌会议嘉宾云集 2019-04-13
  • 《魔兽世界》未解之谜 那些仍不能被玩家们踏足的领地 2019-04-11
  • 钱江潮评 让高质量发展插上人才的翅膀 2019-04-10
  • 打印

    11选5稳中两个算法: [DSP编程] 5509A的FFT例程的结果有点奇怪

    [复制链接]
    403|4
    跳转到指定楼层
    楼主
    根据例程,INPUT的波形如下:



    然后根据CCS软件graph看直接查看INPUT的FFT的结果:



    graph的FFT的结果显示频率大致为0.025Hz。这是正确的(根据时域图可知,波形的周期大约40s)。


    然而例程中的结果也就是DATA数组的波形为:



    DATA的频率变成了2.9Hz。


    所以,这是什么意思???难道在程序中最终FFT的结果不是DATA数组么?因为对FFT不太懂程序,望请教。


    附上FFT例程的代码部分:
    1. #include "myapp.h"
    2. #include "csedu.h"
    3. #include "scancode.h"
    4. #include <math.h>

    5. #define PI 3.1415926
    6. #define SAMPLENUMBER 128

    7. void InitForFFT();
    8. void MakeWave();

    9. int INPUT[SAMPLENUMBER],DATA[SAMPLENUMBER];
    10. float fWaveR[SAMPLENUMBER],fWaveI[SAMPLENUMBER],w[SAMPLENUMBER];
    11. float sin_tab[SAMPLENUMBER],cos_tab[SAMPLENUMBER];

    12. main()
    13. {
    14.         int i;
    15.        
    16.         InitForFFT();
    17.         MakeWave();
    18.         for ( i=0;i<SAMPLENUMBER;i++ )
    19.         {
    20.                 fWaveR[i]=INPUT[i];
    21.                 fWaveI[i]=0.0f;
    22.                 w[i]=0.0f;
    23.         }
    24.         FFT(fWaveR,fWaveI);
    25.         for ( i=0;i<SAMPLENUMBER;i++ )
    26.         {
    27.                 DATA[i]=w[i];
    28.         }
    29.         while ( 1 );        // break point
    30. }

    31. void FFT(float dataR[SAMPLENUMBER],float dataI[SAMPLENUMBER])
    32. {
    33.         int x0,x1,x2,x3,x4,x5,x6,xx;
    34.         int i,j,k,b,p,L;
    35.         float TR,TI,temp;
    36.        
    37.         /********** following code invert sequence ************/
    38.         for ( i=0;i<SAMPLENUMBER;i++ )
    39.         {
    40.                 x0=x1=x2=x3=x4=x5=x6=0;
    41.                 x0=i&0x01; x1=(i/2)&0x01; x2=(i/4)&0x01; x3=(i/8)&0x01;x4=(i/16)&0x01; x5=(i/32)&0x01; x6=(i/64)&0x01;
    42.                 xx=x0*64+x1*32+x2*16+x3*8+x4*4+x5*2+x6;
    43.                 dataI[xx]=dataR[i];
    44.         }
    45.         for ( i=0;i<SAMPLENUMBER;i++ )
    46.         {
    47.                 dataR[i]=dataI[i]; dataI[i]=0;
    48.         }

    49.         /************** following code FFT *******************/
    50.         for ( L=1;L<=7;L++ )
    51.         { /* for(1) */
    52.                 b=1; i=L-1;
    53.                 while ( i>0 )
    54.                 {
    55.                         b=b*2; i--;
    56.                 } /* b= 2^(L-1) */
    57.                 for ( j=0;j<=b-1;j++ ) /* for (2) */
    58.                 {
    59.                         p=1; i=7-L;
    60.                         while ( i>0 ) /* p=pow(2,7-L)*j; */
    61.                         {
    62.                                 p=p*2; i--;
    63.                         }
    64.                         p=p*j;
    65.                         for ( k=j;k<128;k=k+2*b ) /* for (3) */
    66.                         {
    67.                                 TR=dataR[k]; TI=dataI[k]; temp=dataR[k+b];
    68.                                 dataR[k]=dataR[k]+dataR[k+b]*cos_tab[p]+dataI[k+b]*sin_tab[p];
    69.                                 dataI[k]=dataI[k]-dataR[k+b]*sin_tab[p]+dataI[k+b]*cos_tab[p];
    70.                                 dataR[k+b]=TR-dataR[k+b]*cos_tab[p]-dataI[k+b]*sin_tab[p];
    71.                                 dataI[k+b]=TI+temp*sin_tab[p]-dataI[k+b]*cos_tab[p];
    72.                         } /* END for (3) */
    73.                 } /* END for (2) */
    74.         } /* END for (1) */
    75.         for ( i=0;i<SAMPLENUMBER/2;i++ )
    76.         {
    77.                 w[i]=sqrt(dataR[i]*dataR[i]+dataI[i]*dataI[i]);
    78.         }
    79. } /* END FFT */


    80. void InitForFFT()
    81. {
    82.         int i;
    83.        
    84.         for ( i=0;i<SAMPLENUMBER;i++ )
    85.         {
    86.                 sin_tab[i]=sin(PI*2*i/SAMPLENUMBER);
    87.                 cos_tab[i]=cos(PI*2*i/SAMPLENUMBER);
    88.         }
    89. }

    90. void MakeWave()
    91. {
    92.         int i;
    93.        
    94.         for ( i=0;i<SAMPLENUMBER;i++ )
    95.         {
    96.                 INPUT[i]=sin(PI*2*i/SAMPLENUMBER*3)*1024;
    97.         }
    98. }

    复制代码




    沙发
    | 2018-11-6 20:54 | 只看该作者
    FFT 的算法不会有错, 会不会采样时间间隔弄错了

    评论

    三关白马 2018-11-7 15:14 回复TA
    请问一下,这个FFT例程只能适用于128个点的变换??? 
    三关白马 2018-11-7 09:46 回复TA
    采样时间间隔,应该不会改变频率吧,我觉得。 
    板凳
    | 2018-11-7 18:54 | 只看该作者
    这里用
    #define SAMPLENUMBER 128

    定义了 128 点, 应该是测试用的这个值. 其他值比如 512 应该也可以吧, 你可以该一下试试
    扫描二维码,随时随地手机跟帖
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    我要发帖 投诉建议 创建版块 申请版主

    快速回复

    您需要登录后才可以回帖
    登录 | 注册
    高级模式

    论坛热帖

    关闭

    热门推荐上一条 /4 下一条

    快速回复 广东快乐十分走势图 返回列表
  • 孙杨、宁泽涛扛起中国游泳大旗 2019-05-11
  • 那你就写个帖子驳呀!告诉你,我是普通老百姓就是错了也没有啥,你可是论坛“真理化身”呀! 2019-05-11
  • 守住青山不放松 护好绿水不辞难——在渝全国人大代表聚焦“共抓大保护、不搞大开发”专题调研记略 2019-05-10
  • 巴厘岛被爆众多海域垃圾成灾 潜水常被塑料瓶塑料袋包围 2019-05-04
  • 两个不能否定,怎会是停留在30年前? 2019-04-27
  • 国际在线:向世界报道中国,向中国报道世界 2019-04-27
  • 【理上网来喜迎十九大】外媒记者:稳步推进生态文明建设的中国更加美丽 2019-04-26
  • 董卿海霞李思思康辉撒贝宁 央视名嘴童年照谁更萌 2019-04-26
  • “龙江二号”微卫星传回地月合影 2019-04-18
  • 中国驻泰国大使考察合艾国光中学孔子课堂 2019-04-18
  • 【清园20年】感恩大回馈!半价!半价! 2019-04-18
  • 中央纪委紧盯节点 4年来通报曝光问题近4000起 2019-04-13
  • 2013环球企业领袖宁夏圆桌会议嘉宾云集 2019-04-13
  • 《魔兽世界》未解之谜 那些仍不能被玩家们踏足的领地 2019-04-11
  • 钱江潮评 让高质量发展插上人才的翅膀 2019-04-10
  • 香港白小姐玄机图 快乐10分走势图 海南飞鱼彩票网上购买 北京赛车平投稳胜法 快乐8走势图360 街机麻将 今日体彩p3试机号 天津时时彩在哪买 快乐飞艇怎么计划 任选9场开奖规则 彩票走势图-重庆欢乐生肖-五星 中国体育彩票排列5 彩票网 极速飞艇内购破解版 2017双色球北京大奖 重庆时时彩的正规网址