用于ARM上的FFT与IFFT源代码（C语言，不依赖特定平台）

发布时间：2020-12-14 05:27:24 所属栏目：百科来源：网络整理

导读：http://blog.csdn.net/syrchina/article/details/6670517 代码在2011年全国电子大赛结束后（2011年9月3日）发布，多个版本，注释详细。 [cpp] ? view plain copy /******************************************************************************* ? **?程

http://blog.csdn.net/syrchina/article/details/6670517

代码在2011年全国电子大赛结束后（2011年9月3日）发布，多个版本，注释详细。

[cpp]? 
   view plain 
   copy 
    
   
 /*******************************************************************************?
 **?程序名称：快速傅里叶变换(FFT)??
**?程序描述：本程序实现快速傅里叶变换??
**?程序作者：宋元瑞??
**?最后修改：2011年4月5日??
*******************************************************************************/??
 #include?<stdio.h>??
 #include?<math.h>??
 ??
 #define?PI?3.141592653589???//圆周率，12位小数???
 #define?N?8?????????????????//傅里叶变换的点数???
 #define?M?3?????????????????//蝶形运算的级数，N?=?2^M???
 typedef?double?ElemType;????//原始数据序列的数据类型,可以在这里设置??
 ??
typedef?struct??????????????//定义复数结构体???
 {??
 ????ElemType?real,imag;??
 }complex;??
 complex?data[N];????????????//定义存储单元，原始数据与负数结果均使用之???
 ElemType?result[N];?????????//存储FFT后复数结果的模??
//变址???
 void?ChangeSeat(complex?*DataInput)??
 {??
 ????int?nextValue,nextM,i,k,j=0;??
 ????complex?temp;??
 ??????
 ????nextValue=N/2;??????????????????//变址运算，即把自然顺序变成倒位序，采用雷德算法??
 ????nextM=N-1;??
 ????for?(i=0;i<nextM;i++)??
 ????{??
 ????????if?(i<j)?????????????????//如果i<j,即进行变址??
 ????????{??
 ????????????temp=DataInput[j];??
 ????????????DataInput[j]=DataInput[i];??
 ????????????DataInput[i]=temp;??
 ????????}??
 ????????k=nextValue;????????????????//求j的下一个倒位序??
 ????????while?(k<=j)?????????????//如果k<=j,表示j的最高位为1??
 ????????{??
 ????????????j=j-k;??????????????????//把最高位变成0??
 ????????????k=k/2;??????????????????//k/2，比较次高位，依次类推，逐个比较，直到某个位为0??
 ????????j=j+k;??????????????????????//把0改为1??
 ????}?????????????????????????????????????????
 }??
/*?
//变址??
void?ChangeSeat(complex?*DataInput)?
{?
????complex?Temp[N];?
????int?i,n,New_seat;?
????for(i=0;?i<N;?i++)?
????{?
????????Temp[i].real?=?DataInput[i].real;?
????????Temp[i].imag?=?DataInput[i].imag;?
????}?
????for(i=0;?i<N;?i++)?
????{?
????????New_seat?=?0;?
????????for(n=0;n<M;n++)?
????????{?
????????????New_seat?=?New_seat?|?(((i>>n)?&?0x01)?<<?(M-n-1));?
????????}?
????????DataInput[New_seat].real?=?Temp[i].real;?
????????DataInput[New_seat].imag?=?Temp[i].imag;?
}?
*/??
//复数乘法???
 complex?XX_complex(complex?a,?complex?b)??
 ????complex?temp;??
 ??????
 ????temp.real?=?a.real?*?b.real-a.imag*b.imag;??
 ????temp.imag?=?b.imag*a.real?+?a.imag*b.real;??
return?temp;??
 }??
//FFT??
void?FFT(void)??
 ????int?L=0,B=0,J=0,K=0;??
int?step=0;??
 ????ElemType?P=0,T=0;??
 ????complex?W,Temp_XX;??
 ????//ElemType?TempResult[N];??
 ????ChangeSeat(data);??
 ????for(L=1;?L<=M;?L++)??
 ????{??
 ????????B?=?1<<(L-1);//B=2^(L-1)??
for(J=0;?J<=B-1;?J++)??
 ????????????P?=?(1<<(M-L))*J;//P=2^(M-L)?*J???
 ????????????step?=?1<<L;//2^L??
 ????????????for(K=J;?K<=N-1;?K=K+step)??
 ????????????{??
 ????????????????W.real?=??cos(2*PI*P/N);??
 ????????????????W.imag?=?-sin(2*PI*P/N);??
 ??????????????????
 ????????????????Temp_XX?=?XX_complex(data[K+B],W);??
 ????????????????data[K+B].real?=?data[K].real?-?Temp_XX.real;??
 ????????????????data[K+B].imag?=?data[K].imag?-?Temp_XX.imag;??
 ????????????????data[K].real?=?data[K].real?+?Temp_XX.real;??
 ????????????????data[K].imag?=?data[K].imag?+?Temp_XX.imag;??
 ????????????}??
 ????????}??
 ????}??
