一种用于高保真阵列处理的宽带时延估计方法与流程

文档序号:17758245发布日期:2019-05-24 21:28阅读:197来源:国知局
一种用于高保真阵列处理的宽带时延估计方法与流程

本发明属于声纳信号处理领域,尤其涉及一种用于高保真阵列处理的宽带时延估计方法。



背景技术:

水声目标分类识别是水声信号处理的重要课题之一。对于被动声纳目标分类识别,关键是目标辐射噪声的特征提取。为了提高特征提取的性能,波束形成处理是从传感器阵列信号中重构源信号,在空间上抗噪声和分离多目标干扰的有效手段和常用方法。但在水声环境中,由于阵元位置的误差、平面波的假设以及目标方位角的估计误差等,都可能造成常规波束形成方法跟踪波束输出时出现时延失配,从而造成信号在功率谱和解调谱等特征方面的畸变,从而影响目标分类识别的效果。所以对于目标识别这一特殊的声纳使命任务,高保真重构源信号及特征的宽带阵处理方法是必须重视的研究课题。

目前国内外学者提出了许多高保真宽带阵处理方法,主要有:(1)阵形校正:通过估计误差参数来实现对阵形的校正,该类方法首先对阵列误差进行建模,将阵列误差校正转化为参数估计的问题。该类阵形校正方法通常可以分为有源校正类和自校正类。对于有源校正方法,该方法对辅助信号源有着较高的精确方位信息的要求,所以当辅助信号源的方位信息有偏差时,这类校正算法会带来较大的偏差。自校正算法由于阵元位置误差与方位参数之间的耦合和某些病态的阵列结构,参数估计的唯一辨识往往无法保证,更为重要的是参数联合估计对应的高维、多模非线性优化问题带来了庞大的运算量,估计的全局收敛性往往无法保证。(2)聚焦波束形成:根据测量区内不同方位和距离对时延差进行球面波补偿,获得测量区内目标声源的分布图,以实现近场目标的高精度被动定位,根据估计的目标距离和方位采用常规球面波波束形成高保真重构源信号,该技术原理简单、易于实现,能够有效解决近场源条件下平面波假设带来的阵元间时延失配问题,但运算量大且不能解决阵元位置误差带来的时延失配(3)esprit(借助旋转不变性估计信号参数)超分辨测向:利用阵列的平移不变性来进行参数估计,可直接通过特征值获得信号参数,该方法无需进行计算量极大的谱峰搜索,无需存储阵列流行矩阵,但要求阵形精确已知(4)恒定束宽波束形成:选择在宽带信号带宽内各频率分量上波束形成的加权向量,使得基阵系统对不同频率的输入信号具有相同的空间滤波响应,该技术可消除目标方位角估计误差造成的跟踪波束输出信号的非线性失真,但运算量大且不能解决由于阵元位置的误差、平面波的假设以及目标方位角的估计误差等造成的时延失配(5)零导向波束形成:调整权值,使阵响应在信号方向为1,在干扰方向为0,从而达到抑制干扰的目的,它实际上是使波束图在干扰方向形成零点,该技术对消除定向干扰有重要作用,但需要对信号和干扰方位以及阵形信息精确已知。



技术实现要素:

发明目的:针对现有技术中存在的问题,本发明公开了一种用于高保真阵列处理的宽带时延估计方法,该方法基于目标信号不同频率的功率谱线谱具有相同的阵元间时延这一特性,根据线谱的频率和频域信噪比对基于不同频率线谱估计的阵元间相位差进行加权融合,计算更为精确的目标信号宽带时延矢量,获取保真增强的重构源信号。

技术方案:为了实现上述发明目的,本发明采用以下技术方案:一种用于高保真阵列处理的宽带时延估计方法,包括如下步骤:

(1)基于目标信号远场平面波入射和阵形无畸变两个假设,计算基阵采集数据不同方位的预成波束数据,求取波束能量图,估计目标方位

(2)基于目标信号远场平面波入射和阵形无畸变两个假设,从基阵采集数据中按估计的目标方位对源信号进行预重构,检测预重构源信号s1[n],n=1,2,...,n的功率谱线谱位置fk,k=1,2,...,k,计算线谱频率的精测值其中k为检测出的频率小于v/(2d)的功率谱线谱根数,v为声波在水中的传播速度,d为均匀线阵阵元间距;

