关于GD32F450快速傅里叶变换的应用

风吹过后

一、快速傅里叶变换
快速傅里叶变换 (fast Fourier transform), 即利用计算机计算离散傅里叶变换（DFT)的高效、快速计算方法的统称，简称FFT。快速傅里叶变换是1965年由J.W.库利和T.W.图基提出的。采用这种算法能使计算机计算离散傅里叶变换所需要的乘法次数大为减少，特别是被变换的抽样点数N越多，FFT算法计算量的节省就越显著。
4 g# A5 d7 V0 G1 J( P6 x例程以
4 x* A' ?2 U5 Q+ ky=1.6+sin(100πt)+0.3sin(300πt)+0.5sin(700πt)
为测试序列，采样率为6400Hz，基波频率为工频50Hz，采样点数为128点
二、CMSIS-DSP

CMSIS-DSP库是ARM针对各种Arm Cortex-M处理器（例如Cortex-M4，Cortex-M7，Cortex-M33，Cortex-M35和Cortex-M55处理器）优化的丰富的DSP函数的集合。Arm Developer网站包含有关这些处理器的更多信息和支持资源。

CMSIS-DSP在行业中得到了广泛使用，并能够通过各种第三方工具优化C代码生成。

例程所使用的是dsp开发包中的TransformFunctions实数快速傅里叶变换。

三、关键代码
2 }6 I( N3 @( F! Y

: J2 J( r7 m' C8 Q" o; ]

• #include "RFft_test.h"
• arm_rfft_instance_f32 S;
• arm_cfft_radix4_instance_f32 S_CFFT;
• static float32_t pSrc[128]={0.0};
• static float32_t pDst[256]={0.0};
• static float32_t out[64]={0};
• void rfft_test(void)
• {
•         u8 i=0,j=0;
•         for(i=0;i<128;i++)
•         {
•                 pSrc=1.6f+arm_sin_f32(100*PI*i/6400)+0.3f*arm_sin_f32(300*PI*i/6400)+0.5f*arm_sin_f32(700*PI*i/6400);
•         }
•         arm_rfft_init_f32(&S,&S_CFFT,128,0,1);
•         arm_rfft_f32(&S,pSrc,pDst);
•         for(j=0;j<63;j++)
•         {
•                 arm_sqrt_f32(pDst[2*j]*pDst[2*j]+pDst[2*j+1]*pDst[2*j+1],&out[j]);
•                 if(j==0)
•                 {
•                         out[j]=out[j]/128;//直流分量需要特殊处理
•                 }
•                 else
•                 {
•                         out[j]=out[j]/64;//交流分量
•                 }
•         }
•         for(j=0;j<64;j++)
•         {
•                 printf("%f
• ",out[j]);
•         }
• }
• 
  

- S* L' O: N  P- {8 r- h
3 c" Z6 ?: w5 V+ v8 H
由于arm_rfft_f32函数的输出结果分为实部和虚部，在arm_rfft_init_f32函数采用位反转时，输出结果数组为：
7 O$ m, k/ Q/ M) B: J直流实部、直流虚部、基波实部、基波虚部、2次谐波实部、2次谐波虚部...........
arm_rfft_init_f32函数不采用位反转时直接打印输出结果如下：

3 y9 K1 a# L1 v& ^2 ^* Darm_rfft_init_f32函数采用位反转时直接打印输出结果如下：


) w7 c% i7 ~8 {1 g3 f
由于快速傅里叶变换的真实幅值需要特殊处理，需要将各次分量的实部和虚部进行求模运算，即
2 G/ a. t  t9 @2 t! I/ c# G采用arm_sqrt_f32函数实现对实部和虚部的平方和开跟，即可获得真实的幅值结果如图（未加直流分量）

% c5 H! J& O- g) U$ \可见与我们设定的测试序列y=sin(100πt)+0.3sin(300πt)+0.5sin(700πt)
8 L8 V7 h0 ^; s$ F/ u$ R吻合，基波幅值为1，3次谐波分量幅值为0.3，7次谐波分量幅值为0.5

& y5 h/ L+ S- e0 y2 N为完整测试，加入1.6的直流分量，即测试序列为：
y=1.6+sin(100πt)+0.3sin(300πt)+0.5sin(700πt)
测试结果为：
( _* P& X1 W: n# j' B8 k
可以看到直流分量变为了真实的两倍，原因是傅里叶变换幅值还原需要除上采样点数的1半，
2 S& r; }" L! V" K) x变换结果分布在正频和负频，负频在物理上无意义，但其确实含有能量，物理上能量一般与
/ ?; \1 `) e( N/ @' D幅值对应，因此傅里叶变换幅值还原需要除上采样点数的1半，而直流分量既不属于正频也
! X9 X0 i/ V' `5 c不属于负频，也不区分正直流和负直流，因此直流分量幅值还原需除上采样点数，
- p3 v9 h* U0 ^3 @( l2 j处理结果如图：
) A1 N# `! G7 d! @5 ?& @, X3 J
至此完整的实数傅里叶变换完成
2 s, _7 o% d$ j: M1 A) F( P! O1 {四、结果
采用gd32F450ZIT6和CMSIS-DSP开发包对实数傅里叶变换进行测试，并与matlab测试结果进行对比
实际结果如下：
/ d( Q$ |3 X/ L; Q# w! F" p



  z" o# a0 A) n. X3 v3 J
. _$ X: w/ s; ~1 `- X: Y% t
: w8 Q3 Y& V) L: y  ~

" {/ g. S0 a/ |6 r

CRAZY_argentina

可以可以，看一看，学习学习。

Ele_insect

看一看，学习学习。

big_gun

很细致，学习到了，感谢楼主

Terran

很细致，学习到了，感谢楼主