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

呱呱

全国电子设计大赛一等奖作品中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 { xx=0;
for(j=0;j x[j]=0;
for(j=0;j {x[j]=(i/count[j])&0x01;}
for(j=0;j {xx=xx+x[j]*count[n-j-1];}
dataI[xx]=dataR;
}
for(i=0;i {
dataR=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 /* for (3) */
{
TR=dataR[k];
TI=dataI[k];
temp=dataR[k+b];
//QQ=2*pi*p/count[n];
qq=p*count[11-n];
//////////////////////////////查表/////////////////////////////////////
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 {
dataR=sqrt(dataR*dataR+dataI*dataI);
dataR=dataR/count[n-1];
}
dataR[0]=dataR[0]/2;
}

随锋飘舞

顶一下啊！！好拽啊

 

呱呱

一般一般,世界第三. 大家尽管放马过来~

 

classv

顶一个

 

zwdeng123

Thanks!

 

qianhng

丁页！

fang8_2

感谢楼主的无私奉献，慢慢学习中。。

iamabirl

Thanks!

sunlight1888

这类资源很稀少,多谢楼主无私分享~

1133807753

学习了！！！！

mao285909621

FFT的算法总算可以见识一下了

j.w

算法类的好像接触得挺少的，先收藏了

twlk8634109

谢谢  分享

xl0512

学习了:'( :'( :'( :'(

tangyvhfhv56

顶起！！！！！！！！！！！11

平常心

一般一般,世界第三. 大家尽管放马过来

lunlong

感谢楼主，有个注释的话就更好了

r20061474

学习中 、、、、

kebiding

今年很可能还会车fft 看看

hqm1921

好啊！！
我就需要这个