基于特征提取的光纤传感振动信号的异常判断方法和系统

1.本发明属于信号检测技术领域，具体涉及基于特征提取的光纤传感振动信号的异常判断方法和系统。


背景技术：

2.随着光纤传感技术的不断发展，光纤传感系统应用越来越普及，而分布式光纤传感系统以其结构简单、灵敏度高、监测范围广、响应速度快等优点在各个振动传感领域得到了广泛的应用。它在长距离油气管道输送、机场警卫、大型结构安全监测、周边安全等多个领域得到成功运用并迅速发展。不同应用领域中不同场景下对光纤扰动事件的模式识别性能与实时性问题，是对于光纤传感技术发展的核心问题。
3.现有的光纤周界系统的振动信号特征提取方法主要有四种，分别是基于小波分解的方法、基于其他分解模型的方法、提取信号波形图统计特征的方法以及基于梅尔频率倒谱系数的特征提取方法。基于小波分解的方法由于信号的峰值、功率谱和频带能量等特征容易受到高斯噪声的影响，在小波方法处理中存在模糊提取和线性稳态的缺陷，并且当数据量比较大时这种处理方法的计算量较大。基于其他分解模型的方法由于其分解模型的固有缺陷，对信号依靠特征分类产生影响，降低其分类准确率。提取信号波形图统计特征的方法对于类别不同波形相似的信号会造成统计特征的相似，无法作为有效的分类判别依据。基于梅尔频率倒谱系数的特征提取方法常用于声音信号的特征提取，也可用于振动信号特征提取，但算法普适性造成了性能下降，影响分类性能。
4.上述方法在一定程度上解决了光纤传感振动信号模式识别的有效性与实时性问题，但是仍存在着一些限制。梅尔频率倒谱系数特征提取方法是在倒谱技术上依据人耳的感知特性进行改进而提出的一种特征提取算法，它的稳定性与高识别率使得其在语音识别领域得到了广泛应用。然而语音信号与光纤传感振动信号的频率范围存在差异，直接将梅尔频率倒谱系数特征提取方法应用到光纤传感振动信号中并不能完全发挥梅尔频率倒谱系数的优势，造成光纤传感系统的识别率低等问题。因此需要对梅尔频率倒谱系数特征提取方法进行改进，提高光纤传感系统的检测实时性及识别率，实现对异常信号的快速、准确识别。


技术实现要素：

