基于二阶统计量的SOBI盲源分离系统

资 源 简 介

该项目实现了一种基于信号二阶统计量的典型盲源分离方法，即二阶盲辨识（SOBI）算法。 其核心功能是在源信号相互独立且具有时间相关性的前提下，仅根据多路传感器的观测混合信号，通过计算不同时延下的多个协方差矩阵并进行联合近似对角化，从而实现源信号的盲恢复与辨识。 项目首先对采集到的多路观测信号进行中心化和白化预处理，以消除信号间的线性相关性并简化混合矩阵的估计过程。 随后，系统提取信号在多个时间延迟点上的自协方差矩阵，这些矩阵包含了源信号丰富的时域结构信息。 通过调用JADE或同样的联合对角化算法，寻找一个正

详 情 说 明

基于信号二阶统计量的SOBI盲源分离系统

项目详细介绍

本项目实现了一种经典的盲源分离（BSS）算法——二阶盲辨识（SOBI）。与传统的独立成分分析（ICA）不同，SOBI不依赖于信号的高阶统计特性，而是利用源信号在时间上的相关性。该系统通过联合对角化多个滞后时间点的协方差矩阵，从多路传感器的观测混合信号中恢复出相互独立的源信号。该项目提供了一个从信号模拟、混合、分离到性能评估的完整闭环仿真环境。

核心功能特性

• 时域结构特征提取：利用源信号的时间相关性，即使在信号不符合非高斯分布的情况下，只要其功率谱形态不同，即可实现有效分离。
• 预处理流水线：集成中心化与白化处理，有效消除传感器通道间的线性相关性，并将估计混合矩阵的问题简化为寻找一个正交矩阵。
• 联合近似对角化：采用基于Jacobi旋转的迭代算法，通过搜索全局最优的正交变换，使多个时延协方差矩阵同时达到近似对角状态。
• 自动化性能评估：自动处理盲源分离中固有的排列和幅度二义性，计算相关系数和信号干扰比（SIR）以量化分离精度。
• 高保真可视化：通过时域波形对比图和残差分析图，直观展示分离前后的信号吻合程度。

系统逻辑与流程说明

系统的核心逻辑严格遵循SOBI算法的标准流程，具体步骤如下：

1. 仿真源信号生成 系统构造了三种具有显著时间相关性的平稳信号：5Hz正弦波、经过低通滤波的8Hz方波（增加相关性程度）以及通过二阶巴特沃斯滤波器生成的彩色噪声。所有源信号均进行了标准化处理，确保单位方差。

2. 瞬时混合模型构建 通过随机生成的混合矩阵将三路源信号线性组合，模拟传感器采集过程。模型中加入了高斯白噪声以模拟真实的物理监测环境。

3. 信号预白化处理 首先对混合信号进行去均值操作。接着通过对零延迟协方差矩阵进行特征值分解（EVD），构建白化矩阵并作用于原始信号。此步骤确保了白化后信号的各分量在统计上互不相关，且方差为1，极大地减小了后续分离算法的搜索空间。

4. 滞后协方差矩阵计算 系统在指定的时延集合（0到15个采样间隔）上计算白化信号的互协方差矩阵。为了保证数值计算的稳定性及特征值的实数性，对计算出的矩阵进行了对称化处理。这些矩阵序列包含了源信号的自相关规律。

5. 联合近似对角化（JAD） 这是算法的核心部分。系统通过循环调用Jacobi旋转，在每一对信号维度上寻找最优的旋转角度。目标是使所有滞后协方差矩阵的非对角元素平方和最小化。通过不断迭代更新，最终得到一个正交变换矩阵，该矩阵代表了白化信号空间到独立源信号空间的转动量。

6. 信号恢复与对齐 利用得到的正交矩阵与白化矩阵合成最终的分离矩阵，进而恢复出原始信号的估计值。由于盲分离存在次序随机和相位翻转的特性，系统利用互相关准则自动匹配恢复信号与原始信号，并进行幅度和符号修正。

关键算法与实现细节分析

Jacobi旋转实现细节：在联合对角化过程中，系统针对每一对索引分量，提取所有滞后矩阵中对应的$2 times 2$子矩阵信息，通过对这些小矩阵的组合进行特征分解，精确计算出旋转余弦和正弦参数，从而逐步逼近全局对角化。

白化（Whitening）的重要性：通过特征值分解将观测信号投影到其主成分空间，这不仅降低了噪声干扰，更保证了后续只需通过旋转矩阵（正交阵）即可恢复源信号，将复杂参数估计问题简化为空间旋转问题。

滞后时延选择：代码中定义了16个不同的时延点，这种多点观测策略增加了算法的鲁棒性，能够更全面地捕获信号的频谱特征差异，从而在复杂环境下维持较高的分离精度。

互相关系数计算函数：系统内置了自定义的互相关计算逻辑，通过分子计算协方差、分母计算方差乘积根号的方式，准确评估信号间的线性相关程度。

系统要求

• 软件环境：MATLAB R2016b 或更高版本。
• 数学工具箱：需要基本的信号处理函数支持（如butter, filter, square等）。
• 硬件要求：普通PC即可，算法复杂度与信号长度和通道数成线性关系，运行效率高。

使用方法

1. 打开环境并定位到脚本所在目录。
2. 直接运行主函数。
3. 系统将自动生成源信号并模拟混合过程。
4. 在命令行窗口（Command Window）实时输出各通道的分离性能指标（相关系数及SIR数值）。
5. 自动弹出绘图窗口，对比展示原始信号、混合信号、分离信号以及恢复误差。