【大窗除强信号，小窗清残留】基于双尺度广义交叉验证阈值的地震信号自适应剥离和噪声提取方法（MATLAB）

背景知识

在环境噪声层析成像等研究中，我们需要的是纯粹的“噪声”记录，而不是被地震信号“污染”的波形。传统方法是人工剔除含事件的时间段，或者用时间域归一化压制信号，但这些方法要么主观，要么难以彻底去除能量较强的信号。广义交叉验证是一种不依赖噪声先验信息的自适应阈值选取方法，它通过最小化预测误差来找到最优阈值。同步压缩变换则能将信号能量在时频域中高度聚焦，使信号和噪声更容易区分。本算法结合了这两者，专门用于从地震记录中“剥离”有效信号，保留背景噪声，服务于环境噪声研究。

算法总结

这个算法的核心任务是从含噪地震记录中移除有效信号，提取纯粹的背景噪声。它先把信号分成多个滑动窗口，对每个窗口做连续小波变换，得到时频系数。首先，它利用峰度这个统计量做预处理，剔除那些明显呈高斯分布（即纯噪声）的尺度，因为这些尺度不含信号信息。接着，对保留下来的系数，用斐波那契搜索法快速找到使广义交叉验证函数最小的最优阈值，并对系数进行一种特殊的“硬阈值”处理——将大于阈值的系数压缩为阈值的一半，而不是直接置零，这样既能削弱信号又能保留噪声成分。然后，将预处理时剔除的噪声尺度系数“恢复”回来，确保噪声信息完整。最后逆变换回时域，拼接到完整信号中。由于一次处理可能不够彻底，建议用两个不同大小的窗口先后处理两遍：先用大窗口去除主体信号，再用小窗口精雕细琢，把残留的微弱信号也清理掉，最终得到纯净的噪声记录。

算法步骤

读取数据与参数设置

输入含噪地震信号，设定采样间隔、小波类型、每个八度的声音数、高斯修正因子、窗口大小（秒）等参数。

滑动窗口划分

根据设定的窗口大小（秒）计算对应的采样点数，按此宽度将长信号分成若干无重叠的窗口，逐窗口处理。

窗口数据归一化

取当前窗口数据，找到最大值，将整个窗口归一化到 [0,1] 区间，方便后续阈值计算。

连续小波变换

对归一化后的窗口数据做CWT，得到时频系数矩阵（行：尺度，列：时间）。

预处理：基于峰度剔除高斯噪声尺度

对每个尺度，计算其实部系数的峰度值。峰度衡量数据偏离高斯分布的程度。

如果峰度绝对值小于某阈值（由高斯修正因子控制），说明该尺度为高斯噪声，将其系数全部置零，并记录该尺度索引。

广义交叉验证最优阈值搜索

对每个未被置零的尺度，用斐波那契搜索法在系数绝对值范围内寻找使广义交叉验证函数最小的最优阈值

特殊硬阈值处理

将该尺度系数中绝对值大于最优阈值的部分，全部替换为阈值的一半；其余系数保持不变。这种处理方式能有效削弱信号能量，同时保留噪声。

恢复被剔除的噪声尺度

将步骤5中置零的尺度系数，用原始 CWT 系数替换回来。这样做可以保证噪声成分完整，不被误删。

逆连续小波变换

对处理后的系数矩阵做逆变换，重建该窗口的时域信号，并乘以步骤3的归一化因子恢复幅度。

窗口拼接与双重处理

将所有窗口处理后的结果按顺序拼接，得到第一次处理后的信号。

以该信号作为新的“含噪”输入，换一个更小的窗口（比如从500秒换成55秒），重复步骤1-9，进一步去除残留的短时信号，得到最终的纯净噪声记录。

结果输出与绘图

保存最终信号（提取的噪声），并绘制原始信号与提取噪声的对比图。

如果你对信号滤波/降噪，机器学习/深度学习，时间序列预分析/预测，设备故障诊断/缺陷检测/异常检测有疑问，或者需要论文思路上的建议，欢迎咨询

担任《MSSP》《中国电机工程学报》《宇航学报》《控制与决策》等期刊审稿专家，擅长领域：信号滤波/降噪，机器学习/深度学习，时间序列预分析/预测，设备故障诊断/缺陷检测/异常检测