(3)计算参考阵元h1采集数据的功率谱在位置处线谱的频域信噪比snrk,k=1,2,...,k;

(4)计算基阵采集数据的频谱在位置处线谱的阵元间相位差其中,m为基阵阵元数目;

(5)根据线谱频率和频域信噪比对基于不同频率线谱估计的阵元间相位差进行加权融合,估计目标信号宽带时延矢量

(6)根据估计的目标信号宽带时延矢量从基阵采集数据中对源信号进行二次重构,获取保真增强的重构源信号s2[n],n=1,2,...,n。

优选的,步骤(1)中,采用如下方法计算基阵采集数据的波束能量图,估计目标方位,具体包括如下步骤:

(1-1)计算按等余弦划分的预成波束角度:

其中,θi为第i个预成波束的角度,以基阵第一个阵元为坐标原点,建立平面直角坐标系,基阵第m个阵元的坐标为[-(m-1)d,0],m=1,2,...,m,θi角度定义为原点和目标坐标的连线与x轴的夹角,θi∈[0,π],表示计算的反余弦,i为预成波束个数,d为均匀线阵阵元间距;

(1-2)基于目标信号远场平面波入射和阵形无畸变两个假设,计算均匀线阵相邻阵元采集信号在预成波束角度θi下的时延差τi:

其中,cos(θi)表示计算θi的余弦值,v为声波在水中的传播速度;

(1-3)按如下波束形成算法;以基阵阵首第一个阵元h1为参考阵元,计算预成波束角度θi对应的预成波束数据

其中,m为基阵阵元数目,xm[n]表示第m个阵元第n个采样数据,n为预成波束数据离散点数,fs为基阵采集数据时的采样频率;

(1-4)计算基阵采集数据在预成波束角度θi下的波束能量e[i]:

(1-5)估计目标方位

搜索e[i]最大值所对应的预成波束索引imax

其中,argmax1≤i≤i{e[i]}表示在1≤i≤i范围内搜索e[i]最大值所对应的预成波束索引imax。若1<imax<i,将分别记作z-1、z0、z1,则基于抛物线三点插值估计的目标方位为:

优选的,在步骤(2)中,采用如下方法检测预重构源信号的功率谱线谱位置fk,计算线谱频率的精测值具体包括如下步骤:

(2-1)按如下波束形成算法,以基阵阵首第一个阵元h1为参考阵元,从基阵采集数据中按估计的目标方位对源信号进行预重构,获取预重构源信号s1[n]:

(2-2)计算预重构源信号s1[n]的离散频谱:

则预重构源信号s1[n]的功率谱为:

其中,l为s1[l]和p1[l]的离散频率索引,j表示虚数单位,即||代表取模值运算;

(2-3)将功率谱p1[l]通过q阶中值滤波器得到功率谱p1[l]的趋势项c1[l],l=0,1,...,n/2,则功率谱p1[l]与趋势项c1[l]的差值谱d1[l]为:

d1[l]=p1[l]-c1[l],l=0,1,...,n/2

d1[l]的规范化差值谱为:

其中,std(d1[l])是求d1[l]的标准差;

(2-4)设定幅度门限g提取功率谱线谱,若d2(l)满足如下条件则判为功率谱线谱:

假设共提取出k根频率小于v/2d的功率谱线谱,第k根线谱在d2[l]中的索引为jk,则第k根线谱的频率为:

fk=jkδf,k=1,2,...,k

其中,δf为长度为n的离散傅里叶变换的频率分辨率,δf=fs/n;

(2-5)取s1[l]在离散频率索引(jk-1),jk,(jk+1)处的模值,分别记作ak1,ak2,ak3,即ak1=|s1[jk-1]|,ak2=|s1[jk]|,ak3=|s1[jk+1]|,利用所述模值ak1,ak2,ak3计算第k根功率谱线谱频率插值的相对偏差δk,即:

则第k根功率谱线谱的频率精测值为:

