联合变分模态分解和排列熵的GNSS变形监测去噪方法技术

本发明专利技术公开了联合变分模态分解和排列熵的GNSS变形监测去噪方法，针对模态分解去噪方法中的高频噪声判别问题，首次提出通过计算各模态分量排列熵值的方法判别低频有效分量和高频噪声。此方法可将各分量的时间复杂程度量化，定量的表示各分量的特性，根据其数值大小判别高频噪声，对任意GNSS监测序列均具有自适应性。该方法可以有效的缓解EMD模态混叠、局部信息丢失以及小波去噪受小波基等外部选择影响大、不具有自适应性等问题。通过实测数据实验证明，该方法去噪精度较EMD、小波去噪有明显提高，验证了该方法具备有更好的精度和可靠性。在GNSS导航定位精度提高和空间环境的减灾和防灾能力具有重要的应用价值。和防灾能力具有重要的应用价值。和防灾能力具有重要的应用价值。

【技术实现步骤摘要】
联合变分模态分解和排列熵的GNSS变形监测去噪方法


[0001]本专利技术涉及模态分解类去噪
，具体是联合变分模态分解和排列熵的GNSS变形监测去噪方法。

技术介绍

[0002]GNSS观测具有实时获取连续三维坐标观测序列的优点，已经被广泛的应用于工程监测、大地形变、地震预测等诸多领域中。然而由于测量环境的随机性、变形过程的复杂性，监测对象往往受到各种复杂因素的复合影响，产生非平稳、高噪声的观测数据，从而导致GNSS时间序列出现相应的误差，难以准确得到监测对象实际的形变位移。使用包含噪声的数据进行变形分析预测难以达到理想的效果，甚至得出错误的结论，灾害发生时不能及时预警从而危及人民群众的生命财产安全。因此，如何有效降低测量噪声对GNSS观测数据的影响意义重大。
[0003]目前，诸多学者主要利用Kalman滤波、小波分析、经验模态分解(empirical mode decomposition,EMD)等方法对GNSS变形监测数据进行去噪处理。其中Kalman滤波在时域范围内可对监测序列进行很好的分析，但在频域范围无法对非平稳监测序列的特征做出表达，存在很大的局限性。小波变换可以同时对GNSS时间序列的频域和时域进行分析，具有良好的时频特性以及多分辨率分析的特征，在数据去噪方面具有良好的效果，但其去噪效果在很大程度上受到分解尺度、小波基选取等的影响，使之不具有自适应性，分解后的各分量也不具备物理意义。EMD方法适用于非线性、非平稳信号，应用在GNSS变形监测去噪中可有效实现随机噪声的分离并且获得较高的信噪比，近些年被广泛的应用在数据去噪中。然而通过EMD进行变换，存在模态混叠等问题，可能出现无法解释的、缺乏实际意义的频率分量。并且如何区分有效信息和噪声分量难以进行有效判别，可能会致使有效信息的丢失。
[0004]现有的模态分解类去噪方法是将GNSS变形序列看作连续的一维信号，认为所有复杂的信号均由非正弦函数的、简单的、相异的分量信号组成。其去噪理论首先将GNSS监测序列分解为不同频率的模态分量(intrinsic mode function,IMF)和趋势项，然后依据一定方法分离高频噪声，保留低频分量重构以达到去噪的目的。模态分量应该满足两个条件：一是特征模态函数曲线过零点的次数与极大、极小值点的个数必须相同或相差一个；二是由极大、极小值点构成的上下包络线均值为零，即上下包络线关于x轴对称。经验模态分解去噪方法的主要步骤如下：
[0005](1)设GNSS变形监测时间序列为x(t)，找出输入时间序列中所有的极值点；
[0006](2)采用三次样条法对时间序列中的极大值和极小值分别进行差值计算，构成上下包络线m
max
和m
min
；
[0007](3)计算上下包络线的均值m(t)，并求出其与时间序列的差值h
j
(t)；
[0008](4)判断h(t)是否为固有模态函数分量，如h(t)满足两个必要条件，则IMF
i
(t)＝h
j
(t)，然后原始序列减去IMF
i
(t)得到对残余项继续采用上述方法分解直到信号单调为止；
[0009](5)判别各模态分量的频率进行高通、低通或带通滤波得到最终去噪结果。
[0010]这种模态分解方法适用于非线性、非平稳信号，应用在GNSS变形监测去噪中，可以有效实现随机噪声的分离并且获得较高的信噪比，近些年被广泛的应用在数据去噪中。此类方法目前主要依靠相关系数法或标准化模量累计均值法分离噪声，这些方法在不同模态分解的情况下呈现出不同的结果。相关系数法以第一次出现的极小值判别高频噪声，然而在各模态分量中第一次极小值点出现的位置具有不确定性，并不能准确区分高频噪声。计算各模态分量的标准化模量累计均值，通过其累计均值突变区别高频噪声，但此方法计算得到的各模态分量数值量级变化较大，无法准确判断突变模态，仍不能准确确定高频噪声。

