STM32F407实现1Hz分辨率FFT：频率、谐波、THD与相位差测量

简介：基于STM32F407的FFT频谱分析工程，实现1Hz频率分辨率，适用于信号频率测量、谐波分析、失真度评估与双波形相位差测量，面向嵌入式开发、音频处理与电力监测等方向的工程师和进阶学习者。工程利用Cortex-M4硬件浮点单元完成数字信号处理，可覆盖0~2kHz频段内的频谱分析，有效识别主导频率、谐波成分与总谐波失真。压缩包共316个文件，以C源文件、头文件、工程配置及编译输出为主，包含STM32标准外设库、LCD显示、UCOS系统相关代码，并附带hex、axf等烧录文件及PDF说明文档，便于对照工程结构进行二次开发。包体大小约3.71MB，结构紧凑完整。该资源已有10448人学习下载，适合需要直接参考或移植FFT应用的开发者使用。通过该工程可快速获得完整算法实现与硬件适配方案，节省从零搭建的时间。

1. 一块 STM32F407 能把 FFT 分辨率做到 1Hz，顺便把频率、谐波、失真度、相位差一次算完

接到一个活儿：手头只有一块 STM32F407 开发板，要测量一台音频功放的输出信号——基频多少、谐波到哪几次、失真度多大、左右声道相差多少相位。没有台式音频分析仪，那就用 F407 的 ADC 采样，跑一段 FFT，分辨率为 1Hz，把频率、谐波、失真度、相位差全部算出来。这个方向对做音频、电力监控、振动分析和电机驱动的从业者相当实用：一块十几块的芯片就能顶一台信号分析仪的核心功能。整套方案从原理到落地并不复杂，真正的功夫在采样链路的配置和几个容易被忽略的边界条件。

2. 分辨率 1Hz 的硬约束：采样率、点数与 FFT 输出频率的关系

2.1 1Hz 分辨率是怎么来的：Δf = fs / N

FFT 输出的相邻频点间隔由采样率 fs 和点数 N 唯一决定，公式是 Δf = fs / N。想要 1Hz 分辨率，必须保证 fs / N = 1Hz，也就是采样率数值上等于点数。最常用的组合就是 fs = 4096Hz，N = 4096；或者 fs = 2048Hz，N = 2048，同样得到 1Hz。注意这不是软件能“算出来”的，而是采样参数事先定死的。

这里有个概念要区分：分辨率和测量精度不是一回事。分辨率 1Hz 意味着你能区分 100.0Hz 和 101.0Hz 这两个相邻频点；而测量精度指的是读数接近真实值的程度。通过插值算法，频率测量精度可以做到 0.01Hz 量级，但分辨两个靠近的信号仍然受 1Hz 限制。很多人在这一步就搞混，导致后面做谐波分析时对频谱峰值的解读出错。

另外要意识到一个时间约束：Δf = 1 / T，T 是时域采集窗口长度。1Hz 分辨率意味着每次测量至少要采集 1 秒时域数据。无论点数怎么加，采集时间都不会变短。这是实时性设计的第一道红线。

2.2 为什么 STM32F407 做这件事是顺手的

F407 是 Cortex-M4 内核，带硬件 FPU，主频 168MHz。CMSIS-DSP 库提供了现成的 arm_rfft_fast_f32 函数，4096 点实数 FFT 一次耗时才几毫秒，相对 1 秒的采集窗口来说实时性非常充裕。相比 M3 或 M0，省掉了一整轮 FFT 手写优化工作。

它的 RAM 足够跑浮点 FFT。4096 点实数 FFT 的输入数组 16KB，输出复数数组 32KB，加上库内部临时缓冲，几十 KB 的占用对 F407 的 192KB SRAM 来说很宽裕。再往上到 8192 点也能跑，但内存明显紧张，所以 4096 点是这个场景下最舒服的配置。

第三点最关键：F407 支持双 ADC 规则同步采样，也就是 ADC1 和 ADC2 可以同时启动转换。两个波形要测相位差，就必须保证两路信号在同一个采样时刻被采到。用单 ADC 轮流扫描两个通道，通道间天然存在一个转换时间差，相位差结果会引入系统性偏差。双 ADC 硬件同步把这个隐患从根上去掉了。

