stm32-Timer-ADC-DMA-FFT注意要点

总结：利用stm32上的定时器Timer来触发ADC采样，利用DMA搬运采样得到的AD值，最后用DSP库里的有关FFT运算的函数进行各次谐波幅值的获取。

一、设定数据

被采信号频率：f
被采信号周期：T = 1/f
采样频率：fs
采样周期：Ts = 1/fs

采样总点数：NPT
采样总时间：t = NPT * Ts

频谱图频率分辨率：f0 = fs/NPT

采到的被采信号周期数：NT = t/T 
一个被采信号周期的采样点数：fs/f

二、代入数据

假设我们要采集的信号频率为1kHZ，我们以32kHz的采样频率进行采集，一共采集256个点。

被采信号频率：f = 1000Hz
被采信号周期：T = 1/f  = 0.001s
采样频率：fs = 32 000Hz
采样周期：Ts = 1/fs =0.000 031 25s

采样总点数：NPT = 256
采样总时间：t = NPT * Ts = 0.008s

频谱图频率分辨率：f0 = fs/NPT = 32 000Hz/256 = 125Hz

采到的被采信号周期数：NT = t/T = 0.008s/0.001s = 8
一个被采信号周期的被采到的点数：fs/f = 32 000Hz/1000Hz = 32

如何确定采样频率和采样点数？

分辨率就是我们后面fft运算后，出现在频谱图上，横轴是频率，然后以分辨率倍数的方式出现，分辨率f0=125hz。假设输入信号里存在125hz的倍数的频率（比如600hz），那再频谱图上横轴为600hz的点就会有幅值，不过因为有正负半轴，所以真正的幅值应该是频谱图上显示幅值的两倍。

频谱图解析：

FFT后的幅度谱的横坐标是频率，并且是离散的。是以频率f0的n倍展开的。

故横坐标为… , -nf0, …-2f0, -f0, 0, f0, 2f0, … ,nf0, …

纵坐标为对应频点的幅值信息。

本次是设定的f0 = 125Hz, 理想情况下FFT后的幅度谱看起来应该会是下面这个样子：

其中横坐标0刻度的地方即频率为0的点，为直流分量 。如果被采信号没有偏置的话这一点幅度应该为0。而横坐标其它点皆为125Hz的倍数。
由于被采信号的频率为1000Hz，故横坐标的正半轴和负半轴的1000Hz处都会有“擎天柱”。但在运用中我们只会取正半轴部分，毕竟正半轴知道了也就知道了负半轴，所以负半轴可以说是没什么用的。
值得注意的是：对于直流分量来说，纵坐标对应的值就是直流分量的幅度，而其他频率对应的纵坐标值只是其实际幅值的1/2
关于这一点有一个简单的理解方式，就是它的幅度平均分到了正负半轴，导致只有一半。

代码

注意！

ADC+DMA部分： ADC注意配置成外部时钟触发，不连续转换，单通道不扫描。 DMA需要外设不自增，内存自增的方式来存储采到的连续的256个点，非循环模式。

ADC配置

1. 外部时钟触发（External Clock Trigger）：

   • 意义：使用外部时钟触发ADC开始采样。这通常用于同步采样，确保ADC在特定的时间点采样信号。

   • 原因：这样可以与其他系统时序保持一致，避免由于内部时钟的漂移引起的误差。

   • 举例：使用定时器生成触发信号，因为采样周期为32khz，那么我们可以配置定时器产生相应的周期性中断或者触发信号。32khz——0.00003125s，一个周期出发一次中断。然后将定时器的输出配置为ADC的出发源。

   • // 配置定时器，使其每31.25微秒产生一次中断（对应32kHz）
     TIM_HandleTypeDef htim;
     htim.Init.Period = (SystemCoreClock / 32000) - 1;
     htim.Init.Prescaler = 0;
     HAL_TIM_Base_Init(&htim);
     HAL_TIM_Base_Start_IT(&htim);
     
     
     公式推导：
     
     1.目标中断频率：
     
         需要每31.25微秒产生一次中断。
         31.25微秒 = 1 / 32000 秒。
     
     
     // 配置ADC，使其被定时器触发
     ADC_HandleTypeDef hadc;
     hadc.Init.ExternalTrigConv = ADC_EXTERNALTRIGCONV_Tx_TRGO; // 根据具体定时器选择触发源
     HAL_ADC_Init(&hadc);
     

     // 定时器配置（假设使用TIM2）
     void Timer2_Init(void)
     {
         __HAL_RCC_TIM2_CLK_ENABLE();
     
         TIM_HandleTypeDef htim2; //定时器句柄
         htim2.Instance = TIM2;
         htim2.Init.Prescaler = 0; //预分频器，这里是不分频的意思，即定时器时钟频率=系统时钟频率
         htim2.Init.CounterMode = TI