void?IFFT(void)??
int?step=0;??
 ????ElemType?P=0,T=0;??
 ????complex?W,Temp_XX;??
 ????//ElemType?TempResult[N];??
 ????ChangeSeat(data);??
for(L=1;?L<=M;?L++)??
 ????????B?=?1<<(L-1);//B=2^(L-1)??
for(J=0;?J<=B-1;?J++)??
 ????????????P?=?(1<<(M-L))*J;//P=2^(M-L)?*J???
 ????????????step?=?1<<L;//2^L??
 ????????????for(K=J;?K<=N-1;?K=K+step)??
 ????????????{??
 ????????????????W.real?=??cos(2*PI*P/N);??
 ????????????????W.imag?=??sin(2*PI*P/N);//逆运算，这里跟FFT符号相反???
 ??????????????????
 ????????????????Temp_XX?=?XX_complex(data[K+B],W);??
 ????????????????data[K+B].real?=?data[K].real?-?Temp_XX.real;??
 ????????????????data[K+B].imag?=?data[K].imag?-?Temp_XX.imag;??
 ????????????????data[K].real?=?data[K].real?+?Temp_XX.real;??
 ????????????????data[K].imag?=?data[K].imag?+?Temp_XX.imag;??
 ????????????}??
 ????}??
 int?main(int?argc,?char?*argv[])??
int?i?=?0;??
for(i=0;?i<N;?i++)//制造输入序列???
 ????????data[i].real?=?sin(2*PI*i/N);??
 ????????printf("%lf?",data[i]);??
 ????printf("nn");??
 ????FFT();//进行FFT计算???
for(i=0;?i<N;?i++)??
 ????????{printf("%lf?",sqrt(data[i].real*data[i].real+data[i].imag*data[i].imag));}??
 ????IFFT();//进行FFT计算???
 ????printf("nn");??
for(i=0;?i<N;?i++)??
 ????????{printf("%lf?",data[i].real/N);}??
 ????printf("n");??