优选的,在步骤(3)中,采用如下方法计算参考阵元h1采集数据功率谱位置处线谱的频域信噪比snrk,具体包括如下步骤:

(3-1)计算参考阵元h1采集数据x1[n]的功率谱:

其中,l为的离散频率索引;

(3-2)估计参考阵元h1采集数据x1[n]的功率谱在位置处噪声的功率值

其中,wl为噪声功率值估计器的长度,取偶数,且wl≤2(jk-1);

(3-3)计算参考阵元h1采集数据x1[n]的功率谱在位置处信号加噪声的功率值

(3-4)参考阵元h1采集数据x1[n]的功率谱在位置处线谱的频域信噪比snrk为:

优选的,在步骤(4)中,采用如下方法计算基阵采集数据的频谱在位置处的阵元间相位差具体包括如下步骤:

(4-1)计算第m个阵元hm采集数据的频谱在位置处的值:

(4-2)计算各阵元采集数据的频谱在位置处的相位矢量

其中,为求取相位的运算;

(4-3)计算相邻阵元采集数据的频谱在位置处的相位之差

其中,为矢量的第m个元素,m=1,2,...,m;

优选的,在步骤(5)中,采用如下方法对基于不同频率线谱估计的目标信号阵元间相位差进行加权融合,估计目标信号宽带时延矢量具体包括如下步骤:

(5-1)根据线谱频率和频域信噪比snrk,计算第k根线谱的阵元间相位差加权系数wk:

(5-2)根据阵元间相位差加权系数wk对阵元间相位差进行加权融合,估计目标信号宽带时延矢量

其中

优选的,在步骤(6)中,采用如下方法对源信号进行二次重构:根据估计的目标信号宽带时延矢量从基阵采集数据xm[n]中对源信号进行二次重构,获取保真增强的重构源信号s2[n]:

有益效果:与现有的方法相比,本发明的技术方案具有以下有益技术效果:

(1)本发明公开的方法基于目标信号强线谱的阵元间相位差来估计目标信号宽带时延矢量,实现了阵形畸变、近场源条件下数据驱动的、自适应的保真增强源信号重构;

(2)本发明的方法基于目标信号不同频率的功率谱线谱具有相同的阵元间时延这一特性,根据功率谱线谱的频率和频域信噪比对基于不同频率线谱估计的阵元间相位差进行加权融合,综合利用目标各线谱的阵元间相位信息,计算更为精确的目标信号宽带时延矢量,获取保真增强的重构源信号。

附图说明

图1为本发明方法的流程图;

图2为实施例中畸变阵与无畸变阵的阵元坐标位置对比图;

图3为实施例中目标及各阵元位置示意图;

图4为实施例中阵形及目标距离已知时的波束能量图与基于目标信号远场平面波入射及阵形无畸变假设的波束能量图;

图5为实施例中参考阵元采集数据的功率谱;

图6为实施例中以不同方式计算的目标信号阵元间数字时延对比图;

图7为实施例中重构源信号的功率谱线谱特征对比图;

图8为实施例中重构源信号的解调谱线谱特征对比图。

具体实施方式

下面结合附图和实施例对本发明的技术方案作进一步的说明。

如图1所示,一种用于高保真阵列处理的宽带时延估计方法,包括如下步骤:

(1)基于目标信号远场平面波入射和阵形无畸变两个假设,计算基阵采集数据不同方位的预成波束数据,求取波束能量图,估计目标方位

(2)基于目标信号远场平面波入射和阵形无畸变两个假设,从基阵采集数据中按估计的目标方位对源信号进行预重构,检测预重构源信号s1[n],n=1,2,...,n的功率谱线谱位置fk,k=1,2,...,k,计算线谱频率的精测值其中k为检测出的频率小于v/(2d)的功率谱线谱根数,v为声波在水中的传播速度,d为均匀线阵阵元间距;

(3)计算参考阵元h1采集数据的功率谱在位置处线谱的频域信噪比snrk,k=1,2,...,k;

(4)计算基阵采集数据的频谱在位置处线谱的阵元间相位差其中,m为基阵阵元数目;

