从自相关函数看噪声本质：平稳与非平稳噪声的实战识别与处理-拓冰建站

1. 噪声分析的基础：自相关函数是什么？

第一次接触自相关函数这个概念时，我也被这个拗口的名字吓到了。但后来发现，它其实就是信号处理领域的"照妖镜"——能帮我们看清噪声的真面目。简单来说，自相关函数描述的是信号在不同时间点的自我相似程度。想象一下你在KTV唱歌时听到的回声：如果回声和原声完全一致（相关系数为1），说明房间的混响特性非常稳定；如果回声忽大忽小，就说明房间的声学特性在变化。自相关函数就是用来量化这种相似性的数学工具。

在MATLAB中计算自相关函数特别简单，用autocorr()函数就能搞定。比如我们生成一个简单的白噪声信号：

noise = randn(1000,1); % 生成长度1000的高斯白噪声 autocorr(noise); % 绘制自相关函数图

运行后会看到两条红线——那是95%置信区间。理想的白噪声应该只有lag=0时有峰值，其他lag位置都落在置信区间内。但实际工程中遇到的噪声往往没这么"完美"，这时候就需要我们像老中医看脉象一样，从自相关函数图的形态中解读出噪声的特性。

2. 平稳噪声的指纹特征

去年调试工业传感器时，我遇到过典型的平稳噪声。当时采集的振动信号看起来杂乱无章，但自相关函数却暴露了它的本质特征。平稳噪声有三大典型特征：

均值稳定：就像老式挂钟的滴答声，不管什么时候听，响度都差不多
方差恒定：波动幅度不会忽大忽小
自相关函数时不变：不同时间段计算的自相关函数图形状基本一致

用Python的statsmodels库可以直观展示这一点：

import numpy as np import statsmodels.api as sm np.random.seed(42) stationary_noise = np.random.normal(0, 1, 1000) sm.graphics.tsa.plot_acf(stationary_noise, lags=40) plt.show()

这个例子中，除了lag=0处的峰值，其他位置的自相关系数都在蓝色置信带内随机波动，这就是典型的平稳白噪声。但要注意，平稳噪声不一定都是白噪声——比如下面这个有色噪声的例子：

% 生成有色噪声 b = [1 0.5]; a = [1 -0.8]; colored_noise = filter(b, a, randn(1000,1)); autocorr(colored_noise);

这时会看到自相关函数呈现指数衰减，但整体形态在不同时间段保持稳定，这仍然属于平稳噪声范畴。

3. 非平稳噪声的识别技巧

上个月分析一组环境温度数据时，我踩了个坑——一开始误判为非平稳噪声，后来才发现是传感器故障。这个教训让我总结出非平稳噪声的三个报警信号：

自相关函数时变：上午和下午测得的自相关函数图形状明显不同
峰值位置漂移：周期性信号的峰值lag值会随时间变化
衰减模式突变：相关性的衰减速度在不同时段不一致

用MATLAB模拟一个典型的非平稳噪声案例：

t = 0:0.01:10; noise = zeros(size(t)); % 前5秒是平稳噪声 noise(1:500) = randn(1,500); % 后5秒加入时变成分 noise(501:end) = sin(2*pi*0.5*t(501:end)).*randn(1,500); figure subplot(2,1,1) plot(t,noise) title('时变噪声信号') subplot(2,1,2) autocorr(noise(1:500)) % 前5秒自相关 hold on autocorr(noise(501:end)) % 后5秒自相关 legend('前5秒','后5秒')

这个例子中，前后两段信号的自相关函数明显不同——前段符合平稳噪声特征，后段则显示出周期性成分。在实际工程中，这种突变往往意味着系统工作状态发生了变化，比如设备开始出现周期性故障。

4. 实战中的混合噪声处理

真实的工程信号往往更复杂。去年参与风电项目时，我们采集的齿轮箱振动信号就同时包含：

平稳部分：来自背景电磁噪声
非平稳部分：齿轮磨损导致的周期性冲击

处理这类混合噪声时，我常用的方法是分段自相关分析。具体操作步骤：

数据分帧：用滑动窗口将信号分成若干段
逐帧计算：对每帧数据计算自相关函数
动态观察：比较各帧自相关函数的演变

Python实现示例：

from scipy import signal def sliding_autocorr(x, window_size=100, overlap=50): n_frames = len(x) // (window_size - overlap) fig, ax = plt.subplots(n_frames, 2, figsize=(12, 2*n_frames)) for i in range(n_frames): start = i * (window_size - overlap) frame = x[start:start+window_size] # 绘制时域波形 ax[i,0].plot(frame) ax[i,0].set_title(f'Frame {i+1} Time Domain') # 绘制自相关函数 sm.graphics.tsa.plot_acf(frame, ax=ax[i,1], lags=20) ax[i,1].set_title(f'Frame {i+1} Autocorrelation') plt.tight_layout() return fig

通过这种方法，我们成功捕捉到了齿轮损伤导致的周期性冲击成分——在某些特定帧的自相关函数中出现了异常的周期性峰值。