STM32 FFT精准测频测相：采样定理、窗函数与相位解卷绕实战-尧图手机网站定制

简介：本资源是一套基于STM32F407标准库实现的FFT频域分析完整工程，面向嵌入式开发初学者与中级工程师，解决正弦波信号在资源受限MCU上实时、精准提取幅值、频率及相位差的核心难题。项目覆盖ADC采样、汉宁窗预处理、标准库FFT计算、频谱峰值检测与相位解算全流程，适用于电机控制、音频监测、电力谐波分析等需时频联合分析的工业场景。压缩包共283个文件，含46个C源码（如stm32f4xx_rcc.c、stm32f4xx_tim.c）、60个头文件（h）、50个编译中间文件（o）及Keil工程配置文件（uvprojx、uvoptx、axf、hex等），结构完整，可直接编译下载运行。目前已有1666人学习下载，提供可复现的端到端代码框架、DMA+定时器协同采样配置、频谱幅值归一化处理逻辑及相位差计算注释详尽的参考实现，助读者快速掌握嵌入式FFT工程化落地的关键细节。

1. STM32上实现FFT算法精准测量正弦波信号的幅值、频率和相位差：不是调个库就完事，而是要啃透采样定理、窗函数选择和频谱泄露这三块硬骨头

你手头有一块STM32F103C8T6（俗称“蓝 pill”），接了个运放调理后的正弦波传感器信号，想用它实时测出信号的幅值、频率、相位差——不是靠示波器目测，也不是靠过零检测这种粗粒度方法，而是要达到±0.5%幅值误差、±0.1Hz频率分辨率、±2°相位精度。这时候翻官方HAL库，发现arm_math.h里确实有arm_cfft_f32()，但直接套用你会发现：频谱主瓣宽得像摊煎饼、旁瓣高得能盖楼、两个相近频率根本分不开、相位跳变毫无规律……这不是FFT不灵，是你没给它配对的采样策略、没选对窗函数、没做相位校准。本项目完全基于标准库（非HAL），用纯C+CMSIS-DSP实现，不依赖任何第三方GUI或上位机，所有计算在MCU本地闭环完成。适合嵌入式电能质量分析、电机振动监测、音频前端预处理等对实时性与确定性要求高的场景。如果你正在为FFT结果“玄学漂移”抓狂，或者被“为什么理论频率点上没峰值”反复暴击，这篇就是为你写的血泪复现笔记。

2. 从时域到频域：为什么STM32必须手撕FFT流程，而不是直接调arm_cfft_f32()

2.1 为什么不能“拿来即用”：CMSIS-DSP的FFT接口本质是黑匣子

CMSIS-DSP提供的arm_cfft_f32()函数，其输入要求是复数数组（实部存原始采样值，虚部全置0），输出也是复数数组。但问题来了：

它不做任何采样控制，你得自己配置ADC触发、DMA搬运、缓冲区管理；
它不处理频谱泄露，输出的模值直接拿去算幅值会系统性偏高；
它不提供相位解包逻辑，atan2(imag, real)出来的结果是[-π, π]区间，而真实相位差需跨周期连续解卷绕；
它不校准频点偏移，由于实际采样率存在晶振误差，FFT bin中心频率≠理论值，导致频率测量偏差。

提示：arm_cfft_f32()只是DFT计算引擎，不是完整频谱分析仪。就像给你一台发动机，但没配变速箱、离合器和仪表盘——你得自己造。

2.2 完整信号链拆解：6个不可跳过的环节

一个能落地的FFT测量系统，必须串起以下6个环节，缺一不可：

环节	关键动作	为什么必须自控
1. 采样控制	ADC连续扫描+DMA双缓冲+精确触发源（TIMx TRGO）	HAL_Delay()或SysTick会导致采样间隔抖动，破坏奈奎斯特条件
2. 数据预处理	直流偏置消除（滑动平均滤波）、增益归一化（ADC满量程映射为±1.0f）	偏置导致频谱泄漏加剧，未归一化使幅值与ADC参考电压强耦合
3. 窗函数应用	实时乘窗（Hanning窗最常用），窗系数预存在const float win[1024]中	不加窗时矩形窗旁瓣衰减仅-13dB，Hanning窗达-31dB，抗干扰能力翻倍
4. FFT执行	调用`arm_cfft_f32()`+`arm_cmplx_mag_f32()`+`arm_cmplx_phase_f32()`	注意：`arm_cmplx_phase_f32()`输出单位是弧度，需转为角度并做解卷绕
5. 频谱分析	寻峰（3点抛物线插值）、频率校准（根据实测采样率修正bin索引）、幅值补偿（窗函数能量损失补偿）	峰值不在整数bin上？不插值误差超5%；不补偿窗损耗？幅值偏低36%（Hanning窗）
6. 相位差解算	双通道同步采样+共用同一组窗系数+相位差 = φ₁ - φ₂ + 2kπ（k由相邻帧相位跳变判定）	单通道无法测相位差；不同步采样引入固定相位偏移

