-
在沉降监测中,为保证基建工程的顺利建设,使用了如液体静力水准仪、裂缝计等传感器实时监测变形体的变化。传感器所测得数据具有数据量大且非平稳等特点,传统的预测方法如灰色模型、时间序列、回归模型、神经网络都不能有效的对这类数据进行分析与预测[1-3]。目前,就这类数据而言,通常使用小波变换、经验模态分解(empirical mode decomposition, EMD)等时频分解方法对数据进行分析与预测。王新洲 等使用了小波分解和支持向量机相结合来对数据进行分析与预测[4],但小波分解不适于处理离散信号且缺乏一定的自适应性,是一种人为主观的分割;王勇等使用经验模态分解和神经网络的组合模型对监测信号进行预测[5],但经验模态分解会出现过包络、端点效应、欠包络、模态重叠等问题,会对预测模型产生较大影响。任超等使用了集合经验模态分解(ensemble empirical mode decomposition, EEMD)和支持向量机的组合模型对沉降监测数据进行分析与预测[6]。集合经验模态分解很好的抑制了模态混叠问题的出现,但对于3次样条差值本身可能出现的端点效应等问题并没有解决。
作为一种新的时频分解方法局部均值分解(local mean decomposition,LMD)很好地解决了经验模态中过包络、欠包络、端点效应的问题[7],并且在沉降监测中已经得到了应用[8-10]。但在使用中无法抑制模态混叠带来的影响,同样会导致预测过程中出现较大偏差,整体适用性上存在不足。
总体局部均值分解(ensemble local mean decomposition, ELMD)法可解决上述问题。该方法通过对原始信号多次加入白噪声并采用LMD分解,取多次分解后的平均值作为分解值,抑制了模态混叠现象,并且兼备LMD方法的优点。本文基于ELMD设计了一种并联式组合模型,通过实验证明其正确性,并以实际监测数据处理结果证明其优势。
-
任意原始信号x(t),其LMD分解过程归纳如下。
1) 确定信号x(t)的所有局部极值点ni,计算两个相邻极值点ni和ni+1的平均值mi和包络估计值ai:
$$ {m_i} = \frac{{{n_i} + {n_{i + 1}}}}{2}, i = 1, \ldots, n $$ (1) $$ {a_i} = \left| {\frac{{{n_i} + {n_{i + 1}}}}{2}} \right|, i = 1, \ldots, n $$ (2) 2) 用折线将所有相邻mi连接起来,用滑动平均法进行处理,得到局部均值函数m11(t)。用同样方法得到包络估计函数a11(t)。
3) 从原始信号中分离出局部均值函数m11(t),得到:
$$ {h_{11}}\left( t \right) = x\left( t \right) - {m_{11}}\left( t \right) $$ (3) 4) 对h11(t)进行解调:
$$ {s_{11}}\left( t \right) = \frac{{{h_{11}}\left( t \right)}}{{{a_{11}}\left( t \right)}} $$ (4) 理想情况下,s11(t)是一个纯调频信号,即它的包络估计函数a12(t)满足a12(t)=1。如果不满足该情况,则需重复上述过程,直到s1n(t)满足1-Δ≤s1n(t)≤1+Δ后停止迭代。有:
$$ \left\{ \begin{array}{l} {h_{11}}\left( t \right) = x\left( t \right) - {m_{11}}\left( t \right)\\ {h_{12}}\left( t \right) = {s_{11}}\left( t \right) - {m_{12}}\left( t \right)\\ \;\;\;\;\; \vdots \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; \vdots \\ {h_{1n}}\left( t \right) = {s_{1\left( {n - 1} \right)}}\left( t \right) - {m_{1n}}\left( t \right) \end{array} \right. $$ (5) $$ \left\{ {\begin{array}{*{20}{c}} {{s_{11}}\left( t \right) = \frac{{{h_{11}}\left( t \right)}}{{{a_{11}}\left( t \right)}}}\\ {{s_{12}}\left( t \right) = \frac{{{h_{12}}\left( t \right)}}{{{a_{12}}\left( t \right)}}}\\ { \vdots \;\;\;\;\;\;\;\;\;\;\;\;\; \vdots }\\ {{s_{1n}}\left( t \right) = \frac{{{h_{1n}}\left( t \right)}}{{{a_{1n}}\left( t \right)}}} \end{array}} \right. $$ (6) 迭代终止条件为$ \mathop {\lim }\limits_{n \to \infty } {a_{1n}}\left( t \right) = 1 $。实际应用中,为了减少运算迭代次数,可设置微小变量Δ,满足1-Δ≤a1n≤1+Δ时,迭代终止。
5) 把迭代过程中所产生的所有包络估计函数相乘得到瞬时幅值函数(包络信号):
$$ {a_1}\left( t \right) = {a_{11}}\left( t \right){a_{12}}\left( t \right) \ldots {a_{1n}}\left( t \right) = \prod\limits_{q = 1}^n {{a_{1q}}\left( t \right)} $$ (7) 6) 将纯频信号s1n(t)和包络信号a1(t)相乘可得到第一个PF分量PF1(t):
$$ P{F_1}\left( t \right) = {s_{1n}}\left( t \right){a_1}\left( t \right) $$ (8) 从原始信号中将PF1(t)分离出来,可得到一个时间序列信号u1(t),将u1(t)作为原始信号并重复上述过程,直到uk(t)为一个单调函数后停止。
$$ \left\{ \begin{array}{l} {u_1}\left( t \right) = x\left( t \right) - P{F_1}\left( t \right)\\ {u_2}\left( t \right) = {u_t}\left( t \right) - P{F_2}\left( t \right)\\ \;\;\;\;\; \vdots \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; \vdots \\ {u_k}\left( t \right) = {u_{k - 1}}\left( t \right) - P{F_k}\left( t \right) \end{array} \right. $$ (9) -
模态混叠是指原始信号中某频段存在断续或者畸变的情况下,不同频率分量被分解到同一PF分量中或同一频率分量被分解到不同的PF分量中,由此提出ELMD方法来抑制模态混叠问题。首先将白噪声加入到待分解的信号中,利用白噪声均匀地“污染”原信号,之后采用LMD的方法对混入噪声的信号进行分解。由于在分解的过程中会自动地将信号分解到白噪声所确定滤波器组结构当中,因此模态混叠问题会得到抑制。解决白噪声的存在引起分量失真的情况时利用白噪声所具有的统计特性,在多次分解后取分解结果的平均值来抵消白噪声在分量中的残留部分,便能保留真实的分解结果。综上所述,ELMD方法的分解过程归纳如下[11]。
1) 向一组原始的时间序列信号{Y(t)}中加入一组已知的白噪声{nm(t)}。得到加躁后的时间序列信号为:
$$ {Z_m}\left( t \right) = Y\left( t \right) + {n_m}\left( t \right) $$ (10) 2) 对加噪后的时间序列信号Z(t)使用LMD进行分解,得到若干乘积函数PF以及一个剩余分量uk(t)。
3) 多次重复步骤1)、2), 每次参杂不同幅值的白噪声信号。
4) 对每次得到的PF分量和剩余分量uk(t)作总体平均,其均值作为最终分解结果。
-
离散系统的卡尔曼滤波模型的状态方程和观测方程分别为[12]:
$$ \left\{ \begin{array}{l} {\mathit{\boldsymbol{X}}_{k + 1}} = {\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_{k + 1, k}}{\mathit{\boldsymbol{X}}_k} + {\mathit{\boldsymbol{\psi }}_{k + 1, k}}{\mathit{\boldsymbol{U}}_{k + 1}} + {\mathit{\boldsymbol{ \boldsymbol{\varOmega} }}_{k + 1}}\\ {\mathit{\boldsymbol{L}}_{k + 1}} = {\mathit{\boldsymbol{B}}_{k + 1}}{\mathit{\boldsymbol{X}}_{k + 1}} + {\mathit{\boldsymbol{ \boldsymbol{\varDelta} }}_{k + 1}} \end{array} \right. $$ (11) 式中,Xk、Uk、Lk分别为tk时刻的状态向量、控制向量以及观测向量;Φk+1, k为tk时刻至tk+1时刻状态转移矩阵;Ψk+1, k为tk+1时刻控制矩阵;Bk+1为tk+1时刻观测矩阵;Ωk+1和Δ k+1分别为tk+1时刻的动态噪声和观测噪声。
最终可推出其预报方程为:
$$ \left\{ \begin{array}{l} \mathit{\boldsymbol{\hat X}}\left( {k + i/k} \right) = {\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_{k + i,k}}\mathit{\boldsymbol{\hat X}}\left( {k/k} \right) + \sum\limits_{r = k + 1}^{k + i} {{\mathit{\boldsymbol{\psi }}_{k + i,r}}} {\mathit{\boldsymbol{U}}_r}\\ {\mathit{\boldsymbol{D}}_X}\left( {k + i/k} \right) = {\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_{k + i,k}}{\mathit{\boldsymbol{D}}_X}\left( {k/k} \right) + \mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_{k + i,k}^T + \\ \;\;\sum\limits_{r = k + 1}^{k + i} {{\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_{k + i,r}}} {\mathit{\boldsymbol{ \boldsymbol{\varGamma} }}_{r,r - 1}}{\mathit{\boldsymbol{D}}_\Omega }\left( {r - 1} \right)\mathit{\boldsymbol{ \boldsymbol{\varGamma} }}_{r,r - 1}^T\mathit{\boldsymbol{ \boldsymbol{\varPhi} }}_{k + i,r}^T \end{array} \right. $$ (12) 式中,Φj, k = Φj, i Φi, k
-
1) 一组原始的时间序列信号{Y(t)},使用ELMD得到n个PF分量和一个剩余分量uk(t):
$$ X\left( t \right) = P{F_1} + P{F_2} + P{F_3} + \cdots + P{F_n} + {u_k}\left( t \right) $$ (13) 2) 使用支持向量机对PF分量进行预测,使用卡尔曼滤波算法对剩余分量uk(t)进行预测, 所得预测时间序列分别为:F1, F2, F3, …, Fn, Fu
3) 预测组合模型如图 1所示,对所得预测值进行重构得到最终预测值W。
$$ W = {F_1} + {F_2} + {F_3} + \cdots + {F_n} + {F_u} $$ (14) -
为了评价所构建的并联式组合模型的预测精度,采用如下统计量对模型预测效果进行评估:
平均相对百分比误差为:
$$ {\rm{MAPE}} = \frac{1}{n}\sum\limits_{t = 1}^n {\frac{{\left| {{X_t} - {{\hat X}_t}} \right|}}{{\left| {{X_t}} \right|}}} $$ (15) 均方误差为:
$$ {\rm{RMSE}} = \sqrt {\frac{1}{n}\sum\limits_{t = 1}^n {{{\left( {{X_t} - {{\hat X}_t}} \right)}^2}} } $$ (16) 式中,Xt为实测值;$ {{{\hat X}_t}} $为预测值。
-
产生一组仿真信号:
$$ \begin{array}{l} x\left( t \right) = {x_1}\left( t \right) + {x_2}\left( t \right) + {x_3}\left( t \right)\\ \;\;\;\;\;\;\;\;\;\;t \in \left( {0, 1200} \right) \end{array} $$ (17) 式中,
$$ \begin{array}{l} {x_1}\left( t \right) = \sin \left( {2{\rm{\pi }}t/1200} \right), \\ {x_2}\left( t \right) = \sin \left( {2{\rm{\pi }}t/600} \right), \\ {x_3}\left( t \right) = \left\{ \begin{array}{l} 0.2\sin \left( {2{\rm{\pi }}t/50} \right), t \in \left[{200, 300} \right] \cup \\ \;\;\left[{600, 700} \right] \cup \left[{1\;000, 1\;100} \right]\\ 0, t \in {\rm{else}} \end{array} \right. \end{array} $$ 上述仿真信号采样率为1 s, 时长为1 200 s。其序列图如图 2所示。
首先对仿真信号分别进行LMD和ELMD分解。LMD分解结果如图 3所示,得到2个PF分量以及一个剩余分量。在PF1分量中,包含两种不同频率的信号,出现了明显的模态混叠。ELMD分解结果如图 4所示,各阶PF分量的频率尺度统一。在分解过程中,加入了幅值0.1的白噪声,进行500次总体平均,将ELMD所求分量进行重构,与原始信号求差得到残差分量如图 5所示,其标准偏差不足0.01。
-
使用液体静力水准仪对下穿高速铁路某桥墩进行监测得到103期沉降时间序列,如图 6所示。其中,仪器标准量程为100 mm, 精度±0.1%F.S, 灵敏度0.025%F.S,采样间隔3 h。所得时间序列中,实际上是由自由振动响应和外荷载引起的强迫振动响应两部分组成。其中,外荷载的强迫振动主要是由风致振动、机械及交通荷载、水压力等因素的综合作用。
本次实验以前93期作为建模数据,后10期作为预测数据。使用ELMD和LMD进行对比实验,分析测量数据某频率信号不连续情况下的模拟和预测情况。在使用ELMD进行分解的过程中,对建模数据分别加入500次均值为零,方差为0.1的高斯白噪声,取500次分解结果的平均值作为最终分解结果。ELMD与LMD分解所得分量情况如图 7所示。
对比图 7(a)、7(b),由于在LMD分解过程中发生了模态混叠,使得ELMD分解结果要优与LMD。对于LMD所得到的剩余分量uk(t), 在图 7(a)中可以发现,剩余分量uk(t)中(红圈标出)残留了PF3分量的部分信号。
分析ELMD方法所得3个PF分量和1个剩余分量uk(t)。对于高频分量PF1认为是测得噪声部分。对于中频分量PF2、PF3具有一定的趋势性以及部分随机性,主要是由于温度、桥墩自身摆动等因素引起的。对于低频分量uk(t),认为是桥墩的实际沉降量,主要是由于基坑开挖以及区域降水等因素对桥墩造成影响。
对LMD以及ELMD所得分量进行预测,由于各PF分量是具有一定周期性的离散非线性数据,SVM算法在计算过程中,采用结构风险最小原则,能够实现非线性向高维线性转变后支持向量回归模型的建立,能够较好预测此类数据[13]。另外,卡尔曼滤波是一种有效的非线性动力学系统跟踪预测算法,因此本文应用常速卡尔曼滤波方法对非线性分量uk(t)进行预测。根据仪器精度指标,观测噪声取DΔ(k)=±0.1 mm,动态噪声主要包含加速度,可取为DΩ(k)=4·Δt-4·DΔ(k),另外,计算时取X(0/0) 变形值为第一期观测值,初使变形速度为零。预测结果以及重构结果如图 8、图 9所示。
对比分析图 8、图 9可以发现:由于发生了模态混叠,对LMD所得各PF分量和剩余分量uk(t)的预测值都产生了不同程度的影响。剩余分量占整体信号比重较大,两种模型的最终预测值主要是由剩余分量uk(t)的预测值决定的。一般认为低频分量最能反映出变形体的实际变化规律,所以当剩余分量中发生了模态混叠的情况会严重干扰整个模型分析预计预测能力。
最终预测值的预测精度评定如表 1所示。
表 1 模型预测精度评定表
Table 1. Comparison of Prediction Errors
期次 ELMD相对误差 LMD相对误差 94 10.069 15.967 95 4.232 6.020 96 4.252 9.919 97 2.755 5.281 98 3.828 6.521 99 8.660 11.131 100 21.497 21.851 101 16.633 18.142 102 8.379 11.569 103 3.361 5.417 MAPE 8.367 11.182 RMSE 0.149 0.213 在发生模态混叠的情况下,ELMD的预测精度高于LMD。该方法能自适应抑制LMD方法中的模态混叠现象的出现,并且不损失LMD方法的分解精度。抑制模态混叠后,本文设计的并联式方案发挥了时频分解以及传统分析模型的优点,因此该方案能更好的适用于沉降监测数据。
-
本文结合仿真分析说明了ELMD方法能够抑制模态混叠现象,并设计一种基于ELMD的并联式组合模型沉降分析方法,该方法可以分析出桥梁振动响应信息,在整体适用性上要优于基于LMD的并联式组合模型。实验结果表明,能以较高的精度对离散非线性信号进行模拟和预报,可为精密监测提供一定的参考。
值得注意的是,白噪声方差过大会被当作高频分量单独出现。另外,ELMD分解后的重构信号与原始信号之间还是存在少量未被抵消的高斯白噪声。因此针对沉降监测数据使用ELMD时白噪声方差大小以及分解次数还需要进一步研究。
The Feasibility Study on Settlement Monitoring of a Parallel Combination Prediction Method Based on ELMD
-
摘要: 时频分解方法局部均值分解(local mean decomposition,LMD)在沉降监测中已经得到了应用,但在使用中会出现模态混叠现象。总体局部均值分解(ensemble local mean decomposition,ELMD)通过添加辅助噪声可以抑制局部均值分解过程中出现的模态混叠现象。提出了一种基于ELMD的并联式组合沉降预测方法,结合高速铁路某桥梁实际监测数据,在对ELMD模型进行仿真分析的基础上,分别使用ELMD和LMD将一组离散非线性信号分解为3个PF分量和1个剩余分量,并利用支持向量机和卡尔曼滤波进行预测验证。结果表明:使用ELMD进行分解的过程中能够很好地抑制LMD方法中出现的模态混叠问题。在预报精度方面,基于ELMD的并联式组合模型的平均相对误差可以达到8.3%,可为沉降监测的预报工作提供参考和借鉴。Abstract: In the field of time-frequency decomposition, the Local Mean Decomposition(LMD) method is applied in settlement monitoring, but the phenomenon of mode mixing can appear during the application, which results in inaccurate deformation signal extraction.The Ensemble Local Mean Decomposition(ELMD) method can be used to improve the mode of mixing the local mean decomposition by adding auxiliary noise to the original signal, and also can use the statistical characteristics of auxiliary noise to remove the mode mixing. This paper uses simulation data to analyze the model error in the ELMD method and presents a parallel combination prediction method based on ELMD. In the case of high speed railway bridge monitoring data, it divides a series of discrete nonlinear and unstable signal into three product function(PF) components and one remaining component. The method takes advantage of the support vector machine and Kalman filter algorithms to predict these components, and analyses the superiority of ELMD in the case of mode mixing and overall feasibility empirically. The results indicate that: the parallel combination model, based on ensemble local mean decomposition (ELMD), can eliminate the mode mixing problem in the local mean decomposition (LMD) method very well and extracts the deformation signal accurately. In terms of prediction precision, the mean relative error can reach 8.3%, and may provide areference for prediction of deformation monitoring.
-
表 1 模型预测精度评定表
Table 1. Comparison of Prediction Errors
期次 ELMD相对误差 LMD相对误差 94 10.069 15.967 95 4.232 6.020 96 4.252 9.919 97 2.755 5.281 98 3.828 6.521 99 8.660 11.131 100 21.497 21.851 101 16.633 18.142 102 8.379 11.569 103 3.361 5.417 MAPE 8.367 11.182 RMSE 0.149 0.213 -
[1] 陈国栋, 黄甜.灰色卡尔曼滤波在沉降监测中的应用[J].测绘科学技术学报, 2014, 31(2):133-135 http://www.cnki.com.cn/Article/CJFDTOTAL-JFJC201402006.htm Chen Guodong. Huang Tian, Application of Gray Kalman Filtering in the Settlement Monitoring[J].Journal of Geomatics Science and Technology. 2014, 31(02):133-135 http://www.cnki.com.cn/Article/CJFDTOTAL-JFJC201402006.htm [2] 李晓杰, 郭睿, 黄金, 等.神经网络在北斗导航卫星轨道预报中的应用[J].武汉大学学报·信息科学版, 2015, 40(9):1 253-1 258 http://ch.whu.edu.cn/CN/abstract/abstract3330.shtml Li Xiaojie, Guo Rui, Huang Jin, et al. Application of Artificial Neural Network to Orbit Prediction of BeiDou Navigation Satellites[J]. Geomatics and Information Science of Wuhan University, 2015, 40(9): 1 253-1 258 http://ch.whu.edu.cn/CN/abstract/abstract3330.shtml [3] 许昌, 岳东杰, 董育烦, 等.基于主成分和半参数的大坝变形监测回归模型[J].岩土力学, 2011, 34(12):3 738-3 742 http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201112032.htm Xu Chang, Yue Dongjie, Dong Yufan, et al. Regression Model for Dam Deformation Based on Principal Component and Semi-parametric Analysis[J].Rock and Soil Mechanics, 2011, 34(12):3738-3742 http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201112032.htm [4] 王新洲, 范千, 许承权, 等.基于小波变换和支持向量机的大坝变形预测[J].武汉大学学报·信息科学版, 2008, 33(5):469-471+507 http://ch.whu.edu.cn/CN/abstract/abstract1592.shtml Wang Xinzhou, Fan Qian, Xu Chengquan, et al. Dam Deformation Prediction Based on Wavelet Transform and Support Vector Machine[J].Geomatics and Information Science of Wuhan University, 2008, 33(5):469-471+507 http://ch.whu.edu.cn/CN/abstract/abstract1592.shtml [5] 王勇, 杨晶, 张立辉, 等.基于经验模态分解与神经网络的信号预测[J].大地测量与地球动力学, 2011, 31(6):121-123+135 http://youxian.cnki.com.cn/yxdetail.aspx?filename=JSCH201705003&dbname=CJFDPREP Wang Yong, Yang Jing, Zhang Lihui, et al. Analysis of Signal Prediction Based on EMD and ANN[J]. Journal of Geodesy and Geodynamics, 2011, 31(6):121-123+135 http://youxian.cnki.com.cn/yxdetail.aspx?filename=JSCH201705003&dbname=CJFDPREP [6] 任超, 梁月吉, 庞光锋, 等.基于EEMD的大坝变形多步预测方法研究[J].大地测量与地球动力学, 2015, 35(5):816-820 http://www.cnki.com.cn/Article/CJFDTOTAL-DKXB201505023.htm Ren Chao, Liang Yueji, Pang Guangfeng, et al. A Novel Multi-Step Prediction for Dam Deformation Prediction Based on EEMD[J]. Journal of Geodesy and Geodynamics, 2015, 35(5):816-820 http://www.cnki.com.cn/Article/CJFDTOTAL-DKXB201505023.htm [7] 张亢. 局部均值分解方法及其在旋转机械故障诊断中的应用研究[D]. 长沙: 湖南大学, 2012 Zhang Kang. Research on Local Mean Decomposition Method and Its Application to Rotating Machinery Fault Diagnosis[D]. Changsha: Hunan University, 2012 [8] 董正坤, 罗亦泳, 蒋斯斯, 等.基于局域均值分解和BP神经网络的大坝变形预测[J].测绘工程, 2015, 24(4):70-73 http://www.cnki.com.cn/Article/CJFDTOTAL-CHGC201504016.htm Dong Zhengkun. Luo Yiyong, Jiang Sisi, et al. Dam Deformation Prediction Based on LMD-BP Network Modeling[J].Engineering of Surveying and Mapping. 2015, 24(4):70-73 http://www.cnki.com.cn/Article/CJFDTOTAL-CHGC201504016.htm [9] 李胜, 韩永亮, 杨宏伟, 等.露天矿边坡变形的LMD-Elman时序滚动预测研究[J].中国安全科学学报, 2015, 25(6):22-28 http://www.cnki.com.cn/Article/CJFDTOTAL-ZAQK201506004.htm Li Sheng, Han Yongliang, Yang Hongwei, et al. Research on LMD-Elman-Based Time-serise Rolling Prediction of Slope Deformation in Open-pitmine[J]. China Safety Science Journal, 2015, 25(6):22-28 http://www.cnki.com.cn/Article/CJFDTOTAL-ZAQK201506004.htm [10] 韩永亮, 李胜, 杨宏伟, 等.基于LMD-BA-ELM的边坡变形时序非线性预测模型研究[J].中国安全科学学报, 2015, 25(9):59-65 http://www.cnki.com.cn/Article/CJFDTOTAL-ZAQK201509011.htm Han Yongliang, Li Sheng, Yang Hongwei, et al. Research on LMD-BA-ELM-Based Model for Nonlinear Prediction of Slope Deformation Time-series. [J], China Safety Science Journal, 2015, 25(9):59-65 http://www.cnki.com.cn/Article/CJFDTOTAL-ZAQK201509011.htm [11] 程军圣, 张亢, 杨宇.基于噪声辅助分析的总体局部均值分解方法[J].机械工程学报, 2011, 47(03):55-62 http://www.cnki.com.cn/Article/CJFDTOTAL-JXXB201103010.htm Cheng Junsheng, Zhang Kang, Yang Yu.Ensemble Local Mean Decomposition Method Based on Noise-assisted Analysis[J].Journal of Mechanical Engineering, 2011, 47(03):55-62 http://www.cnki.com.cn/Article/CJFDTOTAL-JXXB201103010.htm [12] 崔希璋, 於宗俦, 陶本藻, 等, 广义测量平差[M].武汉:武汉大学出版社, 2009 Cui Xizhang, Yu Zongchou, Tao Benzao, et al.Generalize Surveying Adjustment[M] Wuhan: Wu University Press, 2009 [13] 张豪, 许四法.基于经验模态分解和遗传支持向量机的多尺度大坝变形预测[J].岩石力学与工程学报, 2011, 30(S2):3 681-3 688 http://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2011S2042.htm Zhang Hao, Xu Sifa.Multi-scale Dam Deformation Prediction Based on Empirical Mode Decomposition and Genetic Algorithm for Support Vector Machines[J] Chinese Journal of Rock Mechanics and Engineering.2011, 30(S2): 3 681-3 688 http://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2011S2042.htm -