简介本资源是一套面向嵌入式开发初学者与STM32项目实践者的ADC采样软件滤波实战代码包聚焦传感器数据预处理中的噪声抑制问题覆盖工业采集、智能仪表、环境监测等典型应用场景。压缩包含83个文件以37个头文件.h定义寄存器配置与函数接口、35个源文件.c实现十大经典滤波算法限幅、中位值、算术平均、滑动平均、防脉冲干扰平均、一阶滞后、加权递推、消抖及组合滤波等辅以启动文件、外设驱动、串口调试模块及Keil工程配置文件.uvprojx/.uvoptx结构完整开箱即用。资源大小334KB轻量高效适配STM32F10x系列标准库开发环境。已有1629人学习下载提供从原理说明、代码实现到工程集成的全链路参考特别适合理解滤波算法底层逻辑、对比不同方法滤波效果及快速移植至实际项目。1. 十种滤波算法不是“堆代码”而是解决ADC采样中真实干扰的工程选择在STM32项目里ADC读出来的电压值跳变剧烈、数值毛刺多、温度传感器读数忽高忽低——这不是硬件坏了大概率是原始采样数据没过滤。但直接套用“网上搜到的滤波函数”常会翻车限幅滤波设错阈值导致动态响应迟钝中位值滤波用5点却卡在中断里超时滑动平均用16点后发现RAM爆了。这份源码包的价值不在于它打包了十种算法名称而在于它把每种滤波法都实现在真实STM32F10x平台Keil MDK工程结构标准外设库调用链中且全部通过main.c中统一ADC采集流程驱动。它面向的是需要快速验证滤波效果的嵌入式工程师你不必从零写ADC_Init()也不用纠结DMA_BufferSize怎么配所有滤波函数输入都是uint16_t raw_value输出都是uint16_t filtered_value接口干净可插拔替换。尤其适合工业传感器信号调理、电机电流采样抗干扰、电池电压监测等对实时性与稳定性双重要求的场景。2. 滤波算法选型逻辑与STM32资源约束下的实现边界2.1 为什么必须按干扰类型选算法而不是“越复杂越好”ADC采样干扰本质分三类脉冲型单次尖峰、周期型工频耦合、缓变型温漂/电源纹波。不同算法对这三类干扰的抑制能力差异极大且直接影响CPU占用、内存开销和响应延迟。例如限幅滤波法仅需2个uint16_t变量上次值当前值无数组、无排序、无循环计算执行时间稳定在3~5μs72MHz主频但对50Hz工频干扰完全无效中位值滤波法需N点缓存冒泡排序N5时约需80μsN11时升至320μs但能彻底剔除单次3倍标准差的脉冲干扰一阶滞后滤波法本质是IIR低通output α * new (1-α) * oldα0.25时等效3dB截止频率约120Hz对工频干扰有天然衰减但存在相位滞后不适合闭环控制中的电流环采样。提示本工程中所有滤波函数均定义在adc.c入口参数统一为uint16_t adc_value返回uint16_t。函数内部不调用HAL_Delay()或printf()确保可在ADC中断服务程序ISR中安全调用。若需在中断中使用中位值滤波务必确认N值对应的最大执行时间小于ADC采样周期——例如采样周期1ms时N7约180μs是安全上限。2.2 STM32F10x资源约束下的关键实现细节本工程基于STM32F10x标准外设库FWLibADC配置走ADC1通道PA0采用规则通道单次转换模式非扫描/非DMA触发源为软件触发ADC_SoftwareStartConvCmd(ADC1, ENABLE)。这种配置决定了滤波函数的调用时机每次ADC_GetConversionValue(ADC1)读取后立即调用滤波函数。关键约束点如下2.2.1 RAM占用与数组长度硬限制滤波算法缓存数组长度单元素类型总RAM占用典型适用N值中断安全限幅滤波法1uint16_t2B—✅中位值滤波法Nuint16_t2×N B3,5,7,9⚠️N≤7算术平均滤波法Nuint16_t2×N B4,8,16✅递推平均滤波法1uint32_t4B—✅中位值平均滤波法N×Muint16_t2×N×M BN3,M5→30B❌慎用注意adc.c中#define FILTER_BUF_SIZE 10全局定义缓冲区大小所有依赖数组的滤波函数如MidFilter()、AveFilter()均以此宏为界。若修改N值必须同步调整该宏并检查.map文件中RW Data段是否溢出。Keil编译后查看OBJ\ADC.map中ADC.o的Data节确认filter_buf未跨页STM32F10x默认SRAM起始0x20000000大小20KB。2.2.2 时间敏感操作的汇编级保障以限幅滤波法为例其核心逻辑在adc.c第127行uint16_t LimitFilter(uint16_t new_value) { static uint16_t last_value 0; uint16_t diff; if(new_value last_value) diff new_value - last_value; else diff last_value - new_value; if(diff LIMIT_VALUE) { last_value new_value; return new_value; } else { return last_value; } }此处LIMIT_VALUE定义为#define LIMIT_VALUE 200对应ADC 12位满量程4095的约5%但关键在diff计算——避免使用abs()函数。标准库abs()引入stdlib.h且为函数调用增加压栈开销而手动分支判断仅需2条ARM Thumb指令CMPITT实测在Keil v5.37下编译为4周期指令比abs()快3倍。此细节在ADC中断高频调用时如10kHz采样直接决定系统能否稳定运行。2.2.3 外设库调用链的隐含开销本工程ADC_Init()配置中ADC_InitStructure.ADC_ScanConvMode DISABLE;关闭扫描模式ADC_InitStructure.ADC_ContinuousConvMode DISABLE;关闭连续模式这意味着每次采样需手动触发ADC_SoftwareStartConvCmd(ADC1, ENABLE); // 启动转换 while(!ADC_GetFlagStatus(ADC1, ADC_FLAG_EOC)); // 等待转换结束 raw ADC_GetConversionValue(ADC1); // 读取结果 filtered MidFilter(raw); // 调用滤波此处while循环是忙等待若ADC时钟PCLK2为72MHzADCCLK经RCC_ADCCLKConfig(RCC_PCLK2_Div6)分频为12MHz则单次转换时间≈1.5μs12.5个ADCCLK周期while循环耗时可控。但若误开启连续模式EOC标志将频繁置位while可能因中断优先级问题导致漏判——本工程NVIC_Configuration()中ADC_IRQChannel抢占优先级设为1确保能打断其他外设中断。3. 十大滤波算法在STM32上的逐个实现与参数调优指南3.1 限幅滤波法对抗单次脉冲干扰的轻量级方案3.1.1 核心原理与阈值设定方法限幅滤波本质是“变化率门限”公式为|new_value - last_value| ≤ A其中A非固定值需根据传感器特性设定电位器模拟量A ≈ 满量程5%~10%如12位ADC取200~400温度传感器DS18B20ADCA ≈ ±2℃对应ADC值NTC阻值变化非线性需查表换算电流采样0.1Ω电阻A ≈ 额定电流10%对应mV值经放大后ADC值本工程adc.h中#define LIMIT_VALUE 200适用于多数0~3.3V直连传感器若用于±10V输入经运放衰减需重算A (20V / 3.3V) × 200 ≈ 1212。3.1.2 Keil调试验证步骤在main.c中定位while(1)循环内ADC采样段将filtered_val LimitFilter(raw_val);替换为filtered_val raw_val;直通模式使用逻辑分析仪抓取PA0引脚波形注入10μs宽、2V尖峰干扰用信号发生器MOSFET开关观察串口打印的raw_val序列确认尖峰存在如4095→120→4095恢复LimitFilter()调用观察打印值是否稳定在120尖峰被抑制修改LIMIT_VALUE为50再注入相同干扰观察是否出现“过度限幅”正常缓慢变化也被截断。提示若发现滤波后数据响应迟钝说明LIMIT_VALUE过小应结合传感器最大允许变化率计算。例如电机转速反馈ADC每10ms采样一次机械惯性决定转速变化率≤50rpm/s则对应ADC值变化率≤(50/1000)×10ms×KK为rpm-to-ADC系数。3.2 中位值滤波法剔除随机脉冲的可靠选择3.2.1 排序算法优化与N值选择依据中位值滤波需对N个采样值排序取中值。本工程MidFilter()采用冒泡排序优化版adc.c第189行// 仅排序前(N1)/2轮因只需中位值无需全序 for(i0; i(FILTER_NUM1)/2; i) { for(j0; jFILTER_NUM-1-i; j) { if(buf[j] buf[j1]) { temp buf[j]; buf[j] buf[j1]; buf[j1] temp; } } } return buf[FILTER_NUM/2]; // 返回中位索引值FILTER_NUM定义为#define FILTER_NUM 5则仅执行3轮冒泡而非4轮减少比较次数30%。当FILTER_NUM7时理论比较次数从42次降至28次。N值选择原则N3可剔除1个脉冲RAM占用6B执行时间≈60μsN5可剔除2个脉冲RAM占用10B执行时间≈150μsN7可剔除3个脉冲RAM占用14B执行时间≈280μs注意N必须为奇数偶数N如4导致中位值在buf[1]与buf[2]间线性插值本工程未实现插值逻辑强制FILTER_NUM为奇数宏定义。3.2.2 抗干扰能力量化测试在main.c中构造干扰测试序列// 模拟含脉冲的采样流正常值3000±50叠加2次尖峰4095 uint16_t test_seq[] {2980,3020,4095,2990,3010,4095,2970,3030,2990,3010}; // 调用MidFilter()处理预期输出3000左右剔除两个4095编译后通过J-Link RTT Viewer实时打印test_seq各点滤波结果验证中位值是否稳定在3000±20范围内。若出现4095残留说明FILTER_NUM过小或排序逻辑错误。3.3 算术平均滤波法与递推平均滤波法的资源权衡3.3.1 算术平均滤波法简单但吃内存AveFilter()函数adc.c第245行使用静态数组缓存N点static uint16_t filter_buf[FILTER_NUM]; static uint8_t index 0; uint32_t sum 0; filter_buf[index] new_value; if(index FILTER_NUM) index 0; for(i0; iFILTER_NUM; i) sum filter_buf[i]; return (uint16_t)(sum / FILTER_NUM);致命陷阱sum为uint32_t但FILTER_NUM16时最大sum16×409565520未超uint16_t上限但若误用uint16_t sumsum filter_buf[i]将溢出。本工程已修正但用户自定义N值时必须验算N×4095 65535→N ≤ 15。3.3.2 递推平均滤波法内存友好型替代方案ImproveAveFilter()adc.c第278行采用IIR结构static uint32_t sum 0; sum sum - (sum FILTER_SHIFT) new_value; // 等效α1/(2^SHIFT) return (uint16_t)(sum FILTER_SHIFT);FILTER_SHIFT定义为#define FILTER_SHIFT 4即α1/16等效16点算术平均但RAM仅占4Bsum变量。其频响特性为一阶低通3dB带宽f_c f_s / (2π×15)f_s为采样率对1kHz采样f_c≈10.6Hz适合抑制高频噪声但保留低频信号。提示两种平均法不可混用。算术平均对缓变信号有相位滞后N/2个采样周期递推平均滞后更小但存在稳态误差。若需精确测量直流分量优先选算术平均若需快速跟踪变化选递推平均并增大FILTER_SHIFT降低α。3.4 高阶组合滤波法的工程落地要点3.4.1 中位值平均滤波法防脉冲干扰平均滤波法此算法先取N组数据每组M点做中位值再对N个中位值求平均。本工程MidAveFilter()adc.c第312行中N3M5// 3组每组5点中位值 for(i0; i3; i) { for(j0; j5; j) group_buf[j] GetADCValue(); // 实际中应加延时防耦合 mid_val[i] MidFilter(group_buf); // 复用中位值函数 } return AveFilter(mid_val); // 复用算术平均函数关键约束GetADCValue()需保证组内5点采样间隔足够长干扰周期否则脉冲可能同时污染整组。例如工频干扰50Hz组内采样间隔应20ms。本工程未加延时实际使用时需在for(j0; j5; j)内插入Delay_ms(25)。3.4.2 限幅消抖滤波法应对机械触点抖动专为按键、继电器反馈等设计结合限幅抑制尖峰与消抖延时确认#define DEBOUNCE_TIME 20 // ms static uint16_t last_valid 0; static uint32_t last_time 0; if(abs(new_value - last_valid) LIMIT_VALUE) { if(GetSysTimeMs() - last_time DEBOUNCE_TIME) { last_valid new_value; last_time GetSysTimeMs(); } return last_valid; } else { last_time GetSysTimeMs(); // 重置消抖计时 return last_valid; }GetSysTimeMs()需由SysTick中断维护本工程delay.c中SysTick_Handler()已实现毫秒计数器uwTick可直接调用。4. 滤波效果验证与ADC采样周期协同调优实战4.1 基于串口实时波形的滤波效果可视化本工程usart.c已配置115200bps串口main.c中通过printf(RAW:%d FIL:%d\r\n, raw, filtered)输出。为直观对比需构建测试环境4.1.1 硬件信号注入方案信号源任意波形发生器输出100Hz正弦波0~3.3V叠加1kHz方波干扰幅值0.5V连接信号源→PA0GND共地工具PC端使用XCOM或SSCOM串口助手设置“自动换行”“时间戳”保存日志为CSV数据处理Excel导入CSV作RAW与FIL列折线图观察限幅滤波方波干扰是否被削平正弦波峰值是否失真中位值滤波100Hz正弦波相位是否滞后理论滞后(N-1)/2个采样点递推平均正弦波幅度是否衰减计算20×log10(幅值比)得衰减量4.1.2 关键指标量化表算法干扰抑制比方波正弦波衰减100Hz相位滞后100Hz最大安全采样率限幅滤波法20dB0dB0°100kHz中位值滤波法(N5)30dB0.5dB20ms5kHz递推平均(α1/8)15dB-1.2dB12ms10kHz注意相位滞后影响闭环控制稳定性。若用于PID电流环采样滞后100μs可能导致系统振荡此时必须选限幅或一阶滞后滤波。4.2 ADC采样周期与滤波算法的协同设计ADC采样周期T_s由三部分构成T_s T_setup T_conv T_filter其中T_setup为采样时间ADC_SampleTime_XXX配置T_conv为转换时间固定12.5周期T_filter为滤波函数执行时间。4.2.1 STM32F10x ADC时序计算实例PCLK2 72MHzRCC_ADCCLKConfig(RCC_PCLK2_Div6)→ADCCLK 12MHzADC_SampleTime_239Cycles最长采样时间→T_setup 239 / 12MHz 19.9μsT_conv 12.5 / 12MHz 1.04μsT_filter查表中位值N5约150μs则最小T_s 19.9 1.04 150 ≈ 171μs→ 最大采样率≈5.8kHz。若需求10kHz采样则必须缩短T_setup至ADC_SampleTime_13Cycles1.08μs改用LimitFilter()5μs替代中位值或启用DMA传输避免CPU等待。本工程未启用DMA故T_filter是瓶颈。若需高频采样应在ADC_Init()中开启DMA并将滤波移至DMA传输完成中断中执行。4.2.2 滤波参数在线调整机制为适配不同传感器工程预留参数修改接口。在adc.h中// 可通过串口命令动态修改需扩展usart.c解析逻辑 extern uint16_t LIMIT_VALUE; // 限幅阈值 extern uint8_t FILTER_NUM; // 中位值点数 extern uint8_t FILTER_SHIFT; // 递推平均移位数实际项目中可添加AT指令解析ATLIMIT300→LIMIT_VALUE 300ATNUM7→FILTER_NUM 7ATSHIFT5→FILTER_SHIFT 5修改后需调用ResetFilterBuffer()清空历史缓存否则新参数首周期输出异常。4.3 常见故障排查与性能瓶颈突破4.3.1 “ADC值全为0或满幅”的硬件层检查清单现象检查项操作指引raw_value0PA0是否悬空万用表测PA0对地电压应为0~3.3V若为0V检查上拉/下拉电阻及信号源连接raw_value4095PA0是否短路到VDD断电测PA0对VDD电阻应100kΩ若接近0Ω检查PCB短路或芯片击穿raw_value乱跳ADC参考电压是否稳定测VREF引脚通常接VDD纹波应10mV若超标加10μF钽电容滤波滤波后仍跳变是否未初始化ADC时钟检查RCC_APB2PeriphClockCmd(RCC_APB2Periph_ADC1, ENABLE)是否执行4.3.2 Keil编译报错error: no STM32 target found!的精准修复此错误与滤波算法无关纯属开发环境配置问题打开Project → Options for Target → Device确认选择STM32F10x High-density非Medium-densityTarget页中Xtal(MHz)填8外部晶振频率Use MicroLIB勾选减小printf体积Debug页中Settings → Flash Download勾选Reset and Run若用ST-LinkSettings → SW Device中Connect选Under Reset点击Add添加STM32F10x High-densityFlash算法。提示本工程ADC.uvprojx已配置正确若复制文件到新工程需重新执行上述步骤。切勿直接修改startup_stm32f10x_hd.s中的STACK_TOP地址该值由Keil自动生成。4.3.3 滤波函数执行超时的底层优化当MidFilter()在10kHz采样下超时可实施三级优化编译器级KeilOptions → C/C → Optimization设为Level 3勾选One ELF Section per Function算法级将冒泡排序改为insertion sort对小数组更快adc.c中替换排序循环for(i1; iFILTER_NUM; i) { key buf[i]; j i-1; while(j0 buf[j]key) { buf[j1] buf[j]; j--; } buf[j1] key; }硬件级启用ADC_AnalogWatchdog监控异常值仅对AWD标志置位的采样值启用中位值滤波其余用限幅——大幅降低平均计算量。最终在STM32F103C8T672MHz上经三级优化后MidFilter(N5)执行时间从150μs降至65μs满足10kHz采样需求。本文还有配套的精品资源点击获取