/*for(i=0;?i<N;?i++)?
????????{printf("%lf?",data[i].imag/N);}?
????printf("n");*/??
/*for(i=0;?i<N;?i++)?
??
return?0;??
 }??

copy 
    
   

**?程序版本：V6.5?
**?最后修改：2011年5月16日??
#define?PI??3.14159265358979323846264338327950288419716939937510????//圆周率，50位小数???
#define?PI2?6.28318530717958647692528676655900576839433879875021??
#define?N?1024??????????????//傅里叶变换的点数???
#define?M?10????????????????//蝶形运算的级数，N?=?2^M???
#define?Npart2?512??????????//创建正弦函数表时取PI的1/2???
#define?Npart4?256??????????//创建正弦函数表时取PI的1/4???
#define?FFT_RESULT(x)???(sqrt(data[x].real*data[x].real+data[x].imag*data[x].imag))??
#define?IFFT_RESULT(x)??(data[x].real/N)??
float?ElemType;???????
//定义复数结构体???
????ElemType?real,imag;??
}complex;??
complex?data[N];????????????//定义存储单元，原始数据与负数结果均使用之???
ElemType?SIN_TABLE[Npart4+1];??
//产生模拟原始数据输入???
void?InputData(int?i;??
????????data[i].real?=?sin(2*PI*i/N);???//正弦波???
????????data[i].imag?=?0;??
//创建正弦函数表???
void?CREATE_SIN_TABLE(int?i=0;???
for(i=0;?i<=Npart4;?i++)??
????????SIN_TABLE[i]?=?sin(PI*i/Npart2);//SIN_TABLE[i]?=?sin(PI2*i/N);??
ElemType?Sin_find(ElemType?x)??
int?i?=?(int)(N*x);??
????i?=?i>>1;??
if(i>Npart4)//注意：i已经转化为0~N之间的整数了！???
????{//不会超过N/2???
????????i?=?Npart2?-?i;//i?=?i?-?2*(i-Npart4);??
????}???
return?SIN_TABLE[i];??
ElemType?Cos_find(ElemType?x)??
if(i<Npart4) ????????//i?=?Npart4?-?i;??
return?SIN_TABLE[Npart4?-?i];??
????}???
else//i>Npart4?&&?i<N/2??
????????//i?=?i?-?Npart4;??
return?-SIN_TABLE[i?-?Npart4];??
//变址???
void?ChangeSeat(complex?*DataInput)??
????nextValue=N/2;??????????????????//变址运算，即把自然顺序变成倒位序，采用雷德算法??
????nextM=N-1;??
for?(i=0;i<nextM;i++)??
??
????????????temp=DataInput[j];??
????????????DataInput[j]=DataInput[i];??
????????????DataInput[i]=temp;??
????????k=nextValue;????????????????//求j的下一个倒位序??
??
????????????j=j-k;??????????????????//把最高位变成0??
????????????k=k/2;??????????????????//k/2，比较次高位，依次类推，逐个比较，直到某个位为0??
????????j=j+k;??????????????????????//把0改为1??
????}?????????????????????????????????????????
}?????????????????????????????????????????????
//复数乘法???
/*complex?XX_complex(complex?a,?complex?b)?
{?
????complex?temp;?
?????
????temp.real?=?a.real?*?b.real-a.imag*b.imag;?
????temp.imag?=?b.imag*a.real?+?a.imag*b.real;?
?????
????return?temp;?
}*/??
//FFT运算函数???
int?step=0,?KB=0;??
//ElemType?P=0;??
????ElemType?angle;??
????ChangeSeat(data);//CREATE_SIN_TABLE();??
????????step?=?1<<L; ????????B?=?step>>1;for(J=0;?J<B;?J++)??
????????????//P?=?(1<<(M-L))*J;//P=2^(M-L)?*J???
????????????angle?=?(double)J/B;????????????//这里还可以优化???
????????????W.imag?=??-Sin_find(angle);?????//用C++该函数课声明为inline???
????????????W.real?=???Cos_find(angle);?????//用C++该函数课声明为inline???
//W.real?=??cos(angle*PI);??
????????????//W.imag?=?-sin(angle*PI);??
for(K=J;?K<N;?K=K+step)??
????????????????KB?=?K?+?B;??
????????????????//Temp_XX?=?XX_complex(data[KB],W);??
????????????????//用下面两行直接计算复数乘法，省去函数调用开销???
????????????????Temp_XX.real?=?data[KB].real?*?W.real-data[KB].imag*W.imag;??
????????????????Temp_XX.imag?=?W.imag*data[KB].real?+?data[KB].imag*W.real;??
????????????????data[KB].real?=?data[K].real?-?Temp_XX.real;??
????????????????data[KB].imag?=?data[K].imag?-?Temp_XX.imag;??
//IFFT运算函数???
??????????????
????????????W.imag?=???Sin_find(angle);????? ????????????W.real?=???Cos_find(angle);?????//W.real?=??cos(angle*PI);??
//W.imag?=?-sin(angle*PI);??
for(K=J;?K<N;?K=K+step)??
????????????????KB?=?K?+?B;??
??
//用下面两行直接计算复数乘法，省去函数调用开销???
????????????????Temp_XX.real?=?data[KB].real?*?W.real-data[KB].imag*W.imag;??
????????????????Temp_XX.imag?=?W.imag*data[KB].real?+?data[KB].imag*W.real;??
????????????????data[KB].real?=?data[K].real?-?Temp_XX.real;??
????????????????data[KB].imag?=?data[K].imag?-?Temp_XX.imag;??
//主函数???
????CREATE_SIN_TABLE();?//创建正弦函数表?,这句只需在程序开始时执行一次???
????InputData();????????//输入原始数据?，此处用公式模拟；实际应用时为AD采样数据???
????FFT();??????????????//进行?FFT计算???
????????{printf("%f?",FFT_RESULT(i));}/**/??
//进行?IFFT计算???
copy
1. **?程序描述：本程序实现快速傅里叶变换及其逆变换?
2. **?性能提升：修正了FFT的bug,变量重新命名,并将?N_FFT改为动态值?
3. **?程序版本：V6.6?
4. **?程序作者：宋元瑞?
5. **?最后修改：2011年5月16日?
6. #include?<stdlib.h>??
7. #include?<math.h>??
8. //#define?OUTPRINT?printf??
9. //#define?DEL?/##/??
10. #define?RESULT(x)?sqrt(data_of_N_FFT[x].real*data_of_N_FFT[x].real+data_of_N_FFT[x].imag*data_of_N_FFT[x].imag)??
11. #define?PI??3.14159265358979323846264338327950288419716939937510????//圆周率，50位小数??
12. #define?PI2?6.28318530717958647692528676655900576839433879875021??
13. int?N_FFT=0;????????????????//傅里叶变换的点数??
14. int?M_of_N_FFT=0;???????????//蝶形运算的级数，N?=?2^M??
15. int?Npart2_of_N_FFT=0;??????//创建正弦函数表时取PI的1/2??
16. int?Npart4_of_N_FFT=0;??????//创建正弦函数表时取PI的1/4??
17. //定义复数结构体??
18. }complex_of_N_FFT,*ptr_complex_of_N_FFT;??
19. ptr_complex_of_N_FFT?data_of_N_FFT=NULL;//开辟存储单元，原始数据与负数结果均使用之??
20. ElemType*?SIN_TABLE_of_N_FFT=NULL;??
21. //产生模拟原始数据输入??
22. for?(i=0;?i<N_FFT;?i++)//制造输入序列??
23. ????????data_of_N_FFT[i].real?=?sin(2*PI*i/N_FFT);??//正弦波??
24. ????????data_of_N_FFT[i].imag?=?0;??
25. ????????printf("%f?",data_of_N_FFT[i].real);??
26. //创建正弦函数表??
27. int?i=0;??
28. for?(i=0;?i<=Npart4_of_N_FFT;?i++)??
29. ????????SIN_TABLE_of_N_FFT[i]?=?sin(PI*i/Npart2_of_N_FFT); ElemType?Sin_find(ElemType?x)??
30. int)(N_FFT*x);??
31. ????i?=?i>>1;??
32. if?(i>Npart4_of_N_FFT)//注意：i已经转化为0~N之间的整数了！??
33. //不会超过N/2??
34. ????????i?=?Npart2_of_N_FFT?-?i;return?SIN_TABLE_of_N_FFT[i];??
35. ElemType?Cos_find(ElemType?x)??
36. if?(i<Npart4_of_N_FFT)return?SIN_TABLE_of_N_FFT[Npart4_of_N_FFT?-?i];??
37. return?-SIN_TABLE_of_N_FFT[i?-?Npart4_of_N_FFT];??
38. //变址??
39. void?ChangeSeat(complex_of_N_FFT?*DataInput)??
40. ????complex_of_N_FFT?temp;??
41. ????nextValue=N_FFT/2;?????????????????? ????nextM=N_FFT-1;??
42. //复数乘法??
43. ?
44. ?
45. //FFT运算函数??
46. ????complex_of_N_FFT?W,248); line-height:18px"> ????ChangeSeat(data_of_N_FFT);for?(L=1;?L<=M_of_N_FFT;?L++)??
47. for?(J=0;?J<B;?J++)??
48. //P?=?(1<<(M-L))*J;//P=2^(M-L)?*J??
49. //这里还可以优化??
50. //用C++该函数课声明为inline??
51. //用C++该函数课声明为inline??
52. for?(K=J;?K<N_FFT;?K=K+step)??
53. //用下面两行直接计算复数乘法，省去函数调用开销??
54. ????????????????Temp_XX.real?=?data_of_N_FFT[KB].real?*?W.real-data_of_N_FFT[KB].imag*W.imag;??
55. ????????????????Temp_XX.imag?=?W.imag*data_of_N_FFT[KB].real?+?data_of_N_FFT[KB].imag*W.real;??
56. ????????????????data_of_N_FFT[KB].real?=?data_of_N_FFT[K].real?-?Temp_XX.real;??
57. ????????????????data_of_N_FFT[KB].imag?=?data_of_N_FFT[K].imag?-?Temp_XX.imag;??
58. ????????????????data_of_N_FFT[K].real?=?data_of_N_FFT[K].real?+?Temp_XX.real;??
59. ????????????????data_of_N_FFT[K].imag?=?data_of_N_FFT[K].imag?+?Temp_XX.imag;??
60. //IFFT运算函数??
61. //ElemType?P=0;??
62. ????ElemType?angle;??
63. ????complex_of_N_FFT?W,108); list-style:decimal-leading-zero outside; color:inherit; line-height:18px"> ????ChangeSeat(data_of_N_FFT);//变址??
64. //CREATE_SIN_TABLE();??
65. for?(L=1;?L<=M_of_N_FFT;?L++)??
66. ????????step?=?1<<L; ????????B?=?step>>1;for?(J=0;?J<B;?J++)??
67. //P?=?(1<<(M-L))*J;//P=2^(M-L)?*J??
68. ????????????angle?=?(//这里还可以优化??
69. ????????????W.imag?=???Sin_find(angle);?????//初始化FFT程序??
70. //N_FFT是?FFT的点数，必须是2的次方??
71. void?Init_FFT(int?N_of_FFT)??
72. int?i=0;??
73. int?temp_N_FFT=1;??
74. ????N_FFT?=?N_of_FFT;???????????????????//傅里叶变换的点数?，必须是?2的次方??
75. ????M_of_N_FFT?=?0;?????????????????//蝶形运算的级数，N?=?2^M??
76. for?(i=0;?temp_N_FFT<N_FFT;?i++)??
77. ????????temp_N_FFT?=?2*temp_N_FFT;??
78. ????????M_of_N_FFT++;??
79. ????printf("n%dn",M_of_N_FFT);??
80. ????Npart2_of_N_FFT?=?N_FFT/2;??????//创建正弦函数表时取PI的1/2??
81. ????Npart4_of_N_FFT?=?N_FFT/4;??????//创建正弦函数表时取PI的1/4??
82. //ptr_complex_of_N_FFT?data_of_N_FFT=NULL;//开辟存储单元，原始数据与负数结果均使用之??
83. ????data_of_N_FFT?=?(ptr_complex_of_N_FFT)malloc(N_FFT?*?sizeof(complex_of_N_FFT));??
84. //ptr_complex_of_N_FFT?SIN_TABLE_of_N_FFT=NULL;??
85. ????SIN_TABLE_of_N_FFT?=?(ElemType?*)malloc((Npart4_of_N_FFT+1)?*?sizeof(ElemType));??
86. ????CREATE_SIN_TABLE();?????????????//结束?FFT运算，释放相关内存??
87. void?Close_FFT(if?(data_of_N_FFT?!=?NULL)??
88. ????????free(data_of_N_FFT);??????????//释放内存??
89. ????????data_of_N_FFT?=?NULL;??
90. if?(SIN_TABLE_of_N_FFT?!=?NULL)??
91. ????????free(SIN_TABLE_of_N_FFT);?????//释放内存??
92. ????????SIN_TABLE_of_N_FFT?=?NULL;??
93. //主函数??
94. ????Init_FFT(1024);?????//初始化各项参数，此函数只需执行一次??
95. //输入原始数据??
96. //进行?FFT计算??
97. for?(i=0;?i<N_FFT;?i++)??
98. ????????printf("%f?",RESULT(i));??
99. ????IFFT();//进行?IFFT计算??
100. for?(i=0;?i<N_FFT;?i++)??
101. ????Close_FFT();????????
  copy
  1. **?模块描述：本程序实现快速傅里叶变换??
  2. **?性能提升：已达到网上同类程序最高速度??
  3. **?模块版本：V6.0?
  4. **?模块作者：宋元瑞??
  5. **?最后修改：2011年5月6日?
  6. **??程序说明：?
  7. ????FFT运算输入参数?source_Data(i)?是一个N大小的数组（注意是小括号）?
  8. ????FFT运算输出结果?result_Data(i)?是一个N大小的数组（注意是小括号）?
  9. **??调用举例：?
  10. int?main(int?argc,?char?*argv[])?
  11. ????int?i?=?0;?
  12. ????///以下为FFT运算的调用方式??
  13. ????FFT_prepare();??????//为FFT做好准备，此函数只需程序开始时执行一次即可，切勿写在循环中??
  14. ????while(1)?
  15. ????????for(i=0;i<N_FFT;i++)?//输入原始数据??
  16. ????????????{source_Data(i)?=?sin(2*PI*i/N_FFT);}//注意inputData后面是小括号??
  17. ????????FFT();??????????????//进行FFT计算??
  18. ????????//输出结果：XXX?=??result_Data(i);?
  19. ????}?
  20. ????return?0;?
  21. *******************************************************************************/???
  22. #ifndef?SYRFFT_H??
  23. //#pragma?once??
  24. //#include?<stdio.h>??
  25. #define?FFT_prepare()?CREATE_SIN_TABLE();???//创建正弦函数表???
  26. #define?source_Data(i)?data_FFT[i].imag?????//接收输入数据的数组，大小为N???
  27. #define?result_Data(i)?sqrt(data_FFT[i].real*data_FFT[i].real+data_FFT[i].imag*data_FFT[i].imag)//FFT结果???
  28. #define?PI??3.14159265358979323846264338327950288419716939937510????//圆周率，50位小数???
  29. #define?N_FFT?1024??????????????//傅里叶变换的点数???
  30. #define?M_of_N_FFT?10???????????????//蝶形运算的级数，N?=?2^M???
  31. #define?Npart2_of_N_FFT?512?????????//创建正弦函数表时取PI的1/2???
  32. #define?Npart4_of_N_FFT?256?????????//创建正弦函数表时取PI的1/4???
  33. }complex_of_N_FFT;??
  34. complex_of_N_FFT?data_FFT[N_FFT];???????//存放要进行FFT运输的数据，运算结果也存放这里???
  35. ElemType?SIN_TABLE[Npart4_of_N_FFT+1];??
  36. /*?
  37. void?InputData(void)?
  38. ????int?i;?
  39. ????for(i=0;?i<N_FFT;?i++)???????//制造输入序列??
  40. ????????source_Data(i)?=?sin(2*PI*i/N_FFT);?//模拟输入正弦波??
  41. ????????//data_FFT[i].imag?=?sin(2*PI*i/N);?//正弦波??
  42. ????????//printf("%f?",data_FFT[i].imag);?
  43. for(i=0;?i<=Npart4_of_N_FFT;?i++)??
  44. ????????SIN_TABLE[i]?=?sin(PI*i/Npart2_of_N_FFT);int)(N_FFT*x);??
  45. if(i>Npart4_of_N_FFT)if(i<Npart4_of_N_FFT)return?SIN_TABLE[Npart4_of_N_FFT?-?i];??
  46. return?-SIN_TABLE[i?-?Npart4_of_N_FFT];??
  47. //变址前data_FFT[i].imag存储了输入数据，变址后data_FFT[i].real存储了输入数据??
  48. int?i,New_seat;??
  49. for(i=0;?i<N_FFT;?i++)??
  50. ????????New_seat?=?0;??
  51. for(n=0;n<M_of_N_FFT;n++)??
  52. ????????????New_seat?=?New_seat?|?(((i>>n)?&?0x01)?<<?(M_of_N_FFT-n-1));??
  53. ????????DataInput[New_seat].real?=?DataInput[i].imag;??
  54. for(i=0;?i<N_FFT;?i++)??
  55. ????????DataInput[i].imag?=?0;??
  56. }????????????????????????????????????????????????????????
  57. complex_of_N_FFT?XX_complex(complex_of_N_FFT?a,?complex_of_N_FFT?b)?
  58. ????complex_of_N_FFT?temp;?
  59. ????ChangeSeat(data_FFT);???for(L=1;?L<=M_of_N_FFT;?L++)??
  60. for(K=J;?K<N_FFT;?K=K+step)??
  61. //Temp_XX?=?XX_complex(data_FFT[KB],108); list-style:decimal-leading-zero outside; color:inherit; line-height:18px"> ????????????????Temp_XX.real?=?data_FFT[KB].real?*?W.real-data_FFT[KB].imag*W.imag;??
  62. ????????????????Temp_XX.imag?=?W.imag*data_FFT[KB].real?+?data_FFT[KB].imag*W.real;??
  63. ????????????????data_FFT[KB].real?=?data_FFT[K].real?-?Temp_XX.real;??
  64. ????????????????data_FFT[KB].imag?=?data_FFT[K].imag?-?Temp_XX.imag;??
  65. ????????????????data_FFT[K].real?=?data_FFT[K].real?+?Temp_XX.real;??
  66. ????????????????data_FFT[K].imag?=?data_FFT[K].imag?+?Temp_XX.imag;??
  67. #define?SYRFFT_H??
  68. #endif??
  ?
  
  copy
  1. **?程序作者：宋元瑞??
  2. **?最后修改：2011年5月6日?
  3. *******************************************************************************/??
  4. #include?"syrFFT_H.h"???//包含FFT运算头文件???
  5. //以下3句为FFT运算的调用函数???
  6. ????FFT_prepare();??????//为FFT做好准备，此函数只需程序开始时执行一次即可，切勿写在循环中???
  7. //while(1)??
  8. for(i=0;i<N_FFT;i++)//模拟输入???
  9. ????????????{source_Data(i)?=?sin(2*PI*i/N_FFT);}//注意inputData后面是小括号???
  10. ????????FFT();??
  11. for(i=0;?i<N_FFT;?i++)//输出结果???
  12. ????????????{printf("%lf?",result_Data(i));}??
  13. ?
    
    copy
    
    #ifndef?FFT_H??
    #include?<stdlib.h>??
    
    //傅里叶变换的点数???
    
    //蝶形运算的级数，N?=?2^M???
    
    //创建正弦函数表时取PI的1/2???
    
    //创建正弦函数表时取PI的1/4???
    
    //开辟存储单元，原始数据与负数结果均使用之???
    
    //创建正弦函数表???
    
    int?i=0;???
    
    for(i=0;?i<=Npart4_of_N_FFT;?i++)??
    
    ????????SIN_TABLE_of_N_FFT[i]?=?sin(PI*i/Npart2_of_N_FFT);//SIN_TABLE[i]?=?sin(PI2*i/N);??
    
    //注意：i已经转化为0~N之间的整数了！???
    
    ????{//不会超过N/2???
    
    ????????i?=?Npart2_of_N_FFT?-?i;//i?=?i?-?2*(i-Npart4);??
    
    return?SIN_TABLE_of_N_FFT[i];??
    
    //i?=?Npart4?-?i;??
    
    return?SIN_TABLE_of_N_FFT[Npart4_of_N_FFT?-?i];??
    
    //i>Npart4?&&?i<N/2??
    
    //i?=?i?-?Npart4;??
    
    return?-SIN_TABLE_of_N_FFT[i?-?Npart4_of_N_FFT];??
    
    void?ChangeSeat(complex_of_N_FFT?*DataInput)??
    
    ????complex_of_N_FFT?temp;??
    
    ????nextValue=N_FFT/2;?????????????????? ????nextM=N_FFT-1;??
    
    }?????????????????????????????????????????????
    
    ????complex?temp;?
    
    ????temp.real?=?a.real?*?b.real-a.imag*b.imag;?
    
    ????temp.imag?=?b.imag*a.real?+?a.imag*b.real;?
    
    ????return?temp;?
    
    }*/??
    
    //FFT运算函数???
    
    for(L=1;?L<=M_of_N_FFT;?L++)??
    
    for(J=0;?J<B;?J++)??
    
    //P?=?(1<<(M-L))*J;//P=2^(M-L)?*J???
    
    //这里还可以优化???
    
    ????????????W.imag?=??-Sin_find(angle);?????for(K=J;?K<N_FFT;?K=K+step)??
    
    ????????????????Temp_XX.real?=?data_of_N_FFT[KB].real?*?W.real-data_of_N_FFT[KB].imag*W.imag;??
    
    ????????????????Temp_XX.imag?=?W.imag*data_of_N_FFT[KB].real?+?data_of_N_FFT[KB].imag*W.real;??
    
    ????????????????data_of_N_FFT[KB].real?=?data_of_N_FFT[K].real?-?Temp_XX.real;??
    
    ????????????????data_of_N_FFT[KB].imag?=?data_of_N_FFT[K].imag?-?Temp_XX.imag;??
    
    ????????????????data_of_N_FFT[K].real?=?data_of_N_FFT[K].real?+?Temp_XX.real;??
    
    ????????????????data_of_N_FFT[K].imag?=?data_of_N_FFT[K].imag?+?Temp_XX.imag;??
    
    //IFFT运算函数???
    
    ??????????????
    
    //初始化FFT程序???
    
    //N_FFT是?FFT的点数，必须是2的次方???
    
    //傅里叶变换的点数?，必须是?2的次方???
    
    //蝶形运算的级数，N?=?2^M???
    
    for(i=0;?temp_N_FFT<N_FFT;?i++)??
    
    //printf("n%dn",M_of_N_FFT);??
    
    //创建正弦函数表时取PI的1/2???
    
    //创建正弦函数表时取PI的1/4???
    
    //ptr_complex_of_N_FFT?data_of_N_FFT=NULL;//开辟存储单元，原始数据与负数结果均使用之???
    
    //结束?FFT运算，释放相关内存???
    
    if(data_of_N_FFT?!=?NULL)??
    
    ????????free(data_of_N_FFT);?????????? ????????data_of_N_FFT?=?NULL;??
    
    if(SIN_TABLE_of_N_FFT?!=?NULL)??
    
    ????????free(SIN_TABLE_of_N_FFT);????? ????????SIN_TABLE_of_N_FFT?=?NULL;??
    
    #define?FFT_H??
    
    #endif??
    
    copy
    
    **?程序描述：本程序实现快速傅里叶变换及其逆变换??
    **?性能提升：修正了FFT的bug??
    
    #include?"jxust_fft6_6.h"??
    
    //产生模拟原始数据输入?,在实际应用时替换为AD采样数据???
    
    //输入采样数据???
    
    //主函数?，示例如何调用FFT运算???
    
    char?*argv[])??
    
    int?i?=?0;??
    
    ????Init_FFT(1024);?????//①初始化各项参数，此函数只需执行一次?;?参数为FFT的点数，必须为2的次方???
    
    ????InputData();????????//②输入原始数据?,此处在实际应用时替换为AD采样数据???
    
    ????FFT();??????????????//③进行?FFT计算???
    
    //for(i=0;?i<N_FFT;?i++)//④输出FFT频谱结果??sqrt(data_of_N_FFT[i].real*data_of_N_FFT[i].real+data_of_N_FFT[i].imag*data_of_N_FFT[i].imag)??
    
    //{printf("%f?",RESULT(i));}??
    
    //for(i=0;?i<N_FFT;?i++)//(可选步骤)⑤输出?IFFT结果?，滤波时会用到；暂时不用???
    
    ??
    
    ????Close_FFT();????//结束?FFT运算，释放相关内存?;此函数在彻底结束FFT运算后调用，???
    
    return?0;??
    
    }??
    
    copy
    
    #ifndef?SYRFFT_6_55_H??
    #define?FFT_RESULT(x)???(sqrt(data_of_N_FFT[x].real*data_of_N_FFT[x].real+data_of_N_FFT[x].imag*data_of_N_FFT[x].imag))??
    
    #define?IFFT_RESULT(x)??(data_of_N_FFT[x].real/N_FFT)??
    
    #define?N_FFT???????????1024????????//傅里叶变换的点数???
    
    #define?M_of_N_FFT??????10??????????//蝶形运算的级数，N?=?2^M???
    
    #define?Npart2_of_N_FFT?512?????????//创建正弦函数表时取PI的1/2???
    
    #define?Npart4_of_N_FFT?256?????????//创建正弦函数表时取PI的1/4???
    
    }complex_of_N_FFT,*ptr_complex_of_N_FFT;??
    
    complex_of_N_FFT?data_of_N_FFT[N_FFT];?????????? ElemType?SIN_TABLE_of_N_FFT[Npart4_of_N_FFT+1];??
    
    int)(N_FFT*x); ????i?=?i>>1;//?i?=?i?/?2;??
    
    if(i>Npart4_of_N_FFT)??
    
    //根据FFT相关公式，sin()参数不会超过PI，?即i不会超过N/2???
    
    if(i?<?Npart4_of_N_FFT?)???
    
    ????{?else?//i?>?Npart4?&&?i?<?N/2???
    
    //初始化FFT程序???
    
    void?Init_FFT()??
    
    //结束?FFT运算，释放相关内存???
    
    #define?SYRFFT_6_55_H??
    
    #endif??
    
    copy
    
    **?程序版本：V6.55?
    **?最后修改：2011年5月22日??
    
    #include?"syrFFT_6_55.h"???
    
    ??
    
    int?i;??
    
    //制造输入序列???
    
    ????????data_of_N_FFT[i].real?=?sin(2*PI*i/N_FFT);??//输入采样数据???
    
    ????????data_of_N_FFT[i].imag?=?0;??
    
    //主函数?，示例如何调用FFT运算???
    
    ????Init_FFT();?????????//①初始化各项参数，此函数只需执行一次?;?修改FFT的点数去头文件的宏定义处修改???
    
    ??
    
    //③进行?FFT计算???
    
    //for(i=0;?i<N_FFT;?i++)//④输出FFT频谱结果??sqrt(data_of_N_FFT[i].real*data_of_N_FFT[i].real+data_of_N_FFT[i].imag*data_of_N_FFT[i].imag)??
    
    ??
    
    //进行?IFFT计算???
    
    //for(i=0;?i<N_FFT;?i++)//(可选步骤)⑤输出?IFFT结果?，滤波时会用到；暂时不用???
    
    ??
    
    //结束?FFT运算，此版本此句无用，可不写???
    
    }??
    
    （编辑：李大同）
    
    【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容!