1.一种爆破盲炮检测方法,其特征在于:它包括以下步骤:⑴对原始信号进行EMD分解;
⑵计算原始信号与各IMF分量的互相关函数,具体计算公式如下:计算原始信号y(t)与各分量xi(t)的互相关函数,作归一化处理得到互相关函数rxy(τ)为:式中:t为时间,τ为延迟时间,x(t)为各分量信号,T为信号周期时间,rxy(τ)表示信号x(t)和y(t)的互相关函数值。
⑶定义各分量反映原始信号的敏感度,具体是取互相关函数rxy(τ)的峰值点作为各分量反映爆破信号盲炮特征的敏感度ai:ai=max(rxy(τ))
式中:max(rxy(τ))为互相关函数rxy(τ)的最大值,ai为第i个分量反映爆破信号盲炮特征的敏感度。
⑷计算相邻敏感度的衰减率,具体是将敏感度ai从大到小进行排序,计算相邻敏感度的衰减率kj→i为:kj→i=(aj‑ai)/ai
式中:kj→i为敏感度从大到小排序后第j个分量的敏感度aj相对于第i个分量敏感度ai的衰减率。
⑸选择主IMF分量组合,具体是设衰减率的第一个极大值为km,则取排序后的前m个敏感度所对应的IMF分量作为主分量c(t);
⑹提取主IMF分量包络线;
⑺提取包络线峰值点进行盲炮检测。
2.根据权利要求1所述的爆破盲炮检测方法,其特征在于:所述的对原始信号进行分解,具体是将爆破振动信号y(t)通过EMD分解成一系列IMF分量xi(t)。
3.根据权利要求1所述的爆破盲炮检测方法,其特征在于:所述的提取主IMF分量包络线,具体计算公式如下:对所选取的主分量组合进行希尔伯特变换,提取包络线a(t)为:jφ(t)
z(t)=c(t)+jH[c(t)]=a(t)e式中:pv代表柯西主值,c(t)为选取的主分量,H[c(t)]为主分量c(t)的希尔伯特变换值,j为虚数单位,φ(t)为主分量c(t)的瞬时相位,a(t)为解析信号z(t)的幅值,也称为信号的包络。
4.根据权利要求1所述的爆破盲炮检测方法,其特征在于:所述的盲炮检测是对各主分量包络线峰值点进行汇总,通过对比各炮孔的设计延时与包络线峰值点的实际延时,可对盲炮进行检测。