[发明专利]用于转子系统的振动信号处理的信号重构方法

摘要

本发明公开了一种非均匀欠采样条件下谐波组合信号的重构方法，其中方法包括以下步骤：步骤1、获取N个大于Nyquist采样周期的欠采样周期内j个信号数值的时间序列。步骤2、根据采样数据对应的时间构造测量矩阵。步骤3、构造含测量向量、测量矩阵、傅里叶基矩阵、重构一维稀疏向量的稀疏重构方程。步骤4、对重构表达式进行正交化预操作，构造等效表达式。步骤5、计算重构方程的观测矩阵互相关系数，调整采样时间间隔，重新构造观测矩阵，使观测矩阵互相关系数最小。步骤6、利用凸优化方法、贪婪算法等典型重构算法，求解稀疏向量，获得谐波组合信号的所有频率成分，实现信号重构。本发明提供的方法重构概率大、精度高，可以在无先验信息条件下实现谐波组合信号的盲重构。

主权项

1.一种用于转子系统的振动信号处理的基于稀疏理论的谐波组合信号非均匀欠采样盲重构方法，其特征在于包括以下步骤：步骤S100：获取N个大于Nyquist采样周期的欠采样周期T_n内j个信号数值的非均匀欠采样时间序列；其中，谐波组合信号非均匀欠采样时间序列的获取；包括以下步骤：步骤S101：确定一个虚拟采样间隔T_n/L；其中：L为对每个欠采样周期T_n进行等分的份数，使得等分后的时间间隔T_n/L满足Nyquist采样定理，即虚拟采样频率不低于两倍的待测信号最高频率且各个采样数据的采样时间均对应某个等分点的时刻；步骤S102：将谐波组合信号的非均匀欠采样模型表示为公式：其中r(t)为真实信号，y[t_j]为采样结果的时间序列，δ(t)为狄里克莱函数，L为对欠采样周期T_n的切分倍数，ξ_i为一个小于L的随机正整数，决定每个周期采样的非均匀性，I为每个采样周期T_n内的采样个数，N为采样周期数；步骤S200：根据采样数据对应的时间构造测量矩阵Φ_p×q其具体包括以下步骤：步骤S201：确定测量矩阵的阶数；测量矩阵Φ_p×q的行、列数选取满足以下公式：其中p为测量矩阵的行数，与测量向量维度相同；q为测量矩阵的列数；步骤S202：测量矩阵按以下原则构造：测量矩阵Φ_p×q满足以下形式：其中，矩阵每行均有且只有一个元素为1，即其他元素均为0；该非零元素的坐标可表示为且1元素对应的列坐标k满足其中t_j为步骤S102所提公式的第j,j≤q个采样数据的采样时刻；步骤S300：构造含测量向量y、测量矩阵Φ、傅里叶基矩阵Ψ_FFT、重构一维稀疏向量θ的稀疏重构方程其中y∈R^p×1为时间序列y[t_j]的向量形式、Φ∈Z^p×q为测量矩阵、Ψ_FFT∈C^q×q为傅里叶基矩阵、θ(f)∈C^q为重构的一维稀疏向量，该向量等效可等效成真实信号的频谱；步骤S400：对S300中的稀疏重构方程进行正交化预操作，构造等效表达式z＝Qθ(f)，其中，z∈R^p和Q∈C^p×q分别是y∈R^p和的正交化预变换表示，其中：Q为观测矩阵；步骤S401：构造测量向量的正交化预变换矩阵其中，为R的伪逆矩阵，Q＝orth(R^T)^T步骤S402：构造正交化预处理后的重构方程z＝Qθ(f)，其中，z＝Ty；步骤S500:计算重构方程的观测矩阵互相关系数μ(Q)，调整采样时间间隔，重新构造等效观测矩阵Q，使观测矩阵Q互相关系数最小；其中，q_i、q_j是矩阵Q中的列向量；步骤S600:利用包括凸优化方法、贪婪算法在内的典型重构算法，求解稀疏向量θ，获得谐波组合信号的所有频率成分，实现信号重构。