2.3 参数组合怎么选：拿 4096 点做模板

下面这张表列出常见组合，全部保证 Δf = 1Hz。

最大可测频率受奈奎斯特限制，有效频谱范围是 0 到 fs/2。谐波分析的次数上限也随之确定：基频附近能看的谐波次数最多到 floor((fs/2) / f_base)。以 4096Hz 采样、被测基频 50Hz 为例，理想情况下能数到 40 次谐波。实际工程里我一般看 2 到 10 次谐波就足够评估失真度，后面代码里也会按这个范围写。

选型建议：被测信号最高频率在 2kHz 以内、RAM 预算充足，用 fs = 4096Hz、N = 4096。这个组合既留出了 2kHz 的分析带宽，又不会让 F407 的内存吃紧，后面的代码都按这个配置展开。

3. 采样链路搭建：定时器触发 + 双 ADC 同步采样 + DMA 搬运

3.1 定时器触发频率配到 4096Hz：参数怎么算

FFT 对采样点的时间均匀性极其敏感。用 ADC 连续转换加软件延时的方式，中断响应和代码分支会让采样间隔产生抖动，反映到频谱上就是噪声底抬高。正确的做法是用定时器更新事件触发 ADC 启动转换，让硬件决定采样时刻。

以 TIM8 为例，它挂在 APB2 总线上，时钟 84MHz。目标是输出 4096Hz 触发频率，配置 PSC = 0，ARR = 20507，实际频率 = 84MHz / (20507 + 1) ≈ 4095.97Hz，偏差大约 0.0007%，这个误差对 1Hz 分辨率来说已经足够用了。不是说一定要配出精确的 4096.000Hz，后面坑里会讲怎么校准这个微小偏差。

/* TIM8 定时触发 ADC，目标是 4096Hz 采样率 */ TIM_TimeBaseInitTypeDef timerInit; timerInit.TIM_Prescaler = 0; timerInit.TIM_CounterMode = TIM_CounterMode_Up; timerInit.TIM_Period = 20507; /* 84MHz / (20507+1) ≈ 4095.97 Hz */ timerInit.TIM_ClockDivision = TIM_CKD_DIV1; TIM_TimeBaseInit(TIM8, &timerInit); /* 用更新事件做 ADC 触发源 */ TIM_SelectOutputTrigger(TIM8, TIM_TRGOSource_Update); TIM_Cmd(TIM8, ENABLE);

这里的 PSC 和 ARR 是分频核心参数：定时器时钟先除以 (PSC+1)，再数 (ARR+1) 个脉冲输出一次事件。TIM8 的时钟源来自 APB2，若 APB2 分频系数不为 1，定时器时钟会翻倍，配参数前先确认时钟树。PSC = 0 表示不分频，ARR 才是决定频率的主参量。

3.2 双 ADC 规则同步采样：ADC1 做主，ADC2 从动

相位差测量的基础是两路信号同时采样。F407 的 ADC1 和 ADC2 可以组成规则同步模式：由同一个外部触发信号（TIM8 的 TRGO）同时启动两路转换，转换结果先后落寄存器。这种模式下只需要给 ADC1 配置 DMA 请求，数据会按先 ADC1 后 ADC2 的顺序排进缓冲区，省掉一路 DMA 通道。

ADC 时钟从 APB2 分频得到，84MHz 分 4 得到 21MHz，稳定且满足 F407 最高 36MHz 的约束。采样时间我一般设成 15 周期，既能保证信号源阻抗较高时采样电容充分充电，又不会拖慢整体节奏。每通道单次转换约 (15+12) / 21MHz ≈ 1.29us，对 4096Hz 采样率绰绰有余。

