STM32实战：双滤波算法在传感器数据处理中的应用（附源码）
1. 为什么传感器数据需要双重滤波？
在嵌入式开发中，传感器数据就像个调皮的孩子，总爱给我们制造惊喜（或者说惊吓）。我做过一个电子秤项目，当有人走过秤体附近时，ADC采集到的原始数据会出现明显的毛刺。这时候如果直接用这些数据做计算，结果就会像过山车一样刺激。
中值滤波就像个严厉的班主任，专门对付那些突然冒头的异常值。它的原理很简单：在一个滑动窗口内对数据进行排序，取中间值作为输出。实测下来，对于消除突发性干扰特别有效。比如电子秤上突然落下一滴水，或者电流检测时遇到电机启动的浪涌，都能被它完美化解。
但中值滤波也有软肋。当遇到持续的小幅度波动时（比如电子秤受到风扇气流影响），它就力不从心了。这时候低通滤波就该登场了——它像个和事佬，专门平滑处理高频抖动。我常用的是一阶IIR滤波器，计算量小到连STM32F103这种入门级MCU都能轻松应对。
2. 双滤波算法的具体实现
2.1 中值滤波的优化技巧
原始代码里用了冒泡排序，这在窗口大小为50时确实有点吃力。后来我优化成插入排序，速度直接提升30%。对于资源紧张的STM32F0系列，还可以考虑这些技巧：
// 快速中值查找算法（适用于窗口大小固定为奇数） int quickMedian(int arr[], int size) { int i, j, pos = size/2; for(i = 0; i <= pos; i++) { int min_idx = i; for(j = i+1; j<size; j++) if(arr[j] < arr[min_idx]) min_idx = j; swap(&arr[i], &arr[min_idx]); } return arr[pos]; }窗口大小的选择也很有讲究：
- 电子秤推荐用15-25点（响应速度约0.5秒）
- 电流检测建议5-11点（保证实时性）
- 温度检测可以到31点（信号变化慢）
2.2 低通滤波的参数调校
很多人不知道，低通滤波的效果主要取决于这个α系数：
#define ALPHA_SHIFT 3 // α=1/8 int lowPassFilter(int input) { static int prev = 0; prev += (input - prev) >> ALPHA_SHIFT; return prev; }通过实测发现：
- 位移数3-4（α=1/8~1/16）适合大多数场景
- 位移数2（α=1/4）适合快速变化的信号
- 位移数5（α=1/32）适合超稳定需求
有个小技巧：可以用PWM输出模拟信号源，用逻辑分析仪观察滤波效果。我在调试电流传感器时，就这样找到了最优参数。
3. STM32上的实战应用
3.1 ADC采集的完整流程
以STM32CubeIDE环境为例，完整的数据处理流程应该是这样的：
// 在ADC中断回调函数中 void HAL_ADC_ConvCpltCallback(ADC_HandleTypeDef* hadc) { uint32_t raw = HAL_ADC_GetValue(hadc); int med = medianFilter(&adcFilter, raw); // 先中值滤波 currentValue = lowPassFilter(med); // 再低通滤波 // 后续处理... }注意两个细节：
- 中值滤波器的缓冲区要放在SRAM高速区域
- 如果采样率很高（>1kHz），建议使用DMA+双缓冲
3.2 内存优化方案
当需要同时处理多路传感器时，内存可能吃紧。我的解决方案是：
- 使用共用缓冲区（适用于分时采样）
- 改用更紧凑的数据类型（如int16_t）
- 动态调整窗口大小：
typedef struct { int16_t *buffer; // 改为指针动态分配 uint8_t size; // 运行时确定窗口大小 // ...其他字段 } MedianFilter; void medianInit(MedianFilter* f, uint8_t size) { f->buffer = malloc(size * sizeof(int16_t)); // ...初始化 }4. 常见问题与性能对比
4.1 滤波算法的选择困境
很多工程师会纠结到底该用哪种组合，根据我的项目经验：
| 应用场景 | 推荐配置 | 处理时间(us)@72MHz |
|---|---|---|
| 电子秤 | 中值25点 + 低通α=1/16 | 142 |
| 电池电压检测 | 中值5点 + 低通α=1/8 | 38 |
| 电机转速检测 | 中值11点 + 低通α=1/4 | 67 |
4.2 那些年我踩过的坑
采样率与窗口大小的关系：曾经有个项目，采样率1kHz却设了50点窗口，导致数据延迟明显。后来总结出经验公式：窗口点数 ≤ 采样周期/信号周期
整数运算的溢出问题：低通滤波中的累加操作可能溢出，特别是用int16_t时。解决方法要么用int32_t，要么定期复位基准值：
// 每1000次采样重置基准 static int counter = 0; if(++counter >= 1000) { prev_output = input; counter = 0; }- 多任务环境下的竞争：在RTOS中，滤波器变量可能被多个任务访问。最简单的保护方法是关中断：
uint32_t primask = __get_PRIMASK(); __disable_irq(); int result = medianFilter(&filter, newValue); __set_PRIMASK(primask);源码已经打包好，包含STM32CubeMX配置文件和测试用例。下载后可以直接用J-Link调试观察变量变化，建议先用杜邦线给ADC输入可调电压，亲眼看看滤波效果。