• 中英航母大战谁能赢歼151枚炸弹就能炸瘫女王 2019-03-22
  • 独特“光捕捉器”可控制分子生化性质 2019-03-22
  • 女性之声——全国妇联 2019-03-19
  • 【理上网来·辉煌十九大】德国北威州经济部长:十九大报告表明中国有能力为世界做出更多贡献 2019-03-11
  • 全国小麦收获进度过八成 2019-03-07
  • 【人事】中共临汾市委组织部公示3名拟任职干部 2019-03-07
  • 这泼猴,物理。古文,哲学,基础语文,那一项没出丑?看着就想笑 2019-02-20
  • 云南理发店老板涉嫌杀害女演员因办卡纠纷起杀心 2019-02-20
  • 1至5月我区为企业和社会减负4200余万元 2019-02-09
  • 【学习时刻】参会专家盛小云:坚定文化自信,用高尚的作品引领社会风尚 2019-01-28
  • 绘就生态文明新画卷—甘南州造林绿化工作纪实 2019-01-28
  • 某些人总是用“客观事实”来为自己壮胆,实际上根本不懂什么是“客观事实。 2019-01-12
  • 圣洁的雪山文章中国国家地理网 2018-11-22
  • 崇尚英雄 精忠报国——我们家的报国故事 2018-11-22
  • 618史上最壕“买家”现身 Google以 5.5亿美元投资京东 2018-11-21
  • 打印

    时时彩赚钱技巧: [DSP编程] 5509A的FFT例程的结果有点奇怪

    [复制链接]
    330|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. }

    复制代码




    本帖子中包含更多资源

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

    x
    沙发
    | 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 应该也可以吧, 你可以该一下试试
    扫描二维码,随时随地手机跟帖
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

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

    快速回复

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

    论坛热帖

    关闭

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

    快速回复 广东快乐十分走势图 返回列表
  • 中英航母大战谁能赢歼151枚炸弹就能炸瘫女王 2019-03-22
  • 独特“光捕捉器”可控制分子生化性质 2019-03-22
  • 女性之声——全国妇联 2019-03-19
  • 【理上网来·辉煌十九大】德国北威州经济部长:十九大报告表明中国有能力为世界做出更多贡献 2019-03-11
  • 全国小麦收获进度过八成 2019-03-07
  • 【人事】中共临汾市委组织部公示3名拟任职干部 2019-03-07
  • 这泼猴,物理。古文,哲学,基础语文,那一项没出丑?看着就想笑 2019-02-20
  • 云南理发店老板涉嫌杀害女演员因办卡纠纷起杀心 2019-02-20
  • 1至5月我区为企业和社会减负4200余万元 2019-02-09
  • 【学习时刻】参会专家盛小云:坚定文化自信,用高尚的作品引领社会风尚 2019-01-28
  • 绘就生态文明新画卷—甘南州造林绿化工作纪实 2019-01-28
  • 某些人总是用“客观事实”来为自己壮胆,实际上根本不懂什么是“客观事实。 2019-01-12
  • 圣洁的雪山文章中国国家地理网 2018-11-22
  • 崇尚英雄 精忠报国——我们家的报国故事 2018-11-22
  • 618史上最壕“买家”现身 Google以 5.5亿美元投资京东 2018-11-21