技术实现思路

[0011]本专利技术要解决的技术问题就是克服以上的技术缺陷，提供联合变分模态分解和排列熵的GNSS变形监测去噪方法，针对模态分解类去噪方法存在噪声判别标准不统一，无法准确区分高频噪声和有效变形信息，致使去噪不彻底或过度去噪的问题。本专利技术将各模态分量的复杂程度量化，可根据实际情况定量的准确的判别高频噪声，无需复杂的参数设置即可达到精准去噪的目的。
[0012]为了解决上述问题，本专利技术的技术方案为：联合变分模态分解和排列熵的GNSS变形监测去噪方法，具体步骤如下：
[0013](1)利用VMD算法将GNSS监测序列分解不同频率的模态分量，VMD算法中约束变分模型分为等式、不等式两种，等式约束指各IMF分量叠加后的信号依旧等于原始信号，不等式约束指各IMF分量相应的估计带宽和应为最小，具体模型如下所示：
[0014][0015]式中，δ
t
为狄拉克函数；*为卷积运算符；u
k
为通过VMD分解出的K个IMF分量，即{u
k
}＝{u1,u2,u3,
…
,u
k
}；{w
k
}为对应IMF分量的中心频率，即{w
k
}＝{w1,w2,w3,
…
,w
k
}；
[0016](2)构造增广Lagrange函数，引入惩罚因子α以及Lagrange乘子α，将约束性变分问题转化为非约束性进行求解；
[0017][0018](3)利用交替方向乘法算子法解决以上无约束变分问题，通过公式交替更新λ
n+1
，迭代进行解算，寻找Lagrange表达式鞍点，当满足一定的收敛条件时输出K个模态分量，其公式如下所示：
[0019][0020][0021]式中为迭代次数；为u(t)、λ(t)、X(t)的Fourier变换；
[0022](4)计算经VMD算法分解后的各模态分量的排列熵值，设有各模态分量时间序列为{u
t
}＝{u1,u2,u3,...,u
n
}，通过相空间重构，得到对应的重构矩阵如下：
[0023][0024]式中j＝1,2,...k；t为时间序列；m为嵌入维数；K为重构分量的个数；
[0025](5)以行向量的形式分为K个重构分量，对每个重构分量按大小进行升序排序，如若其值相等则按索引值排列：
[0026]U(i)＝{u(i+(j1‑
1)t),u(i+(j2‑
1)t),...,u(i+(j
m
‑
1)t)}
[0027]u(i+(j1‑
1)t)≤u(i+(j2‑
1)t)≤
…
≤u(i+(j
m
‑
1)t)
[0028]则获得一组新的符号序列：
[0029]S(i)＝(j1,j2,...j
m
)
本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.联合变分模态分解和排列熵的GNSS变形监测去噪方法，其特征在于，具体步骤如下：(1)利用VMD算法将GNSS监测序列分解不同频率的模态分量，VMD算法中约束变分模型分为等式、不等式两种，等式约束指各IMF分量叠加后的信号依旧等于原始信号，不等式约束指各IMF分量相应的估计带宽和应为最小，具体模型如下所示：式中，δ
t
为狄拉克函数；*为卷积运算符；u
k
为通过VMD分解出的K个IMF分量，即{u
k
}＝{u1,u2,u3,
…
,u
k
}；{w
k
}为对应IMF分量的中心频率，即{w
k
}＝{w1,w2,w3,
…
,w
k
}；(2)构造增广Lagrange函数，引入惩罚因子α以及Lagrange乘子α，将约束性变分问题转化为非约束性进行求解；(3)利用交替方向乘法算子法解决以上无约束变分问题，通过公式交替更新λ
n+1
，迭代进行解算，寻找Lagrange表达式鞍点，当满足一定的收敛条件时输出K个模态分量，其公式如下所示：所示：式中为迭代次数；为u(t)、λ(t)、X(t)的Fourier变换；(4)计算经VMD算法分解后的各模态分量的排列熵值，设有各模态分量时间序列为{u
t
}＝{u1,u2,u3,
…
,u
n
}，通过相空间重构，得到对应的重构矩阵如下：式中j＝1,2,
…
k；t为时间序列；m为嵌入维数；K为重构分量的个数；(5)以行向量的形式分为K个重构分量，对每个重构分量按大小进...

【专利技术属性】
技术研发人员：汤俊，李垠健，林海飞，董晓燕，
申请(专利权)人：华东交通大学，
类型：发明
国别省市：