(5)根据线谱频率和频域信噪比对基于不同频率线谱估计的阵元间相位差进行加权融合,估计目标信号宽带时延矢量

(6)根据估计的目标信号宽带时延矢量从基阵采集数据中对源信号进行二次重构,获取保真增强的重构源信号s2[n],n=1,2,...,n。

步骤(1)中,采用如下方法计算基阵采集数据的波束能量图,估计目标方位,具体包括如下步骤:

(1-1)计算按等余弦划分的预成波束角度:

其中,θi为第i个预成波束的角度,以基阵第一个阵元为坐标原点,建立平面直角坐标系,基阵第m个阵元的坐标为[-(m-1)d,0],m=1,2,...,m,目标角度定义为原点和目标坐标的连线与x轴的夹角,θi∈[0,π],表示计算的反余弦,i为预成波束个数,d为均匀线阵阵元间距;

(1-2)基于目标信号远场平面波入射和阵形无畸变两个假设,计算均匀线阵相邻阵元采集信号在预成波束角度θi下的时延差τi:

其中,cos(θi)表示计算θi的余弦值,v为声波在水中的传播速度;

(1-3)按如下波束形成算法;以基阵阵首第一个阵元h1为参考阵元,计算预成波束角度θi对应的预成波束数据

其中,m为基阵阵元数目,xm[n]表示第m个阵元第n个采样数据,n为预成波束数据离散点数,fs为基阵采集数据时的采样频率;

(1-4)计算基阵采集数据在预成波束角度θi下的波束能量e[i]:

(1-5)估计目标方位

搜索e[i]最大值所对应的预成波束索引imax

其中,argmax1≤i≤i{e[i]}表示在1≤i≤i范围内搜索e[i]最大值所对应的预成波束索引imax。若1<imax<i,将分别记作z-1、z0、z1,则基于抛物线三点插值估计的目标方位为:

在步骤(2)中,采用如下方法检测预重构源信号的功率谱线谱位置fk,计算线谱频率的精测值具体包括如下步骤:

(2-1)按如下波束形成算法,以基阵阵首第一个阵元h1为参考阵元,从基阵采集数据中按估计的目标方位对源信号进行预重构,获取预重构源信号s1[n]:

(2-2)计算预重构源信号s1[n]的离散频谱:

则预重构源信号s1[n]的功率谱为:

其中,l为s1[l]和p1[l]的离散频率索引,j表示虚数单位,即||代表取模值运算;

(2-3)将功率谱p1[l]通过q阶中值滤波器得到功率谱p1[l]的趋势项c1[l],l=0,1,...,n/2,则功率谱p1[l]与趋势项c1[l]的差值谱d1[l]为:

d1[l]=p1[l]-c1[l],l=0,1,...,n/2

d1[l]的规范化差值谱为:

其中,std(d1[l])是求d1[l]的标准差;

(2-4)设定幅度门限g提取功率谱线谱,若d2(l)满足如下条件则判为功率谱线谱:

假设共提取出k根频率小于v/2d的功率谱线谱,第k根线谱在d2[l]中的索引为jk,则第k根线谱的频率为:

fk=jkδf,k=1,2,...,k

其中,δf为长度为n的离散傅里叶变换的频率分辨率,δf=fs/n;

(2-5)取s1[l]在离散频率索引(jk-1),jk,(jk+1)处的模值,分别记作ak1,ak2,ak3,即ak1=|s1[jk-1]|,ak2=|s1[jk]|,ak3=|s1[jk+1]|,利用所述模值ak1,ak2,ak3计算第k根功率谱线谱频率插值的相对偏差δk,即:

则第k根功率谱线谱的频率精测值为:

在步骤(3)中,采用如下方法计算参考阵元h1采集数据功率谱位置处线谱的频域信噪比snrk,具体包括如下步骤:

(3-1)计算参考阵元h1采集数据x1[n]的功率谱:

其中,l为的离散频率索引;

(3-2)估计参考阵元h1采集数据x1[n]的功率谱在位置处噪声的功率值

其中,wl为噪声功率值估计器的长度,取偶数,且wl≤2(jk-1);

(3-3)计算参考阵元h1采集数据x1[n]的功率谱在位置处信号加噪声的功率值