/* 双 ADC 规则同步模式初始化 */ ADC_CommonInitTypeDef adcCommonInit; adcCommonInit.ADC_Mode = ADC_Mode_RegSimult; /* 规则同步 */ adcCommonInit.ADC_Prescaler = ADC_Prescaler_Div4; /* 84MHz / 4 = 21MHz */ adcCommonInit.ADC_DMAAccessMode = ADC_DMAAccessMode_1; /* 双 ADC 共用 ADC1 的 DMA，先读 ADC1 后读 ADC2 */ adcCommonInit.ADC_TwoSamplingDelay = ADC_TwoSamplingDelay_5Cycles; ADC_CommonInit(&adcCommonInit); /* ADC1 主通道：PA0 */ ADC_InitTypeDef adcInit; adcInit.ADC_Resolution = ADC_Resolution_12b; adcInit.ADC_ScanConvMode = DISABLE; /* 每路只采一个通道 */ adcInit.ADC_ContinuousConvMode = DISABLE; /* 由外部触发启动 */ adcInit.ADC_ExternalTrigConvEdge = ADC_ExternalTrigConvEdge_Rising; adcInit.ADC_ExternalTrigConv = ADC_ExternalTrigConv_T8_TRGO; /* TIM8 触发 */ ADC_Init(ADC1, &adcInit); /* ADC2 从通道：PA1，触发源不用再配，由同步模式接管 */ adcInit.ADC_ExternalTrigConv = 0; ADC_Init(ADC2, &adcInit); ADC_ExternalTrigConvCmd(ADC1, ENABLE); ADC_DMACmd(ADC1, ENABLE); ADC_Cmd(ADC1, ENABLE); ADC_Cmd(ADC2, ENABLE);

两个 ADC 的通道配置要放在外面单独完成，这里略去 ADC_RegularChannelConfig 的调用。关键点是 ADC2 不要再独立配外部触发，否则同步时序会被打破。ADC_DMAAccessMode_1 决定 DMA 读出顺序，这个字段在不同固件库版本里名字可能略有差异，但含义不变。

3.3 DMA 循环模式与双缓冲切换

采样数据的搬运交给 DMA，CPU 只处理中断标志。用循环模式让 DMA 不停地把 ADC 结果写进内存，每写满半块缓冲区触发一次半传输中断，每写满整块触发传输完成中断。CPU 在中断里只做标志位登记，主循环检测到标志后处理已填充的数据块，这样不会丢采样点。

以 4096 点为例，缓冲区要容纳两个通道各 4096 个 16 位 ADC 原始数据，所以是 8192 个 uint16_t。为什么不用 float 数组直接给 FFT？因为 DMA 搬运的是 ADC 外设寄存器里的 12 位原始值，先以 uint16_t 存放，主循环里再统一转换成 float，职责清晰。

#define N_POINTS 4096 volatile uint16_t adc_buffer[2 * N_POINTS]; /* 8192 个 u16，按 [ADC1][ADC2] 交替排列 */ volatile uint8_t dma_half_done = 0; volatile uint8_t dma_done = 0; void DMA2_Stream0_IRQHandler(void) { if (DMA_GetITStatus(DMA2_Stream0, DMA_IT_HTIF)) { dma_half_done = 1; /* 前半块填满 */ DMA_ClearITPendingBit(DMA2_Stream0, DMA_IT_HTIF); } if (DMA_GetITStatus(DMA2_Stream0, DMA_IT_TCIF)) { dma_done = 1; /* 整块填满 */ DMA_ClearITPendingBit(DMA2_Stream0, DMA_IT_TCIF); } }

DMA 配置时要指定外设地址、存储器地址、传输长度和循环模式，其中传输长度就是 8192 个半字。外设地址是 ADC1 的 DR 寄存器，存储器地址指向 adc_buffer。注意 DMA_IT_HTIF 和 DMA_IT_TCIF 的清除顺序，先清标志再处理业务逻辑，避免中断重入。两个标志位在主循环里处理完当前数据后必须手动清零。

主循环处理逻辑是这样：检测到 dma_half_done 时处理 adc_buffer 的前半段（前 4096 个 u16），检测到 dma_done 时处理后 4096 个 u16，处理完清零对应标志。这样数据仍在 DMA 连续写入，不会互相踩踏。F407 的 DMA2 流 0 连到 ADC1，中断优先级建议配到较高档位，保证采样不丢。

