嵌入式数据采集中的数字滤波算法对比与实践

1. 滤波技术在数据采集中的核心价值

在工业自动化、环境监测、医疗设备等实时数据采集场景中，传感器原始信号往往混杂着各种噪声干扰。上周调试一个温湿度监测项目时，就遇到了典型的信号毛刺问题——当空调突然启动时，温度传感器的ADC读数会出现剧烈跳变，导致控制算法误判。这种场景下，选择合适的数字滤波算法就成了保证数据可靠性的关键。

一阶低通滤波（First-order Low-pass Filter）和去极值平均滤波（Outlier-removing Average Filter）是嵌入式领域最常用的两种预处理方案。前者通过频域处理平滑信号，后者在时域上剔除异常值后计算均值。我在多个工业级数据采集项目中实测发现：对于采样率1kHz以下的系统，这两种算法在STM32F103上的执行时间差异不超过50μs，但滤波效果却大相径庭。

2. 算法原理深度解析

2.1 一阶低通滤波的数学本质

一阶低通滤波的差分方程表示为：

code复制y(n) = α * x(n) + (1-α) * y(n-1)

其中α=Δt/(RC+Δt)，Δt为采样间隔。在STM32CubeMonitor实测中，当α取0.2时，对50Hz工频干扰的抑制效果可达-15dB。但要注意截止频率fc与α的关系：

code复制fc = α / (2πΔt)

去年在电机振动监测项目中，就曾因错误设置α=0.5（对应fc≈160Hz）导致关键的100Hz谐波成分被滤除。经验表明，对于多数工业信号，α取值0.05-0.3较为合适。

2.2 去极值平均滤波的工程实现

该算法的典型实现流程：

1. 缓存最近N个采样值（通常N取4-16）
2. 排序后剔除最大和最小的K个值（K通常为N/4）
3. 计算剩余数据的算术平均值

在电化学传感器项目中验证发现：当N=8、K=2时，能有效抑制±20%的突发干扰。但需注意内存消耗——对于100通道的采集系统，每个通道缓存8个float类型数据将占用3.2KB RAM。

3. 对比实验设计与结果分析

3.1 测试环境搭建

使用标准信号发生器输出叠加噪声的测试信号：

• 基波：10Hz正弦波（模拟真实信号）
• 干扰：50Hz工频噪声（幅值30%）
• 脉冲：随机出现的2倍幅值尖峰（模拟电磁干扰）

通过STM32H743的16位ADC采样，分别应用两种滤波算法处理。用J-Scope实时观测波形，记录关键参数：

3.2 典型场景适用性分析

一阶低通滤波优势场景：

• 高频噪声为主的信号（如开关电源干扰）
• 对相位延迟不敏感的控制系统
• 资源受限的MCU（RAM<1KB）

去极值平均滤波优势场景：

• 存在突发性大幅干扰（如继电器动作）
• 需要保持相位特性的测量系统
• 具有充足RAM的处理器

在无人机气压计数据处理中，就曾混合使用两种算法：先用去极值平均滤除异常跳变，再用一阶低通平滑高频波动，组合后高度数据方差降低62%。

4. 嵌入式实现中的优化技巧

4.1 定点数优化方案

对于没有FPU的MCU（如STM32F030），可将一阶低通滤波改为定点运算：

c复制#define ALPHA_Q8 51  // 0.2 in Q8 format
int16_t filter_lp(int16_t x_new, int16_t y_prev) {
    int32_t tmp = (ALPHA_Q8 * x_new) + ((256 - ALPHA_Q8) * y_prev);
    return (int16_t)(tmp >> 8);
}

实测在Cortex-M0上，该实现比浮点版本快5倍。

4.2 去极值算法的快速实现

避免每次全排序的优化方法：

c复制int16_t filter_avg(int16_t *buf, uint8_t N) {
    int16_t min = INT16_MAX, max = INT16_MIN;
    int32_t sum = 0;
    
    for(uint8_t i=0; i<N; i++) {
        if(buf[i] < min) min = buf[i];
        if(buf[i] > max) max = buf[i];
        sum += buf[i];
    }
    return (sum - min - max) / (N - 2);
}

此方法时间复杂度从O(NlogN)降至O(N)，在N=12时执行时间从56μs降至18μs。

5. 实际项目中的选择建议

根据多年项目经验总结的决策流程图：

1. 明确干扰特征：

   • 用示波器观察噪声频谱
   • 记录异常幅值的发生频率

2. 评估系统资源：

   • 可用RAM是否大于N通道数2字节
   • CPU是否具备DSP指令或FPU

3. 验证相位影响：

   • 对阶跃响应的延迟要求
   • 控制系统的稳定性边际

4. 混合方案考量：

   • 级联使用时的顺序影响
   • 参数协调（如先平均后低通时α需减小）

在最近的光伏逆变器项目中，最终选用α=0.15的一阶低通配合硬件RC滤波，在满足-3dB截止频率100Hz的同时，将THD控制在0.8%以下。而农业大棚监测系统则采用N=5的去极值平均，成功滤除灌溉水泵启动造成的温度数据跳变。