(3-4)参考阵元h1采集数据x1[n]的功率谱在位置处线谱的频域信噪比snrk为:

优选的,在步骤(4)中,采用如下方法计算基阵采集数据的频谱在位置处的阵元间相位差具体包括如下步骤:

(4-1)计算第m个阵元hm采集数据的频谱在位置处的值:

(4-2)计算各阵元采集数据的频谱在位置处的相位矢量

其中,为求取相位的运算;

(4-3)计算相邻阵元采集数据的频谱在位置处的相位之差

其中,为矢量的第m个元素,m=1,2,...,m;

在步骤(5)中,采用如下方法对基于不同频率线谱估计的目标信号阵元间相位差进行加权融合,估计目标信号宽带时延矢量具体包括如下步骤:

(5-1)根据线谱频率和频域信噪比snrk,计算第k根线谱的阵元间相位差加权系数wk:

(5-2)根据阵元间相位差加权系数wk对阵元间相位差进行加权融合,估计目标信号宽带时延矢量

其中

在步骤(6)中,采用如下方法对源信号进行二次重构:根据估计的目标信号宽带时延矢量从基阵采集数据xm[n]中对源信号进行二次重构,获取保真增强的重构源信号s2[n]:

实施例:

仿真信号参数分别设置为:均匀线阵阵元数目m=32,阵元间距d=0.375m,畸变阵与无畸变阵阵元坐标位置对比如图2所示。目标及无畸变阵阵元坐标位置示意如图3所示,以阵首第一个阵元为坐标原点,目标与阵首方向的夹角θ=120°,目标与参考阵元的距离为r=3000m。采样频率fs=10khz,声波在水中的传播速v=1500m/s。单阵元接收数据信噪比为-9.62db,包含100.08hz、300.1hz、500.2hz、800.3hz、1000.41hz、1200.7hz、1500.5hz、1801.5hz、1999.5hz共9根线谱,时间分析窗长为8s时,对应的频域信噪比分别为-3db、12db、0db、16db、18db、-3db、-3db、-3db、10db;设定2hz、4hz、6hz、8hz共4根调制谱线谱,调制深度分别为0.04、0.04、0.04、0.06。

在第(1)步中,设置预成波束个数i=65,基于目标信号远场平面波入射和阵形无畸变两个假设,求取的波束能量图如图4所示,估计的目标方位为

依据第(2)步,设置中值滤波器的阶数为q=51,规范化差值谱线谱提取幅度门限g=20。检测出的功率谱线谱根数为k=4,4根功率谱线谱的频率分别为

基于插值估计的4根功率谱线谱的频率精测值分别为

依据第(3)步,参考阵元采集数据的功率谱如图5所示,计算线谱频域信噪比时,参与平均的频点个数wl=80,考阵元h1采集数据的功率谱对应频率线谱的频域信噪比分别为

依据第(5)步,4根功率谱线谱频率的阵元间相位差加权系数分别为

根据加权系数对4根线谱的阵元间相位差进行加权融合,估计目标信号宽带时延矢量以不同方式计算的目标信号阵元间数字时延如图6所示,可以看出相对于基于目标信号远场平面波入射和阵形无畸变假设估计的阵元间数字时延,基于目标信号多线谱阵元间相位差加权融合估计的阵元间数字时延与目标信号阵元间数字时延真值具有更好的一致性。

依据第(6)步,根据估计的目标信号宽带时延矢量从基阵采集数据中对源信号进行二次重构,获取保真增强的重构源信号。图7为预重构源信号功率谱与二次重构源信号功率谱线谱特征对比图,可以看出,相较于预重构源信号,二次重构源信号的功率谱线谱特征保真度明显增强。对于可以从预重构源信号检测出的功率谱线谱,从二次重构源信号检测出的功率谱线谱具有更高的信噪比,且可从二次重构源信号检测出无法从预重构源信号检测出的功率谱线谱。图8为预重构源信号解调谱与二次重构源信号解调谱线谱特征对比图,可以看出,相较于预重构源信号,二次重构源信号的解调谱结构更清晰、线谱特征信噪比更高。

当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1