多特征值分解的稀疏混沌信号盲源分离算法研究

周双红，王玲玲

与传统的通信方式相比，激光混沌通信具有很多优点。由于其对初始值的灵敏度具有长期不可预测性，非常适合安全通信应用；通过使用不同的混沌吸引子或相同吸引子的不同初始值和参数产生的混沌信号可以认为是不相关的，因此容易产生大量相互正交的混沌信号，这使得激光混沌信号在多用户通信中具有广泛的应用前景[1]。然而，在激光通信和信号处理等领域中混沌技术的应用也可能面临一些挑战。对于系统现有的信道间干扰，时变衰落或混合多用户载波，盲源分离是一个必须解决的问题。一些研究人员使用混沌信号的动态属性进行盲源分离[2-4]，这种方法只能在源信号动态方程的条件下应用,独立分量分析(ICA)方法也用于分离混沌信号[5-7]。这种方法假设每个源之间是互相统计独立，并通过使用高阶统计性质分离混合信号。然而，激光混沌流信号本质上是确定单一的，并且具有对初始值和宽带光谱非常敏感的特性，因此仅通过使用统计特性难以有效地进行盲源分离。本文使用相空间重构理论描述盲源分离问题,提出了一种针对上述问题的新型混沌信号盲源分离方法。

1 混沌信号的盲源分离方法

表示为矩阵形式：

式中：

图1 混沌流信号的盲源分离模型Fig. 1 Blind source separation model of chaotic stream signal

2 基于相空间的盲源分离理论分析

相空间重建广泛应用于激光混沌序列的分析，如预测、诊断和激光混沌信号的李亚普诺夫指数的计算等。文献[9]通过相空间重建，定义指数来测量流模式混沌序列相空间的生长速率——生长指数。本节将使用生长指数来描述混沌信号的盲源分离问题。

文献[9]证明，当重构的观测序列具有零均值时，具有以下性质。

3 盲源分离算法的过程

前一节中生长指数的性质1是在观测序列为零均值时成立的，因此首先要对观测信号进行均值去除：

生长指数的性质2还要求信号序列彼此不相关，这可以通过预白化观测[信号来实现。假设自相关矩阵特征分解为，其中， Q 是正交矩阵，是对角矩阵，那么称为白化矩阵，在线性变换后，的每个分量彼此不相关。

这种方法可以减少要优化的参数，从而显著提高算法的收敛速度和鲁棒性。本节后续部分采用以下形式的参数化矩阵来表示2阶和3阶正交矩阵：

对于非约束优化问题，粒子群优化算法是一种非常有效的方法。将参数向量视为粒子位置，将目标函数式(14)视为适应度函数，可以使用粒子群优化算法估计最优分离矩阵，以重构每个信道源信号。整个算法过程为：

4) 对于每个粒子，根据式 (11)、(12)、(13)计算分离矩阵B，根据式(3)计算分离信号，根据式(5)实现相位空间重构；根据式(7)～(9)计算适度函数值；

5) 应用粒子群算法优化式(14)直到满足终止条件，记录优化粒子位置θopt；

4 仿真实验

本节将通过仿真实验来评估盲源分离算法的性能。源信号通过式(1)产生混合信号，混合矩A的元素在每次仿真中通过服从[-1,1]独立均匀分布随机数产生。盲源分离算法的精度将通过性能指标PI来测量：

式中gij为的元素，PI越小越好，当G满足式(6)时，PI得到最小值0。

仿真实验中使用的4个源信号分别由Rossler吸引子式(16)、洛伦兹吸引子式(17)、Duffing吸引子式(18)和Mackey-Glass吸引子式(19)产生，前3个信号通过4阶龙格库塔法积分得到，积分步长分别为0.05、0.05、0.01，信号的长度为500 s，如图2所示。利用4个信号组成如表1所示的三通道混合信号和双通道混合信号，用于测试盲源分离效果的。

图2 源信号波形Fig. 2 source signal waveform

表1 测试信号组Table 1 Test signal grouping

4.1 混沌流无噪声信号的盲源分离仿真结果

对于表1中的三通道源信号，利用本文提出的方法进行4次盲源分离实验，性能指标PI随迭代步骤变化，如图3所示。随着迭代步长的增加，PI快速减小，在所有实验中，算法可以通过几十次迭代收敛。表2给出了无噪声环境下，对表1中两组源信号进行100次蒙特卡罗模拟之后的实验结果，可以看出，对于组1中的三通道混合信号和组2中的双通道混合信号的分离，提出的算法的分离精度优于快速独立分量分析(FastICA)、蝙蝠算法(BA)和差分进化算法(DEA)。

图3 无噪声条件下的收敛条件Fig. 3 The convergence condition of the algorithm at no noise

表2 无噪声条件下盲源分离算法的性能Table 2 The blind source separation performance at no noise

4.2 噪声环境下混沌流信号的盲源分离仿真结果

本节将评估基于蚁群技术的多特征值分解盲源分离算法与基于FastICA、最小互信息和极大似然估计法对噪声污染混合信号的分离性能。对于表1中的两组源信号的混合信号加入高斯白噪声。让观测噪声从20 dB变化到50 dB，以5 dB的间隔测试点，在每个测试点重复100次蒙特卡罗模拟，各分离算法的分离性能指数如图4所示。

图4 不同盲源分离算法在不同噪声强度下的性能Fig. 4 The performance of blind source separation algorithm under different noise intensities

5 结束语

本文描述了使用相空间重构理论的盲源分离，并为激光混沌流信号提出了一种新的盲源分离算法。该算法通过构建分离信号相位空间中的目标函数，将混沌流信号的盲源分离转换为无约束优化问题，并通过粒子群优化算法求解。同时，它采用分离矩阵上的正交矩阵的参数表示，有效地减少了优化问题的维数，因此该算法可以快速收敛。仿真结果表明，该算法不仅具有快速收敛的特点，其在各种SNR下的精度明显优于FastICA、最小互信息和极大似然算法。