发明内容
本发明的目的是提供一种可分辨多个同时到达的回波信号的传感器阵列回波方向估计方法,同时提供一种应用该方法的高分辨力的多波束回波测深底检测方法。
为实现上述目的,本发明所采取的技术方案是:其传感器阵列回波方向估计方法包括以下步骤:
(1)对传感器阵列的接收数据进行波束扫描及回波检测,确定回波对应的空间区域,在该空间区域内选取两个以上期望角度;
(2)执行方案一或方案二:
方案一:对所述传感器阵列进行划分得到三个以上子阵,并得到各子阵之间的延时关系;对每个子阵的接收数据分别进行步骤(1)所述各期望角度上的波束驾驶,得到每个子阵在所述各期望角度上的输出波束;
方案二:对所述传感器阵列的接收数据进行步骤(1)所述各期望角度上的波束驾驶,得到所述传感器阵列在所述各期望角度上的输出波束;对所述传感器阵列在所述各期望角度上的输出波束进行虚拟变换而形成三个以上虚拟子阵,并得到各虚拟子阵在所述各期望角度上的输出波束以及各虚拟子阵之间的延时关系;
(3)根据步骤(2)所述方案一的各子阵或方案二的各虚拟子阵之间的延时关系,对步骤(2)中得到的各子阵或各虚拟子阵在所述各期望角度上的输出波束的数据进行排列,采用基于多子阵分解的高分辨力信号到达方向角估计算法获取各回波的方位估计。
进一步地,本发明所述步骤(1)中,所述各期望角度或者在所述回波对应的空间区域内均匀分布,或者在所述回波对应的空间区域内随机分布。
进一步地,本发明所述步骤(2)中,所述方案一的各子阵或者无重叠相邻分布,或者无重叠间隔分布,或者重叠分布。
本发明的多波束回波测深底检测方法包括以下步骤:
(1)对传感器阵列的接收数据进行波束扫描及回波检测,估计每个扫描波束内的回波个数,并执行方案三或方案四:
方案三:如果所述回波个数为一个,则估计所述回波的到达时刻和方位,并执行步骤(5);
方案四:如果所述回波个数大于一个,则执行步骤(2);
(2)根据步骤(1)所述扫描波束的输出数据,确定回波对应的空间区域,并在该空间区域内选取三个以上期望角度;
(3)执行方案一或方案二:
方案一:对所述传感器阵列进行划分得到三个以上子阵,并得到各子阵之间的延时关系;对每个子阵的接收数据分别进行步骤(1)所述各期望角度上的波束驾驶,得到每个子阵在所述各期望角度上的输出波束;
方案二:对所述传感器阵列的接收数据进行步骤(1)所述各期望角度上的波束驾驶,得到所述传感器阵列在所述各期望角度上的输出波束;对所述传感器阵列在所述各期望角度上的输出波束进行虚拟变换而形成三个以上虚拟子阵,并得到各虚拟子阵在所述各期望角度上的输出波束以及各虚拟子阵之间的延时关系;
(4)根据步骤(3)所述方案一的各子阵或方案二的各虚拟子阵之间的延时关系,在回波到达时刻对步骤(3)中得到的各子阵或各虚拟子阵在所述各期望角度上的输出波束的数据进行排列,采用基于多子阵分解的高分辨力信号到达方向角估计算法获取各回波的方位估计;
(5)对所述回波到达时刻和各回波的方位估计结果进行运动姿态修正和声速修正,生成水底各测量点的位置与深度信息。
进一步地,本发明所述步骤(2)中,所述回波对应的空间区域是扫描方向上的输出波束时间维上最大值对应角度的组合,或者是扫描方向上的输出波束时间维上最大值对应角度的平均值的邻近区域。
进一步地,本发明所述步骤(2)中,所述各期望角度或者在所述回波对应的空间区域内均匀分布,或者在所述回波对应的空间区域内随机分布,或者是所述扫描方向上的输出波束在回波对应的空间区域内的角度。
进一步地,本发明所述步骤(3)中,所述方案一的各子阵或者无重叠相邻分布,或者无重叠间隔分布,或者重叠分布。
进一步地,本发明所述步骤(4)中,所述回波到达时刻之间的间隔为传感器阵列接收数据的采样时间间隔的p倍,p≥1且p为整数。
与现有技术相比,本发明方法具有如下有益效果:
(1)将基于多子阵分解的高分辨力信号到达方向角估计算法与常规波束形成技术相结合,能够提高处理信噪比,减少有效阵元个数,进而降低高分辨力处理算法对于数据样本点数的要求及计算要求;与现有处理技术相比,本发明中的方法适用性更强,应用于多子阵的输出波束的高分辨力处理算法的操作使得本发明方法在较少样本点、较低信噪比的情况下,能得到更精确、更稳定的回波方位估计,同时计算速度快,对传感器阵列的要求低,便于工程应用。
(2)将高分辨力的回波方向估计方法应用于多波束回波测深系统中,可以突破阵列孔径对分辨力的限制,分辨出在单个波束宽度内同时到达的多个方向的回波信号,提高尤其是边缘波束的测深精度;
(3)对单个形成波束中的回波个数进行估计,并根据估计结果进行常规或高分辨的相应处理,在提高测深精度的同时保证计算效率,有利于数据的实时处理。
具体实施方式
本发明的基本构思为:提供一种高分辨力的传感器阵列回波方向估计方法,此方法结合多子阵处理技术、虚拟变换技术和波束驾驶技术,能获得更高精度的多个方向回波的方位估计,同时降低计算的复杂性,有利于实时应用。同时将此方法应用于多波束回波测深系统信号处理机中,得到一种能分辨多个方向回波的底检测方法,实现对底部回波信号的高空间分辨力处理。多波束系统的典型组成如图1所示。信号处理机的主要功能是实现传感器阵列接收信号的实时采集与处理,将采集数据和处理结果通过高速通讯口传给显示/控制平台,同时协调显示/控制平台与系统其它设备之间的操作。
下面结合附图以及具体实施方式对本发明做进一步的详细描述:
如图4所示,本发明提供的传感器阵列回波方向估计方法的处理模块包括:
(1)波束形成单元,用于得到传感器阵列扫描方位上的输出波束;
(2)区域划分单元,用于根据波束形成单元的输出数据,通过回波检测得到回波对应的空间区域,并选取N个期望角度,且N≥2;
(3)子阵划分单元,用于对传感器阵列进行划分,得到M个子阵,且M≥3;
(4)子阵构造单元,用于对传感器阵列进行虚拟变换,得到M个虚拟子阵,且M≥3;
(5)波束驾驶单元,用于得到各子阵或各虚拟子阵在各期望角度上的输出波束;
(6)高分辨力信号到达方向估计算法单元,用于对各子阵或各虚拟子阵在各期望角度上的输出波束进行处理,得到回波到达方向的估计。
其中,所述区域划分单元包括:
1)回波空间区域单元,用于根据波束形成单元的输出波束进行回波检测,得到回波对应的空间区域;
2)期望角度选取单元,用于在回波对应的波束方向及其邻域内选取N个期望角度,且N≥2。
所述子阵划分单元包括:
1)子阵间延时关系单元,用于产生各子阵之间的延时关系;
2)子阵接收数据单元,用于产生各子阵的接收数据。
本发明所述子阵构造单元包括:
1)虚拟子阵间延时关系单元,用于产生各虚拟子阵之间的延时关系;
2)虚拟转换单元,用于将传感器阵列在所述期望角度上的输出波束转换为各虚拟子阵在所述期望角度上的输出波束。
所述高分辨力信号到达方向估计算法单元包括:
1)数据合成单元,用于将输入数据进行排列,得到所述矩阵束;
2)方向个数单元,用于根据数据合成单元输出的矩阵束,估计回波方向的个数;
3)相位差空间单元,用于根据数据合成单元输出的矩阵束,结合方向个数单元输出的回波方向的个数,得到所述多个方向回波的相位差空间;
4)回波方向角估计单元,用于根据相位差空间单元的输出结果,结合所述各子阵或各虚拟子阵之间的延时关系,得到各回波方向的角度值。
本发明提供的传感器阵列回波方向估计方法的处理流程包括以下步骤:
(1)通过波束形成单元,对传感器阵列的接收数据进行波束扫描及回波检测,通过区域划分单元确定回波对应的空间区域,在该空间区域内选取两个以上期望角度;
在本步骤中,回波对应的空间区域是从扫描方位上的输出波束中得到,回波对应的空间区域有多种获得方式,如每个时刻的扫描方位上的输出波束最大值对应角度的组合、或者扫描方位上的输出波束最大值对应角度的时间均值的临近区域。本发明优选所述回波对应的空间区域为回波对应的波束方向及其邻域,但对于其他的获得方式,本发明同样适用。
在本步骤中,所述期望角度有多种选取方式,例如:在所述回波对应的空间区域内均匀选取、在所述回波对应的空间区域内随机选取或者从所述扫描方位上的输出波束在回波对应的空间区域内的对应角度中选取。期望角度个数的增加,能提高回波方位的估计精度,但同时也会增加计算量。本发明优选两个以上的期望角度,但对于其他的期望角度个数,本发明同样适用。
(2)执行方案一或方案二:
方案一:通过子阵划分单元,对传感器阵列进行划分得到三个以上子阵,并得到各子阵之间的延时关系;通过波束驾驶单元对每个子阵的接收数据分别进行步骤(1)所述各期望角度上的波束驾驶,得到每个子阵在各期望角度上的输出波束;
在本步骤中,子阵的划分方式可在图3所示的三种方式中任选一种,其中,(a)为不重叠的相邻划分方式,此方式计算量较小,适用于较短的线列阵;(b)为重叠的划分方式,此方式能增加阵元的利用率,但是会增加计算量;(c)为不重叠的间隔划分方式。
在本步骤中,子阵的个数可以是大于二的任意数。子阵的个数越大,计算量也越大;但子阵的个数较少时,分辨率可能无法达到要求。本发明中优选子阵的个数为三个以上,但对于其他个数,本发明也是适用的。
方案二:通过波束驾驶单元,对传感器阵列的接收数据进行步骤(1)所述各期望角度上的波束驾驶,得到传感器阵列在各期望角度上的输出波束;通过子阵构造单元,对传感器阵列在各期望角度上的输出波束进行虚拟变换,形成三个以上虚拟子阵,并得到各虚拟子阵在各期望角度上的输出波束以及各虚拟子阵之间的延时关系;
在本步骤中,虚拟子阵的个数可以是大于二的任意数。虚拟子阵的个数越大,计算量也越大;但虚拟子阵的个数较少时,分辨率可能无法达到要求。本发明中优选虚拟子阵的个数为三个以上,但对于其他个数,本发明也是适用的。
(3)根据步骤(2)所述方案一的各子阵或方案二的各虚拟子阵之间的延时关系,通过高分辨力信号到达方向估计算法单元,对步骤(2)中得到的各子阵或各虚拟子阵在各期望角度上的输出波束的数据进行排列,采用基于多子阵分解的高分辨力信号到达方向角估计算法获取各回波的方位估计。
在本步骤中,高分辨力信号到达方向估计算法为基于多子阵分解的高分辨力信号到达方向角估计算法,本发明中优选基于旋转不变技术的信号参数估计方法(estimation of signal parameters via rotational invariance techniques,简写为ESPRIT)。
本发明的多波束回波测深底检测方法是以上本发明传感器阵列回波方向估计方法在多波束回波测深系统中的具体应用。如图5所示,本发明多波束回波测深底检测方法的处理模块包括:
(1)波束形成单元,用于得到传感器阵列扫描方位上的输出波束;
(2)回波判定单元,用于根据波束形成单元的输出数据,得到每个扫描波束内回波个数的估计;
(3)区域划分单元,用于根据波束形成单元的输出数据,通过回波检测得到回波对应的空间区域,并选取N个期望角度,且N≥3;
(4)子阵划分单元,用于对传感器阵列进行划分,得到M个子阵,且M≥3;
(5)子阵构造单元,用于对传感器阵列进行虚拟变换,得到M个虚拟子阵,且M≥3;
(6)波束驾驶单元,用于得到各子阵或各虚拟子阵在各期望角度上的输出波束;
(7)高分辨力信号到达方向估计算法单元,用于对各子阵或各虚拟子阵在各期望角度上的输出波束进行处理,得到回波到达方向的估计。
其中,所述区域划分单元包括:
1)回波空间区域单元,用于根据波束形成单元的输出波束进行回波检测,得到回波对应的空间区域;
2)期望角度选取单元,在回波对应的波束方向及其邻域内选取N个期望角度,且N≥3。
所述子阵划分单元包括:
1)子阵间延时关系单元,用于产生各子阵之间的延时关系;
2)子阵接收数据单元,用于产生各子阵的接收数据。
所述子阵构造单元包括:
1)虚拟子阵间延时关系单元,用于产生各虚拟子阵之间的延时关系;
2)虚拟转换单元,用于将传感器阵列在所述期望角度上的输出波束转换为各虚拟子阵在所述期望角度上的输出波束。
所述高分辨力信号到达方向估计算法单元包括:
1)数据合成单元,用于将输入数据进行排列,得到所述矩阵束;
2)方向个数单元,用于根据数据合成单元输出的矩阵束,估计回波方向的个数;
3)相位差空间单元,用于根据数据合成单元输出的矩阵束,结合方向个数单元输出的回波方向的个数,得到所述多个方向回波的相位差空间;
4)回波方向角估计单元,用于根据相位差空间单元的输出结果,结合所述各子阵或各虚拟子阵之间的延时关系,得到各回波方向的角度值。
本发明多波束回波测深底检测方法具体包括以下步骤:
(1)通过波束形成单元,对传感器阵列的接收数据进行波束扫描及回波检测;通过回波判定单元,估计每个扫描波束内的回波个数,并执行方案三或方案四:
方案三:如果回波个数为一个,则进行回波的到达时刻和方位的估计,并执行步骤(5);
方案四:如果回波个数为两个以上,则执行步骤(2);
(2)根据步骤(1)所述扫描波束的输出数据,通过区域划分单元,确定回波对应的空间区域,并在该空间区域内选取三个以上期望角度;
在本步骤中,回波对应的空间区域是从扫描方位上的输出波束中得到,所述回波对应的空间区域有多种获得方式,如每个时刻的扫描方位上的输出波束最大值对应角度的组合、或者扫描方位上的输出波束最大值对应角度的时间均值的临近区域。本发明优选回波对应的空间区域为回波对应的波束方向及其邻域,但对于其他的获得方式,本发明同样适用。
在本步骤中,期望角度有多种选取方式,例如:在回波对应的空间区域内均匀选取、在回波对应的空间区域内随机选取或者从扫描方位上的输出波束在回波对应的空间区域内的对应角度中选取。期望角度个数的增加,能提高回波方位的估计精度,但同时也会增加计算量。本发明优选三个以上的期望角度,但对于其他的期望角度个数,本发明同样适用。
(3)执行方案一或方案二:
方案一:通过子阵划分单元,对传感器阵列进行划分得到三个以上子阵,并得到各子阵之间的延时关系;通过波束驾驶单元,对每个子阵的接收数据分别进行步骤(1)所述各期望角度上的波束驾驶,得到每个子阵在各期望角度上的输出波束;
在本步骤中,子阵的划分方式主要有三种,如图3所示:(a)为不重叠的相邻划分方式,此方式计算量较小,适用于较短的线列阵;(b)为重叠的划分方式,此方式能增加阵元的利用率,但是会增加计算量;(c)为不重叠的间隔划分方式。
在本步骤中,子阵的个数可以是大于二的任意数。子阵的个数越大,计算量也越大;但子阵的个数较少时,分辨率可能无法达到要求。本发明中优选子阵的个数为三个以上,但对于其他个数,本发明也是适用的。
方案二:通过波束驾驶单元,对传感器阵列的接收数据进行步骤(1)所述各期望角度上的波束驾驶,得到传感器阵列在各期望角度上的输出波束;通过子阵构造单元,对传感器阵列在各期望角度上的输出波束进行虚拟变换,形成三个以上虚拟子阵,并得到各虚拟子阵在各期望角度上的输出波束以及各虚拟子阵之间的延时关系;
在本步骤中,虚拟子阵的个数可以是大于二的任意数。虚拟子阵的个数越大,计算量也越大;但虚拟子阵的个数较少时,分辨率可能无法达到要求。本发明中优选虚拟子阵的个数为三个以上,但对于其他个数,本发明也是适用的。
(4)根据步骤(3)所述方案一的各子阵或方案二的各虚拟子阵之间的延时关系,通过高分辨力信号到达方向估计算法单元,在回波到达时刻对步骤(3)中得到的各子阵或各虚拟子阵在各期望角度上的输出波束的数据进行排列,采用基于多子阵分解的高分辨力信号到达方向角估计算法获取各回波的方位估计;
在本步骤中,高分辨力信号到达方向估计算法为基于多子阵分解的高分辨力信号到达方向角估计算法,本发明中优选ESPRIT高分辨算法。
(5)根据传感器接收到的姿态以及声速数据,对回波到达时刻和各回波到达方向的估计结果进行运动姿态修正和声速修正,生成水底各测量点的位置与深度信息,并将各测量点的位置与深度信息进行存储。
在本步骤中,回波到达时刻之间的间隔为传感器阵列接收数据的采样时间间隔的p倍,p≥1且p为整数。因此不需要对传感器阵列的接收数据进行插值处理,简化操作流程,减少计算量。
如图6所示,本实施例是本发明多波束回波测深底检测方法在阵元均匀分布的线列阵上的具体应用,假定存在两个以上回波情况时,其具体步骤如下:
(1)传感器阵列为阵元均匀分布的线列阵,得到线列阵各阵元上的时域数据;
本步骤中,阵元均匀分布的线列阵如图2所示。本实施例中选用的传感器阵列为阵元均匀分布的线列阵,但本发明可应用于任意形状的平面阵,如阵元非均匀分布的线列阵、圆弧阵、矩形阵等。
(2)对该线列阵的接收数据进行波束扫描,确定回波对应的空间区域,并在回波对应的波束方向及其邻域内选取N个期望角度,且N≥3;将该线列阵的接收数据作为下个步骤的输入;
本实施例中选用扫描方位上的输出波束最大值对应角度的时间均值的临近区域作为所述回波对应的空间区域,但对于其他的获得方式,本发明同样适用。
本实施例中采用期望角度是从所述扫描方位上的输出波束在回波对应的空间区域内的对应角度中选取的方式,但对于其他的选取方式,本发明同样适用。本实施例中选取N的值为4,可以衡量在期望角度个数较小时本发明的估计精度。
(3)将该线列阵划分为M个子阵,且M≥3,得到各子阵之间的延时关系和各子阵的接收数据;将各子阵的接收数据作为下个步骤的输入;当采用方案二时,该步骤跳过;
本实施例中采用(b)重叠的划分方式,对于其他的划分方式,本发明同样适用。本实施例中选取M的值为3,两个相邻子阵的阵中心之间的距离为一个阵元间距,对于其他的设定本发明同样适用。
(4)对输入数据进行所述各期望角度上的波束驾驶处理,得到在所各述期望角度上的输出波束,并将该输出波束数据进行存储;
(5)得到虚拟变换矩阵。当采用方案一时,该步骤跳过;
具体地说,本步骤中,虚拟变换矩阵由以下式子定义:
这里
为传感器阵列的阵列流型矩阵,
是传感器阵列在
方向上的驾驶向量;
为虚拟子阵的阵列流型矩阵,
是虚拟子阵在
方向上的驾驶向量,
表示虚拟变换矩阵,符号
代表共轭转置,
表示第n个所述期望角度。根据所需虚拟子阵的个数,确定形成的虚拟变换矩阵的个数。
(6)根据由步骤(5)所得的虚拟变换矩阵,对存储的输出波束数据进行变换,得到各虚拟子阵在所述期望角度上的输出波束,并将变换的结果替换掉原存储数据进行存储。当采用方案一时,该步骤跳过;
具体地说,本步骤中,变换方式由以下式子定义:
这里
为传感器阵列的输出波束,
为虚拟子阵的输出波束,
表示虚拟变换矩阵,符号
代表共轭转置。
(7)读取存储的数据,根据数据对应的样本时刻得到回波的到达时刻并存储;在每个回波到达时刻将存储的波束域数据进行排列,构造为矩阵束,并计算该矩阵束的协方差矩阵;
具体地说,本步骤中,矩阵束及其协方差矩阵的表达式由以下式子定义:
这里
为第m个子阵或虚拟子阵的输出波束,
为所述矩阵束,
为该矩阵束的协方差矩阵,符号
表示求期望,符号
代表共轭转置。
(8)对步骤(7)所得的协方差矩阵进行特征值分解,估计回波的个数,并获得每个子阵的信号子空间的估计;根据不同的准则,得到旋转变换矩阵,通过对该旋转变换矩阵进行特征值分解,得出特征值的估计值;
本步骤是高分辨力信号到达方向估计算法单元的具体实现,本实施例中采用的是ESPRIT算法,这只是给出了其中的一种,其它基于多子阵分解的算法也是可行的。ESPRIT算法的准则主要有两种:最小二乘(least squares,简写为LS)准则和总体最小二乘(total least squares,简写为TLS)准则。本实施例中采用的是LS准则。
(9)获得各子阵或各虚拟子阵之间的延时关系的表达式;
本步骤中,对于真实子阵和虚拟子阵的处理方式相同。首先选定一个参考阵列,计算每个子阵与该参考阵列之间,在所述各期望角度下的延时关系,该延时关系即为所述各子阵之间的延时关系。通常选取一个子阵为参考阵列,则各子阵之间的延时关系对应为其它子阵与该子阵之间,在所述各期望角度下的延时。本实施例中选取第一个子阵为参考阵列。
具体地说,本步骤中,各子阵之间的延时关系由以下式子定义:
这里选第一个子阵为参考,
表示第m+1个子阵与该参考阵列之间的延时关系,
为第m+1个子阵的首阵元与该参考阵列的首阵元之间的距离,d是阵元间距,
表示第n个所述期望角度,
表示声速,
表示以输入向量为对角元素构造矩阵。
(10)根据由步骤(9)所得的各子阵或各虚拟子阵之间的延时关系的表达式,结合由步骤(8)所得的特征值,得到回波方向的相位差空间的估计;
(11)计算回波到达方向的估计,并存储估计结果;
(12)读取存储的回波到达时刻和到达方向,结合相对应的运动姿态数据以及声速数据,将其转换为各测量点的位置以及深度信息,并根据需求进行存储或显示。
仿真条件为:线列阵的阵元数为64,阵元间距为半波长,划分的子阵个数为3,每个子阵的阵元数为62,两个相邻子阵之间的距离为一个阵元间距。在25.6°和26.7°方向有回波同时到达线列阵,不同方向的回波信号是相互独立的,均与噪声不相关,且信噪比相同。回波信号的样本点数为10,选取的期望角度个数为4。图7是在信噪比较高为35dB情况下常规波束形成的处理结果,波束主瓣对应着回波的方向,图中只有一个主瓣,因此常规波束形成算法无法分辨多个方向同时到达的回波信号。 图8是ESPRIT高分辨算法的方位估计结果,可以看到在小信噪比(SNR)的情况下,估计的偏差较大,随着SNR的增大,估计的偏差减小,但与仿真设定仍有一定的差距。图9是本发明实施例的传感器阵列高分辨力回波信号到达方向估计方法的处理结果,可以看到随着信噪比的增加,本发明的估计结果接近仿真设定并吻合。将图8与图9进行比较可以发现:在较小SNR时本发明的估计偏差较小;在较大SNR时本发明能提供回波方向的精确估计而ESPRIT算法的估计结果有偏差。因此本发明与ESPRIT算法相比,在信噪比较低以及样本点数较少的情况下,能提供更好的分辨性能,对回波方位的估计精度更高。
需要说明的是,本发明中的子阵可以是真实子阵,也可以是虚拟子阵,不同的子阵形式对应不同的处理方案。方案一为真实子阵的情况,处理步骤中不包括步骤(5)和步骤(6),其余步骤不变;而方案二为虚拟子阵的情况,处理步骤中不包括步骤(3),其余步骤不变。
最后应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管实施例对本发明进行了详细说明,本领域的普通技术人员都应当了解,对本发明的技术方案进行修改或者等价替换,都不脱离本发明技术方案的精神和范围,均涵盖于本发明的权利要求范围内。