2.3 标准库下的ADC+DMA配置实操（以STM32F103为例）

// ADC初始化：关键参数必须锁定 void ADC1_Init(void) { ADC_InitTypeDef ADC_InitStructure; DMA_InitTypeDef DMA_InitStructure; RCC_APB2PeriphClockCmd(RCC_APB2Periph_ADC1 | RCC_APB2Periph_GPIOA, ENABLE); RCC_AHBPeriphClockCmd(RCC_AHBPeriph_DMA1, ENABLE); // PA0 作为 ADC1_IN0，模拟输入 GPIO_InitTypeDef GPIO_InitStructure; GPIO_InitStructure.GPIO_Pin = GPIO_Pin_0; GPIO_InitStructure.GPIO_Mode = GPIO_Mode_AIN; GPIO_Init(GPIOA, &GPIO_InitStructure); // ADC配置：12位、连续转换、右对齐、不扫描 ADC_DeInit(ADC1); ADC_InitStructure.ADC_Mode = ADC_Mode_Independent; ADC_InitStructure.ADC_ScanConvMode = DISABLE; // 关闭扫描，单通道 ADC_InitStructure.ADC_ContinuousConvMode = ENABLE; // 连续模式 ADC_InitStructure.ADC_ExternalTrigConv = ADC_ExternalTrigConv_T1_CC1; // TIM1 CH1触发 ADC_InitStructure.ADC_DataAlign = ADC_DataAlign_Right; ADC_InitStructure.ADC_NbrOfChannel = 1; ADC_Init(ADC1, &ADC_InitStructure); // ADC通道0：采样时间设为239.5周期（最高精度档） ADC_RegularChannelConfig(ADC1, ADC_Channel_0, 1, ADC_SampleTime_239Cycles5); // DMA配置：双缓冲，循环模式，每次传输1024点 DMA_DeInit(DMA1_Channel1); DMA_InitStructure.DMA_PeripheralBaseAddr = (uint32_t)&ADC1->DR; DMA_InitStructure.DMA_MemoryBaseAddr = (uint32_t)adc_buffer_a; // 缓冲区A DMA_InitStructure.DMA_DIR = DMA_DIR_PeripheralSRC; DMA_InitStructure.DMA_BufferSize = SAMPLE_NUM; // SAMPLE_NUM = 1024 DMA_InitStructure.DMA_PeripheralInc = DMA_PeripheralInc_Disable; DMA_InitStructure.DMA_MemoryInc = DMA_MemoryInc_Enable; DMA_InitStructure.DMA_PeripheralDataSize = DMA_PeripheralDataSize_HalfWord; DMA_InitStructure.DMA_MemoryDataSize = DMA_MemoryDataSize_HalfWord; DMA_InitStructure.DMA_Mode = DMA_Mode_Circular; // 循环模式，自动翻页 DMA_InitStructure.DMA_Priority = DMA_Priority_High; DMA_InitStructure.DMA_M2M = DMA_M2M_Disable; DMA_Init(DMA1_Channel1, &DMA_InitStructure); // 启用DMA半传输/全传输中断，用于双缓冲切换 DMA_ITConfig(DMA1_Channel1, DMA_IT_TC | DMA_IT_HT, ENABLE); NVIC_EnableIRQ(DMA1_Channel1_IRQn); ADC_DMACmd(ADC1, ENABLE); // ADC启动DMA ADC_Cmd(ADC1, ENABLE); ADC_ResetCalibration(ADC1); while(ADC_GetResetCalibrationStatus(ADC1)); ADC_StartCalibration(ADC1); while(ADC_GetCalibrationStatus(ADC1)); }

参数说明：

