血氧的滤波算法代码怎么写

血氧饱和度测量中，滤波算法通常是为了减少噪声影响和提高读数稳定性。常用的滤波技术有简单平均滤波、滑动平均滤波以及数字低通滤波等。以下是使用简单滑动平均滤波的一个基本示例：

#include <stdio.h>
#include <math.h>

#define FILTER_SIZE 5 // 滤波窗口大小

// 滤波器结构体
typedef struct {
    int samples[FILTER_SIZE];
    int index;
} MovingAverageFilter;

// 初始化滤波器
void filter_init(MovingAverageFilter* filter) {
    for (int i = 0; i < FILTER_SIZE; ++i) {
        filter->samples[i] = 0;
    }
    filter->index = 0;
}

// 添加新的血氧值到滤波器
void add_sample_to_filter(int sample, MovingAverageFilter* filter) {
    filter->samples[filter->index] = sample;
    filter->index = (filter->index + 1) % FILTER_SIZE; // 循环数组
}

// 计算和返回滤波后的血氧值
int filtered_spo2(MovingAverageFilter* filter) {
    int sum = 0;
    for (int i = 0; i < FILTER_SIZE; ++i) {
        sum += filter->samples[i];
    }
    return sum / FILTER_SIZE;
}

int main() {
    MovingAverageFilter filter;
    filter_init(&filter);

    // 获取连续的血氧样本...
    for (int i = 0; i < N_SAMPLES; ++i) {
        int new_sample = get_spo2(); // 假设get_spo2()为获取最新血氧值函数
        add_sample_to_filter(new_sample, &filter);
        int filtered = filtered_spo2(&filter);
        printf("滤波后血氧: %d%%\n", filtered);
        // 更新显示或其他操作
    }
    return 0;
}

这个例子中，每次获取新的血氧值后，将其加入滤波器，然后计算前`FILTER_SIZE`个值的平均作为结果。滤波效果取决于窗口大小的选择，较大的窗口可以提供更平滑的结果但响应会变慢。