全国电子设计大赛一等奖作品中FFT算法的源码

georg

全国电子设计大赛一等奖作品中FFT算法的源码 [复制链接]

第八届全国大学生电子设计大赛终于落下帷幕了.经过半年的努力,总算有了回报,拿了个全国一等奖.经过这段时间的历练,懂得了很多事和某些人,我觉得这是远远超过技术本身的收获. 我们做的是A题音频信号分析仪.主要核心是快速傅立叶算法(FFT)的实现.以下是我们写的实现1,2,4,8,16,32,64,128,256,512,1024,2048点的浮点FFT计算代码,利用查表的方法： ///////////////////////////////////////////////////////////////////////////////////////////////////////////// //Writed by levension @SCUT // //Date:2007.09.04 // ////////////////////////////////////////////////////////////////////////////////////////////////////////////// #include"table.h" #include intcount[]={1,2,4,8,16,32,64,128,256,512,1024,2048}; /****************************************************** **基2碟形算法 **dataR--采样点的实部的指针 **dataI--采样点的虚部的指针 **n--点数的幂指数 **0,1,2,3,4,5,6,7,8,9,10,11分别对应者点数 **1,2,4,8,16,32,64,128,256,512,1024,2048 ******************************************************/ void FFT(float *dataR,float *dataI,int n) { int i,L,j,k,b,p,xx,qq; int x[11]={0}; float TR,TI,temp; /********** following code invert sequence************/ for(i=0;i; } for(i=0;i=dataI; dataI=0; } /************** following code FFT*******************/ for(L=1;L<=n;L++) { /* for(1) */ b=1; i=L-1; while(i>0) { b=b*2; i--; } /* b= 2^(L-1) */ for(j=0;j<=b-1;j++) /* for (2)*/ { p=1; i=n-L; while(i>0) /* p=pow(2,7-L)*j; */ { p=p*2; i--;} p=p*j; for(k=j;k查表///////////////////////////////////// dataR[k]=dataR[k]+dataR[k+b]*cos_tab[qq]+dataI[k+b]*sin_tab[qq]; dataI[k]=dataI[k]-dataR[k+b]*sin_tab[qq]+dataI[k+b]*cos_tab[qq]; dataR[k+b]=TR-dataR[k+b]*cos_tab[qq]-dataI[k+b]*sin_tab[qq]; dataI[k+b]=TI+temp*sin_tab[qq]-dataI[k+b]*cos_tab[qq]; //////////////////////////////////////////////////////////////// } /* END for (3) */ } /* END for (2) */ } for(i=0;i=sqrt(dataR*dataR+dataI*dataI); dataR=dataR/count[n-1]; } dataR[0]=dataR[0]/2; } 下面附上TABLE.C文件

 

huangliang112

hao xia kan kan

 

ganice

kankan yixia

 

zou_tian

楼主 怎么没有附件

 

xinzai

收录到资源整理贴

 

yfshangbs

看看 学习

zcy8311

楼主 附件呢？？？

renchunyan

好难啊 看不懂

huanyi414

用什么平台做的啊

rong_dianzi

就算看懂了也看花眼了

676748500

好东西 好东西 我喜欢 我爱楼主

cugtan

晕了，没有附件！

guoxiongmao

好乱啊。。。

guoxiongmao

好乱！好乱！

小志

好了，呵呵

duxiao

thank you very much

ccssa0k

楼主，下来看看，学习学习，顶！！！！！！！！！！！！！！

xiaofeng09

很好呀 感谢能共享呀。。。。。

xujm2008

感谢能共享。。。。。
不过得慢慢看。。。。。

shuining

怎么没有附件呀？