4. 从频谱到测量结果：频率、谐波、失真度和相位差的计算

4.1 峰值搜索与频率读数：1Hz 粗读加抛物线插值

拿到时域数据先转 float，然后调用 CMSIS-DSP 的实数 FFT 函数把时域信号变成频谱。4096 点输入经过 arm_rfft_fast_f32 后得到复数频谱，实部和虚部交替存放，长度 2N=8192 个 float，其中前 2048 个复数对应 0 到 2048Hz 的有效频点。注意跳过第 0 个频点，那是直流分量，幅度通常远大于交流信号，不跳过会让峰值搜索误判。

/* 实数 FFT：4096 点时域 -> 8192 个浮点复数输出 */ arm_rfft_fast_instance_f32 fftInst; arm_rfft_fast_init_f32(&fftInst, 4096); arm_rfft_fast_f32(&fftInst, timeData, fftOutput, 0); /* 0 表示正向变换 */ /* 计算幅度谱：2048 个有效频点 */ float32_t mag[2048]; arm_cmplx_mag_f32(fftOutput, mag, 2048); /* 跳过直流，找最大幅度频点 */ uint32_t kMax = 1; float32_t vMax = 0; for (uint32_t k = 1; k < 2048; k++) { if (mag[k] > vMax) { vMax = mag[k]; kMax = k; } } /* Δf = 1Hz 时，bin 编号就是频率粗略值 */ float freqCoarse = (float)kMax;

arm_rfft_fast_f32 的第三个参数 0 是正向变换标志。输出排布是实部虚部交替，比如频点 k 的实部在 fftOutput[2k]，虚部在 fftOutput[2k+1]。这个布局后面算相位时会再次用到。峰值搜索得到的 kMax 直接就是 1Hz 分辨率下的频率读数，但想要更高精度，可以在 kMax 附近做抛物线插值。

float y0 = mag[kMax - 1]; float y1 = mag[kMax]; float y2 = mag[kMax + 1]; float delta = 0.5f * (y0 - y2) / (y0 - 2.0f * y1 + y2); float freqFine = (float)kMax + delta;

插值公式的数学含义是：把峰值附近三个幅度点拟合成一条抛物线，delta 是抛物线顶点的偏移量。信噪比良好、信号较纯时，freqFine 的精度能到 0.01Hz 量级。但要注意两个边界：kMax 若等于 1，y0 取了直流点；kMax 若等于 2047，y2 越界。实际代码里要先做边界保护，否则数组越界会返回一个垃圾值。

4.2 谐波分析与 THD 失真度：幅度平方和开根号

失真度也叫 THD，定义是所有谐波功率之和除以基波功率，开根号后换算成百分比。谐波在频谱上的位置是基波频点的整数倍：基波在 kMax，二次谐波在 2×kMax，三次在 3×kMax。用这个关系直接查频谱幅度，不需要额外搜索峰值。

/* 从 2 次到 10 次谐波累加幅度平方，算 THD 百分比 */ float harmSum = 0.0f; for (uint32_t n = 2; n <= 10; n++) { uint32_t k = n * kMax; if (k >= 2048) break; /* 超出奈奎斯特范围就停 */ harmSum += mag[k] * mag[k]; } float thdPercent = sqrtf(harmSum) / mag[kMax] * 100.0f;

为什么用幅度而不是有效值？因为 THD 公式里分子分母的有效值都带同一个 1/√2 系数，约掉之后直接用幅度比。如果你的目的是对比不同设备，统一用幅度比就好。谐波次数上限取 10 次是一种折中：太高次的谐波幅度通常已经淹没在噪声底里，算进去反而拉高 THD。若被测基波频率较高，循环里要判断 k 是否超出 2048，超出就 break。

这一个计算块有个前提：kMax 必须是基波的真实位置。如果信号本身失真严重，或者加窗不合适，频谱峰值可能偏离实际基频，谐波索引链全部错位。这也是为什么谱峰搜索要带插值，同时被测信号进来之前最好有一个带通滤波器，把带外噪声和直流偏置压一压。