ADC_ExternalTrigConv_T1_CC1：使用TIM1通道1的PWM边沿作为触发源，确保采样时刻绝对稳定；
DMA_Mode_Circular：循环模式让DMA自动在adc_buffer_a和adc_buffer_b间切换，避免数据覆盖；
DMA_IT_TC | DMA_IT_HT：全传输（TC）和半传输（HT）中断分别标记1024点采满和512点采满，用于双缓冲乒乓操作；
ADC_SampleTime_239Cycles5：最大采样时间，牺牲速度换精度，对低频信号尤其关键。

3. 窗函数、频谱泄露与幅值补偿：为什么你的FFT幅值总比真值小36%

3.1 矩形窗 vs Hanning窗：一张图看懂频谱泄露的本质

当采集N点有限长信号时，相当于原始无限长信号乘以一个长度为N的矩形窗。矩形窗的频域响应是sinc函数，主瓣宽、旁瓣高。若信号频率不恰好落在FFT bin中心（即f₀ ≠ k·fs/N），能量就会泄漏到邻近bin，导致：

主峰幅度降低；
旁瓣抬高，淹没弱信号；
频率估计偏移。

Hanning窗定义为：
$$ w(n) = 0.5 - 0.5 \cos\left(\frac{2\pi n}{N-1}\right), \quad n = 0,1,\dots,N-1 $$
其频域主瓣宽为矩形窗的2倍，但旁瓣衰减达-31dB，且第一旁瓣已低于-40dB。这意味着：即使信号频率偏移半个bin，主峰幅度损失也<0.1%，远优于矩形窗的>20%。

3.2 Hanning窗系数生成与内存优化技巧

// 预生成Hanning窗系数（float型，1024点） const float hanning_win[SAMPLE_NUM] = { 0.000000f, 0.000006f, 0.000024f, 0.000054f, 0.000096f, 0.000150f, // ... 中间省略990项（实际代码中应完整列出或用脚本生成） 0.000150f, 0.000096f, 0.000054f, 0.000024f, 0.000006f, 0.000000f }; // 应用窗函数（定点优化版：避免浮点除法） void apply_hanning_window(float *pSrc, uint16_t len) { uint16_t i; for(i = 0; i < len; i++) { pSrc[i] *= hanning_win[i]; // 直接查表乘，无计算开销 } }

关键细节：

系数声明为const并放在.rodata段，不占RAM；
不用sin/cos实时计算，查表速度提升10倍以上；
系数精度用float足够，double在Cortex-M3上无硬件加速，反而拖慢。

3.3 幅值补偿系数推导：为什么Hanning窗后要×2.0

窗函数会衰减信号总能量。Hanning窗的能量衰减因子（Coherent Gain）定义为：
$$ CG = \frac{1}{N}\sum_{n=0}^{N-1} w(n) $$
对Hanning窗，CG ≈ 0.5。这意味着：若原始信号幅值为A，加窗后时域均方根值降为0.5A。而FFT输出的幅值谱（经arm_cmplx_mag_f32()后）对应的是单边谱，需乘以2才能还原真实幅值（因双边谱对称，单边谱只取一半能量）。故总补偿系数 = 2 / CG = 2 / 0.5 =4.0。

但注意：CMSIS-DSP的arm_cfft_f32()输出是未归一化DFT，其模值需再除以N才得真实频谱密度。因此最终幅值公式为：
$$ A_{real} = \frac{2}{N} \times \frac{1}{CG} \times |X[k]| = \frac{2}{N} \times 4.0 \times |X[k]| = \frac{8.0}{N} \times |X[k]| $$

对N=1024，补偿系数 = 8.0 / 1024 =0.0078125。这就是你在代码里看到mag_val * 0.0078125f的由来。

3.4 频率精修：3点抛物线插值法实战

FFT bin是离散的，真实频率往往落在bin之间。3点抛物线插值利用峰值bin及其左右邻点，拟合抛物线求顶点：
设峰值在bink，幅值为Y[k-1],Y[k],Y[k+1]，则频率偏移量：
$$ \delta = \frac{Y[k+1] - Y[k-1]}{2 \times (2Y[k] - Y[k-1] - Y[k+1])} $$
真实频率：
$$ f = (k + \delta) \times \frac{f_s}{N} $$

