实时中值滤波 + 低通滤波 示例程序(STM32环境)

一、功能概述

本示例实现两个滤波器:

  • 中值滤波器(Median Filter):对短期异常值(如尖峰噪声)有良好的抑制能力;
  • 低通滤波器(Low-Pass Filter):对数据进行平滑处理,适合去除高频抖动。

适用于如 ADC传感器数据处理、电子秤、电流电压读取等场景。


二、代码结构说明

✅ 头文件 filter.h

c 复制代码
#ifndef FILTER_H
#define FILTER_H

#define WINDOW_SIZE 50

typedef struct {
    int buffer[WINDOW_SIZE]; // 数据缓冲区
    int index;               // 当前写入索引(环形)
    int count;               // 当前已存数据个数
} MedianFilter;

extern MedianFilter myMedianFilter;

void medianInit(MedianFilter* filter);
int  medianFilter(MedianFilter* filter, int new_value);
int  lowPassFilter(int input);

#endif

✅ 源文件 filter.c

🔹 中值滤波器实现
c 复制代码
#include "filter.h"
#include <stdlib.h>

// 冒泡排序用于求中值(可替换为更高效算法)
static void sort_array(int* src, int* dst, int size) {
    for (int i = 0; i < size; i++) dst[i] = src[i];
    for (int i = 0; i < size - 1; i++) {
        for (int j = i + 1; j < size; j++) {
            if (dst[i] > dst[j]) {
                int temp = dst[i];
                dst[i] = dst[j];
                dst[j] = temp;
            }
        }
    }
}

// 初始化滤波器状态
void medianInit(MedianFilter* filter) {
    for (int i = 0; i < WINDOW_SIZE; i++)
        filter->buffer[i] = 0;
    filter->index = 0;
    filter->count = 0;
}

// 中值滤波主函数
int medianFilter(MedianFilter* filter, int new_value) {
    filter->buffer[filter->index] = new_value;
    filter->index = (filter->index + 1) % WINDOW_SIZE;
    if (filter->count < WINDOW_SIZE)
        filter->count++;

    int sorted[WINDOW_SIZE];
    sort_array(filter->buffer, sorted, filter->count);

    if (filter->count % 2 == 1)
        return sorted[filter->count / 2];
    else
        return (sorted[filter->count / 2 - 1] + sorted[filter->count / 2]) / 2;
}
🔹 低通滤波器实现
c 复制代码
int lowPassFilter(int input) {
    static int prev_output = 0;
    uint8_t alpha_shift = 3;  // α = 1/8,越大越平滑
    prev_output = prev_output + ((input - prev_output) >> alpha_shift);
    return prev_output;
}

✅ 主函数测试示例

c 复制代码
#include <stdio.h>
#include "filter.h"

MedianFilter myMedianFilter;

int mainTest() {
    medianInit(&myMedianFilter);

    int data[] = {100, 102, 98, 500, 101, 99, 97, 103};
    int len = sizeof(data) / sizeof(data[0]);

    for (int i = 0; i < len; i++) {
        int medFiltered = medianFilter(&myMedianFilter, data[i]);
        int lpFiltered = lowPassFilter(data[i]);

        printf("Raw: %d\tMedian: %d\tLowPass: %d\n", data[i], medFiltered, lpFiltered);
    }

    return 0;
}

三、使用说明

  • 可移植至 STM32 项目,适合 ADC 滤波等场景;
  • 中值滤波适合抗冲击干扰,低通滤波适合抑制小幅高频抖动;
  • WINDOW_SIZE 可根据实际采样频率和响应速度调整;
  • 建议中值窗口不宜过大,推荐 5~51 范围内奇数;
  • 实际效果图:
相关推荐
半杯咖啡半行码1 分钟前
STM32开发之IIC通信详解:从协议基础到FM24CL04实战
stm32·单片机·嵌入式硬件
硅基手札1 小时前
【I2C 技术系列 05】I2C 多主仲裁与时钟同步
驱动开发·stm32·单片机·嵌入式硬件·mcu·计算机外设
szxinmai主板定制专家1 小时前
面向虚拟电厂协调控制:RK3588+FPGA+AI边缘控制器设计思路与算力评估
人工智能·嵌入式硬件·fpga开发·机器人·jetson·zynq·振动采集
Zw-awa1 小时前
DO 和 AO 到底怎么选:比较器替你决定了什么
stm32
硅基手札1 小时前
【I2C 技术系列 01】I2C 电气特性与物理层
驱动开发·stm32·单片机·嵌入式硬件·mcu·计算机外设
comyccccc1 小时前
TB6612电机驱动板原理图与PCB设计学习笔记
单片机·嵌入式硬件
gscsded2 小时前
MCUboot裸机移植实战
单片机·嵌入式硬件·arm
北京盛世宏博3 小时前
多协议以太网温湿度变送器:Modbus TCP/UDP/SNMP 协议应用场景详解
单片机·tcp/ip·udp·传感器·温湿度
HRTOS3 小时前
HRTOS 内核 BUG 修复:任务等待对象未初始化导致的异常唤醒问题
单片机·嵌入式硬件·算法·51单片机
@MMiL4 小时前
SImulink: PowerGui下的信号FFT分析
单片机