4.3 相位差测量：两个通道同帧 FFT，取基波 bin 相位做差

相位差的实现思路是：双通道数据分别做 FFT，在两个频谱的基波频点处各自取相位角，然后相减。取相位的函数用 atan2f，它能根据实部和虚部的正负号正确判断象限，得到的角度范围是 -180° 到 180°。atanf 只有一个参数，遇到实部为负的分母会出错，所以这里必须用 atan2f。

/* 通道 A 和通道 B 的基波相位角，单位弧度 */ float reA = fftOutA[2 * kMax]; float imA = fftOutA[2 * kMax + 1]; float reB = fftOutB[2 * kMax]; float imB = fftOutB[2 * kMax + 1]; float phA = atan2f(imA, reA); float phB = atan2f(imB, reB); /* 相位差并折算到 -180° ~ 180° */ float diff = phA - phB; while (diff > 3.14159265f) diff -= 2 * 3.14159265f; while (diff < -3.14159265f) diff += 2 * 3.14159265f; float diffDeg = diff * 180.0f / 3.14159265f;

需要注意这里的正负含义：diffDeg 为正表示通道 A 相位超前通道 B，为负表示滞后。如果你习惯以通道 B 为参考，就取相反数。另一个坑是相位差公式里必须用同一帧数据里的两个通道，程序里要确保 fftOutA 和 fftOutB 是从同一次双 ADC 同步触发采集得到的样本算出来的，不是两帧拼接。

相位测量精度依赖两个条件：一是基波频点 kMax 找得准，二是两通道采样时刻严格同步。前者用前面的抛物线插值解决不了的——相位插值比幅度插值复杂得多，所以工程上更常用的是把采样率做好、把触发做好，让整周期采样成立；后者就是双 ADC 同步模式。如果这两点没保证，相位差读数会无规律抖动，下文会在踩坑章节里展开说。

5. 相位差测量的三个典型坑与排查顺序

5.1 先定位现象，再动手改代码

遇到读数不对的时候，我习惯按这个顺序排查：先看时域波形，确认 ADC 有没有采到干净的数据；再看频谱形态，判断是不是泄漏或混叠；最后才怀疑算法本身。很多相位差问题其实出在频率测量不准上，基波频点偏了一个 bin，相位结果自然全乱。把这个思路理顺，很多“玄学”问题都能在十分钟内定位。

5.2 坑一：频率读数整体偏大或偏小，偏差固定

现象：信号发生器输出精确的 1000Hz 正弦，FFT 读数是 1000.4Hz 或 999.6Hz，而且偏差比例固定。

原因：定时器分频配不出精确的 4096Hz，晶振本身也存在 ppm 级误差。ARR = 20507 只是把频率凑到接近 4096Hz，实际采样率与标称值之间存在系统性偏差，这个偏差会线性地影响频率读数。

解决：用已知频率的参考信号做一次校准。测出一个频率读数误差比例，记到 Flash 里，后续所有频率读数乘以修正系数。校准流程我一般写成一个独立函数：输入参考频率 fRef，测量读数 fMeas，计算修正系数 fRef / fMeas，下次测量时freqTrue = freqCoarse * coeff。注意校准要在信号稳定、插值算法生效的条件下进行。

5.3 坑二：频谱在信号附近“变胖”，基波旁边全是裙边，谐波幅度失真

现象：输入一个干净的 1000Hz 正弦，幅度谱上 1000Hz 附近出现一大片衰减拖尾，而不是一个窄峰。谐波分析的结果里多了很多不存在的频率分量。

原因：非整周期采样造成的频谱泄漏。4096 点窗口刚好采集 1 秒，如果被测信号的周期不能整除 1 秒，窗口首尾就不是完整的整数个周期，FFT 默认认为时域信号首尾相连，于是产生剧烈跳变，能量泄漏到所有频点。

