[发明专利]一种基于量化信息熵的一维构件应力波信号奇异点检测的方法

说明书

技术领域

本发明涉及信号特征识别领域，具体涉及的是一种基于量化信息熵的一维构件应力波信号奇异点检测的方法。

背景技术

信号的奇异点检测是应用非常广泛的技术，特别是在结构体的完整性检测上，例如一维构件的完整性检测。

目前在一维构件完整性检测上应用最广泛的方法为应力波法，即：在一维构件的同一端安置应力波激发器和信号接收传感器，激发应力波之后，应力波传感器会接收到完整的应力波信号波形，波形中包含丰富的信息，通过分析信号波形特征量，可以得到一维构件的完整性信息。目前常见的应力波特征量包括信号的能量(功率谱密度)，方差等，提取特征量的方式有在信号整个时程上提取特征量，还有先将信号分段，提取每段信号的特征量。无论是在信号整个时程上，还是分段提取信号的特征量，都难免会造成奇异点的遗漏和检测分辨率的降低，如2014102188871所公开的以中国焊接缺陷提取方法及焊接缺陷检测方法。

本发明提出量化信息熵的概念，并将量化信息熵作为信号的特征量；构造时间窗为工具，在信号上移动提取特征量；将信号特征量降维，分析应力波信号的奇异点位置。与已知的基于模极大值法等奇异点检测方法相比，本发明实现了特征量的降维，在信号存在多处奇异点的情况下可以得到优异的检测结果，提高了检测分辨率，并且抗噪性能优良，不占有过多的时间资源。

发明内容

本发明的目的在于提供一种基于量化信息熵的一维构件应力波信号奇异点检测的方法。

本发明的目的是这样实现的：

(1)将原始应力波信号进行预处理，进行n层小波包分解，得到不同频率范围的2ⁿ个子波信号；

(2)在信息熵的基础上提出量化信息熵，作为一维构件应力波信号的特征量：

信号的量化信息熵为：

其中DH(X)即为信号的量化信息熵，D为幅值区间的个数，Num_i表示第i个幅值区间内采样点个数，p_i为第i个幅值区间内采样点数的统计概率

信号采样幅值集合为X＝{x₁，x₂,…,x_n}，X的概率分布为p_i＝P(x_i)(i＝1,2,…,n)，变量X的最大值与最小值，分别记为x_max和x_min，幅值由x_min到x_max平均分成D个量化幅值区间，则每段幅值的取值长度为：

每段幅值的量化取值范围为：

(x_min+iΔd,x_min+(i+1)Δd,)(i＝0,1,2,…,D-1)；

第i段幅值区间内的采样点个数为Num_i，第i段幅值区的统计概率为：

应力波经小波包k层分解后，得到m＝2^k个信号子波，时间窗步进数为x，提取特征量，即量化信息熵之后，每个子波熵值向量为：

H_m＝{DH_m1,DH_m2,DH_m3,…,DH_mn}(m＝1,2,…,2^k)；

应力波量化信息熵矩阵为：

其中每行数据代表一个子波的量化信息熵向量，m个子波构成m×n特征量矩阵

(3)构造固定宽度W的时间窗，在信号的时程上按步长B移动提取信号特征量，构造信号的特征量矩阵；

时间窗口宽度为W，步进长度为B，时间窗口宽度W的取值覆盖至少一个完整的奇异点，不超过两个奇异点；步进长度B的取值满足W/8<B<W/4；

(4)对信号的特征量矩阵做降维处理，分析判定信号奇异点的位置信息：

所有子波对应同一位置的特征值的均值，得到一维的信号特征量，求得对应同一段信号的时间窗的量化信息熵均值：

处理得到量化信息熵均值向量为：

若信号的特征值左右均收敛至最小值，则对应的信号部分存在奇异点，的第i个元素在左右对称的特定范围内均收敛为最小值，则：

((i-1)B,(i-1)B+W)(i∈[1,n])；

该段区域为应力波信号奇异点位置，其中B为步进长度，W为时间窗宽度。

本发明的有益效果在于：