一种高维的振动信号的降维处理方法技术

本发明专利技术公开了一种高维的振动信号的降维处理方法，其通过计算信号与信号之间的欧氏距离获得信号的近邻矩阵，再根据信号的近邻矩阵，并利用稀疏约束条件获得信号的重构权值矩阵，最后利用信号的重构权值矩阵获得降维后的振动信号，降维过程简单；其在利用稀疏约束条件获得信号的重构权值矩阵的过程中，在稀疏约束条件中引入了L1范数，使得重构权值矩阵具有很好地稀疏性，有效地剔除了噪声点的影响，提高了抗噪声能力，从而保障了本发明专利技术方法的鲁棒性；其在最后获取降维后的振动信号时，是求解一个稀疏、对称、半正定的矩阵的特征向量，因此可以降低本发明专利技术方法的计算复杂度。

【技术实现步骤摘要】

【技术保护点】
一种高维的振动信号的降维处理方法，其特征在于包括以下步骤：①假定待处理的高维的振动信号为x，以矩阵形式表示为x=[x1?x2?…?xi?…?xN]，其中，x∈Rq×N，R表示全体实数集，此处N表示x中包含的振动信号的条数，q表示x中的每条振动信号的维数，1≤i≤N，x1表示x中的第1条振动信号，x2表示x中的第2条振动信号，xi表示x中的第i条振动信号，xi∈Rq×1，xN表示x中的第N条振动信号；②将x中当前待处理的第i条振动信号xi定义为当前信号；③计算当前信号与x中的每条振动信号之间的欧氏距离，将当前信号与x中的第j条振动信号xj之间的欧氏距离记为dij，其中，1≤j≤N；④从当前信号与x中的每条振动信号之间的欧氏距离中，找出除当前信号与自身之间的欧氏距离外的K个最小的欧氏距离，将这K个最小的欧氏距离各自对应的振动信号均作为当前信号的近邻信号，其中，2≤K≤min(d+20,N×15%)，min()为取最小值函数，d表示x待降到的维数，N表示x中包含的振动信号的条数；⑤根据当前信号的K条近邻信号，构建当前信号的近邻矩阵，记为Zi，Zi=[xi(1)?xi(2)?…?xi(k)?…?xi(K)]，其中，Zi∈Rq×K，1≤k≤K，xi(1)表示当前信号的第1条近邻信号，xi(2)表示当前信号的第2条近邻信号，xi(k)表示当前信号的第k条近邻信号，xi(k)∈Rq×1，xi(K)表示当前信号的第K条近邻信号，xi(1)对应的欧氏距离≤xi(2)对应的欧氏距离≤……≤xi(k)对应的欧氏距离≤……≤xi(K)对应的欧氏距离；⑥根据当前信号及当前信号的近邻矩阵Zi，计算当前信号的重构权值矩阵，记为wi，wi为满足如下稀疏约束条件的最小解：且wi满足：wi中的所有重构权值的和为1，其中，wi∈RK×1，wi的维数为K，||xi?Ziwi||2表示xi?Ziwi的L2范数，||wi||1表示wi的L1范数，λ为正则化参数，λ≥0；⑦令i“=i+1,i=i“，将x中下一条待处理的振动信号作为当前信号，然后返回步骤③继续执行，直至x中的N条振动信号均处理完毕，得到x中的N条振动信号各自的重构权值矩阵，然后将这N个重构权值矩阵合并成x的重构权值矩阵，记为W，W=[w1?w2?…?wi?…?wN]，其中，i“的初始值为0，i“=i+1,i=i“中的“=”为赋值符号，W∈RK×N，w1表示x中的第1条振动信号x1的重构权值矩阵，w2表示x中的第2条振动信号x2的重构权值矩阵，wi表示x中的第i条振动信号xi的重构权值矩阵，wN表示x中的第N条振动信号xN的重构权值矩阵；⑧令M=(I?W)T(I?W)，然后利用MATLAB中的eigs函数获取M的所有特征值及每个特征值对应的特征向量，接着从M的所有特征值中提取出所有非零特征值，并对这些非零特征值按从小到大的顺序进行排列，再选取出值最小的d个非零特征值，最后由选取出的d个非零特征值对应的特征向量构成x降维处理后的低维振动信号Y，其中，I为d阶单位矩阵，(I?W)T为(I?W)的转置矩阵，符号“”表示向下取整符号。FDA00003185534400011.jpg,FDA00003185534400012.jpg,FDA00003185534400022.jpg,FDA00003185534400023.jpg...

【技术特征摘要】

【专利技术属性】
技术研发人员：叶庆卫，孙洋，周宇，王晓东，
申请(专利权)人：宁波大学，
类型：发明
国别省市：