很感谢@STMCU  把论坛bug修复了，自从换了手机号，好几年的时间都不能参加兑换了，一直是手机号格式不合法，虽然以前也反映过多次，不过bug一直存在。这次不但可以参加了，还兑换到了一块G4开发板，原本是想兑换NUCLEO_STM32F446的
G4系列还是第一次接触，看官网介绍，G4系列，处于主流行产品定位。模拟外设资源丰富，更像F3的升级产品。nucleo-G4板最直观的变化就是板载的仿真器从V2版升级到V3版了，之前主控是STM32F103C8，现在居然是高大上的F7了





因为是第一次使用，电脑还需要安装V3版的驱动


不论是keil还是cubeMX,在安装的PCK里都有DSP的运算库，可以直接调用





直接添加相关的文件到工程中去，因为这里是要测试arm_cfft_f32与arm_rfft_fast_f32两个函数，因此添加如下文件

其中，在lib目录下，有四个文件

其后缀分别是Big endian，Floating Point Unit，Little endian，大端还是小端存储模式，是否支持浮点，因此加入了arm_cortexM4lf_math.lib。

计算FFT要涉及大量的浮点运算，所以首先要打开浮点运算单元，实际在现在的keil版本中，默认都是打开的


省去了以前的开启宏定义：__FPU_USED。
要计算FFT，还有定义几个另外的宏
ARM_MATH_CM4,__CC_ARM,ARM_MATH_MATRIX_CHECK,ARM_MATH_ROUNDING,__FPU_PRESENT



最后包含头文件

1. #include "arm_math.h"
   
2. #include "arm_const_structs.h"
   
3. #include "stdio.h"
   
4. #include "stdbool.h"
   
5. #include <stdlib.h>

复制代码

这里省去了FFT的验证内容，就是生成一个合成波数据，对标准数据进行FFT,验证代码的准确性。验证无误后，进行主要的目的，时间测算。ADC采集数据，进行1024点FFT计算

1. HAL_GPIO_WritePin(Ceshi_GPIO_Port, Ceshi_Pin, GPIO_PIN_RESET);
   
2. HAL_Delay(50);
   
3. HAL_GPIO_WritePin(Ceshi_GPIO_Port, Ceshi_Pin, GPIO_PIN_SET);
   
4. HAL_Delay(50);
   
5. HAL_GPIO_WritePin(Ceshi_GPIO_Port, Ceshi_Pin, GPIO_PIN_RESET);
   
6. HAL_TIM_Base_Start_IT(&htim3);
   
7. HAL_ADC_Start_DMA(&hadc1, (uint32_t*)&ADC_Value_Buff, SIZE);
   
8. HAL_TIM_Base_Start(&htim8);
   
9. 
   
10. while(Buf_OK!=true);
    
11. 
    
12. 
    
13. HAL_GPIO_WritePin(Ceshi_GPIO_Port, Ceshi_Pin, GPIO_PIN_SET);
    
14. 
    
15. 
    
16. Buf_OK=false;
    
17. 
    
18. arm_rfft_fast_init_f32(&sfft, FFT_LENGTH);
    
19. for(i=0;i<FFT_LENGTH;i++)
    
20. {
    
21. fft_inputbuf[i]=(float)ADC_Value_Buff[j]*33/40960;
    
22. j=j+2;
    
23. }
    
24. 
    
25. HAL_GPIO_WritePin(Ceshi_GPIO_Port, Ceshi_Pin, GPIO_PIN_RESET);
    
26. 
    
27. arm_rfft_fast_f32(&sfft,fft_inputbuf,fft_outputbuf,0);
    
28. 
    
29. HAL_GPIO_WritePin(Ceshi_GPIO_Port, Ceshi_Pin, GPIO_PIN_SET);
    
30. 
    
31. arm_cmplx_mag_f32(fft_outputbuf,fft_outputbuf2,FFT_LENGTH/2);
    
32. 
    
33. 
    
34. HAL_GPIO_WritePin(Ceshi_GPIO_Port, Ceshi_Pin, GPIO_PIN_RESET);
    
35. 
    
36. 
    
37. for(int i=0;i<FFT_LENGTH/2;i++)
    
38. {
    
39. if(i!=0)
    
40. {
    
41. printf("%d: %.2f\r\n",i,fft_outputbuf2[i]*2/1024);
    
42. }
    
43. else
    
44. {
    
45. printf("%d: %.2f\r\n",i,fft_outputbuf2[i]/1024);
    
46. }        
    
47. }
    
48. 
    
49. HAL_GPIO_WritePin(Ceshi_GPIO_Port, Ceshi_Pin, GPIO_PIN_SET);
    
50. 
    
51. 
    
52. for(int i=0;i<FFT_LENGTH;i++)
    
53. {
    
54. fft_inputbuf[2*i]=(float)ADC_Value_Buff[j]*33/40960;
    
55. 
    
56. fft_inputbuf[2*i+1]=0;
    
57. j+=2;
    
58. }
    
59. 
    
60. HAL_GPIO_WritePin(Ceshi_GPIO_Port, Ceshi_Pin, GPIO_PIN_RESET);
    
61. 
    
62. 
    
63. arm_cfft_f32(&arm_cfft_sR_f32_len1024,fft_inputbuf,0,1);
    
64. 
    
65. HAL_GPIO_WritePin(Ceshi_GPIO_Port, Ceshi_Pin, GPIO_PIN_SET);
    
66. 
    
67. 
    
68. 
    
69. arm_cmplx_mag_f32(fft_inputbuf,fft_outputbuf,FFT_LENGTH);
    
70. 
    
71. 
    
72. HAL_GPIO_WritePin(Ceshi_GPIO_Port, Ceshi_Pin, GPIO_PIN_RESET);

复制代码

arm_cfft_f32与arm_rfft_fast_f32全部用来计算实数，虚数部分程序里做赋0处理。


50ms的启动信号：


ADC采集时间：


数组赋值：


arm_rfft_fast_f32的计算时间，mag_32时间：


printf输出：



数组赋值：


arm_cfft_f32计算时间：


mag_32计算时间：


缩小验证：

示波器测量：






用了好几种方法来测试时间，相互验证，printf函数、RTC函数、定时器、示波器、逻辑分析仪，最后结果统一，说明没问题。可以看出arm_rfft_fast_f32确实要比arm_cfft_f32快很多，对于1024点的计算，一个是434us，一个是710us。