解决：两个方向。一是加窗函数，汉宁窗或者布莱克曼窗能把泄漏压低几个数量级，代价是主瓣变宽，1Hz 分辨率下两个相距很近的频率分量更难分开。二是尽量让采样点和信号周期匹配，比如被测基波固定为 50Hz 时，把采样率调成 3200Hz 且点数仍为 4096，这样 50Hz 落在整周期上。两者权衡下来，通用性最强的还是加窗结合插值，但要明确：加窗之后 THD 公式里的幅度需要乘以窗函数的相干增益修正系数，否则基波和谐波幅度都比真实值小。

5.4 坑三：相位差读数忽大忽小，但频率读数很稳

现象：同一对信号重复测量，相位差今天 35°，明天 42°，测量期间仪器和接线都没动过；频率读数却始终稳定。

原因：两路信号不是同一时刻采的，或者触发时刻在抖动。单 ADC 轮询两个通道时，两路之间天然隔了一个通道的转换时间，这个时间差换算成相位偏转会随频率变化；触发信号带噪声时，每次采样启动时刻漂移，相位差自然跟着跳。

解决：改用双 ADC 同步采样，也就是第 3 章描述的那种硬件同步结构。再做一次硬件整形：被测信号进 ADC 之前，先用比较器把正弦整成方波，再用方波沿触发定时器同步，触发抖动可以从几十微秒压到几百纳秒。如果只是低频测量、要求不高，多采几帧求平均值也能缓解，但治标不治本。

5.5 坑四：FFT 结果偶尔整帧丢失，或数据被 DMA 覆盖

现象：程序跑一段时间，串口输出的频谱突然变了一个样，幅度整体不对劲，过一会儿又恢复；或者某次测量结果里混进了上一帧的数据。

原因：DMA 在循环写的同时，CPU 处理数据不及时，半传输中断和传输完成中断里的标志位被主循环遗漏；下一次 DMA 半帧到来时，CPU 还在跑上一次的 FFT，缓冲区被覆盖。

解决：把标志位处理改成“缓冲区分块锁”的思路：中断里切换当前可读写缓冲区的索引，主循环检测到新索引变化才处理对应的数据块，处理完标记释放。同时给 DMA 中断配较高的抢占优先级，防止被其他外设中断长时间打断。实际上 4096 点 FFT 在 F407 上只要几毫秒，远快于 4096 点采集所需的 1 秒，正常情况下 CPU 有大量空闲，出现覆盖基本可以断定是中断优先级或标志位管理出了逻辑问题。

6. 闭环验收：三组测试数据和一个进阶技巧

整套系统写完，先别急着接真实信号，用信号发生器和示波器做三组闭环验证。

第一组验证频率读数：分别输出 100Hz、101Hz、1000Hz 正弦，确认读数对应 100.0、101.0、1000.0 左右，误差在插值后小于 0.05Hz。再把两路信号输出一个 100Hz 和一个 101Hz 混合，确认频谱上能清晰看到两个峰。第二组验证谐波和 THD：输出一个 1kHz 正弦加 3kHz 正弦（幅度为基波的 1/3），理论 THD 约 33.3%，程序算出来应该非常接近；如果信号本身是标准方波，THD 理论值约 48.3%，这也是一个不错的参考点。第三组验证相位差：用信号发生器双通道输出同频、相位差 30°，程序读取应该在 29.5° 到 30.5° 之间；再把相位差调成 -90° 和 179°，确认正负和边界情况没被角度折算搞乱。

一个比较实用的进阶技巧：如果你后续只关心基波附近一小段频谱，可以不用跑完整 FFT，改用 Goertzel 算法只算几个目标频点，内存占用和计算开销都大幅下降。但像这个标题里的完整需求——频率、谐波、失真度、相位差全都算——整帧 FFT 的通用性远胜过后缀算法，因为谐波分析需要把整个有效频段的幅度谱都拿出来用。

最后说一个我的习惯：每次改完采样率或点数参数，我都会先用信号发生器扫一遍全频段，确认频率读数和 THD 都稳定了才接真实前端。这套方案做下来最深的感受是，FFT 本身不难，难的是把采样触发和两通道同步这种硬件细节压稳。参数链路理顺，后面的算法都是顺水推舟。希望帮到你。

本文还有配套的精品资源，点击获取