CN107037418B - 一种低频超指向性多极子矢量阵的有界空间校正方法 - Google Patents
一种低频超指向性多极子矢量阵的有界空间校正方法 Download PDFInfo
- Publication number
- CN107037418B CN107037418B CN201610945407.0A CN201610945407A CN107037418B CN 107037418 B CN107037418 B CN 107037418B CN 201610945407 A CN201610945407 A CN 201610945407A CN 107037418 B CN107037418 B CN 107037418B
- Authority
- CN
- China
- Prior art keywords
- array
- directivity
- super
- vector
- sound source
- 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
- 230000005405 multipole Effects 0.000 title claims abstract description 19
- 238000012937 correction Methods 0.000 claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 11
- 238000012360 testing method Methods 0.000 claims abstract description 7
- 238000010586 diagram Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000003491 array Methods 0.000 abstract description 2
- 239000003643 water by type Substances 0.000 abstract 1
- 238000013461 design Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 101100042630 Caenorhabditis elegans sin-3 gene Proteins 0.000 description 2
- 101100234408 Danio rerio kif7 gene Proteins 0.000 description 2
- 101100221620 Drosophila melanogaster cos gene Proteins 0.000 description 2
- 101100398237 Xenopus tropicalis kif11 gene Proteins 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52004—Means for monitoring or calibrating
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明属于水下声矢量阵校正领域,具体涉及一种20‑1000Hz低频范围内超指向性多极子矢量阵的校正的低频超指向性多极子矢量阵的有界空间校正方法。本发明包括:建立超指向性多极子矢量阵包括0号‑4号阵元;调节信号发生器产生CW脉冲信号,经功率放大器加载到声源上形成测试声场;在0‑360°范围内旋转超指向性多极子矢量阵,采集与存储0号~4号基元的接收信号,每个阵元包括一路声压信号和两路振速通道信号等。本发明克服了目前小尺度基阵的校正需要大尺度自由场空间的难题,提供一种在有界空间水域、可操作性强、方便实用的校准方法,可广泛应用于低频超指向性矢量阵的校准,因此本发明可以广泛应用于低频水声计量各领域。
Description
技术领域
本发明属于水下声矢量阵校正领域,具体涉及一种20-1000Hz低频范围内超指向性多极子矢量阵的校正的低频超指向性多极子矢量阵的有界空间校正方法。
背景技术
在低频甚低频工作段,传统声呐设计的瑞利限严重制约了水下声基阵的应用及其方位估计性能。而超指向性基阵研究工作的突破为本问题提供了解决思路。近几年,基于小孔径的超指向性基阵波束形成方法成为研究热点,大量的设计分析及计算机仿真表明了本类型基阵的优势,但目前只有有限的成品***被应用于实际工作中,其主要原因之一便是本类型基阵较强的敏感性。
在利用基阵进行方位估计时,需要已知基阵的阵列流形。理论上说,阵列流形在设计时便已给出,但实际中传感器灵敏度误差、相位不一致、阵元间耦合等因素都会导致二者之间的差异。现有的大量研究及发明工作都对这些因素进行了补偿校正,但都只是针对常规基阵。超指向性基阵的校正工作更为复杂,由于超指向性基阵的基元之间排列紧密,因此容易引起互散射等效应,而互散射通常与频率有关,而且安装误差等会影响基阵指向性,这些都难以建模并消除。水声领域超指向性基阵的校正工作更为缺乏。由于实际工作场地的原因,能在水池环境进行校正的频率被限定的非常严格,低频校正难度大大增加;由于超指向性基阵高幅相误差灵敏度等因素的存在,必须考虑安装误差及传感器后续电路不完全一致等***误差,使得超指向性基阵整体校正迫在眉睫,而其基阵尺寸对于校正空间的不可忽略性使得问题更为复杂。综合以上因素,建立超指向性矢量多极子阵的有限空间校正模型及其关键技术突破是超指向性基阵应用中必须解决的问题。
发明内容
本发明的目的在于针对水下有界空间的影响,提供一种操作简单、方便实用的低频超指向性多极子矢量阵的有界空间校正方法。
本发明的目的是这样实现的:
(1)建立超指向性多极子矢量阵包括0号-4号阵元,构成一个十字型,十字原点处为0 号阵元,1号~4号阵元与0号阵元间距d相等,且阵元间距d与波数k之间满足kd<1;声源的主轴方向对准超指向性多极子矢量阵的正横方向;
(2)调节信号发生器产生CW脉冲信号,经功率放大器加载到声源上形成测试声场;
(3)在0-360°范围内旋转超指向性多极子矢量阵,采集与存储0号~4号基元的接收信号,每个阵元包括一路声压信号和两路振速通道信号;
(4)将工作环境设置参数代入到球面波声场中的矢量阵列模型中,计算理论阵列流形;
声源距基阵几何中心距离为r,任一初始位置位于与X正轴夹角θ0位置的基元,在任意声源入射方向θ时,距声源位置为
其矢量通道指向性图为
有界空间中的阵列流形矢量为:
其中,表示第i号基元的X通道指向性;表示第i号基元的Y通道指向性, i=0,1,…,4;
(5)对存储的接收信号进行分析,将结果与理论模型进行对比:对于每个通道,画出实测数据的幅度随声源入射角度变化图,即指向性图,并与理论指向性图对比;画出实测数据的相位随声源入射角度变化图,并与理论情况对比;若实测情况与本发明中的理论计算结果相符,则表明阵列误差与方位无关,阵列仅存在通道幅相误差;
(6)建立基阵阵列流形校正的加权最小二乘算法;根据所建模型及实测数据分析,设定校正矩阵为对角阵
和分别代表第n个测量角度上的实测阵列流形和理论阵列流形,An中的元素表示对应位置上通道的幅相不一致量;而基阵的整体校正算法为:
或
利用参数αn对数据可信度进行调整;加权最小二乘算法的最优解为
选择数据均方幅值作为加权值:
本发明的有益效果在于:本发明考虑了在有界空间中小尺度基阵校准的主要矛盾,对于有界空间的影响进行了建模并去除。本发明克服了目前小尺度基阵的校正需要大尺度自由场空间的难题,提供一种在有界空间水域、可操作性强、方便实用的校准方法,可广泛应用于低频超指向性矢量阵的校准,因此本发明可以广泛应用于低频水声计量各领域。
附图说明
图1基阵校正模型;
图2发明内容程序流程图;
图3 1X通道幅度特性;
图4 1X通道相位特性;
图5(a)为水池试验所得校正因子实部;
图5(b)为水池试验所得校正因子虚部;
图6(a)为cos2θ基阵高阶指向性;
图6(b)为sin2θ基阵高阶指向性;
图6(c)为cos3θ基阵高阶指向性;
图6(d)为sin3θ基阵高阶指向性;
图7(a)cos2θ校正后基阵高阶指向性;
图7(b)sin2θ校正后基阵高阶指向性;
图7(c)cos3θ校正后基阵高阶指向性;
图7(d)sin3θ校正后基阵高阶指向性。
具体实施方式
下面结合附图对发明做更详细地描述:
(1)超指向性多极子矢量阵包括0号~4号阵元,它们构成一个十字型,十字原点处为0 号阵元,1号~4号阵元与0号阵元间距d相等,且阵元间距d与波数k之间满足kd<<1。将基阵通过转接法兰安装在测量水池的行车上,放入水池深度约1/2处,其中,0号基元位于阵中心,且吊放在机械旋转装置的正下方。同时将测量用辅助声源吊放到同一深度,声源距基阵几何中心距离为r满足声学远场条件,声源的主轴方向对准超指向性多极子矢量阵的正横方向。
(2)调节信号发射器产生需要频率的CW脉冲信号,经功率放大器加载到声源输入端,形成测试声场。
(3)在试验过程中,保持声源不动,利用机械装置在0°~360°范围内旋转基阵。利用信号采集***对信号进行采集,存储,其中信号包括0号~4号基元的接收信号,每个阵元包括一路声压信号和两路振速通道信号。
(4)根据工作环境,建立球面波声场中的矢量阵列接收模型。由于水池空间有限,且在低频工作段,为使直达波信号尽量长,并尽量提高信噪比,基阵中心与发射器间距离有限,且与基阵孔径可比拟,因此必须考虑声源按照球面波方式传播。如图1中所示,设声源在红色圆点处,声源距基阵几何中心距离为r。根据几何关系,任一初始位置位于与X正轴夹角θ0位置(逆时针为正)的基元,在声源旋转一周的过程中,当声源从任意方向θ入射时,它距声源位置为
其中幅度随传播距离衰减,且与距离减成反比关系,即对于声压通道,其信号可以表示为
而对于各传感器的矢量通道而言,由于其投影轴的关系,声源相对各基元坐标轴的入射角度与原点处基元并不相同,因此各基元矢量通道的指向性图需要考虑图1中θ'的变化情况。根据几何关系,任一初始位置位于与X正轴夹角θ0位置(逆时针为正)的基元矢量通道的指向性可以表示为
由此,可建立球面波声场中的矢量阵列流形矢量:
尽管此矢量是在五元阵的情形下得出,但对于其他阵型本分析方法同样适用。
(5)对采集到的不同角度入射的15路信号进行后置处理:将各通道接收的所有角度上的信号进行截取,获得含有直达波的有效信号段,然后将数据进行FFT处理。保存所有通道所有角度上FFT处理后的频域信号,作为后续分析的主要目标。
对于特定通道,分析接收信号随声源入射角度变化情况。此时,取对应位置上所保存数据的幅度值,将幅度值随角度的变化情况在极坐标图中画出,并与建模情况进行对比,1X通道的幅度对比情况如图3所示。取对应位置上所保存数据的相位值,将相位随角度的变化情况在直角坐标系中画出,并与建模情况进行对比,1X通道相位的对比情况如图4所示。其它通道的对比情况可同理画出。分析对比结果,在结果基本一致的前提下,认为实测情况符合所建模型,可进行后续分析。
通过相位随角度的变化情况可以看出,基元的实测相位变化规律与自身站位造成的相位(即理论相位)变化规律一致,仅相差一个常数,此常数即是表明基元间的初始相位不一致的量。由以上分析可知,阵列误差与方位基本无关,阵列仅存在通道幅相误差。因此校正矩阵可设定为对角阵。
(6)对于特定入射角度(第n个测量角度),取对应位置上的数据。选取工作频点处的数值,作为第n个测量角度上的实测阵列流形而理论阵列流形由实际参数代入到下式中即可获得。
最后对公式进行求解
选择数据均方幅值作为加权值:
由此获得校正因子,获得的校准因子矩阵为复对角阵,其实部和虚部如图5(a)和5(b)所示。从校正因子可以看出,X与Y通道基本保持一致,幅相误差很小,而它们与P通道之间差别则较大,且1基元的P通道与其它基元P之间反相,这可以从实验数据波形中得到验证。通过校正因子的质量可以看出,本基阵中各基元一致性较好,无须校正即可实现较理想的二阶指向性图。将校准前基阵的高阶指向性图(图(6))与校准后基阵的高阶指向性图(图(7)) 对比,可见二阶指向性得到了优化,三阶指向性得到了改善,尤其是指向性的凹点部分改善明显。
Claims (1)
1.一种低频超指向性多极子矢量阵的有界空间校正方法,其特征在于:
(1)建立超指向性多极子矢量阵包括0号-4号阵元,构成一个十字型,十字原点处为0号阵元,1号~4号阵元与0号阵元间距d相等,且阵元间距d与波数k之间满足kd<1;声源的主轴方向对准超指向性多极子矢量阵的正横方向;
(2)调节信号发生器产生CW脉冲信号,经功率放大器加载到声源上形成测试声场;
(3)在0-360°范围内旋转超指向性多极子矢量阵,采集与存储0号~4号基元的接收信号,每个阵元包括一路声压信号和两路振速通道信号;
(4)将工作环境设置参数代入到球面波声场中的矢量阵列模型中,计算理论阵列流形;
声源距基阵几何中心距离为r,任一初始位置位于与X正轴夹角θ0位置的基元,在任意声源入射方向θ时,距声源位置为
其矢量通道指向性图为
有界空间中的阵列流形矢量为:
其中,表示第i号基元的X通道指向性;表示第i号基元的Y通道指向性,i=0,1,…,4;
(5)对存储的接收信号进行分析,将结果与理论模型进行对比:对于每个通道,画出实测数据的幅度随声源入射角度变化图,即指向性图,并与理论指向性图对比;画出实测数据的相位随声源入射角度变化图,并与理论情况对比;若实测情况与本发明中的理论计算结果相符,则表明阵列误差与方位无关,阵列仅存在通道幅相误差;
(6)建立基阵阵列流形校正的加权最小二乘算法;根据所建模型及实测数据分析,设定校正矩阵为对角阵
和分别代表第n个测量角度上的实测阵列流形和理论阵列流形,An中的元素表示对应位置上通道的幅相不一致量;而基阵的整体校正算法为:
或
利用参数αn对数据可信度进行调整;加权最小二乘算法的最优解为
选择数据均方幅值作为加权值:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610945407.0A CN107037418B (zh) | 2016-11-02 | 2016-11-02 | 一种低频超指向性多极子矢量阵的有界空间校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610945407.0A CN107037418B (zh) | 2016-11-02 | 2016-11-02 | 一种低频超指向性多极子矢量阵的有界空间校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107037418A CN107037418A (zh) | 2017-08-11 |
CN107037418B true CN107037418B (zh) | 2019-10-18 |
Family
ID=59533125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610945407.0A Active CN107037418B (zh) | 2016-11-02 | 2016-11-02 | 一种低频超指向性多极子矢量阵的有界空间校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107037418B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763740B (zh) * | 2018-05-28 | 2019-12-27 | 西北工业大学 | 一种基于双振速传感器声探头灵活指向性图的设计方法 |
CN109283511B (zh) * | 2018-09-01 | 2022-07-29 | 哈尔滨工程大学 | 一种宽覆盖多波束接收基阵校准方法 |
CN109374107B (zh) * | 2018-10-23 | 2021-01-05 | 哈尔滨工程大学 | 一种多途环境中利用极化特征截取直达波的方法 |
CN109375197B (zh) * | 2018-10-23 | 2023-02-10 | 哈尔滨工程大学 | 一种小尺寸矢量阵低频散射校正方法 |
CN109471113B (zh) * | 2018-11-06 | 2023-06-02 | 哈尔滨工程大学 | 基于相位法的多波束声呐海底地形测量质量实时评估方法 |
CN110133627B (zh) * | 2019-05-21 | 2022-06-14 | 哈尔滨工程大学 | 水声定位导航***阵元位置校准测量点间距优化方法 |
CN114280618A (zh) * | 2021-12-08 | 2022-04-05 | 哈尔滨工程大学 | 一种基于声矢量水听器的小尺寸体积阵三维波束形成方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1580813A (zh) * | 2004-05-17 | 2005-02-16 | 中国人民解放军空军雷达学院 | 用高精度辅助阵元进行阵列校正与信源测向的方法 |
CN101436890A (zh) * | 2008-12-10 | 2009-05-20 | ***通信集团广东有限公司 | 一种天线阵列误差参数估计方法及其设备 |
CN102901965A (zh) * | 2012-09-26 | 2013-01-30 | 同济大学 | 一种海洋杂波环境中电大尺寸目标雷达成像方法 |
CN103983338A (zh) * | 2014-05-09 | 2014-08-13 | 哈尔滨工程大学 | 多极子矢量接收阵校准方法 |
-
2016
- 2016-11-02 CN CN201610945407.0A patent/CN107037418B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1580813A (zh) * | 2004-05-17 | 2005-02-16 | 中国人民解放军空军雷达学院 | 用高精度辅助阵元进行阵列校正与信源测向的方法 |
CN101436890A (zh) * | 2008-12-10 | 2009-05-20 | ***通信集团广东有限公司 | 一种天线阵列误差参数估计方法及其设备 |
CN102901965A (zh) * | 2012-09-26 | 2013-01-30 | 同济大学 | 一种海洋杂波环境中电大尺寸目标雷达成像方法 |
CN103983338A (zh) * | 2014-05-09 | 2014-08-13 | 哈尔滨工程大学 | 多极子矢量接收阵校准方法 |
Non-Patent Citations (2)
Title |
---|
一种新型低频声矢量场采集***的涉及与实现;刑传玺 等;《仪器仪表学报》;20131231;第34卷(第12期);2800-2805 * |
基于超指向性多极子矢量阵的水下低频声源方位估计方法研究;郭俊媛 等;《物理学报》;20160530;第65卷(第13期);1-14 * |
Also Published As
Publication number | Publication date |
---|---|
CN107037418A (zh) | 2017-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107037418B (zh) | 一种低频超指向性多极子矢量阵的有界空间校正方法 | |
CN106443587B (zh) | 一种高分辨率的快速反卷积声源成像算法 | |
CN111948615A (zh) | 一种星载全极化sar数据的极化定标方法及装置 | |
CN102997988B (zh) | 大型潜标矢量水听器低频声学指向性的水池测试方法 | |
CN107884741A (zh) | 一种多球阵列多宽带声源快速定向方法 | |
CN102944872B (zh) | 雷达散射截面近场到近场的变换方法 | |
CN107783086B (zh) | 用于诊断天线阵口径幅相场的畸变位置的方法 | |
CN104200110B (zh) | 基于去耦合的二维波束扫描测角方法 | |
CN107976671B (zh) | 一种适用于稀疏阵列天线的雷达目标角度计算方法 | |
Sun et al. | Array geometry calibration for underwater compact arrays | |
CN111220942B (zh) | 一种接收换能器阵列幅相一致性近场校准方法 | |
CN105277932B (zh) | 一种基于下变频波束形成中的多普勒频移校正方法 | |
CN109375197B (zh) | 一种小尺寸矢量阵低频散射校正方法 | |
CN107202975B (zh) | 一种二维矢量阵阵元姿态误差校正方法 | |
Yuan et al. | Correction method for magnitude and phase variations in acoustic arrays based on focused beamforming | |
Luo et al. | Accurate underwater localization through phase difference | |
CN104020465A (zh) | 基于八单元小孔径圆阵天线的外辐射源雷达测角方法 | |
CN111965602B (zh) | 一种相控阵雷达幅相一致性检测方法和*** | |
CN113821907B (zh) | 针对大型平面天线阵列***的幅度和相位自动校准方法 | |
CN116203505B (zh) | 基于块稀疏贝叶斯的正交匹配追踪声源识别方法及装置 | |
CN113311403B (zh) | 基于时间反演技术的雷达远场目标定位方法 | |
Hu et al. | Amplitude and phase errors analysis on spherical multi-probe antenna near-field measurement | |
Filik et al. | Uniform and nonuniform V-shaped planar arrays for 2-D direction-of-arrival estimation | |
CN108957411A (zh) | 一种高精度阵列信号生成方法 | |
CN114355280B (zh) | 多传感器复合阵列天线布阵和多信息融合分选测角方法 |
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 |