// 寻峰+插值（输入：mag_spectrum[512]，输出：*freq_out, *amp_out） void find_peak_and_interpolate(const float *mag_spectrum, uint16_t len, float *freq_out, float *amp_out) { uint16_t k_max = 1; // 跳过DC bin（bin 0） float max_mag = mag_spectrum[1]; // 粗寻峰：找最大值索引（限于1~255 bin，避开高频噪声） for(uint16_t k = 2; k < 256; k++) { if(mag_spectrum[k] > max_mag) { max_mag = mag_spectrum[k]; k_max = k; } } // 三邻点插值（确保不越界） if(k_max == 0 || k_max == len-1) { *freq_out = k_max * SAMPLING_FREQ / len; *amp_out = max_mag * 0.0078125f; // 补偿系数 return; } float y0 = mag_spectrum[k_max-1]; float y1 = mag_spectrum[k_max]; float y2 = mag_spectrum[k_max+1]; // 抛物线插值公式（δ ∈ [-0.5, 0.5]） float delta = (y2 - y0) / (2.0f * (2.0f*y1 - y0 - y2)); *freq_out = (k_max + delta) * SAMPLING_FREQ / len; *amp_out = (y1 + 0.25f*(y2 - y0)*delta) * 0.0078125f; // 插值后幅值补偿 }

参数说明：

SAMPLING_FREQ：实测采样率（非理论值！见第4章校准）；
len=1024→ 单边谱长度512，故k_max搜索范围为1~255（对应0.5~125Hz，若fs=250Hz）；
y1 + 0.25*(y2-y0)*delta是插值后幅值的快速近似，比直接用y1精度高一个数量级。

4. 相位差测量的致命陷阱与双通道同步方案

4.1 为什么单通道FFT无法测相位差？相位的物理意义是什么

FFT输出的相位φ(k)是该频率分量相对于采样窗口起始时刻的初相。若两路信号分别用独立ADC+DMA采集，即使触发源相同，DMA搬运延迟、缓存对齐差异也会导致两路数据在时间轴上存在亚微秒级偏移。此时计算φ₁ - φ₂得到的是系统延迟引入的固定偏移，而非真实信号相位差。必须保证：

两路ADC共用同一触发源；
两路DMA使用同一时钟域、同一DMA控制器（如都走DMA1_Channel1/2）；
两路缓冲区严格对齐（起始地址、长度、双缓冲翻页时刻完全一致）。

4.2 双通道同步采集硬件连接与寄存器配置

// 双通道ADC配置（PA0=CH0, PA1=CH1） void ADC1_Init_Dual(void) { // ... GPIO初始化：PA0, PA1 均设为模拟输入 ... ADC_DeInit(ADC1); ADC_InitStructure.ADC_Mode = ADC_Mode_Independent; ADC_InitStructure.ADC_ScanConvMode = ENABLE; // 必须开启扫描！ ADC_InitStructure.ADC_ContinuousConvMode = ENABLE; ADC_InitStructure.ADC_ExternalTrigConv = ADC_ExternalTrigConv_T1_CC1; ADC_InitStructure.ADC_DataAlign = ADC_DataAlign_Right; ADC_InitStructure.ADC_NbrOfChannel = 2; // 2通道 ADC_Init(ADC1, &ADC_InitStructure); // 通道0（PA0）排第1，采样时间239.5周期 ADC_RegularChannelConfig(ADC1, ADC_Channel_0, 1, ADC_SampleTime_239Cycles5); // 通道1（PA1）排第2，采样时间239.5周期 ADC_RegularChannelConfig(ADC1, ADC_Channel_1, 2, ADC_SampleTime_239Cycles5); // DMA配置：传输2048字（1024次×2通道），内存地址按通道交错排列 DMA_InitStructure.DMA_PeripheralBaseAddr = (uint32_t)&ADC1->DR; DMA_InitStructure.DMA_MemoryBaseAddr = (uint32_t)adc_dual_buffer; // [ch0_0,ch1_0,ch0_1,ch1_1,...] DMA_InitStructure.DMA_BufferSize = SAMPLE_NUM * 2; // 2048 half-word // ... 其余DMA参数同单通道 ... }

关键设计：

ADC_ScanConvMode = ENABLE：扫描模式下，ADC自动在CH0→CH1→CH0…间切换，保证严格时间对齐；
DMA_BufferSize = 1024*2：DMA一次搬2048个16位数据，内存布局为[ch0_0, ch1_0, ch0_1, ch1_1, ..., ch0_1023, ch1_1023]；
解包时，ch0_data[i] = adc_dual_buffer[2*i]，ch1_data[i] = adc_dual_buffer[2*i+1]，天然对齐。

4.3 相位解卷绕（Phase Unwrapping）：如何让-179°平滑过渡到+175°

