[发明专利]一种基于非负张量分解的多源信号分离方法

权利要求书

1.一种基于非负张量分解的多源信号分离方法，其特征在于，包括以下步骤：

步骤一：提取多源振动信号A^m×n，其中m为每个源信号数据点个数，n为信号个数；

步骤二：对于多源振动信号A^m×n，选择不同窗长下的汉宁窗对A^m×n进行短时傅里叶变换构建时频数据，得到不同窗长下的时频数据B^m×m×n，根据时频数据B^m×m×n求解时域方向和频域方向的平均信息熵，绘制不同窗长下的平均信息熵图形，选择时域方向和频域方向的平均信息熵的交点所对应的窗长作为最优窗长；

步骤三：根据最优窗长，对多源振动信号A^m×n进行短时傅里叶变换得到时频张量数据X^t×f×n，其中t、f、n分别代表时间、频率和信号个数；

步骤四：对时频张量数据X^t×f×n进行非负张量分解，根据分解过程中的收敛误差、迭代步数和核一致性，采用基于收敛误差及迭代步数的评价方法和基于核一致性的源信号个数的优化估计方法，其中基于收敛误差及迭代步数的评价方法是通过非负张量分解过程中算法的迭代步数和相邻层之间的相对误差来衡量的，而基于核一致性的源信号个数的优化估计方法通过将非负张量分解和同等规格的非负Tucker分解的内核数据结构间的相互关系进行评估，核一致性指标用公式表示为：

上式中：F—非负Tucker分解的核中所有元素的平方和；g_def∈G—非负张量分解的核元素；t_def∈T—非负Tucker分解的核元素，最终根据得到的收敛误差及迭代步数和核一致性的变化规律曲线，找到最优源信号个数r；

步骤五：根据最优源信号个数r，对时频张量数据X^t×f×n进行非负张量分解得到X₁∈R^t×r，X₂∈R^f×r，X₃∈R^n×r，由公式：

X_j＝X₁(:,i)×X₂(:,i)' (2)

其中i＝1:1:r，j＝1:1:r，X_j代表第j个源信号的重构时频矩阵特征分量，得到了r个源信号的重构时频矩阵特征分量；

步骤六：选择步骤三中得到的时频张量数据X^t×f×n中任一个信号的时频数据X^t×f×k，进行源信号的时频矩阵重构，即用公式：

Y_j＝X_j·X^t×f×k (3)

其中k＝1:1:n,j＝1:1:r，Y_j代表第j个源信号的时频矩阵，得到r个重构源信号时频矩阵，然后利用短时傅里叶逆变换公式：

其中f代表频率，t代表时间，得到r个源信号，成功地从多源振动信号中提取出源信号。