5.为解决上述问题，本发明提出了基于特征提取的光纤传感振动信号的异常判断方法和系统，采用标准差局部线性嵌入对特征向量进行降维，解决了在梅尔频率倒谱系数特征提取过程中因语音信号与光纤传感振动信号频率差异而导致的识别率低的问题，提高了光纤传感系统的检测实时性及识别率，实现对异常信号的快速、准确识别。
6.为达到上述目的，本发明采用如下技术方案：
7.第一方面，本发明提供基于特征提取的光纤传感振动信号的异常判断方法，具体步骤如下：
8.采集光缆中的振动信号；
9.对振动信号进行预处理；
10.将预处理后的振动信号通过基于频谱能量分析的梅尔频率倒谱系数特征提取方法提取特征向量矩阵；
11.将特征向量矩阵通过标准差局部线性嵌入方法进行特征降维；
12.根据降维后的特征向量判断信号是否为异常信号；
13.若信号为异常信号，则进行报警。
14.优选地，将预处理后的振动信号通过基于频谱能量分析的梅尔频率倒谱系数特征提取方法提取特征向量矩阵的步骤包括：
15.获取振动信号的功率谱；
16.基于频谱能量分析的梅尔滤波器组；
17.根据梅尔滤波器组与功率谱得到对数能量；
18.对对数能量进行离散余弦变换，获取基于频谱能量分析的梅尔频率倒谱系数；
19.对基于频谱能量分析的梅尔频率倒谱系数进行动态显著性差分处理，得到差分系数；
20.组合基于频谱能量分析的梅尔频率倒谱系数与动态差分系数，得到光纤传感振动信号的联合特征向量。
21.优选地，获取基于频谱能量分析的梅尔滤波器组包括如下步骤：
22.根据不同频段能量占比，通过阈值判断得到频率范围，确定截止频率；
23.根据截止频率，调整梅尔转换公式；
24.根据调整的梅尔转换公式，得到基于频谱能量分析的梅尔滤波器组。
25.优选地，根据不同频段能量占比，通过阈值判断得到频率范围，确定截止频率，包括如下步骤：
26.对光纤振动传感信号进行小波包分解，获得不同频段能量占比；通过对信号进行i层小波包分解产生j个不同的尺度，则第i层分解的输出系数记为c
i,j
，其中j＝0,1,2,...,2
i-1，i为分解层次，j为相应的尺度。
27.每个频段对应的能量为ej,每个频段的能量占比为pj，根据不同频段能量占比，通过阈值l判断，得到频率范围，确定截止频率；选择能量占比之和大于等于l的频段，确定截止频率为f
cut-off
,具体实现过程如下:
28.每个频段对应的能量为ej:
[0029][0030]
总能量为e
total
:
[0031]etotal
＝∑ej[0032]
每个频段的能量占比为pj:
[0033]
pj＝ej/e
total
[0034]
频段范围应该满足判断公式，即
[0035]
[0036]
其中，l为0.95；
[0037]
根据频段集合确定频段范围f
cut-off
：
[0038][0039]
其中fs，为每个频带的频率范围，j为频段集合，j＝{1,
…
,j},
[0040]
优选地，所述调整前的调整梅尔-频率转换公式为：
[0041][0042]
其中f是频率，m是信号频率f对应的梅尔尺度，c为梅尔-频率转换公式系数。
[0043]
优选地，依据截止频率设计用于调整梅尔-频率转换公式的调整参数为：
[0044][0045]
其中，k1,k2为调整系数。
[0046]
优选地，依据调整参数调整梅尔-频率转换公式，调整后的梅尔-频率转换公式为：
[0047][0048]
其中n是补偿拟合常数，为n＝(k
1-k2)cln2，调整后梅尔-频率转换公式在截止频率内保持线性函数关系并增幅线性关系，在超过截止频率后保持对数函数关系并平缓对数关系。
[0049]
优选地，将特征向量矩阵通过标准差局部线性嵌入进行特征降维，包括如下步骤：
[0050]
对于特征向量矩阵，对于各个维度的每个样本点求相邻点的标准化欧式距离，根据标准化欧式距离选择最近邻域点；
[0051]
对每个样本点与它的邻域集进行重构权值计算，获得重构权值矩阵；
[0052]
对损失函数φ(y)进行极小化运算，获得低维嵌入矩阵y：
[0053][0054]
其中，y为所得低维嵌入矩阵，w为权值矩阵，wi,ii分别是单位矩阵w,i的第i列，t为转置。
[0055]
优选地，对振动信号进行预处理的过程包括去噪及端点检测。
[0056]
第二方面，本发明提供了基于特征提取的光纤传感振动信号的异常判断系统，包括：存储介质和处理器；
[0057]
所述存储介质用于存储指令；
[0058]
所述处理器用于根据所述指令进行操作以执行根据第一方面中中任一项所述的方法
[0059]
本发明所达到的有益效果：本发明提供基于特征提取的光纤传感振动信号的异常
判断方法和系统，通过对梅尔频率倒谱系数特征提取方法的改进，解决了在梅尔频率倒谱系数特征提取过程中因语音信号与光纤传感振动信号频率差异而导致的识别率低的问题，提高了光纤传感系统的检测实时性及识别率，实现对异常信号的快速、准确识别。本发明通过基于频谱能量分析的梅尔频率倒谱系数特征提取方法提取特征向量，使用标准差局部线性嵌入方法进行特征降维，进而确定光纤振动信号是否存在异常，可用于各种基于光纤振动信号的安防系统，具有较强的准确率和适应性。
附图说明
[0060]
图1是本发明实施例一提供的基于特征提取的光纤传感振动信号的异常判断方法的流程示意图；
[0061]
图2是本发明实施例一中基于频谱能量分析的梅尔频率倒谱系数特征提取方法的流程示意图。
[0062]
图3是标准差局部线性嵌入方法的流程示意图。
[0063]
图4是本发明实施例一中的频段能量占比图。
[0064]
图5是本发明实施例一中的改进梅尔-频率转换关系图
具体实施方式
[0065]
下面根据附图并结合实施例对本发明的技术方案作进一步阐述。
[0066]
实施例一：
[0067]
图1是本发明实施例的基于特征提取的光纤传感振动信号的异常判断方法的流程示意图，包括以下步骤：
[0068]
s1：采集光缆中的振动信号，对振动信号进行预处理：
[0069]
振动信号包括因非法入侵产生的异常信号、恶意破坏产生的异常信号、风吹产生的异常信号和正常走动等产生的干扰信号；
[0070]
对振动信号进行预处理的过程包括去噪以及端点检测：去噪的方法包括小波去噪；使用峰值估计以及小波包络截取法进行端点检修以及截取。
[0071]
s2：将预处理后的振动信号通过基于频谱能量分析的梅尔频率倒谱系数特征提取方法提取特征向量矩阵，如图2所示，包括如下步骤：
[0072]
21)获取振动信号的功率谱：
[0073]
212)对预处理后的振动信号进行预加重，获得加重信号；
[0074]
预加重的过程包括：对预处理后的振动信号的高频部分进行预加重处理，根据如下公式进行预加重处理：y[n]＝x[n]-αx[n-1]。其中，y[n]是经过预加重处理后的输出信号，x[n]是光纤振动信号输入，x[n-1]是光纤振动信号下一个输入，α是滤波器系数，α的值是0.97；
[0075]
212)对加重信号进行分帧，获得分帧信号；
[0076]
本实施例中信号采样频率预设10khz，计算对应的采样长度，选定帧长为25.6ms，采样频率乘以帧长为采样点，即每一帧包括256个采样点；为了避免相邻两帧之间的变化过大，因此让两相邻帧之间有一段重叠区域即为帧移，此重叠区域的采样点数一般为每帧采样点数的1/3。
[0077]
213)对每一帧分帧信号加汉明窗，得到加窗信号；
[0078]
214)对每一帧加窗信号进行快速傅里叶变换，得到每一帧加窗信号的频谱；
[0079]
215)对频谱进行计算得到功率谱。
[0080]
22)获取基于频谱能量分析的梅尔滤波器组：
[0081]
221)对光纤振动传感信号进行小波包分解，获得不同频段能量占比；通过对信号进行i层小波包分解产生j个不同的尺度，则第i层分解的输出系数记为c
i,j
，其中j＝0,1,2,...,2
i-1，i为分解层次，j为相应的尺度。
[0082]
222)每个频段对应的能量为ej,每个频段的能量占比为pj，根据不同频段能量占比，通过阈值l判断，得到频率范围，确定截止频率；选择能量占比之和大于等于l的频段，确定截止频率为f
cut-off
，具体实现过程如下:
[0083]
每个频段对应的能量为ej：
[0084][0085]
总能量为e
total
：
[0086]etotal
＝∑ej[0087]
每个频段的能量占比为pj：
[0088]
pj＝ej/e
total
[0089]
频段范围应该满足判断公式，即
[0090][0091]
其中，l为0.95；
[0092]fs
为每个频带的频率范围，j为频段集合，j＝{1,
…
,j},根据频段集合确定截止频率f
cut-off
：
[0093][0094]
图4为本实施例的频段能量占比图，根据频段能量占比的不同，第一类经过计算判定为第一频段与第二频段，频率范围为0-1250hz，截止频率为1250hz。
[0095]
223)根据截止频率，调整梅尔-频率转换公式系数c：
[0096]
调整前的梅尔-频率转换公式为：
[0097][0098]
其中f是频率；m是信号频率f对应的梅尔尺度；c为梅尔-频率转换公式系数；f
cut-off
为梅尔-频率转换公式中截止频率。依据截止频率设计用于调整梅尔-频率转换公式的调整参数：
[0099][0100]
其中，k1,k2为调整系数。
[0101]
依据调整参数调整梅尔-频率转换公式，调整后梅尔-频率转换公式为：
[0102][0103]
其中n是补偿拟合常数，为n＝(k
1-k2)cln2，调整后梅尔-频率转换公式在截止频率内保持线性函数关系并增幅线性关系，在超过截止频率后保持对数函数关系并平缓对数关系。在本实施例中，f
cut-off
为1250hz，k1为2,k2为1/2。图5为本发明实施例提供的改进梅尔-频率转换关系图。
[0104]
224)根据调整后的梅尔-频率转换公式，得到基于频谱能量分析的梅尔滤波器组；
[0105]
23)根据梅尔滤波器组与功率谱得到对数能量s(m)：
[0106][0107]
其中p(k)是信号的功率谱，hm(k)是滤波器组，k为频域中第k个谱线，在此实施例中有n-1个谱线。
[0108]
24)对对数能量进行离散余弦变换，获取基于频谱能量分析的梅尔频率倒谱系数：
[0109][0110]
其中，s(m)是对数能量，得到c(n)为基于频谱能量分析的梅尔频率倒谱系数。
[0111]
25)对基于频谱能量分析的梅尔频率倒谱系数进行动态显著性差分处理，得到差分系数；
[0112]
26)组合基于频谱能量分析的梅尔频率倒谱系数与差分系数，得到联合特征向量。
[0113]
s3：将联合特征向量通过标准差局部线性嵌入进行特征降维，如图3所示：
[0114]
31)对每个样本点求相邻点的标准化欧式距离，根据标准化欧式距离选择最近邻域点。
[0115]
32)对每个样本点与它的邻域集进行重构权值计算，获得重构权值矩阵。
[0116]
33)对损失函数φ(y)进行极小化运算，获得低维嵌入矩阵y：
[0117][0118]
其中，y为所得低维嵌入矩阵，w为权值矩阵，wi,ii分别是单位矩阵w,i的第i列，t为转置。
[0119]
s4：通过降维后的特征向量矩阵判断信号是否为异常信号：若判断异常信号为扰动信号，则进入步骤s5；若判断异常信号不是扰动信号，则执行步骤s1；
[0120]
降维后基于频谱能量分析的梅尔频率倒谱系数特征向量矩阵判断信号是否为异常信号的方法为：信号经过特征提取得到特征向量矩阵，然后进行降维处理后，送入预先训练好的支持向量机判别是否为异常情况，为何种异常行为。
[0121]
s5：对异常信号根据不同情况进行报警。
[0122]
实施例二：
[0123]
基于与实施例一相同的发明构思，本发明提供了基于特征提取的光纤传感振动信号的异常判断系统，包括：存储介质和处理器；
[0124]
存储介质用于存储指令；
[0125]
处理器用于根据所述指令进行操作以执行根据实施例一中任一项所述的方法。
[0126]
以上所述仅为本发明的优选实施方式，用以说明本发明的技术方案，而非对其限制；应当指出：对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。