本发明专利技术的目的在于提供一种基于可变带宽滤波器的多波束测深数据处理方法，包括以下步骤：向水下发射信号，对海底反向散射信号进行带通采样，经过正交变换和低通滤波后，得到通道数为M的多波束复信号，根据深度值计算波束形成角度的各波束的带宽；设计一组阶数相同、截止频率不同的滤波器，作为可变带宽滤波器的子滤波器；选取滤波器的截止频率，产生一组滤波器选择系数；采用以上的滤波器系数和选择系数，对各通道信号进行滤波处理，之后对信号进行波束形成，并通过测深算法得到深度值用于下一个周期的深度初始值。本发明专利技术既能够弥补固定带宽滤波器原理上的不足，又具有工程可行性的滤波器可以进一步提高多波束测深系统的精度。
【技术实现步骤摘要】
本专利技术涉及的是一种水下的测深方法。
技术介绍
在多波束测深系统中，海底回波的信噪比直接影响多波束测深结果的精度。为了提高多波束测深精度必须提高信噪比，对信号进行滤波处理是一种有效手段。目前国内外的多波束测深系统广泛采用固定带宽滤波处理，其方法是根据发射脉冲的长度选取滤波器的带宽，其不足之处是由于带宽依赖于波束宽度、海底深度、波束和垂直方向夹角、脉冲宽度等很多因素，来自条带上不同波束方向的海底回波信噪比并不相同，因此固定带宽滤波器的效果并不是最优的。2004 ^ OCEANS ^^ Variable Bandwidth Filter for Multibeam Echo-sounding Bottom Detection 一文禾口 2010 年 ICIA 会议 Estimation of Multibeam Phase Difference Using Variable Bandwidth Filter—文提出了采用可变带宽滤波器， 其方法是在波束形成之后，针对不同波束具有不同的带宽，采用可变带宽滤波器对相位差序列进行自适应滤波处理，其不足之处是由于滤波在波束形成之后进行，计算量很大，难以实时实现并应用于工程实际。
技术实现思路
本专利技术的目的在于提供一种既能够弥补固定带宽滤波器原理上的不足，又具有工程可行性的滤波器可以进一步提高多波束测深系统的精度的。本专利技术的目的是这样实现的本专利技术，其特征是(1)向水下发射信号，并对海底反向散射信号进行带通采样，经过正交变换和低通滤波后，得到通道数为M的多波束复信号，在下述信号处理步骤中将对信号进行波束形成， 波束形成的角度为θ2， ，θ η，在测量的第一个周期设置默认深度初始值H，根据深度值H计算波束形成角度为θ” θ2， ，θ η的各波束的带宽；(2)由步骤(1)各波束的带宽设计一组阶数相同、截止频率不同的滤波器，作为可变带宽滤波器的子滤波器；(3)对于每一采样时刻回波到达的波束，根据其中的最大带宽选取滤波器的截止频率，产生一组滤波器选择系数，用来选择不同时刻采用的子滤波器； (4)采用以上的滤波器系数和选择系数，对各通道信号进行滤波处理，之后对信号进行波束形成，并通过测深算法得到深度值用于下一个周期的深度初始值。本专利技术还可以包括1、所述的波束角度θ2， ，θ η的选择满足θ2， ，θ η对波束形成的角度范围全面覆盖。2、所述的根据深度值H计算波束形成角度为θ ” θ 2， ，θ n的各波束的带宽的方法为c为声速，τ为脉冲宽度，λ为信号的波长，对于阵列长度为L的均勻直线阵，H为深度，β i为波束指向为θ i时的-3dB半波束宽度2 β i ^ 2sin_1 (0. 443 λ /L) sec θ j (i = 1,2, · , η),、为_3dB波束宽度内的回波持续时间tj = 2 (cos-1 ( θ j+ β j) -cos-1 ( θ j- β j)) H/c+ τ (i = 1,2, · , η),Qi方向波束输出信号的-3dB带宽Afi为Afi = 1/、 (i = l，2，*，n)。3、所述的设计一组阶数相同、截止频率不同的滤波器的方法为采用HR低通滤波器，子滤波器的截止频率f。i根据各波束的带宽选取4 = kX Af^k为常数，k的范围为 1到2之间。4、所述的产生一组滤波器选择系数的方法为当前周期的深度近似等于前一周期得到的深度值，对波束角度为Qi得到其-3dB波束宽度，该波束的回波角度范围为 ；回波最先到达的时间为、=2^0^( Qi-β J/c，最后到达的时间为t2 =2HCOS-1 ( θ # i)/c+ τ，采样频率为fs，则回波持续的采样点为H1 = Vfs到n2 = t2/fs ； 在同一采样时刻的波束有效回波中大的带宽选择滤波器的截止频率，按照时间顺序保存采用的滤波器的编号。本专利技术的优势在于考虑到各波束的回波信号来自于不同的角度，回波的到达时间有先后顺序，比如正下方最先收到回波，而外侧波束方向较晚收到回波，因此在波束形成之前，在时间轴上对不同波束方向的回波到达时间进行划分，对某个时间段采用这个波束方向回波的带宽作为滤波器的截止频率、对另一个时间段采用另一个波束方向回波的带宽作为滤波器的截止频率，…，以此类推，这样既能够达到可变带宽滤波的效果，又可以极大程度的减小计算量。附图说明图1为本专利技术的不同波束角度回波信号的带宽；图2不同时刻滤波器的截止频率；图3本专利技术流程图；图4信号的频谱；图5相位差曲线对比。具体实施例方式下面结合附图举例对本专利技术做更详细地描述结合图1 5，为实现本专利技术所述目的，通过以下技术方案实现1.多波束测深系统向水下发射信号，并对海底反向散射信号进行带通采样，经过正交变换和低通滤波后，得到通道数为M的多波束复信号，在后续信号处理中将对信号进行波束形成，波束形成的角度为θ2， ，θ η，在测量的第一个周期设置默认深度初始值H，其后的每一个周期，H由前一周期测量得到，根据深度值H计算波束形成角度为Θ” θ2， ，θ η的各波束的带宽；2.由各波束的带宽设计一组阶数相同，截止频率不同的滤波器，作为可变带宽滤波器的子滤波器；3.对于每一采样时刻，有若干个波束的回波到达，根据其中的最大带宽选取滤波器的截止频率，产生一组滤波器选择系数，用来选择不同时刻采用的子滤波器；4.采用以上的滤波器系数和选择系数，对各通道信号进行滤波处理，之后对信号进行波束形成，并通过测深算法得到深度值H用于下一个周期的处理。本专利技术还可以包括(1)步骤1中的各波束角度θ” θ2， ，θ 是根据波束形成的角度范围和进行滤波的精度选择，要满足θ2， ，θ η对波束形成的角度范围全面覆盖。(2)步骤1中根据前一周期测得的深度值，求当前周期各波束θ ” θ 2， ，θ η的带宽方法如下设c为声速，τ为脉冲宽度，λ为信号的波长，对于阵列长度为L的均勻直线阵， H为深度，β i为波束指向为θ i时的-3dB半波束宽度2 β i ^ 2sin_1 (0. 443 λ /L) sec θ j (i = 1,2, · , η)、为_3dB波束宽度内的回波持续时间tj = 2 (cos-1 ( θ j+ β j) -cos-1 ( θ j- β j)) H/c+ τ (i = 1,2, · , η)Qi方向波束输出信号的-3dB带宽Afi为Afi = 1/、 (i = l，2，*，n)。(3)步骤2中滤波器的设计方法是采用HR低通滤波器，子滤波器的截止频率f。i 根据各波束的带宽选取f。i = kX Afi，其中k根据环境情况选取。(4)步骤3中滤波器选择系数按如下步骤计算(a)假设深度为渐变的，当前周期的深度近似等于前一周期得到的深度值为H，对波束角度为Qi，由权利要求2中的计算方法可以得到其-3dB波束宽度β”该波束的回波角度范围为[θ「β" θ +β」；(b)回波最先到达的时间为、=？！！⑶^⑶厂^八，最后到达的时间为、= 2^08^(0 i+β i)/c+ τ，采样频率为fs，则回波持续的采样点为H1 = Vfs到n2 = t2/fs ；(c)在同一采样时刻可能有若干个波束都有有效回波，根据其中较大的带宽选择滤波器的截止频率，按照时间顺序保存采用的滤波器的编本文档来自技高网...
【技术保护点】
１．一种基于可变带宽滤波器的多波束测深数据处理方法，其特征是：（１）向水下发射信号，并对海底反向散射信号进行带通采样，经过正交变换和低通滤波后，得到通道数为Ｍ的多波束复信号，在下述信号处理步骤中将对信号进行波束形成，波束形成的角度为θ１，θ２，·，θｎ，在测量的第一个周期设置默认深度初始值Ｈ，根据深度值Ｈ计算波束形成角度为θ１，θ２，·，θｎ的各波束的带宽；（２）由步骤（１）各波束的带宽设计一组阶数相同、截止频率不同的滤波器，作为可变带宽滤波器的子滤波器；（３）对于每一采样时刻回波到达的波束，根据其中的最大带宽选取滤波器的截止频率，产生一组滤波器选择系数，用来选择不同时刻采用的子滤波器；（４）采用以上的滤波器系数和选择系数，对各通道信号进行滤波处理，之后对信号进行波束形成，并通过测深算法得到深度值用于下一个周期的深度初始值。
【技术特征摘要】
【专利技术属性】
技术研发人员：李海森，魏玉阔，周天，李珊，陈宝伟，李若，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：93
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。