[发明专利]基于结构化稀疏子空间聚类的欠定盲源分离方法和设备

权利要求书

1.一种基于结构化稀疏子空间聚类的欠定盲源分离方法，其特征在于，所述方法包括以下步骤：

将源信号变换至时频域，采用时频区间加窗法去除低能量点，筛选出单源点；

采用结构化稀疏子空间聚类发掘所述单源点中分布的子空间，确定信号源数量以及每个所述单源点所属的簇；

采用改进势函数法求簇中心，基于所述簇中心估计混合矩阵；

基于所述混合矩阵，使用基于源数目估计的空间投影法结合最小范数法对所述源信号进行重构；

使用短时傅里叶逆变换将时频域的重构源信号恢复到时域，得到估计的时域源信号。

2.如权利要求1所述的基于结构化稀疏子空间聚类的欠定盲源分离方法，其特征在于，采用结构化稀疏子空间聚类发掘所述单源点中分布的子空间，确定信号源数量以及每个所述单源点所属的簇，包括步骤：

初始状态下，单源点的分类位置未知，结构化稀疏子空间聚类的模型不含结构化稀疏范数项，结构化稀疏子空间聚类的模型表示为：

s.t.X＝XC+E,diag(C)＝0

其中，X表示单源点组成的数据集，C为单源点的自我表示系数矩阵，E表示噪声矩阵，λ＞0为平衡参数，表示对高斯噪声的正则化约束项，||C||₁表示单源点的自我表示系数矩阵的范数；

根据所述模型，采用交替方向乘子法求得自我表示系数矩阵，基于所述自我表示系数矩阵，使用谱聚类进行初始类别划分，从而确定源信号数量以及每个单源点所属的簇。

3.如权利要求2所述的基于结构化稀疏子空间聚类的欠定盲源分离方法，其特征在于，所述使用谱聚类进行初始类别划分，确定源信号数量的步骤包括：

基于第一数学式求相似度矩阵：

其中，C^T表示自我表示系数矩阵C的转置，W为相似度矩阵；

基于第二数学式求标准化的拉普拉斯矩阵：

其中，L_sym为标准化的拉普拉斯矩阵，I为单位矩阵，D为度矩阵，具体为一个对角矩阵，对角线上的元素基于第四数学式求得：

对矩阵L_sym进行特征分解，并将其特征值按照由小到大的顺序排列，定义特征值间隙为：

Δλ_m＝λ_m+1-λ_m，

其中，λ_m+1和λ_m表示顺序排列的特征值序列中相邻的两个特征值，Δλ_m表示相邻特征值的差，即特征值间隙；

第一个极大的所述特征值间隙所对应的下标的数值即为源信号数量，也为簇的数目。

4.如权利要求2所述的基于结构化稀疏子空间聚类的欠定盲源分离方法，其特征在于，基于所述自我表示系数矩阵，使用谱聚类进行初始类别划分，从而确定源信号数量以及每个单源点所属的簇，包括步骤：

对所述拉普拉斯矩阵，基于所述簇的数目，作为谱聚类的输入，基于规范割集准则对所述自我表示系数矩阵进行切割，得到任一单源点的簇的划分。