第四节：噪声与误差处理

做皮肤阻抗测量，最头疼的是什么？

不是算法不够准，而是信号太脏。我刚开始做这个项目时，满心以为ADC采回来就能直接用。结果一看波形，好家伙，50Hz工频干扰比信号还大，稍微动一下探头就飘到天上去。今天咱们就把这些坑一个一个填上。

4.1 工频50Hz陷波：干掉那个讨厌的嗡嗡声

你想想看，我们生活在一个充满50Hz电磁场的世界里。墙壁里的电线、桌上的充电器、甚至日光灯，都在往外辐射这个频率。皮肤阻抗信号本身就很微弱，很容易被它淹没。

我个人的习惯是，在硬件上先做一级模拟陷波，但别指望它能完全滤干净。模拟滤波器有温漂，元件有误差，所以软件上必须再补一刀。

数字陷波器怎么做？最经典的是双二阶IIR陷波。我直接给代码：

// 50Hz陷波器系数计算（采样率Fs=1000Hz）
// 陷波频率f0=50Hz，Q值取30
float b0, b1, b2, a0, a1, a2;
float w0 = 2 * PI * 50.0 / 1000.0;  // 归一化角频率
float alpha = sin(w0) / (2 * 30.0);  // Q=30

b0 = 1.0;
b1 = -2.0 * cos(w0);
b2 = 1.0;
a0 = 1.0 + alpha;
a1 = -2.0 * cos(w0);
a2 = 1.0 - alpha;

// 归一化系数
b0 /= a0; b1 /= a0; b2 /= a0;
a1 /= a0; a2 /= a0;

// 滤波函数
float notch_filter(float x) {
    static float x1=0, x2=0, y1=0, y2=0;
    float y = b0*x + b1*x1 + b2*x2 - a1*y1 - a2*y2;
    x2 = x1; x1 = x;
    y2 = y1; y1 = y;
    return y;
}

4.2 运动伪迹抑制：别让手抖毁了数据

做美容仪测量时，用户不可能像做心电图那样躺着不动。探头一滑动，接触压力一变，阻抗值就跟着乱跳。这种运动伪迹，频率通常在0.5~5Hz之间，和真正的皮肤阻抗变化混在一起。

我在项目中遇到过最夸张的情况：用户一边敷面膜一边刷手机，探头在脸上画圈圈，测出来的阻抗曲线跟心电图似的。怎么办？

我的做法分两步：

1. 加速度辅助检测：在探头里加一颗三轴加速度计。检测到加速度突变时，标记这段时间的数据为「可疑」。
2. 自适应中值滤波：对于可疑数据段，用滑动中值滤波代替均值滤波。中值滤波对突发性干扰有奇效。

// 运动伪迹检测与抑制
#define WINDOW_SIZE 5
float motion_suppress(float new_sample, float accel_magnitude) {
    static float buffer[WINDOW_SIZE];
    static int idx = 0;
    static int motion_flag = 0;
    
    // 运动检测阈值，根据实际调试
    if (accel_magnitude > 0.5) {  // 单位：g
        motion_flag = 1;
    } else {
        motion_flag = 0;
    }
    
    buffer[idx] = new_sample;
    idx = (idx + 1) % WINDOW_SIZE;
    
    if (motion_flag) {
        // 运动期间用中值滤波
        float sorted[WINDOW_SIZE];
        memcpy(sorted, buffer, sizeof(sorted));
        // 简单冒泡排序找中值
        for (int i=0; i<WINDOW_SIZE-1; i++)
            for (int j=0; j<WINDOW_SIZE-1-i; j++)
                if (sorted[j] > sorted[j+1]) {
                    float t = sorted[j];
                    sorted[j] = sorted[j+1];
                    sorted[j+1] = t;
                }
        return sorted[WINDOW_SIZE/2];
    } else {
        // 正常情况用均值滤波
        float sum = 0;
        for (int i=0; i<WINDOW_SIZE; i++) sum += buffer[i];
        return sum / WINDOW_SIZE;
    }
}

4.3 接触阻抗影响：电极没贴好，一切白搭

皮肤阻抗测量用的是两电极还是四电极？这里有个大坑。

两电极法简单，但电极-皮肤接触阻抗会直接串入测量结果。你以为是皮肤阻抗变了，其实是用户出汗导致接触阻抗下降了。四电极法能消除接触阻抗，但电路复杂，美容仪上很少用。

那怎么办？我分享一个实用技巧：

• 测量前做接触质量检测：在正式测量前，先注入一个已知的小电流，测一下接触阻抗。如果超过阈值（比如10kΩ），提示用户重新贴好电极。
• 动态补偿：测量过程中实时监测接触阻抗的变化。如果发现接触阻抗突变，用查表法修正测量值。

接触阻抗范围	建议操作	补偿系数
< 5kΩ	正常测量	1.0
5kΩ ~ 10kΩ	轻度补偿	1.05
10kΩ ~ 20kΩ	提示用户调整	1.15
> 20kΩ	停止测量	无效

4.4 温度补偿：皮肤也是会「热胀冷缩」的

皮肤阻抗对温度极其敏感。你想想看，温度升高1°C，皮肤血管扩张，血流量增加，阻抗能下降2%~5%。如果用户从空调房走到阳光下，阻抗值能差出20%。

我踩过这个坑。有一次做临床测试，上午测一批，下午测一批，数据对不上。查了半天，原来是上午空调开得足，下午太阳晒进来，室温差了8°C。

温度补偿怎么做？我的经验是：

1. 内置温度传感器：在电极附近放一颗数字温度传感器，精度±0.1°C就够了。
2. 建立补偿模型：在实验室里测出不同温度下的阻抗变化曲线，拟合出一个补偿公式。

// 温度补偿函数
// 基准温度：25°C
// 补偿系数：每°C变化约3.5%（经验值，需实际标定）
float temp_compensate(float impedance, float temperature) {
    const float REF_TEMP = 25.0;
    const float TC = 0.035;  // 温度系数，单位：1/°C
    
    float delta_temp = temperature - REF_TEMP;
    float compensation = 1.0 + TC * delta_temp;
    
    // 补偿后的阻抗值
    return impedance / compensation;
}

嗯，噪声和误差处理就讲到这里。说白了，做嵌入式算法，一半时间在写功能，另一半时间在跟噪声作斗争。但正是这些细节，决定了你的产品是「能用」还是「好用」。

下一节咱们聊聊阻抗谱的频域分析，那又是另一个有意思的话题了。