FFT相位输出范围是[-π, π]，当真实相位连续增加时，会在±π处发生跳变（如178°→-178°）。解卷绕就是检测这种跳变并加减2π修正。双通道相位差解算流程：

对每通道分别计算相位谱phase_ch0[k],phase_ch1[k]；
计算初始相位差delta_phi_raw[k] = phase_ch0[k] - phase_ch1[k]；
对delta_phi_raw序列进行一维解卷绕（非对每个bin单独解）；
取基波bin（k_fund）处的解卷绕后相位差，转为角度制。

// 一维相位解卷绕（输入：delta_phi_raw[512]，输出：delta_phi_unwrap[512]） void phase_unwrap_1d(float *phi, uint16_t len) { for(uint16_t i = 1; i < len; i++) { float diff = phi[i] - phi[i-1]; if(diff > M_PI) phi[i] -= 2.0f * M_PI; if(diff < -M_PI) phi[i] += 2.0f * M_PI; } } // 双通道相位差计算主函数 void calc_phase_diff(float *ch0_buf, float *ch1_buf, float *freq_out, float *phase_diff_out) { // 1. 分别加窗、FFT、求幅值/相位 apply_hanning_window(ch0_buf, SAMPLE_NUM); apply_hanning_window(ch1_buf, SAMPLE_NUM); arm_cfft_f32(&S, ch0_buf, 0, 1); // 正向FFT arm_cfft_f32(&S, ch1_buf, 0, 1); arm_cmplx_mag_f32(ch0_buf, mag_ch0, SAMPLE_NUM); arm_cmplx_mag_f32(ch1_buf, mag_ch1, SAMPLE_NUM); arm_cmplx_phase_f32(ch0_buf, phase_ch0, SAMPLE_NUM); arm_cmplx_phase_f32(ch1_buf, phase_ch1, SAMPLE_NUM); // 2. 寻基波峰（假设单频信号） uint16_t k_fund; find_fundamental_peak(mag_ch0, &k_fund); // 自定义函数，返回峰值bin索引 // 3. 计算初始相位差（单点） float delta_phi_raw = phase_ch0[k_fund] - phase_ch1[k_fund]; // 4. 解卷绕（此处简化为单点，实际应基于多帧历史） static float last_delta = 0.0f; float diff = delta_phi_raw - last_delta; if(diff > M_PI) delta_phi_raw -= 2.0f * M_PI; if(diff < -M_PI) delta_phi_raw += 2.0f * M_PI; last_delta = delta_phi_raw; *phase_diff_out = delta_phi_raw * 180.0f / M_PI; // 弧度转角度 }

注意：单帧解卷绕不可靠，工业级做法是维护一个长度为5~10帧的delta_phi_raw滑动窗口，对整个窗口做解卷绕后再取中位数，抗脉冲噪声能力极强。

5. 避坑指南：STM32 FFT项目中踩过的5个真实大坑

5.1 现象：FFT后频谱主瓣异常宽，两个100Hz/101Hz信号无法分辨

原因：ADC采样率不稳定。使用SysTick作为触发源，但SysTick被其他中断抢占，导致采样间隔抖动，等效于加了随机噪声窗。
解决：改用高级定时器（TIM1/TIM8）的TRGO事件作为ADC外部触发源，并将TIM1配置为主模式：Update事件触发，时基精度由APB2时钟分频决定，抖动<1ns。

5.2 现象：幅值测量结果随输入信号频率变化，50Hz时准确，100Hz时偏低15%

原因：未做窗函数能量补偿。Hanning窗CG=0.5，但代码中补偿系数写成1.0/N（应为4.0/N）。
解决：重新推导补偿系数，确认arm_cmplx_mag_f32()输出的是未归一化模值，最终补偿系数=2.0 * (1/CG) / N = 4.0/N。对N=1024，用0.00390625f（=4/1024）替代错误的0.0009765625f（=1/1024）。

5.3 现象：相位差测量值在±180°附近剧烈跳变，无法稳定跟踪

原因：未做相位解卷绕，且arm_cmplx_phase_f32()输出范围是[-π, π]，当真实相位差为179°时，下一帧变为-179°，差值达-358°。
解决：在calc_phase_diff()中加入状态记忆，比较当前帧与上一帧相位差，若跳变>150°，则加减360°修正。更鲁棒的做法是维护3帧滑动窗口，用中值滤波抑制跳变。

