一种基于量化信息熵的一维构件应力波信号奇异点检测的方法
【技术领域】
[0001] 本发明涉及信号特征识别领域，具体涉及的是一种基于量化信息熵的一维构件应 力波信号奇异点检测的方法。
【背景技术】
[0002] 信号的奇异点检测是应用非常广泛的技术，特别是在结构体的完整性检测上，例 如一维构件的完整性检测。
[0003] 目前在一维构件完整性检测上应用最广泛的方法为应力波法，即：在一维构件的 同一端安置应力波激发器和信号接收传感器，激发应力波之后，应力波传感器会接收到完 整的应力波信号波形，波形中包含丰富的信息，通过分析信号波形特征量，可以得到一维构 件的完整性信息。目前常见的应力波特征量包括信号的能量（功率谱密度），方差等，提取 特征量的方式有在信号整个时程上提取特征量，还有先将信号分段，提取每段信号的特征 量。无论是在信号整个时程上，还是分段提取信号的特征量，都难免会造成奇异点的遗漏和 检测分辨率的降低。
[0004] 本发明提出量化信息熵的概念，并将量化信息熵作为信号的特征量；构造时间窗 为工具，在信号上移动提取特征量；将信号特征量降维，分析应力波信号的奇异点位置。与 已知的基于模极大值法等奇异点检测方法相比，本发明实现了特征量的降维，在信号存在 多处奇异点的情况下可以得到优异的检测结果，提高了检测分辨率，并且抗噪性能优良，不 占有过多的时间资源。
【发明内容】
[0005] 本发明的目的在于提供一种基于量化信息熵的一维构件应力波信号奇异点检测 的方法。
[0006] 本发明的目的是这样实现的：
[0007] (1)将原始应力波信号进行预处理，进行n层小波包分解，得到不同频率范围的2n 个子波信号；
[0008] (2)在信息熵的基础上提出量化信息熵，作为一维构件应力波信号的特征量：
[0009] 信号的量化信息熵为：
[0010]
【主权项】
1. 一种基于量化信息熵的一维构件应力波信号奇异点检测的方法，其特征在于，包括 如下步骤： (1) 将原始应力波信号进行预处理，进行n层小波包分解，得到不同频率范围的2n个子 波信号； (2) 在信息熵的基础上提出量化信息熵，作为一维构件应力波信号的特征量： 信号的量化信息熵为：
其中DH(X)即为信号的量化信息熵，D为幅值区间的个数，Nunii表示第i个幅值区间内 采样点个数，Pi为第i个幅值区间内采样点数的统计概率 信号采样幅值集合为X= {x^x2, . . .，xn)，X的概率分布为Pi=P(x)(i= 1，2,...， n)
变量X的最大值与最小值，分别记为x_^Px_，幅值由Xujljx_平均分成 D个量化幅值区间，则每段幅值的取值长度为：
每段幅值的量化取值范围为： (xmin+iAd，xmin+(i+l)Ad，）（i= 0,1，2, ? ? ?，D-1); 第i段幅值区间内的采样点个数为Nmv第i段幅值区的统计概率为：
应力波经小波包k层分解后，得到m=2k个信号子波，时间窗步进数为X，提取特征量， 即量化信息熵之后，每个子波熵值向量为： Hm={DHml，DHm2，DHm3，? ? ?，DHmn} {m= 1，2，...，2); 应力波量化信息熵矩阵为：
其中每行数据代表一个子波的量化信息摘向量，m个子波构成mXn特征量矩阵； (3) 构造固定宽度W的时间窗，在信号的时程上按步长B移动提取信号特征量，构造信 号的特征量矩阵； 时间窗口宽度为W，步进长度为B，时间窗口宽度W的取值覆盖至少一个完整的奇异点， 不超过两个奇异点；步进长度B的取值满足W/8<B<W/4 ; (4) 对信号的特征量矩阵做降维处理，分析判定信号奇异点的位置信息： 所有子波对应同一位置的特征值的均值，得到一维的信号特征量，求得对应同一段信 号的时间窗的量化信息熵均值：
处理得到量化信息熵均值向量为：
若信号的特征值左右均收敛至最小值，则对应的信号部分存在奇异点，ff的第i个元 素在左右对称的特定范围内均收敛为最小值，则： ((i-l)B, (i-l)B+ff)(iG[l,n])； 该段区域为应力波信号奇异点位置，其中B为步进长度，W为时间窗宽度。
【专利摘要】本发明涉及信号特征识别领域，具体涉及的是一种基于量化信息熵的一维构件应力波信号奇异点检测的方法。本发明包括：(1)将原始应力波信号进行预处理，进行n层小波包分解，得到不同频率范围的2n个子波信号；(2)在信息熵的基础上提出量化信息熵，作为一维构件应力波信号的特征量；(3)构造固定宽度W的时间窗，在信号的时程上按步长B移动提取信号特征量，构造信号的特征量矩阵；(4)对信号的特征量矩阵做降维处理，分析判定信号奇异点的位置信息。与现有算法技术相比，本发明的有益效果是实现了特征量的降维，在信号存在多处奇异点的情况下可以得到优异的检测结果，提高了检测分辨率，并且抗噪性能优良，不占有过多的时间资源。
【IPC分类】G01N29-44
【公开号】CN104698091
【申请号】CN201510112344
【发明人】康维新, 李敬德, 高爽
【申请人】哈尔滨工程大学
【公开日】2015年6月10日
【申请日】2015年3月13日