简介基于STM32F407的FFT频谱分析工程面向嵌入式开发者与信号处理工程师实现1Hz高分辨率频谱分析支持信号频率测量、谐波成分分析、总谐波失真度评估以及双波形相位差检测适用于音频工程、电力系统监测和通信调试等场景。资源共316个文件以C/H源码和MDK工程文件为主同时包含编译生成的axf/hex固件、PDF讲解文档以及LCD显示与UCOS相关代码压缩包大小3.71MB结构完整可直接打开工程参考或烧录测试。已有10448人学习下载在STM32 FFT应用领域具有一定参考价值。项目覆盖0~2KHz分析范围从底层ADC采样到FFT运算再到LCD显示均有对应实现包含FFT计算、谐波提取和相位差测量等关键模块适合需要开展实时频谱分析、失真度检测或频率测量开发的读者作为基础模板便于二次修改与移植。1. 用STM32F407做1Hz分辨率FFT测量这块板子完全扛得住坑都在细节里STM32F407上做FFT不稀奇稀奇的是把分辨率做到1Hz再同时输出信号频率、谐波、失真度和双通道相位差这在168MHz的MCU上并不是默认就能跑好的方案。很多人第一次接这个需求总以为“FFT是库函数调用就有结果”结果一上电发现频率读数飘、谐波幅度对不上、相位差总有固定偏移最后只能怪硬件不稳。实际上1Hz分辨率背后藏着采样时间窗、内存布局、窗函数和通道同步四道坎每一道都有明确的数学约束和调试手段。这篇文章把我自己做过的采集、变换、测算和标定流程完整拆开从采样参数怎么定到相位差为什么翻车按步骤复现一遍新手能跟着跑通老手能直接拿去核对边界和参数。2. 分辨率计算的物理账采样点数、采样率与时间窗怎么凑出1Hz2.1 为什么1Hz分辨率不是调参问题而是内存和时间窗问题FFT的频率分辨率由公式Δf Fs / N决定Fs是采样率N是FFT点数。想要分辨率等于1Hz只要让Fs和N的比值等于1就行比如Fs4096Hz、N4096或者Fs8192Hz、N8192。表面上看数学很简单但落到F407上要同时算清两笔账一笔是采样时间窗T N / Fs不管怎么组合1Hz分辨率都需要至少1秒的连续采样窗口另一笔是内存账FFT缓冲区的体积随着N翻倍。4096点单通道复数FFT需要32KB的float数组8192点直接翻到64KB再把窗函数、ADC DMA缓冲和双通道副本放进去192KB的SRAM很快就紧张了。很多人把采样率提到48kHz又想要1Hz分辨率那需要48000点这已经不是F407直接算的点了。FFT点数必须是2的幂所以采样率也最好选2的幂。常见组合就三种Fs1024、N1024奈奎斯特频率512Hz适合40次工频谐波以内的测量Fs4096、N4096奈奎斯特频率2048Hz测到音频没问题Fs8192、N8192奈奎斯特频率4096Hz内存占用量变高但高频成分更全。选型时还要注意一点1Hz分辨率的含义是“能区分间隔1Hz的两个等幅信号”并不是说测量读数只能显示整数。配合后面讲的插值方法单一频率的读数精度可以远超1Hz。2.2 采样率-点数组合三套配置内存与测量范围选个平衡下面的表格给出一组能直接用于F407的组合内存估算包含双通道FFT缓冲、窗函数和ADC原始缓冲区实际工程还要留出栈、堆和显示缓冲别卡得太死。采样率FsFFT点数N频率分辨率最高分析频率时间窗双通道预估内存1024 Hz10241.0 Hz512 Hz1 s约 20 KB4096 Hz40961.0 Hz2048 Hz1 s约 80 KB8192 Hz81921.0 Hz4096 Hz1 s约 160 KB我一般优先选4096点这套理由是内存余量充足还能放心开LCD缓冲或者跑个简单的UI。如果被测信号是工频50Hz2次谐波100Hz40次谐波2000Hz2048Hz的奈奎斯特频率刚好够用但50Hz的41次谐波就会折返进频谱里造成假的“高次谐波”所以真正做谐波分析时不要只看表上分辨率。8192点那套适合需要看3kHz以上分量但内存又有余量的场景代价是FFT计算耗时接近翻倍F407仍然能跑只是刷新率会掉到每秒几帧。采样率还要配合抗混叠滤波器ADC输入前级必须有一个截止频率在奈奎斯特频率附近的RC或运放低通否则高频折叠进来会让谐波和失真度结果完全失真。2.3 CMSIS-DSP库的实数与复数FFT取舍F407带FPU和DSP指令算FFT靠的是芯片原厂提供的DSP库里面有两类常用接口实数FFTarm_rfft_fast_f32和复数FFTarm_cfft_f32。实数FFT输入是N个float内存占用更小速度更快但输出格式是打包的需要按约定位置“拆”出幅度谱拆错一次结果就全乱。复数FFT输入是N个复数也就是2N个float实部和虚部交错存放填数据时把每个采样点放在偶下标、虚部填0就行输出格式直观直接用arm_cmplx_mag_f32求幅值。对于4096点复数FFT在F407上耗时大约几毫秒完全不影响1秒一次的数据刷新所以我更推荐用复数接口代码容易读也容易排错。直接定义数组并调用变换的核心代码是这样#include arm_math.h #define FFT_SIZE 4096u #define SAMPLE_RATE 4096.0f static float32_t fft_buf[FFT_SIZE * 2]; /* 复数输入缓冲区 */ static float32_t mag_buf[FFT_SIZE / 2]; /* 单边幅值谱 */ static float32_t win_buf[FFT_SIZE]; /* 窗函数 */ void fft_init(void) { /* 周期Hann窗n以FFT_SIZE为周期 */ for (uint32_t n 0; n FFT_SIZE; n) { win_buf[n] 0.5f * (1.0f - cosf(2.0f * PI * n / FFT_SIZE)); } /* 只需初始化一次arm_cfft_f32内部会建旋转因子表 */ arm_cfft_init_f32(arm_cfft_sR_f32_len4096, FFT_SIZE); }这段代码里的窗口函数要特别说明一下很多人用0.5*(1-cos(2PIn/(N-1)))那是信号处理教科书里的对称窗而FFT实际处理的是一帧“周期延拓”的信号所以这里用周期窗更合适让窗函数首尾连续基波落在整数谱线上时泄漏接近于零。arm_cfft_init_f32里的旋转因子表是DC表static数组占一点SRAM但能避免每次变换重复建表。选复数接口的代价是内存比实数接口多一个虚部数组但对4096点这个量级80KB之内能塞下仍然在F407的SRAM承受范围内。3. 从ADC采样到频率输出4096点FFT的完整测频流程3.1 采样触发与DMA缓冲先保证数据不丢频率测量的准确性首先不取决于FFT而取决于采样率是否稳定。F407的ADC建议用定时器触发让每个采样点之间的时间间隔严格等于1/Fs不要用HAL_ADC_Start加延时轮询的方式那会被中断和任务调度打乱产生抖动。用定时器触发配合DMA循环模式是常见做法定时器输出PWM或TRGO信号给ADC外部触发ADC完成转换后由DMA把结果搬进内存整个过程不占CPU。采样率由定时器的PSC和ARR决定计算公式是Fs TIM_CLK / (PSC1) / (ARR1)。实现时要注意TIM_CLK不是168MHz通常挂在APB1上最高84MHz除以分频和重装载值后才是真正的采样率。DMA配置用循环模式每当传输完成一轮4096点就触发一次转换完成回调。下面代码是HAL风格的大致结构具体引脚号和DMA通道需要按实际板子接线调整#define ADC_SAMPLES 4096u static uint16_t adc_buf[ADC_SAMPLES]; void sample_init(TIM_HandleTypeDef* htim) { /* ADC1启用外部触发触发源为定时器TRGO */ ADC_ChannelConfTypeDef sConfig {0}; sConfig.Channel ADC_CHANNEL_4; sConfig.SamplingTime ADC_SAMPLETIME_84CYCLES; HAL_ADC_ConfigChannel(hadc1, sConfig); /* DMA循环模式外设到内存数据宽度半字 */ hdma_adc.Init.Direction DMA_PERIPH_TO_MEMORY; hdma_adc.Init.Mode DMA_CIRCULAR; HAL_DMA_Init(hdma_adc); __HAL_LINKDMA(hadc1, DMA_Handle, hdma_adc); HAL_TIM_Base_Start(htim); HAL_ADC_Start_DMA(hadc1, (uint32_t *)adc_buf, ADC_SAMPLES); } void HAL_ADC_ConvCpltCallback(ADC_HandleTypeDef *hadc) { if (hadc-Instance ADC1) { process_frame(); /* 一轮4096点采集完成 */ } }这里的SamplingTime不是采样周期而是ADC单次转换的采样保持时间84个ADC时钟周期是为了降低信号源内阻带来的采样误差。真正决定Fs的是定时器所以线上调试时必须回读定时器的实际频率而不是想当然地把SAMPLE_RATE写成4096.0f。如果外部晶振不准确或者定时器分频有偏差测频结果会出现比例性偏移比如真实采样率是4080Hz你按4096算所有频率读数都会偏大0.4%。我习惯在调试串口打一帧采样数据用示波器测一下相邻采样点之间的间隔确认偏差小于0.1%再继续往下做。3.2 加窗与对数插值让频率输出稳定在1Hz档位ADC拿到的原始数据在送入FFT之前先要减去直流偏置再乘窗函数。减去2048这个中点值可以防止直流分量占据了0号bin把动态范围吃掉。接着把数据归一化到-1到1之间然后乘上Hann窗。为什么要加窗只要被测信号的频率不是采样率除以点数的整数倍基波能量就会泄漏到相邻bin峰值位置不再精确。Hann窗把泄漏压低到主瓣附近但代价是主瓣变宽峰值位置必须靠插值恢复。频率估计用对数幅值抛物线插值比线性插值稳定得多。原因是Hann窗的频谱主瓣在对数坐标下更接近抛物线插值残差更小。具体步骤先找到幅度谱中最大值所在的bin取这个bin左右各一个bin的对数幅值计算偏移量d。代码如下uint32_t peak 0; float max_mag 0.0f; for (uint32_t i 1; i FFT_SIZE / 2 - 1; i) { float m mag_buf[i]; if (m max_mag) { max_mag m; peak i; } } if (peak 2) { peak 2; /* 远离DC避免低频噪声干扰 */ } float p logf(mag_buf[peak - 1]); float c logf(mag_buf[peak]); float n logf(mag_buf[peak 1]); float d 0.5f * (p - n) / (p - 2.0f * c n); float freq ((float)peak d) * SAMPLE_RATE / FFT_SIZE;这个d的理论范围是-0.5到0.5之间如果算出来的d超过±1说明信号太弱或旁边有个等幅干扰这时freq结果不要信宁可标记为无效帧也不要把跳变值显示出来。为了让读数稳定我还会对连续几帧的freq结果做一次滑动中值滤波只取中间值对抑制偶发的噪声毛刺很有效。这里的1Hz分辨率说的是能分辨两个相邻频率分量而插值让单一频率读数的精度达到了0.01Hz量级两者不要混为一谈。3.3 主循环里跑FFT与峰值搜索一块代码说清整个链路把前面的环节串起来整个处理流程可以放在DMA回调里直接调用也可以把回调里只置标志位主循环里再处理。我一般倾向于回调置标志、主循环处理这样FFT计算不会阻塞DMA采样也不容易在中断里踩堆栈。完整的一帧处理函数大致如下static void process_frame(void) { float vin 3.3f / 4096.0f; /* 12位ADC满量程3.3V */ for (uint32_t i 0; i FFT_SIZE; i) { float x (float)adc_buf[i]; x (x - 2048.0f) * vin; /* 去除直流偏置得到实际电压 */ x * win_buf[i]; /* 施加Hann窗 */ fft_buf[2 * i] x; fft_buf[2 * i 1] 0.0f; } arm_cfft_f32(arm_cfft_sR_f32_len4096, fft_buf, 0, 1); arm_cmplx_mag_f32(fft_buf, mag_buf, FFT_SIZE / 2); float freq estimate_fundamental_freq(); float thd calculate_thd(freq); send_result(freq, thd); }这段代码把整条链路的顺序固定下来了先预处理采样数据再做复数FFT接着求幅值谱最后估频和算谐波。arm_cmplx_mag_f32的输出是双单边谱下标0对应DC下标FFT_SIZE/2对应奈奎斯特频率这里只取前半部分合理。注意预处理中减去的2048是按12位ADC中值算的如果ADC参考电压不是3.3V要换成VREF的实测值。整个处理过程中堆栈开销最大的就是arm_cfft_f32它内部会使用临时缓冲区所以在中断里调用时务必确认栈空间足够否则会出现随机进HardFault的玄学问题。4. 谐波与失真度在基波周围做“局部寻峰”别拿全谱段硬算4.1 谐波搜索的排除带与窗口宽度设定谐波分析不是把FFT幅值谱里每个整数倍频点直接取值就完事因为基波频率本身有微小误差实际谐波不可能正好落在整数bin上。正确做法是先拿到基波频率f0再对每个谐波次数k在k*f0附近开一个局部搜索窗口窗口里找局部峰值。窗口宽度怎么设Hann窗主瓣约4个bin也就是4Hz加上参数估计误差我一般用±4Hz作为搜索半径也就是±4个bin。如果信噪比很差或者信号有频率抖动可以放宽到±6Hz但太宽会把邻近的窄带噪声当成谐波。排除带也很关键。当基波非常强时其泄漏会蔓延到2倍频附近造成二次谐波“虚高”。Hann窗在远离主瓣处衰减比较快但如果基波比二次谐波大40dB以上泄漏仍然可能压过真实谐波。常见做法是在谐波搜索时至少跳开基波位置左右各4个bin同时把参与THD计算的谐波幅度设置一个门限低于门限的直接置零。这里的门限我建议用相对基波的比值比如基波幅值的0.001倍对应-60dB也在ADC的12位分辨能力之内。门限设太低噪声全部算成谐波失真度虚高门限设太高真实谐波被丢掉失真度偏小。第一次调试时可以把所有搜索到的分量都打出来看一眼再定门限。4.2 THD计算与噪声门限为什么算出来的失真度总偏大THD的计算公式是各次谐波幅值平方和开根号再除以基波幅值通常算到10次谐波。IMD和噪声不参与THD但这在测量中是很难严格分开的。用4096点、1Hz分辨率的FFT噪声底在理想条件下大约-70dB到-80dB之间但ADC本身的量化噪声、前级放大器的噪声、还有采样抖动都会把底噪抬起来。如果门限设在-60dB而噪声底已经被抬到-55dB那么每个谐波位置附近都可能抓到噪声峰THD会明显偏大。解决方向有两种一是降低门限但配合幅度一致性检查单独一个点超过门限不算要求连续两个搜索点都超过才算二是对多帧求平均幅值谱平滑后再算THD代价是响应变慢。计算代码里要特别注意谐波搜索范围不要越界同时用int32_t处理bin下标避免unsigned减法回绕#define MAX_HARMONIC 10u #define GATE_RATIO 0.001f /* 相对基波的-60dB门限 */ static float calculate_thd(float fundamental_freq) { float h1 mag_buf[peak_bin]; if (h1 1e-6f) { return 0.0f; } float acc 0.0f; float gate h1 * GATE_RATIO; for (uint32_t h 2; h MAX_HARMONIC; h) { float target_bin fundamental_freq * (float)h * FFT_SIZE / SAMPLE_RATE; int32_t k0 (int32_t)(target_bin 0.5f); int32_t start k0 - 4; int32_t end k0 4; float local_max 0.0f; for (int32_t j start; j end; j) { if (j 1 || j (int32_t)(FFT_SIZE / 2)) { continue; } float m mag_buf[j]; if (m local_max) { local_max m; } } if (local_max gate) { acc local_max * local_max; } } return sqrtf(acc) / h1; }这段代码里h1直接用峰值bin的幅值没有做窗增益补偿但THD是比值计算窗增益在分子分母里抵消所以不影响结果。如果你顺带要显示各次谐波的绝对幅度就必须除以Hann窗的相干增益0.5也就是把幅值乘2。另外当基波频率低于100Hz时10次谐波也才1kHz还在奈奎斯特范围内但搜索窗口下标的起点要避开DC附近排除带。我遇到过测一个50Hz信号二次谐波搜索结果里包含了基波泄漏的尾巴门限一放THD直接跳到5%实际上信号失真只有0.5%这就是排除带没留够。4.3 幅值恢复与校正窗函数增益、ADC基准与分压系数如果仪器要输出绝对的谐波电压而非百分比那么幅值链路上有三处必须校准窗函数增益、ADC参考电压、前端分压比。Hann窗的相干增益是0.5意味着对一个纯正弦波FFT峰值bin的幅值只有真实幅值乘以0.5再乘以N/2的缩放关系直接用arm_cmplx_mag_f32的原始结果会偏低一倍。恢复时需要乘2再乘上采样数据的归一化系数和分压系数。校准流程我一般是这样做的用一个已知幅值的正弦信号接进去频率选在整bin位置比如Fs4096时选512Hz这样128点一个周期正好落在128号bin上。测量到的峰值幅值记为m已知有效值或峰值记为A那么总增益系数就是A/m。调试时把这组系数存进flash开机时读取比理论推算省心很多。分压系数如果是一个固定电阻分压网络还要考虑电阻精度和温度漂移工业级方案会预留软件校准接口让产线用标准源标定一次。还有一个常见的坑ADC的参考电压不是3.3V整。很多板子用LDO输出3.3V误差可能到±1%对频率不影响但谐波幅度会跟着偏所以幅值校准时要用万用表实测VREF代入计算不要直接写3.3。校准之后失真度百分比因为比值关系不受VREF影响但绝对谐波电压会。做产品方案时至少保留这组校准系数的软件入口后续产线调起来会省很多事。5. 相位差测量与同步采集避坑偏角、跳变、反相的三类现场5.1 双ADC同步采样相位差测量的第一步测两个波形的相位差最核心的前提是两路采样必须严格同时而不是先后采。先后采会有固定的时间偏差折算成相位会随信号频率变化频率越高误差越大。F407上有三个ADC可以配置ADC1和ADC2为规则同步模式由同一个定时器触发同时启动转换。配置时重点在TwoSamplingDelay它决定两路ADC启动之间的间隔设置为1个ADC周期可以把通道间的时间偏移压到很小。下面是用HAL库配置双ADC同步模式的核心片段ADC_MultiModeTypeDef multimode {0}; multimode.Mode ADC_DUALMODE_REGSIMULT; /* 规则通道同步模式 */ multimode.DMAAccessMode ADC_DMAACCESSMODE_DISABLED; multimode.TwoSamplingDelay ADC_TWOSAMPLINGDELAY_1CYCLE; HAL_ADCEx_MultiModeConfigChannel(hadc1, multimode);这段配置只负责把ADC1和ADC2“绑定”成同步触发实际信号仍然分别走各自的通道和DMA缓冲区。同步模式建立后ADC1和ADC2在同一个定时器边沿开始转换DMA搬运各自的采样结果到两个独立的数组这样后续FFT处理时两路数据的起点就是对齐的。不要在回调里把两路数据分别拿去做FFT再比较相位只要数据起点是对齐的就没问题。还有一点容易被忽略如果两个输入通道用了不同的抗混叠滤波器比如通道A是二阶RC低通通道B是一阶RC低通那么在分析频率附近会引入不同的相位滞后相位差结果会随频率变化。做相位测量时两个输入通道的前级滤波器必须做到一样的截止频率和阶数否则软件怎么校准都补不齐这个频变误差。5.2 基波相位的提取与相位差计算公式两个通道各做一次FFT之后从各自幅值谱中找到基波峰值bin再取出该bin对应的复数的实部和虚部用atan2f求相位。受窗函数和栅栏效应影响单通道的绝对相位和实际信号过零点并不严格对应但两路信号频率相同、窗函数相同、泄漏模式相同相位差值仍然是可信的。相位差的范围是-π到π直接相减后要归一到这个区间否则会看到0到2π之间跳变的读数。float reA fft_a[2 * peak_a]; float imA fft_a[2 * peak_a 1]; float reB fft_b[2 * peak_b]; float imB fft_b[2 * peak_b 1]; float phase_a atan2f(imA, reA); float phase_b atan2f(imB, reB); float diff phase_a - phase_b; /* 归一到 [-PI, PI] */ while (diff PI) diff - 2.0f * PI; while (diff -PI) diff 2.0f * PI;这段代码里peak_a和peak_b分别是两路基波的峰值bin。正常工作时两个通道应该用同一个bin因为频率理论上一致但实际会有微小差别所以各自取峰值更稳。如果两路信号幅度相差很大比如一路饱和、一路只有几毫伏那么弱信号通道的相位会被噪声淹没diff会乱跳。我一般会加一个有效状态判断两路基波幅值都大于设定门限时才更新相位差输出否则保持上一帧有效值这样仪表不会在弱信号时显示无意义的跳变数字。5.3 现象排查相位差为什么总有一个固定偏角第一个高频踩坑现象是相位差读数始终比真实值偏大或偏小一个固定角度而且频率越高差得越多。原因是两路没有真正同步采集比如用ADC1的注入组轮询采A路再切到B路中间隔了几个微秒。解决办法是改用双ADC同步模式同步之后残留偏差很小。如果同步模式下仍然有固定偏角就可能是两个输入前级的滤波电容参数不一致这时候只能做软件补偿接入同相信号记录初始偏差phase_offset之后每帧直接减掉。补偿代码只有一行diff - phase_offset; if (diff PI) diff - 2.0f * PI; if (diff -PI) diff 2.0f * PI;第二个常见现象是相位差读数在接近±180°时来回跳一会儿显示179°一会儿显示-179°。这不是电路问题而是相位差本来就接近边界atan2f的一个小抖动就会让结果越过±π。解决方法是不要对瞬时帧直接输出而是把连续几帧的diff都归到同一方向后再平均或者记录上一帧的符号发现跳变幅度大于π时按“越过边界”处理。更简单的方法是在后端加一个小的低通滤波但要注意接近边界时不能简单平均必须先把角度展开成连续时间序列再滤波。第三个现象是反相显示180°而且很稳定。常见原因是某个通道的输入极性反了或者前级差分放大器接反也有可能是ADC通道引脚的采样保持电压反了。处理办法是先让两个输入端口短接同一信号源看相位差是否约等于0°如果显示180°先查硬件极性再查软件里两路缓冲区的通道映射关系是否颠倒了。校准逻辑里要保留这类自检模式产线调试时能省很多排查时间。6. 用信号发生器给这套1Hz分辨率方案“验明正身”验证这套系统不能只看内部计算最直接的办法是用一台信号发生器输出标准正弦波把频率设成整数和非整数两种比如50Hz、50.3Hz、997Hz再用高精度频率计做参考。判定标准我一般按0.05Hz以内算合格因为信号发生器本身的频率精度也在0.01Hz量级插值估频的误差主要来自ADC噪声和窗函数0.05Hz是个合理余量。相位差验证则用双通道信号发生器同时输出同频信号一个通道固定0°另一个通道从0°扫到300°每30°记录一次误差在1°以内就可以认为相位差测量可用。下面是一张典型验证记录表参考频率FFT读数误差参考相位差实测相位差50.00 Hz50.002 Hz0.002 Hz0.0°0.08°50.30 Hz50.296 Hz-0.004 Hz90.0°90.15°997.00 Hz997.03 Hz0.03 Hz300.0°299.92°验证无头绪时可以从一个桌面技巧开始不接任何仪表直接把一个已知频率的方波信号接进去。方波的谐波非常丰富如果谐波列表能同时看到3次、5次、7次谐波且幅值比接近1/3、1/5、1/7说明整条链路是通的如果看到的谐波位置不对或者THD算出来离理论值48%差太多问题大概率出在采样率配置或频谱峰值搜索上。这个技巧对新手尤其好用不用高精度仪表也能快速判断方案是否落地。我现在做这类测量仪器有一个固定习惯先信串口打印的第一帧原始数据再信FFT输出最后才信校准后的显示值。任何一个环节的数据“看起来合理”都不是结论必须用已知信号做交叉验证。FFT方案最大的优势是没有模拟鉴相电路的温度漂移最大的风险是数字链路里某个细节错了却表现得很正常。希望帮到你。本文还有配套的精品资源点击获取