[发明专利]一种基于波形动态匹配的子波拉伸校正处理方法

权利要求书

1.一种基于波形动态匹配的子波拉伸校正处理方法，其特征在于，包括以下步骤：

对参考序列与匹配序列进行波形动态匹配，使得两个序列相匹配后之间的误差和达到最小值；

基于误差和最小值对应的时差，回溯找到最短路径对应的时移量序列；

将得到的时移量序列应用到需要匹配的序列上，即得到校正后的处理结果，

其中，对参考序列与匹配序列进行波形动态匹配的步骤包括，由原始数据叠加求取参考序列f[i]，定义参考序列f[i]和匹配序列g[i]间的校正误差函数e[i,l]：

e[i,l]＝(f[i]-g[i+l])²,i＝0,1,...,N-1

则e[i,l]表示对匹配序列g的第i个采样点进行时移l后与参考序列f第i个采样点之间的距离；

计算参考序列f[i]和匹配序列g[i]匹配后的误差和，即总距离函数：

2.根据权利要求1所述的基于波形动态匹配的子波拉伸校正处理方法，其特征在于，计算总距离函数包括定义时间应变约束条件，并依此推导出最优化求解中误差积累递归表达式，假设时间应变约束条件为|u[i_u]-u[i_u-1]|≤1，则累积过程为：

d[0,l]＝e[0,l],

i_u＝1,2,...,N-1,

式中，u[i_u]为时移量序列；

d[i_u,l]为距离函数。

3.根据权利要求2所述的基于波形动态匹配的子波拉伸校正处理方法，其特征在于，所述累积过程中当i_u＝N-1时，计算出所有时差l对应的距离函数，找到最小的距离，即得到最小的总距离函数：

此时最小距离对应的时差l，为u[N-1]的值。

4.根据权利要求3所述的基于波形动态匹配的子波拉伸校正处理方法，其特征在于，回溯找到最短路径对应的时移量序列的步骤包括，回溯求解最短路径对应的时移量序列u[0:N-1]，回溯表达式为：

i_u＝N-1,N-2,...,1,

通过往前递归地回溯，最后找到每一个采样点对应的时差，即所述的时移量序列。

5.根据权利要求4所述的基于波形动态匹配的子波拉伸校正处理方法，其特征在于，将时移量序列应用于匹配序列g[i]得到校正后的结果