CN109407052B - 一种线性阵换能器通道一致性校准方法 - Google Patents
一种线性阵换能器通道一致性校准方法 Download PDFInfo
- Publication number
- CN109407052B CN109407052B CN201811546316.5A CN201811546316A CN109407052B CN 109407052 B CN109407052 B CN 109407052B CN 201811546316 A CN201811546316 A CN 201811546316A CN 109407052 B CN109407052 B CN 109407052B
- Authority
- CN
- China
- Prior art keywords
- phase
- channel
- sound source
- distance
- array
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种线性阵换能器通道一致性校准法。该方法在近场环境下,将多波束接收阵的中心与声源的中心设在同一吃水深度,通过空间扫描的方式获取最优位置,利用最优估计的算法,找到幅度和相位的误差补偿曲线,降低了试验复杂度,也降低了算法处理的复杂度,通用性比较好。得到的幅度和相位补偿曲线补偿到波束指向性上,从而提高主旁瓣增益比。
Description
技术领域
本发明涉及水下定位领域,尤其涉及一种线性阵换能器通道一致性校准方法。
背景技术
多波束接收机的性能对测深结果非常重要,而接收换能器作为接收机的重要组成部分,它影响到信号的输入和采集,影响到后期信号处理的整个过程。且接收换能器是把声压信号转换到电信号,连接模拟电路,然后才进入AD采样前端。为了保证采集信号的重要性,首先必须保证采集信号的正确性,反应到换能器上主要是幅度的一致性和相位的一致性上,通过对换能器的通道一致性的校准,可以把换能器的加工误差和模拟电路的误差同时进行校准,可以把这些误差降到最低。
国内外一般是远场条件下进行换能器的通道校准,这样就对试验条件和试验环境要求比较高,消声水池要求比较大,并且配合NI的控制系统,由于消声水池造价高,国内一般单位很难拥有这样的测试环境。且有的在校准方法中,在测量过程中,对声源和接收阵的吃水深度没有做具体的要求,造成声源和接收阵按照三维空间的几何模型模型,造成算法模型的复杂度。
发明内容
本发明的目的旨在解决上述问题,从而提供一种线性阵换能器通道一致性校准方法。
为实现上述目的,本发明提供了一种线性阵换能器通道一致性校准方法,该方法包括以下步骤:
1)在消声水池中安装声源和多波束接收阵,且保证在同一深度的吃水位置,声源的旋转中心和多波束的旋转中心在同一水平位置,所述声源安装在声源法兰上,所述多波束接收阵安装在多波束法兰上;
2)通过驱动单元驱动声源工作,声源法兰固定不动,通过控制台控制多波束法兰按照[-α:l:α]旋转,其中α为旋转角度,l为步长;且每一个旋转角度位置都利用与接收阵通道数目相同的采集板保存数据,共要保存k=(2α/l+1)组数据,k为保存数据的组数。
3)在声源发射方向正对接收阵法线方向时(如图1所示位置),测量声源到多波束接收阵中心的垂直距离Ver,单位为m,并根据测量的垂直距离Ver计算出声源到多波束接收阵的各个阵元的相位值,并根据各个相位值及对应的阵元获取相位曲线j。
4)利用步骤2)保存的k组数据,根据FFT求相位的方式求出每组数据的接收到的相位值,并根据相位值及对应通道获取相位曲线,从求出的k个位置的所有通道的相位曲线中找出与步骤3)中几何位置计算的相位曲线j重合度最高的一条c;
5)根据步骤4)得出的相位曲线c对应的相位值数据建立近场的物理模型,根据物理模型建立关于角度和距离的目标函数,然后对目标函数求偏导数,获取与实际相符合的最优估计的距离和角度所述获取与实际相符合的最优估计的角度与步骤4)中得出相位曲线c对应的多波束接收阵旋转的角度接近,所述获取与实际相符合的最优估计的距离与步骤3)中距离Ver接近,根据与实际相符合的最优估计的距离和角度得出所有通道相位值,并根据各个相位值及对应的通道获取相位曲线g;
6)将步骤5)估计求出的各通道相位值与和步骤4)中相位曲线c对应的各通道相位值相减,得到各通道相位误差,并根据相位误差和对应通道获取相位补偿曲线x;
7)将步骤4)中相位曲线c对应的采集数据,通过连续多ping的统计分析,让第一通道作为参考,通过希尔伯特方法求取所有通道的幅度,然后与第一通道相除后取倒数,得到所有通道的幅度补偿曲线f。
进一步地,所述步骤2)中驱动单元包括信号源和功率放大器,所述信号源和功率放大器相连,所述功率放大器与所述声源相连。
进一步地,所述步骤5)具体包括:
其中,ψi为测量数据FFT计算的每个通道的相位值,M为阵元总数目,f为发射信号的频率,c为水中的声速为1500m/s,ri为声源到阵元的距离,r为声源到多波束接收阵中心阵元的距离,xi为阵元到中心阵元的距离,θ为多波束接收阵法线方向和声源中心到阵元的距离之间的夹角;
更进一步地,所述步骤6)包括:
进一步地,所述多波束法兰按照[-4:0.1:4]旋转。
进一步地,所述多波束接收阵包括64个阵元,所述采集板的通道数为64。
本方法在近场环境下,仅仅需要保证接收阵的中心与声源的中心在同一个吃水深度,通过空间扫描的方式获取最优位置,利用最优估计的算法,找到幅度和相位的误差补偿曲线,降低了试验复杂度,也降低了算法处理的复杂度,通用性比较好。得到的幅度和相位补偿曲线补偿到波束指向性上,从而提高主旁瓣增益比。
附图说明
图1为本发明实施例提供的一种线性阵换能器通道一致性校准方法的流程图;
图2为本发明实施例提供的声源发射方向正对接收阵法线方向时几何模型图;
图3为本发明实施例提供的近场模型下阵元之间的几何模型图;
图4为本发明实施例提供的几何位置计算、FFT采样数据计算、最优估计三种方式得到的64通道的相位曲线图;
图5为本发明实施例提供的得到的相位补偿曲线图;
图6为本发明实施例提供的得到的幅度补偿曲线图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。需要说明的是,附图仅为示例性说明,并未按照严格比例绘制,而且其中可能有为描述便利而进行的局部放大、缩小,对于公知部分结构亦可能有一定缺省。
图1为本发明提供的一种线性阵换能器通道一致性校准方法的流程图。
在步骤101中,在消声水池中安装声源和多波束接收阵,且保证在同一深度的吃水位置,声源的旋转中心和多波束的旋转中心在同一水平位置,声源安装在声源法兰上,多波束接收阵安装在多波束法兰上。
在步骤102中,通过驱动单元驱动声源工作,声源法兰固定不动,通过控制台控制多波束法兰按照[-α:l:α]旋转,其中α为旋转角度,l为步长;且每一个旋转角度位置都利用与接收阵通道数目相同的采集板保存数据,共要保存k=(2α/l+1)组数据,k为保存数据的组数。
本实施例中驱动单元包括信号源和功率放大器,通过信号源连接功率放大器,然后放大器端接声源,给信号源端加换能器频率的正弦信号,然后通过功率放大器来驱动声源。本实施例中的控制台步长为0.1度,假设按照[-4:0.1:4]度(其中4度是待估计的角度量),且每一个旋转角度位置都利用与接收阵通道数目相同的高速采集板保存数据(本实施例换能器接收阵为64通道,采集卡的通道数为64),也就是多波束接收阵每旋转0.1度,便进行64通道的数据采集,一共要保存81组数据。
在步骤103中,在声源S发射方向正对接收阵法线方向时(如图2所示的声源与多波束接收阵的位置),测量声源S到多波束接收阵中心的垂直距离Ver,单位为m,并根据测量的垂直距离Ver带入公式计算出声源到多波束接收阵的各个阵元xi的相位值,并根据各个相位值及对应的阵元获取相位曲线j。其中f为发射信号的频率,d为多波束接收阵的阵元之间的间距,C为水中的声速为1500m/s,m为多波束接收阵的阵元序号。
在步骤104中,利用步骤102保存的k组数据,根据FFT求相位的方式求出每组数据的接收到的相位值,并根据相位值及对应通道获取相位曲线,从求出的k个位置的所有通道的相位曲线中找出与步骤3)中几何位置计算的相位曲线j重合度最高的一条c。
如按照[-4:0.1:4]度旋转,则k=81,根据FFT求相位的方式求出81组数据的对应的64通道相位值,会得到81条64通道的相位曲线。算出的81个位置的所有通道的相位曲线,由于是按照-4到4度旋转方式扫描覆盖的,所以肯定存在一组数据的相位曲线和与步骤103中几何位置计算的相位曲线的重合度最高,则认为此位置是,声源正对着多波束接收阵的法线上发射,充分利用该组数据进行幅度和相位一致性的分析。
其中,ψi为测量数据FFT计算的每个通道的相位值,M为阵元总数目,f为发射信号的频率,c为水中的声速为1500m/s,ri为声源到阵元的距离,r为声源到多波束接收阵中心阵元的距离,xi为阵元到中心阵元的距离,θ为多波束接收阵法线方向和声源中心到阵元的距离之间的夹角,参考图3。
分别对目标函数Q的角度和距离进行求偏导,获取与实际相符合的最优估计的距离和角度(由于方程可能存在几组最优解,选取与实际相符合的一组最优解),获取与实际相符合的最优估计的角度与步骤104中得出相位曲线c对应的多波束接收阵旋转的角度接近,获取与实际相符合的最优估计的距离与步骤103中距离Ver接近。将最优估计的距离和角度代入公式求取各通道的相位值,并根据各个相位值及对应的通道获取相位曲线g。
在步骤106中,将步骤105估计求出的各通道相位值与和步骤104中相位曲线c对应的各通道相位值相减,得到各通道相位误差,并根据相位误差和对应通道获取相位补偿曲线x,相位补偿曲线x参见图5。
在步骤107中,将步骤104中相位曲线c对应的采集数据,通过连续多ping的统计分析,让第一通道作为参考,通过希尔伯特方法求取所有通道的幅度,然后与第一通道相除后取倒数,得到所有通道的幅度补偿曲线f,幅度补偿曲线x参见图6。
步骤106和107得到的幅度和相位补偿曲线可补偿到波束指向性上,从而提高主旁瓣增益比。步骤103几何位置相位曲线j、步骤104中FFT计算相位曲线c、步骤105最优估算相位曲线g可参见图4,如图4所示,最优估算相位曲线和几何位置相位曲线基本重合,即验证了本发明的最优估算法的正确性。
综上,本发明建立近场的物理模型,降低了对消声水池的要求,且只要保证吃水深度一样,降低了模型的复杂度。由于直接让声源正对换能器法线位置发射的位置不容易控制,通过几何位置,旋转多波束发射,通过与几何位置的重合性,找到正对法线的位置,方法简单明了,降低了测试过程的复杂度,且不需要太发杂的控制平台操作,只利用了一组数据,不用利用繁琐的谱估计方法。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种线性阵换能器通道一致性校准方法,其特征在于,包括:
1)在消声水池中安装声源和多波束接收阵,且保证在同一深度的吃水位置,声源的旋转中心和多波束的旋转中心在同一水平位置,所述声源安装在声源法兰上,所述多波束接收阵安装在多波束法兰上;
2)通过驱动单元驱动声源工作,声源法兰固定不动,通过控制台控制多波束法兰按照[-α:l:α]旋转,其中α为旋转角度,l为步长;且每一个旋转角度位置都利用与接收阵通道数目相同的采集板保存数据,共要保存k=(2α/l+1)组数据,k为保存数据的组数;
3)在声源发射方向正对接收阵法线方向时,测量声源到多波束接收阵中心的垂直距离Ver,单位为m,并根据测量的垂直距离Ver计算出声源到多波束接收阵的各个阵元的相位值,并根据各个相位值及对应的阵元获取相位曲线j;
4)利用步骤2)保存的k组数据,根据FFT求相位的方式求出每组数据的接收到的相位值,并根据相位值及对应通道获取相位曲线,从求出的k个位置的所有通道的相位曲线中找出与步骤3)中几何位置计算的相位曲线j重合度最高的一条c;
5)根据步骤4)得出的相位曲线c对应的相位值数据建立近场的物理模型,根据物理模型建立关于角度和距离的目标函数,然后对目标函数求偏导数,获取与实际相符合的最优估计的距离和角度所述获取与实际相符合的最优估计的角度与步骤4)中得出相位曲线c对应的多波束接收阵旋转的角度接近,所述获取与实际相符合的最优估计的距离与步骤3)中距离Ver接近,根据与实际相符合的最优估计的距离和角度得出所有通道相位值,并根据各个相位值及对应的通道获取相位曲线g;
6)将步骤5)估计求出的各通道相位值与和步骤4)中相位曲线c对应的各通道相位值相减,得到各通道相位误差,并根据相位误差和对应通道获取相位补偿曲线x;
7)将步骤4)中相位曲线c对应的采集数据,通过连续多ping的统计分析,让第一通道作为参考,通过希尔伯特方法求取所有通道的幅度,然后与第一通道相除后取倒数,得到所有通道的幅度补偿曲线f。
2.根据权利要求1所述的一种线性阵换能器通道一致性校准方法,其特征在于,所述步骤2)中驱动单元包括信号源和功率放大器,所述信号源和功率放大器相连,所述功率放大器与所述声源相连。
4.根据权利要求1所述的一种线性阵换能器通道一致性校准方法,其特征在于,所述步骤5)具体包括:
其中,ψi为测量数据FFT计算的每个通道的相位值,M为阵元总数目,f为发射信号的频率,c为水中的声速为1500m/s,ri为声源到阵元的距离,r为声源到多波束接收阵中心阵元的距离,xi为阵元到中心阵元的距离,θ为多波束接收阵法线方向和声源中心到阵元的距离之间的夹角;
6.根据权利要求1所述的一种线性阵换能器通道一致性校准方法,其特征在于,所述多波束法兰按照[-4:0.1:4]旋转。
7.根据权利要求1所述的一种线性阵换能器通道一致性校准方法,其特征在于,所述多波束接收阵包括64个阵元,所述采集板的通道数为64。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811546316.5A CN109407052B (zh) | 2018-12-18 | 2018-12-18 | 一种线性阵换能器通道一致性校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811546316.5A CN109407052B (zh) | 2018-12-18 | 2018-12-18 | 一种线性阵换能器通道一致性校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109407052A CN109407052A (zh) | 2019-03-01 |
CN109407052B true CN109407052B (zh) | 2022-09-13 |
Family
ID=65459836
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811546316.5A Active CN109407052B (zh) | 2018-12-18 | 2018-12-18 | 一种线性阵换能器通道一致性校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109407052B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110418093B (zh) * | 2019-07-30 | 2021-01-29 | 成都芯通软件有限公司 | 一种hfc放大器多路端口输出一致性校准方法 |
CN110703258B (zh) * | 2019-10-22 | 2022-09-09 | 哈尔滨工程大学 | 水下探测基阵通道间相位一致性校准方法 |
CN110703259B (zh) * | 2019-10-22 | 2021-06-18 | 哈尔滨工程大学 | 基于运动声源的水下声学基阵通道间相位一致性校准方法 |
CN111964757B (zh) * | 2020-08-13 | 2022-11-11 | 浙江苍南仪表集团股份有限公司 | 一种基于回波信号特征参数的换能器一致性评价方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644765B (zh) * | 2009-04-23 | 2012-03-14 | 中国科学院声学研究所 | 一种用于水声换能器线阵的幅相误差校正方法 |
CN102890271B (zh) * | 2012-10-25 | 2013-11-27 | 北京理工大学 | 一种外辐射源雷达阵列天线幅相一致性校正方法 |
CN103983946A (zh) * | 2014-05-23 | 2014-08-13 | 北京神州普惠科技股份有限公司 | 一种声源定位过程中的多测量通道信号处理方法 |
CN105649606A (zh) * | 2014-08-07 | 2016-06-08 | 中国石油大学(北京) | 多探头声波幅度测井仪探头一致性的补偿式标定法 |
CN105629216B (zh) * | 2016-03-09 | 2018-03-02 | 无锡市海鹰加科海洋技术有限责任公司 | 一种多波束测深仪接收换能器性能检测装置 |
CN105973447A (zh) * | 2016-05-05 | 2016-09-28 | 哈尔滨工程大学 | 直线水听器阵列幅度和相位一致性测量方法 |
-
2018
- 2018-12-18 CN CN201811546316.5A patent/CN109407052B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109407052A (zh) | 2019-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109407052B (zh) | 一种线性阵换能器通道一致性校准方法 | |
CN108037496B (zh) | 一种自由场水听器复数灵敏度精确测量方法 | |
CN107329125B (zh) | 消除短时突发干扰信号的自闭环校准方法 | |
CN109283511B (zh) | 一种宽覆盖多波束接收基阵校准方法 | |
CN104345303A (zh) | 车辆用雷达校准系统 | |
CN108469615A (zh) | 一种基于蝙蝠仿生声纳的小目标定位装置及其构建方法与应用 | |
CN108535682B (zh) | 一种基于旋转非均匀双l阵的水下二维doa估计方法及装置 | |
CN206161743U (zh) | 一种数字阵列天线的通道特性测试及配相校准系统 | |
CN112014812B (zh) | 一种相控阵风廓线雷达标定系统及方法 | |
US10470009B2 (en) | Test device and test method | |
CN110361720B (zh) | 一种声呐收声板控制方法 | |
CN109581274A (zh) | 基于夹角可调三维阵的非圆信号水下doa估计方法和装置 | |
CN110515082B (zh) | 一种基于超声波的自动测距系统 | |
CN115219804A (zh) | 一种太赫兹天线和天线罩电性能近场测试装置及方法 | |
JP2009236707A (ja) | 信号源位置推定方法 | |
CN103983946A (zh) | 一种声源定位过程中的多测量通道信号处理方法 | |
CN109752705B (zh) | 高频水声阵列性能参数测量方法及系统、设备及存储介质 | |
CN109839543B (zh) | 一种天线的幅相一致性测试系统及测试方法 | |
CN113093107A (zh) | 水声信号采集传输系统和方法 | |
CN107991657B (zh) | 一种用于双波束天馈的波束对准系统 | |
WO2019061173A1 (zh) | 一种校准阵列天线的方法、装置和系统 | |
RU2476899C1 (ru) | Гидроакустический комплекс для измерения азимутального угла и горизонта источника звука в мелком море | |
KR20210044079A (ko) | 지향성 배열 센서 기반의 원격 연안 음향 토모그래피 시스템 | |
RU2284043C1 (ru) | Способ калибровки компьютерно-интерферометрических систем на подвижных платформах | |
RU2591030C1 (ru) | Гидроакустический комплекс для обнаружения движущегося источника звука, измерения азимутального угла на источник и горизонта источника звука в мелком море |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |