专利名称：一种相关测速声纳测量载体对底速度的方法和系统的制作方法
技术领域：
本发明涉及用声纳测速领域，更具体地说，涉及一种相关测速声纳测量载体对底速度的方法和系统。
背景技术：
目前声相关测速声纳的测速方法概括如下(1)S.E.Bradley等人的美国专利5315562“Correlation Sonar System”中介绍了相关声纳，可用于测量流速的剖面和水中载体对底的速度。该发明的主要内容有4个方面(A)发射一复杂信号，它的自相关函数在延时τ＝0和τ＝τc时各有一个幅值不等的峰，避免采用以前的发射两个脉冲的技术，后者在测量流速时会引起流体介质层间的干扰。
(B)以级数形式给出了声纳阵时空相关函数对于流体介质和底介质的理论表式，表式中包括了贝赛尔函数和勒让德函数，以及由经验得出的简易表式。在进行数据处理时，因为前者复杂，所以用的是后者。
(C)基于最大似然原理，采用单纯形方法，将理论的和实验数据获得的声纳阵时空相关函数进行最佳拟合，求得流体的速度和载体对底的速度。
(D)采用匹配滤波器的方法检测海底回波。
(E)声纳阵中部分发射基元与接收基元合置。
(2)S.E.Bradley等人的美国专利5422860“Correlation Sonar System”中详细介绍了相关声纳发射信号产生的方法，发射调相伪随机信号，其自相关函数在延时τ＝0和τ＝τc时各有一个幅值不等的峰。
上述的方法在测量载体对底的速度时，存在明显的不足(1)它的声纳阵时空相关函数的理论表式相当复杂，实际中难以应用，而由经验得出的简易表式缺乏足够的物理基础。这是相关测速声纳系统的最关键技术。(2)用单纯形将理论的和实验的时空相关函数拟合，并不是最好的方法。(3)采用与声纳阵时空相关函数最大值对应的速度作为速度估计的初始值，这不是最好的方法。(4)部分发射基元与接收基元合置，限制了发射波束和接收波束的宽度选择。
本申请人在申请号为“03119665.9”的中国专利申请中提供了一种“相关测速声纳测量载体对底速度的方法及其系统”，该方法主要包括1)对底介质回波信号进行解调滤波，并计算底介质时空相关函数矩阵；2)设置未知参数系综 的搜索范围，其中V1x、V1y分别是x、y方向载体对底的速度均值，γ是宽度因子；3)在未知参数系综 的搜索范围内，基于最大似然原理，采用序列二次规划方法，将底介质相关函数矩阵与理论的底介质声纳阵时空相关函数进行拟合运算获得载体对底的速度均值{V1x，V1y}。在该申请中提出一种新的底介质声纳阵时空相关函数 其中B是常数，θ是声波入射角，τ是延时，d是相邻声纳阵基元间的间距，J0是零阶贝塞尔函数；β2=ω0c((τV‾1x+dx)2+(τV‾1y+dy)2)1/2,]]>其中ω0是发射信号的中心频率，c是声速，dx和dy是d在x和y方向的分量。
该方法中的底介质声纳阵时空相关函数不仅适用远场区，即平面波区，而且适用于夫琅和费区，即球面波区。而且该底介质声纳阵时空相关函数用零阶贝赛尔(Bessel)函数简洁表示，与实验符合良好。但是，在该申请中仍然存在明显的不足(1)它采用最大似然原理估计速度，在实际工作情况中，环境噪声在空间上往往是不均匀的，声纳阵基元的幅度和相位存在不一致性，因此最大似然原理估计方法的性能有明显的退化；最大似然原理计算相当复杂，计算量相当大；(2)现有的信号计算方法采用非局域化的方法，在时空函数的整个区域搜索，包括信噪比低的区域，这使测速精度下降，有时下降明显；(3)它没有指出如何使公式(1)中的B成为常数。此外，该方法中没有给出声纳阵基元的分布。
发明内容
本发明的目的之一是对现有技术的测量载体对底速度方法中的拟合方法进行改进；本发明的另一个目的是对现有技术的测量载体对底速度方法中的矩阵数据处理方法进行改进；本发明的再一个目的是对现有技术中声纳阵基元分布进行改进。
为了实现上述发明目的，本发明提供一种相关测速声纳测量载体对底速度的方法，包括以下步骤(1)选择发射编码，该发射编码在相关延时非零时具有峰值。
(2)根据发射编码发射声脉冲到流体介质中，并接收各流层或底介质反向散射的回波信号。
(3)如步骤(2)中的回波信号包含底介质回波信号，则进行步骤(4)；如未包含，则返回步骤(1)重新开始。
(4)对所述底介质回波信号进行解调滤波，并根据底介质回波信号计算数据底介质时空相关函数矩阵。
(5)将数据底介质时空相关函数矩阵中的所有元素作绝对值运算，获得数据底介质时空相关函数绝对值矩阵，该矩阵中元素的最大值为EMax。
(6)设定一门限值x，其中0＜x≤1，优选0.7＜x≤1；在所述数据底介质时空相关函数绝对值矩阵中，将数值小于xEMax的元素置零，数值大于或等于xEMax的元素保持不变，对所有元素进行此操作后获得局域化的底介质时空相关函数绝对值矩阵。
(7)设置未知参数系综 的搜索范围，其中V1x、V1y分别是x、y方向载体对底的速度均值，γ是宽度因子。
(8)在未知参数系综 的搜索范围内，基于非线性最小二乘原理，采用序列二次规划方法，将步骤(6)中局域化的底介质时空相关函数绝对值矩阵与理论的底介质时空相关函数的绝对值矩阵进行拟合运算。其中，所述理论的底介质声纳阵时空相关函数绝对值为 其中B是常数，θ是声波入射角，τ是延时，d是相邻声纳阵基元间的间距，J0是零阶贝塞尔函数；β2=ω0c((τV‾1x+dx)2+(τV‾1y+dy)2)1/2,]]>其中ω0是发射信号的中心频率，c是声速，dx和dy是d在x和y方向的分量。
(9)由拟合结果得到载体对底的速度均值{V1x，V1y}，并将结果存储。
在一次测量结束后可重复步骤(1)～(9)，进行下一次载体对底速度的测量，在重复步骤(7)时，将上一次或多次测量得到的载体对底速度均值的平均作为未知参数系综 搜索范围的初始值。
本发明还提供一种相关测速声纳系统，包括声纳阵200和电子分机，所述声纳阵200包括发射声纳阵与接收声纳阵，电子分机内包括计算机406，所述计算机406内包括如下模块初始化模块，用于软件和硬件的初始化；发射编码模块，用于选择发射编码，该发射编码在相关延时非零时具有峰值；发射/接收模块，用于根据发射编码向流体介质内发射声脉冲，并接收各流层或底介质反向散射的回波信号；判断模块，用于判断发射/接收模块中接收到的回波信号是否包含底介质回波信号；截取模块，用于从回波信号中截取出底介质回波信号；解调滤波模块，用于对底介质回波信号进行解调滤波；第一计算模块，用于根据解调滤波后的底介质回波信号计算数据底介质时空相关函数矩阵；第二计算模块，用于将数据底介质时空相关函数矩阵作绝对值运算获得数据底介质时空相关函数绝对值矩阵；第三计算模块，用于在数据底介质时空相关函数绝对值矩阵中择取出其最大值EMax并设定一门限值x，其中0＜x≤1；该模块还用于在所述数据底介质时空相关函数绝对值矩阵中，将数值小于xEMax的元素置零，数值大于或等于xEMax的元素保持不变，对所有元素进行此操作后获得局域化的底介质时空相关函数绝对值矩阵；参数模块，用于存储未知参数系综 的搜索范围，其中V1x、V1y分别是x、y方向载体对底的速度均值，γ是宽度因子；参数模块所存储的未知参数系综 搜索范围的初始值优选是上一次或多次测得的载体对底速度均值的平均；第四计算模块，用于在参数模块中的未知参数系综 的搜索范围内，基于非线性最小二乘原理，采用序列二次规划方法，将第三计算模块中计算得到的局域化的底介质时空相关函数绝对值矩阵与理论的底介质时空相关函数绝对值矩阵进行拟合运算；所述理论底介质时空相关函数绝对值为 其中B是常数，θ是声波入射角，τ是延时，d是相邻声纳阵基元间的间距，J0是零阶贝塞尔函数；β2=ω0c((τV‾1x+dx)2+(τV‾1y+dy)2)1/2,]]>其中ω0是发射信号的中心频率，c是声速，dx和dy是d在x和y方向的分量；速度模块，用于存储由第四计算模块的拟合结果得到的载体对底的速度均值{Vx，Vy}。
在本发明的相关测速声纳系统中，由发射换能器组成的发射声纳阵与由接收换能器组成的接收声纳阵分离设置。所述接收声纳阵中的接收换能器的排列使换能器间距矢量差的数最大。所述发射声纳阵中的发射换能器紧密排列。在一个实施例中，所述发射声纳阵由7个发射换能器组成，所述接收声纳阵由8个接收换能器组成。
本发明的信号处理方法，基于非线性最小二乘原理，只处理底介质时空相关函数的绝对值矩阵和此绝对值矩阵幅值较大的部分，因此称它为信号局域化非线性最小二乘方法。
本发明的优点在于(1)在测量载体对底的速度时，本发明基于非线性最小二乘原理估计速度，与最大似然原理相比，它的鲁棒性好，计算量小，特别是在相关测速声纳的实际工作情况下，环境噪声在空间上往往是不均匀的，声纳阵基元的幅度和相位存在不一致性，它们对最小二乘原理的不利影响要明显小于对最大似然原理的影响。
(2)本发明的信号处理对象是底介质时空相关函数的绝对值，这使得公式(1)中的B是常数。相关函数绝对值只与载体的平均水平运动速度Vx和Vy有关，这样获得的速度值精度高，计算量小。对于绝大数实际使用情况，知道Vx和Vy已是足够了。
(3)本发明采用局域化的方法，用数据时空相关函数绝对值矩阵中幅值大的区域计算速度，去除了信噪比低的区域，提高了测量精度。
(4)本发明提出了新型声纳阵，给出了基元的分布，发射基元与接收基元分置，有利于合理选择发射波束宽度和接收波束宽度。
图1是相关测速声纳系统工作示意图；图2是相关测速声纳系统的装置示意图；图3是相关测速声纳系统中声纳阵的基元分布图；图4是相关测速声纳系统的软件流程图；图5是图4的流程图中步骤619的详细流程图；图6是本发明的相关测速声纳系统(ACL)与差分GPS分别测得的载体速度对比，工作水域水深3500米左右；其中，图6a是两设备分别测量的速度幅值701和703；图6b是两设备分别测量的速度方向702和704。
图面说明载体100声纳阵200水下电子分机300干端400终端500具体实施方式
下面结合附图与具体实施方式
对本发明作进一步详细描述。
如图1所示，使用相关测速声纳系统测量载体对底的速度时，该系统装载于载体100上，相关测速声纳系统总体上包括声纳阵200和电子分机。电子分机包括水下电子分机300、干端400和终端500。其中，声纳阵200和水下电子分机300安装在水下，干端400和终端500安装在水面上。声纳阵200中的发射声纳阵向水中发射声脉冲，一个脉宽的声脉冲102在水中传播，与海底相交在环103上，产生海底回波。海底回波被声纳阵200中的接收声纳阵接收，再经电子分机处理后，得载体的速度。
相关测速声纳系统的详细结构如图2所示，声纳阵200中包括接收换能器201和发射换能器202，由接收换能器201组成接收声纳阵，同时由发射换能器202组成发射声纳阵。声纳阵200的基元分布见图3，其中，1到7基元是发射换能器，它们紧密地紧靠排列；8到15基元是接收换能器，它们的排列原则是使得换能器间距矢量差的数最大。
再次参见图2，水下电子分机300包括多个与接收换能器201相连接的前置放大器301，水下电子分机300内还包括温度传感器302、漏水传感器303和姿态传感器304，分别与干端400中的声纳接口控制板407连接。
干端400中的发射机401与发射换能器202连接，多路接收机402与前置放大器301连接，多路同步数据采集板403与多路接收机402连接，DSP处理机404与多路同步数据采集板403连接；干端400还包括计算机406，并通过数据/控制总线405分别与DSP处理机404和多路同步数据采集板403连接；干端400还包括声纳接口控制板407，分别与多路接收机402、发射机401、DSP处理机404和计算机406连接；AC/DC电源408分别与声纳接口控制板407、多路接收机402、发射机401、数据/控制总线405、温度传感器302、漏水传感器303和姿态传感器304连接；干端400还包括与计算机406连接的GPS409和电罗经410。
终端500包括通过网络501与计算机406连接的远程计算机502。
该系统的专用测速程序装载在计算机406的存储器中，该程序包括初始化模块、发射编码模块、发射/接收模块、判断模块、截取模块、解调滤波模块、第一计算模块、第二计算模块、第三计算模块、参数模块、第四计算模块和速度模块，并按照图4所示流程图中的步骤执行。
步骤601是开始步骤，由终端计算机502发出指令，经网络501传输到计算机406，启动计算机406存储器中的程序，使声纳系统处于工作状态。在步骤602和603中，程序中的初始化模块对软件和系统硬件进行初始化。在步骤613中，发射编码模块根据底深度选择发射编码，该发射编码在相关延时非零时具有峰值。在步骤614中，发射/接收模块将计算机406的指令经数据/控制总线405发送至DSP处理机404，使得DSP处理机404产生发射信号驱动发射机401，驱动发射换能器202，发出声脉冲到流体介质中。在步骤615中，发射/接收模块控制接收换能器201接收从流体介质和海底介质的反向散射的回波信号，经前置放大器301馈送给多路接收机402，再经多路同步数据采集板403馈送到DSP处理机404。在步骤616中，判断模块控制DSP处理机404根据回波信号判断其中是否包含底介质回波信号，如果“否”，则返回步骤613；如果“是”，则进行步骤617。在步骤617中，截取模块控制DSP处理机404从回波信号中截取出底介质回波信号。在步骤618中，解调滤波模块控制DSP处理机404对底介质回波信号进行解调滤波。
在步骤619中，由第一、二和三计算模块控制DSP处理机404通过底介质回波信号计算局域化的数据底介质时空相关函数绝对值矩阵，其具体流程如图5所示。在图5中，第一计算模块首先根据解调滤波后的底介质回波信号计算数据底介质时空相关函数矩阵，然后第二计算模块将数据底介质时空相关函数矩阵作绝对值运算获得数据底介质时空相关函数绝对值矩阵，最后第三计算模块将数据底介质时空相关函数绝对值矩阵进行局域化操作。其中，这里的局域化是指在数据底介质时空相关函数绝对值矩阵中择取出其最大值EMax，并设定一门限值x，其中0＜x≤1，然后将绝对值矩阵中数值小于xEMax的元素置零，数值大于或等于xEMax的元素保持不变，对所有元素进行此操作后即可获得局域化的数据底介质时空相关函数绝对值矩阵。在这里，进行的局域化操作只取大于或等于xEMax的元素，即只取信噪比大的区域，舍掉信噪比小的区域，简化了计算，提高了测量精度。在实际应用时，门限x优选在0.7到1之间。
在申请号为“03119665.9”的中国专利申请中提供了一种如公式(1)描述的底介质声纳阵时空相关函数。在本发明中，理论的底介质声纳阵时空相关函数形式为 其中B是f(Vz)的函数，f为某函数，Vz是z方向载体对底速度的均值，d是相邻声纳阵基元间的间距，τ是延时，θ是声波入射角，J0是零阶贝塞尔函数；β2=ω0c((τV‾1x+dx)2+(τV‾1y+dy)2)1/2,]]>其中ω0是发射信号的中心频率，c是声速，dx和dy是d在x和y方向的分量。
由公式(2)可见， 与Vx，Vy和Vz有关，如果同时估计此三个速度，则计算复杂，精度也不会高。将公式(2)所描述的理论的底介质声纳阵时空相关函数取绝对值后，得到
此时，B是常数，将公式(3)描述的理论底介质时空相关函数绝对值构成矩阵，该矩阵称为理论的底介质时空相关函数绝对值矩阵，它只与Vx和Vy有关。这样计算简化，精度提高。在实际情况下，往往只需要Vx和Vy已是足够了，况且Vz可由其它设备测量。公式(1)与公式(3)虽然形式相同，但获得的方法不同。公式(1)是人为设B是常数，而公式(3)是由正确的时空相关函数公式(2)取绝对值得到，物理意义明确，显然比公式(1)好。
在步骤620中，由参数模块设置并存储公式(3)中未知系综 的非线性最小二乘法的搜索范围，其中，在首次测量时，可以将未知系综 的搜索范围设置的尽可能的大，从而将真实的载体对底速度值包含在该范围内。在以后的测量中优选采用前一次或多次的测量结果的均值作为搜索范围的初始值，这样可提高运算速度、减少硬件开销。
在步骤621中，第四计算模块控制DSP处理机404将步骤619中得到的局域化数据底介质时空相关函数绝对值矩阵与本发明理论的底介质时空相关函数绝对值矩阵进行拟合运算，得到载体对底的速度均值。该拟合算法基于非线性最小二乘原理，可以采用单纯形方法或序列二次规划方法，优选为序列二次规划方法。
在步骤622中，速度模块将步骤621得到的拟合结果经数据/控制总线405馈送给计算机406，并存储在存储器中。步骤622之后，可返回步骤613进行下一次测量。
最后还可以将温度传感器302、漏水传感器303和姿态传感器304的数据经声纳接口控制器407馈送给计算机406。计算机406再融合GPS409和电罗经410的数据，最终结果经网络501馈送给终端计算机502。
图6是使用本发明的相关测速声纳系统在中心频率23.5kHz，带宽4.4kHz时与差分GPS分别测得载体100对底的速度图。其中，图6a表示的是在不同时间测得的载体对底的速度，横坐标是时间，纵坐标是速度；图6b表示的是在不同时间测得的载体对底速度的方向，横坐标是时间，纵坐标是角度。曲线701和703是本发明的相关测速声纳系统测得的数据，曲线702和704是差分GPS测得的数据，两者相当符合。
权利要求
1.一种相关测速声纳测量载体对底速度的方法，包括以下步骤(1)选择发射编码；(2)根据发射编码发射声脉冲到流体介质中，并接收各流层或底介质反向散射的回波信号；(3)如步骤(2)中的回波信号包含底介质回波信号，则进行步骤(4)；如未包含，则返回步骤(1)重新开始；(4)对所述底介质回波信号进行解调滤波，并根据底介质回波信号计算数据底介质时空相关函数矩阵；(5)将数据底介质时空相关函数矩阵中的所有元素作绝对值运算，获得数据底介质时空相关函数绝对值矩阵，该矩阵中元素的最大值为EMax；(6)设定一门限值x，其中0＜x≤1；在所述数据底介质时空相关函数绝对值矩阵中，将数值小于xEMax的元素置零，数值大于或等于xEMax的元素保持不变，对所有元素进行此操作后获得局域化的底介质时空相关函数绝对值矩阵；(7)设置未知参数系综 的搜索范围，其中V1x、V1y分别是x、y方向载体对底的速度均值，γ是宽度因子；(8)在未知参数系综 的搜索范围内，基于非线性最小二乘原理，将步骤(6)中局域化的底介质时空相关函数绝对值矩阵与理论的底介质时空相关函数的绝对值矩阵进行拟合运算；所述理论的底介质声纳阵时空相关函数绝对值为 其中B是常数，θ是声波入射角，τ是延时，d是相邻声纳阵基元间的间距，J0是零阶贝塞尔函数；β2=ω0c((τV‾1x+dx)2+(τV‾1y+dy)2)1/2,]]>其中ω0是发射信号的中心频率，c是声速，dx和dy是d在x和y方向的分量；(9)由拟合结果得到载体对底的速度均值{V1x，V1y}，并将结果存储。
2.根据权利要求1所述的相关测速声纳测量载体对底速度的方法，其特征在于，所述门限值0.7＜x≤1。
3.根据权利要求1所述的相关测速声纳测量载体对底速度的方法，其特征在于，重复步骤(1)～(9)，进行下一次载体对底速度的测量。
4.根据权利要求1所述的相关测速声纳测量载体对底速度的方法，其特征在于，步骤(8)中的拟合算法为序列二次规划方法。
5.根据权利要求3所述的相关测速声纳测量载体对底速度的方法，其特征在于，在重复步骤(7)时，将上一次或多次测量得到的载体对底速度均值的平均作为未知参数系综 搜索范围的初始值。
6.根据权利要求1所述的相关测速声纳测量载体对底速度的方法，其特征在于，步骤(1)中的发射编码在相关延时非零时具有峰值。
7.一种实施权利要求1所述方法的相关测速声纳系统，包括声纳阵(200)和电子分机，所述声纳阵(200)包括发射声纳阵与接收声纳阵，电子分机内包括计算机(406)，其特征在于，所述计算机(406)内包括如下模块初始化模块，用于软件和硬件的初始化；发射编码模块，用于选择发射编码；发射/接收模块，用于根据发射编码向流体介质内发射声脉冲，并接收各流层或底介质反向散射的回波信号；判断模块，用于判断发射/接收模块中接收到的回波信号是否包含底介质回波信号；截取模块，用于从回波信号中截取出底介质回波信号；解调滤波模块，用于对底介质回波信号进行解调滤波；第一计算模块，用于根据解调滤波后的底介质回波信号计算数据底介质时空相关函数矩阵；第二计算模块，用于将数据底介质时空相关函数矩阵作绝对值运算获得数据底介质时空相关函数绝对值矩阵；第三计算模块，用于在数据底介质时空相关函数绝对值矩阵中择取出其最大值EMax并设定一门限值x，其中0＜x≤1；该模块还用于在所述数据底介质时空相关函数绝对值矩阵中，将数值小于xEMax的元素置零，数值大于或等于xEMax的元素保持不变，对所有元素进行此操作后获得局域化的底介质时空相关函数绝对值矩阵；参数模块，用于存储未知参数系综 的搜索范围，其中V1x、V1y分别是x、y方向载体对底的速度均值，γ是宽度因子；第四计算模块，用于在参数模块中的未知参数系综 的搜索范围内，基于非线性最小二乘原理，将第三计算模块中计算得到的局域化的底介质时空相关函数绝对值矩阵与理论的底介质时空相关函数绝对值矩阵进行拟合运算；所述理论底介质时空相关函数绝对值为 其中B是常数，θ是声波入射角，τ是延时，d是相邻声纳阵基元间的间距，J0是零阶贝塞尔函数；β2=ω0c((τV‾1x+dx)2+(τV‾1y+dy)2)1/2,]]>其中ω0是发射信号的中心频率，c是声速，dx和dy是d在x和y方向的分量；速度模块，用于存储由第四计算模块的拟合结果得到的载体对底的速度均值{Vx，Vy}。
8.根据权利要求7所述的相关测速声纳系统，其特征在于，第四计算模块是采用序列二次规划方法进行拟合运算的计算模块。
9.根据权利要求7所述的相关测速声纳系统，其特征在于，参数模块所存储的未知参数系综 搜索范围的初始值是上一次或多次测得的载体对底速度均值的平均。
10.根据权利要求7所述的相关测速声纳系统，其特征在于，发射编码模块产生的发射编码在相关延时非零时具有峰值。
11.根据权利要求7所述的相关测速系统，其特征在于，由发射换能器组成的发射声纳阵与由接收换能器组成的接收声纳阵分离设置。
12.根据权利要求11所述的相关测速系统，其特征在于，所述接收声纳阵中的接收换能器的排列使换能器间距矢量差的数最大。
13.根据权利要求11或12所述的相关测速系统，其特征在于，所述发射声纳阵中的发射换能器紧密排列。
14.根据权利要求13所述的相关测速系统，其特征在于，所述发射声纳阵由7个发射换能器组成，所述接收声纳阵由8个接收换能器组成。
全文摘要
本发明公开了一种相关测速声纳测量载体对底速度的方法和系统。本发明采用局域化的方法，用数据底介质时空相关函数绝对值矩阵中幅值大的区域来计算速度，去除了信噪比低的区域。本发明基于非线性最小二乘原理，采用序列二次规划方法，将局域化的数据底介质时空相关函数绝对值矩阵与理论的底介质时空相关函数绝对值矩阵进行拟合运算来估计速度。采用本发明的方法测量载体对底速度时，鲁棒性好，计算量小，测量精度高。在本发明相关测速声纳中，发射基元与接收基元分置，有利于合理选择发射波束宽度和接收波束宽度。
文档编号G01S7/523GK1621857SQ20031011515
公开日2005年6月1日 申请日期2003年11月24日 优先权日2003年11月24日
发明者朱维庆, 王长红 申请人:中国科学院声学研究所