5.4 现象：FFT计算耗时超预期，1024点FFT占用CPU时间>8ms（理论应<2ms）

原因：arm_cfft_f32()要求输入数组长度为2的幂，但CMSIS-DSP的F1系列优化版本（ARM_MATH_CM3）未启用编译器内联优化，且未关闭浮点单元异常。
解决：

在arm_math.h前定义#define ARM_MATH_CM3；
在main.c开头添加#pragma GCC optimize ("O3")；
关闭FPU异常：FPU->FPCCR &= ~FPU_FPCCR_ASPEN_Msk;；
确保链接脚本中.bss段未被意外填充为0（会触发内存访问异常）。

5.5 现象：双通道相位差恒为固定值（如-45°），不随信号变化

原因：两路ADC未真正同步。虽然共用TIM1触发，但CH0和CH1的采样时间设置不同（如CH0=13.5周期，CH1=239.5周期），导致CH1比CH0晚采样数百ns，引入固定延迟。
解决：双通道必须使用完全相同的采样时间（ADC_SampleTime_xxx），且在ADC_RegularChannelConfig()中按顺序配置，确保扫描时序严格一致。

6. 实测验证与工程化技巧：如何用一块示波器+信号发生器30分钟验真结果

6.1 验证工具链搭建：低成本高精度标定方案

不用昂贵的矢量网络分析仪，用以下组合即可完成全指标验证：

信号源：普通DDS信号发生器（如RIGOL DG1022），输出纯净正弦波；
参考测量：数字示波器（如DS1054Z）开启FFT功能，作为黄金标准；
待测设备：你的STM32板，UART输出freq: xx.xx Hz, amp: yy.yyy V, phase: zz.z °；
校准负载：1%精度金属膜电阻+电容构成RC分压网络，产生已知相位差（如1kHz时-45°）。

提示：示波器FFT的窗函数、采样率、分辨率带宽（RBW）必须与STM32设置严格一致，否则对比无意义。

6.2 四步验证法：逐层剥离误差源

验证层级	操作	判定标准	典型问题定位
L1：ADC采样率校准	输入1kHz方波，用示波器测DMA缓冲区首尾时间差	实测采样率 = 1024 × (首尾时间差)⁻¹，误差应<0.1%	若误差>1%，检查TIM1时基配置、DMA传输延迟
L2：幅值线性度	输入0.1V~3.0Vpp正弦波（步进0.1V），记录STM32输出幅值	全量程内非线性度<±0.5%，斜率偏差<±0.2%	若低端偏差大，检查运放输入偏置电流；若高端饱和，检查ADC参考电压
L3：频率分辨率	输入99.9Hz/100.0Hz/100.1Hz信号，观察能否分辨	100Hz信号下，可分辨≤0.1Hz频差（对应1024点@100kHz采样）	若无法分辨，检查窗函数选择（矩形窗不行）、插值是否启用
L4：相位差精度	用RC网络产生-45°相位差（理论值），输入双通道	测量值在-45°±2°内，且10秒内波动<1°	若漂移大，检查电源纹波、PCB地线分割、运放共模抑制比

6.3 工程化技巧：把FFT模块封装成可复用的“黑盒”

为便于移植到不同项目，我将核心逻辑封装为三个接口函数，隐藏全部CMSIS-DSP细节：

// fft_engine.h typedef struct { float freq_hz; // 基波频率 float amp_v; // 基波幅值（V） float phase_deg; // 相位差（°） uint8_t valid_flag; // 1=测量有效，0=信噪比不足 } fft_result_t; // 初始化：传入采样率、缓冲区指针、FFT点数 void fft_engine_init(uint32_t sampling_freq, float *buf_ch0, float *buf_ch1, uint16_t len); // 启动一次测量（非阻塞，内部触发ADC） void fft_engine_start(void); // 获取结果（调用前确保fft_engine_start()已执行且DMA完成） uint8_t fft_engine_get_result(fft_result_t *result); // 使用示例 float ch0_buf[1024], ch1_buf[1024]; fft_result_t res; fft_engine_init(25000, ch0_buf, ch1_buf, 1024); // 25kHz采样，1024点 while(1) { fft_engine_start(); delay_ms(50); // 等待ADC+DMA+FFT完成（实测约42ms） if(fft_engine_get_result(&res)) { printf("F:%.2fHz A:%.3fV P:%.1f°\r\n", res.freq_hz, res.amp_v, res.phase_deg); } }

封装要点：