基于经验模态分解的桥梁应变信号降噪方法
【专利摘要】本发明提出基于经验模态分解法的桥梁应变信号降噪方法，首先通过桥梁关键部位上设置的应变信号采集单元进行数据采集，提取出桥梁关键部位的应变信号；其次是对提取出的应变信号做经验模态分解EMD处理，具体分析经验模态分解得到的本征模函数IMF分量，对IMF分量中含噪较大的分量进行相对应的阈值滤波；最后将通过阈值滤波后的IMF分量与原始IMF分量中含噪较小的IMF分量以及残余分量一起进行EMD信号重构，从而得到降噪后的桥梁应变信号。
【专利说明】
基于经验模态分解的桥梁应变信号降噪方法 一、
技术领域
[0001] 本发明涉及一种信号降噪方法，涉及到桥梁结构健康监测，具体是一种基于经验 模态分解的桥梁应变信号降噪方法。 二、
【背景技术】
[0002] 桥梁结构健康检测在桥梁的管理和养护中发挥着不可替代的作用。结构健康监测 数据是分析桥梁结构损伤以及外界激励作用的重要参考。而获得结构健康监测数据的一个 重要途径是则通过分析桥梁各关键部位的应变数据得到。但是我们通过结构健康监测系统 得到的应变数据中总是不可避免的包含各种各样的噪声，甚至在某些时刻，噪声会淹没桥 梁应变信号的原始特征，如果我们直接使用通过结构健康监测系统采集到的应变数据做结 构健康分析，则很可能导致我们得到的结构健康检测数据不准确，从而造成重大的安全隐 患。而现有的降噪方法在处理桥梁应变信号的时候，总是不能完美地将信号中的噪声去除。
[0003] 采用基于经验模态分解的阈值滤波来对桥梁应变信号进行降噪处理，其能够在保 留原始应变数据基本特征的前提下，将其中的噪声基本去除。我们可以根据不同IMF分量的 时频特性，设计相对应的阈值滤波器对其进行处理。最后将处理后的IMF分量与未处理的 頂F分量以及残余分量一起进行信号重构，得到降噪后的应变信号。该方法对于确保桥梁结 构健康检测的准确性有着积极作用，同时对桥梁的管理与维护也有着不可或缺的作用。 三、
【发明内容】
[0004] 本发明的目的在于，提供一种方便高效的桥梁应变信号降噪方法，该方法可以将 桥梁应变信号中的噪声基本去除，从而确保通过分析桥梁应数据得到的桥梁结构健康监测 数据的准确性。该方法还可广泛用于大型建筑、交通隧道等会常见的重大基础设施的应变 数据的降噪处理。
[0005] 本发明的目的是这样实现的：基于经验模态分解法的桥梁应变信号降噪方法，首 先通过桥梁关键部位上设置的应变信号采集单元进行数据采集，提取出桥梁关键部位的应 变信号;其次是对提取出的应变信号做经验模态分解EMD处理，具体分析经验模态分解得到 的本征模函数IMF分量，对IMF分量中含噪较大的分量进行相对应的阈值滤波;最后将通过 阈值滤波后的IMF分量与原始IMF分量中含噪较小的IMF分量以及残余分量一起进行EMD信 号重构，从而得到降噪后的桥梁应变信号。
[0006] 所述的经验模态分解EMD处理，是指对桥梁应变信号进行从高频段到低频段的分 解，其步骤如下:首先，找出通过桥梁上设置的应变信号采集单元采取的关键部位的原始应 变信号X(t)的所有极大值点和极小值点，原始应变信号通过应变信号采集单元取得、取得 的极大值和极小值点指需要降噪信号的整个时间段内的;并用三次样条函数对所有极大值 点和极小值点进行拟合，分别对应得到上包络线m(t)和下包络线u 2(t)，并求出"仂与仍 (t)的均值&1(〇，即&1(〇 = (1!1(〇+1!2(〇)/2;接着用原始信号乂(〇减去&1(〇后得到一个除 去低频的新信号nu(t)，即nuUhXUhadt);通常情况下nu(t)并不满足本征模函数頂F分 量的要求，即信号的极大值点与极小值点的数目和过零点的数目相等或最多相差一个以及 由局部极大值构成的上包络线和由局部极小值构成的下包络线的平均值为零。因此对Hl1 (t)重复进行上述操作，即找出m(t)的所有极大值与极小值点，然后求出两者平均a2(t)，再 用m(t)减去a 2(t)得到m2(t)，如果m2(t)仍然不满足MF分量的要求，则继续重复上述操作， 直到得到的某一级m n(t)符合頂F分量的要求，这样就得到了第一级MF分量C1(t)，C1(t)表 示原始应变信号中频率最高的部分;接下来用X(t)减去 C1(t)得到n(t)，并对n(t)重复上 述操作得到第二级頂F分量c2(t);不断重复上述步骤，得到其余所有的頂F分量以及残余分 量;桥梁应变原始信号X(t)经过EMD分解后可以用如下的公式表示：
[0007]
[0008] 其中X(t)代表桥梁应变原始信号，Cl(t)代表所有的IMF分量，rn(t)代表残余分量。 取得的极大值和极小值点指需要降噪信号的整个时间段内的。
[0009] 所述的阈值通过如下公式计算：
[0010]
[0011] 其中var为頂F分量的方差，L为頂F分量的长度，ThrValue为每级頂F分量所对应的 阈值。
[0012] 所述的阈值滤波是指对含噪大的IMF分量进行相对应的阈值滤波，最常见的一种 通过如下公式来实现：
[0013]
[0014] 其中，cn(t)为原始各级頂F分量，ThrValue为各级頂F分量对应的阈值，average为 各级IMF分量的平均值，q (t)为阈值滤波后的各级IMF分量。
[0015] 所述的EMD信号重构是指将经过阈值滤波后的MF分量，含噪较小的MF分量以及 残余分量进行信号重构，从而得到降噪后的桥梁应变信号，通过如下公式计算实现：
[0016]
[0017]其中ci(t)代表未经过阈值滤波的MF分量，即原始信号经验模态分解后含噪较小 的頂F分量，％⑴代表经过阈值滤波后的MF分量，rn(t)为残余分量，又⑴代表降噪后的桥 梁应变信号。
[0018] 本发明的有益效果包括：
[0019] (1)采用结构健康监测系统的实测数据，从实际数据中得出桥梁应变数据的独特 性，解决了模拟数据的不精确问题；
[0020] (2)采用经验模态分解法对桥梁应变信号进行处理，只对含噪较大的MF分量进行 处理，解决了一般降噪方法对整个信号进行处理而带来的桥梁应变信号原始特征丢失问 题；
[0021] (3)采用针对各頂F分量设计专门的阈值滤波器对其进行滤波，解决了使用同一种 阈值滤波方法对不同的IMF分量进行滤波而带来的IMF分量原始特征丢失问题；
[0022] (4)通过本发明创造，可以方便高效地对桥梁应变信号进行降噪处理，从而提高桥 梁结构健康监测数据的准确性和可靠性；
[0023] (5)可根据本发明对建筑物以及其他基础设施的应变信号进行降噪处理；
[0024] (6)通过对桥梁应变信号做检验模态分解，结合特定的针对桥梁应变信号的阈值 滤波器，提供一种效率高，降噪效果好的桥梁应变信号降噪方法。该方法可以针对桥梁关键 部位应变信号的特点，对其进行降噪处理，基本上消除桥梁应变原始信号中的噪声。该方法 还可以广泛应用于大型建筑、交通隧道等常见的重大基础设施应变信号的降噪处理。 四、
【附图说明】
[0025] 图1是本发明所述结构健康监测系统采集到的桥梁应变数据；
[0026] 图2是本发明的系统工作流程图；
[0027] 图3是桥梁应变数据进行经验模态分解得到的頂F分量以及残余分量图；
[0028] 图4是对含噪较大MF分量进行相对应阈值滤波前后对比图；a、b分别对含噪较大 IMF分量进行相对应阈值滤波前后对比图。
[0029] 图5是进行EMD信号重构后得到的降噪信号(b)与原始桥梁应变信号对比图(a)。 五、
【具体实施方式】
[0030] 下面结合附图和具体实例对本发明做详细阐述，但本发明应用范围不限于此例。
[0031] 结合图1为某大桥结构健康监测系统测到的大桥箱梁内某结点的应变数据，可以 看到原始应变数据中含有较大的噪声，甚至在某些时刻，噪声淹没了信号的原始特征。
[0032] 结合图2,为本发明的系统工作流程图，其具体过程为先将原始应变信号进行经验 模态分解，找出分解后得到的IMF分量中含噪较大的IMF分量进行相对应的阈值滤波处理， 再用滤波处理后得到的IMF分量与含噪较小的IMF分量以及残余分量进行EMD信号重构，从 而得到降噪后的信号。
[0033] 结合图3,为将原始桥梁应变信号进行经验模态分解之后得到的MF分量以及残余 分量，原可以看出原始信号被分解成从高频到低频段的12个IMF分量以及1个残余分量，如 以下公式所示：
[0034]
[0035] 其中X (t)代表桥梁应变原始信号，Ci (t)代表所有的IMF分量，rn (t)代表残余分量。 通过进一步分析可以发现，噪声基本分布在前6级IMF分量中，即imfl~imf6,所以只需要对 前6级IMF分量进行相对应的阈值滤波就可以将原始信号中的噪声全部去除。
[0036] 结合图4,为将原始信号经验模态分解后的前6级MF分量进行如下式所示的阈值 滤波前后的对比：
[0037]
[0038] 其中，cn⑴为原始各级IMF分量，average为各级IMF分量的平均值，尤⑴为阈值滤 波后的各级頂F分量，ThrValue为各级頂F分量对应的阈值，通过以下公式算出：
[0039]
[0040] 其中var为頂F分量的方差，L为頂F分量的长度，ThrValue为每级頂F分量所对应的 阈值。
[0041]通过对比图4中（a)和(b)可以看出，经过特定设计的阈值滤波后，原始信号前6级 頂F分量中的噪声基本被去除，也就意味着原始桥梁应变信号中的噪声基本被去除。
[0042]结合图5,为用滤波后的前6级IMF分量以与后6级IMF分量以及残余分量进行如下 式所示的信号重构后得到的降噪后的桥梁应变信号与原始信号的对比图：
[0043]
[0044] 其中Cl(t)代表未经过阈值滤波的頂F分量，即原始信号经验模态分解后含噪较小 的IMF分量，Kt)代表经过阈值滤波后的IMF分量，rn(t)为残余分量，X(t)代表降噪后的桥 梁应变信号。
[0045] 通过对比图5中（a)和(b)可以得出，经过该基于经验模态分解法的桥梁应变信号 降噪方法后，原始桥梁应变信号中的噪声基本被去除，且很好的保留了信号的原始特征。 通过该方法可以对桥梁应变信号以及一些大型建筑物的应变信号进行降噪处理，从而使得 通过这些应变信号得出的桥梁结构健康监测数据更加准确可靠，有助于相关部门对桥梁进 行根号的管理和维护。
【主权项】
1. 基于经验模态分解法的桥梁应变信号降噪方法，其特征是首先通过桥梁关键部位上 设置的应变信号采集单元进行数据采集，提取出桥梁关键部位的应变信号；其次是对提取 出的应变信号做经验模态分解EMD处理，具体分析经验模态分解得到的本征模函数IMF分 量，对IMF分量中含噪较大的分量进行相对应的阔值滤波;最后将通过阔值滤波后的IMF分 量与原始IMF分量中含噪较小的IMF分量W及残余分量一起进行EMD信号重构，从而得到降 噪后的桥梁应变信号。2. 根据权利要求1所述的降噪方法，其特征是所述的经验模态分解EMD处理，是指对桥 梁应变信号进行从高频段到低频段的分解，其步骤如下:首先，找出通过桥梁上设置的应变 信号采集单元采取的关键部位的原始应变信号X(t)的所有极大值点和极小值点，原始应变 信号通过应变信号采集单元取得、取得的极大值和极小值点指需要降噪信号的整个时间段 内的；并用Ξ次样条函数对所有极大值点和极小值点进行拟合，分别对应得到上包络线m (t)和下包络线U2(t)，并求出m(t)与U2(t)的均值ai(t)，即ai(t) = (山(t)+U2(t))/2;接着 用原始信号X(t)减去ai(t)后得到一个除去低频的新信号mi(t)，即mi(t)=X(t)-ai(t);通 常情况下mi(t)并不满足本征模函数IMF分量的要求；因此对mi(t)重复进行上述操作，即找 出mi(t)的所有极大值与极小值点，然后求出两者平均a2(t)，再用mi(t)减去a2(t)得到m2 (t)，如果m2(t)仍然不满足IMF分量的要求，则继续重复上述操作，直到得到的某一级mn(t) 符合IMF分量的要求，运样就得到了第一级IMF分量ci(t)，ci(t)表示原始应变信号中频率最 高的部分;接下来用X(t)减去ci(t)得到ri(t)，并对ri(t)重复上述操作得到第二级IMF分量 C2(t);不断重复上述步骤，得到其余所有的IMF分量W及残余分量;桥梁应变原始信号X(t) 经过EMD分解后可W用如下的公式表示：其中X (t)代表桥梁应变原始信号，ci (t)代表所有的IMF分量，rn( t)代表残余分量。3. 根据权利要求1所述的桥梁应变信号降噪方法，其特征是所述的IMF分量阔值，通过 如下公式来计算：其中var为IMF分量的方差，L为IMF分量的长度，ThrValue为每级IMF分量所对应的阔 值。4. 根据权利要求1所述的桥梁应变信号降噪方法，其特征是所述的对含噪IMF分量进行 针对桥梁应变信号设计的阔值滤波，针对不同的IMF分量的性质设计不同的阔值滤波器，从 而将噪声信号滤除，其中常见的一种滤波通过如下公式实现： 其中，Cn(t)为原始各级IMF分量，Thr化lue为各级IMF分量对应的阔值，average为各级 IMF分量的平均值，常(巧为阔值滤波后的各级IMF分量。5. 根据权利要求1所述的桥梁应变信号降噪方法，其特征是将阔值滤波得到的IMF分 量，未经阔值滤波处理的IMF分量W及残余分量进行EMD重构，通过如下公式实现：其中ci(t)代表未经过阔值滤波的IMF分量，即原始信号经验模态分解后含噪较小的IMF 分量，τ',·(〇代表经过阔值滤波后的IMF分量，rn(t)为残余分量，X(t)代表降噪后的桥梁应变 信号。
【文档编号】G06Q50/08GK105938542SQ201610150664
【公开日】2016年9月14日
【申请日】2016年3月16日
【发明人】张晖, 丁华平, 鲁明明, 耿朋, 李锦 , 谢兆金, 陈红梅
【申请人】南京大学, 江苏金晓电子信息股份有限公司