CN109375197B - 一种小尺寸矢量阵低频散射校正方法 - Google Patents
一种小尺寸矢量阵低频散射校正方法 Download PDFInfo
- Publication number
- CN109375197B CN109375197B CN201811236822.4A CN201811236822A CN109375197B CN 109375197 B CN109375197 B CN 109375197B CN 201811236822 A CN201811236822 A CN 201811236822A CN 109375197 B CN109375197 B CN 109375197B
- Authority
- CN
- China
- Prior art keywords
- array
- small
- vector
- size
- size vector
- 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
- 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低频范围内任意对称结构小尺寸声呐基阵的校正领域,具体涉及一种小尺寸矢量阵低频散射校正方法。该方法包含以下步骤:(1)布放基阵与发射器;(2)调节信号发射设备形成测试声场;(3)在水平面等间隔旋转基阵一周,采集与存储所有角度接收信号;(4)利用离散傅里叶变换求取总声场傅里叶展开项系数估计值;(5)将期望指向性输出与基阵实际输出分别展开;(6)利用最小二乘法拟合期望指向性输出与实际输出,求取加权矢量矩阵。本发明直接将基阵输出作为校准依据,避免了基元不一致性、安装误差等因素对声散射校正的影响,使基阵系统声散射校正问题更直接,可操作证性强,避免了误差传递的影响,为声呐基阵后续应用提供保障。
Description
技术领域
本发明属于20-1000Hz低频范围内任意对称结构小尺寸声呐基阵的校正领域,具体涉及一种小尺寸矢量阵低频散射校正方法。
背景技术
在低频甚低频工作段,传统声呐设计的瑞利限严重制约了水下声呐基阵的探测性能。而布阵间隔远小于半波长的小尺寸基阵研究工作的突破为本问题提供了解决思路。近几年,基于小尺寸声呐基阵的超指向性波束形成方法成为研究热点,大量的设计分析及计算机仿真表明了本类型基阵的优势,但目前只有有限的成品系统被应用于实际工作中,其主要原因之一便是基阵周围声散射的影响。虽然小尺寸基阵可以灵活安装在小平台上,但工作平台或基阵载体等结构体,都是声散射体,其散射效应会导致接收声信号畸变。加之这些结构体大多都是复合而成,其材料与形状的参数复杂,使得其影响程度难以精确估计。
水声领域中的目标声散射特性研究工作,重点集中在对远场情况的讨论上。而在近场声散射问题中,大多也只考虑水声换能器自身散射对测量结果的影响,直到近些年才开始关注周围障板散射体的作用,但绝大多数解决方法只考虑了具有解析解的球体、长圆柱体等规则结构体周围的声场,对于复合结构引起的声散射,则难以精确计算。同时,要解决声障板影响下矢量传感器阵的应用,不仅要研究障板散射特性的计算方法,更是要研究相应的声信号畸变规律及解决办法。但到目前为止,特别是对于小尺寸矢量阵阵列而言,仍缺少一种灵活的、且考虑更广泛形状散射体情况的校正方法。
发明内容
本发明的目的在于针对基阵载体声散射的影响,发明一种操作简单、方便实用的小尺寸基阵校正方法,尤其是20-1000Hz低频范围内小尺寸矢量阵的校正。
本发明的目的是这样实现的,该方法包括以下步骤:
(1)布放基阵与发射器;
(2)根据所布放的基阵与发射器,调节信号发射设备形成测试声场;
(3)根据形成的测试声场,在水平面等间隔旋转基阵一周,采集与存储所有角度接收信号;
(4)根据采集与存储的所有角度接收信号,利用离散傅里叶变换求取总声场傅里叶展开项系数估计值;
(5)根据总声场傅里叶展开项系数估计值,将期望指向性输出与基阵实际输出分别展开;
(6)根据期望指向性输出与基阵实际输出的展开式,利用最小二乘法拟合期望指向性输出与实际输出,求取加权矢量矩阵。
所述步骤(1)布放基阵与发射器时以阵元间距d与波数k之间满足kd<<1设计小尺寸矢量阵,将小尺寸矢量阵固定在薄圆板上,并在圆板下方设计短圆柱腔以放置调理电路等电子设备。
所述步骤(2)调节信号发射设备形成测试声场时布放声源与小尺寸矢量阵,二者之间距离满足声学远场条件,声源的主轴方向对准小尺寸矢量阵的正横方向;调节信号发生器产生CW脉冲信号,经功率放大器加载到声源上形成测试声场。
基于结构的对称性,将基阵接收声场矢量进行傅里叶展开:
将小尺寸矢量阵看作一个高阶传感器,将其输出的高阶指向性作为依据进行校正,为此,需要分别将期望指向性输出和基阵实际输出展开,其中,将中心频率ω时各基元加权系数表示为下式:
将信号表达式限定到前M阶,得到基阵输出表达式为下式:
将基阵输出表达式写成矩阵形式为下式:
将步骤(5)所述期望指向性输出展开为下式:
式中,
γ(θ)=diag[e-iMθ,...,1,...,e-iMθ]
b=[b-M,...,b0,...,bM]T
建立小尺寸矢量阵散射校正的加权最小二乘算法,令基阵实际输出与期望指向性输出相等,得下式:
式中,上标H表示共轭转置。
利用此系数矩阵对小尺寸矢量阵进行加权处理,即可修正低频声散射的影响,解决小尺寸矢量阵的低频校正问题。本发明考虑了安装平台或基阵载体的低频散射问题与小尺寸矢量阵波束形成实现之间的主要矛盾,对于低频散射的影响进行测量并去除。
本发明的有益效果在于:将小尺寸矢量阵看作一个高阶传感器,因此可以忽略单基元的校准,直接将高阶指向性输出作为校准依据,解决了载体低频散射问题的同时,可操作性强、方便实用。本发明可以应用于低频水声计量与声散射修正的各个领域。
附图说明
图1小尺寸矢量阵示意图。
图2发明内容程序流程图。
图3(a)0号基元X通道实测指向性图。
图3(b)0号基元Y通道实测指向性图。
图3(c)0号基元P通道实测指向性图。
图3(d)1号基元X通道实测指向性图。
图3(e)1号基元Y通道实测指向性图。
图3(f)1号基元P通道实测指向性图。
图3(g)2号基元X通道实测指向性图。
图3(h)2号基元Y通道实测指向性图。
图3(i)2号基元P通道实测指向性图。
图4(a)二阶实测波束校正前后基阵输出对比图。
图4(b)三阶实测波束校正前后基阵输出对比图。
具体实施方式
下面结合附图对发明做更详细地描述:
本发明是一种针对20-1000Hz范围内小尺寸矢量阵的结构声散射校正方法,它适用于任何对称结构矢量阵系统的声散射校正工作。本发明中以五元小尺寸矢量阵系统为例,它包括五个阵元,它们构成一个十字型,阵元间距d与波数k之间满足kd<<1;载体结构由薄圆板和短圆柱腔组成;声源与小尺寸矢量接收阵之间距离满足声学远场条件。调节信号发生器产生CW脉冲信号,经功率放大器加载到声源上形成测量声场;采集与存储各阵元的接收信号。利用傅里叶求和表示总接收声场,并对实测声场进行离散傅里叶变换以求取展开项系数。忽略多项误差来源影响,不严格区分基元的校准与基阵的校准,而将小尺寸基阵看作一个高阶传感器,将其总体输出的高阶指向性作为校准依据。最后将基阵输出与期望指向性图间拟合,提出利用加权最小二乘算法的校正方法。本发明所提校正方法为水声领域小尺度基阵的声散射校正提供了一种新思路,直接将基阵输出作为校准依据,避免了基元不一致性、安装误差等因素对声散射校正的影响,使基阵系统声散射校正问题更直接,可操作证性强,避免了误差传递的影响,为声呐基阵后续应用提供保障。
(1)以图1所示小尺寸矢量基阵为例。所研究小尺寸矢量阵包括0号~4号阵元,它们构成一个十字型,十字原点处为0号阵元,1号~4号阵元与0号阵元间距d相等,且阵元间距d与波数k之间满足kd<<1。将基阵安装在比基阵尺寸略大的薄圆板上,圆板下方设计短圆柱腔以放置调理电路等电子设备,由此,安装载体与小尺寸矢量阵组成一个可自容式工作的声呐系统。
(2)将基阵吊放到自由场条件相对较好的水域中,同时将测量用辅助声源吊放到同一深度,声源距基阵几何中心距离为r满足声学远场条件,声源的主轴方向对准小尺寸矢量阵的正轴方向。调节信号发生器产生CW脉冲信号,经功率放大器加载到声源输入端,形成测试声场。
(3)在试验过程中,保持声源不动,利用机械装置在0°~360°范围内旋转基阵。利用信号采集系统对信号进行采集,存储,其中信号包括0号~4号基元的接收信号,每个阵元包括一路声压信号和两路振速通道信号。将采集的Q个不同方向(q=1,…,Q)实测导向矢量的集合写成矢量形式:
(4)基于结构的对称性,将基阵接收声场矢量进行傅里叶展开:
其中,由于平台散射影响,声压通道系数项在m取正和取负时相等,而矢量通道系数项在m取正和取负时却不一定相等,即矢量通道信号具有更复杂的形式。
(5)图3给出了不同基元声压与矢量通道的指向性图,由于基阵对称性,3#及4#基元的指向性与1#及2#类似,限于篇幅,文中将只给出0#及1#、2#基元各通道的实测指向性图,其中,每行的三个子图依左至右分别为X、Y、P通道结果。可以看出,0#基元(中心位置处)声压具有完好的全指向性,矢量通道的“8”字形指向性也依然存在;而1#与2#基元的声压通道指向性有明显的偏向,且其偏向方位与其所在位置有关,其矢量通道指向性,尤其是径向方向的通道已不再是“8”字形。可以明显观察到结构散射导致的基元指向性的变化。其中,切向方向矢量通道受到结构的影响是对称的,而不位于中心位置的基元,其径向方向矢量通道受到的影响是不对称的。
在小尺寸基阵的校正中,误差来源复杂,除了载体散射的影响外,还有基元位置误差、幅度和相位误差等因素,难以分别考虑因此将不严格区分基元的校准与基阵的校准,而将小尺寸基阵看作一个高阶传感器,将其总体波束输出的高阶指向性与输出增益作为校准依据。为此,需要分别将期望指向性输出和基阵实际输出展开。其中,将中心频率ω时各基元加权系数表示为
并将信号展开表达式限定到前M阶,得到基阵输出表达式:
写成矩阵形式为
将期望波束表示为
θ为导向角,M称为波束图的阶数,而各系数am项决定了波束图的形状。为了更直观地与接收信号形式对应,可以将波束图写为如下形式:
式中,
γ(θ)=diag[e-iMθ,...,1,...,e-iMθ]
b=[b-M,...,b0,...,bM]T
(6)建立小尺寸矢量阵散射校正的加权最小二乘算法。令基阵实际输出与期望指向性输出相等,可得
因此可以对其求取广义逆来获取加权系数矩阵的估计值:
利用此系数矩阵对小尺寸矢量阵进行加权处理,即可修正低频声散射的影响,解决小尺寸矢量阵的低频校正问题。
在图1所示系统下,由于实体结构复杂,散射系数难以理论计算,因此难以用常规基阵校正与波束形成方法。而根据本专利所述,可以利用实测结果,依据公式(3),对每个基元的各阶展开系数进行求解,并进行小尺寸矢量阵的系统校正。
以具有
系数形式的二阶和三阶波束为例,导向角为90°时基阵输出的结果如图4所示。图中,以方形标志曲线给出本节考虑散射情况下的波束形成算法结果,而不考虑散射时的波束图如三角标记曲线所示,同时用无标志线给出理想波束图作为参考。可以看出,除旁瓣略高外,声散射校正后的二阶波束主瓣及三阶波束主瓣基本与理论值一致,其中二阶波束的旁瓣畸变很小,凹点基本存在,但三阶波束与理论情况相差略大,可见二阶波束相较三阶波束更稳健。而不考虑散射引起的信号畸变时,基阵输出结果基本没有明显的主瓣出现,算法失效。
所述的由于载体或安装平台结构散射的影响,建立的矢量阵列接收模型,其特征在于基阵接收声场矢量可以进行傅里叶展开
式中,声压通道系数项在m取正和取负时相等,而矢量通道系数项在m取正和取负时却不一定相等,即矢量通道信号具有更复杂的形式。
其中第n个通道各展开项的系数估计值由下式计算获得:
将小尺寸矢量阵看作一个高阶传感器,将其输出的高阶指向性作为依据进行校正。
而期望指向性图
可展开为:
其中,θ为导向角,M称为波束图的阶数,各系数am项决定了波束图的形状,各矢量意义如下
令基阵实际输出与期望指向性输出相等,建立小尺寸矢量阵散射校正的加权最小二乘算法;加权系数矢量可由下式估计:
利用此系数矩阵对小尺寸矢量阵进行加权处理,即可修正低频声散射的影响,解决小尺寸矢量阵的低频校正问题。
Claims (5)
1.一种小尺寸矢量阵低频散射校正方法,其特征在于,该方法包括以下步骤:
(1)布放基阵与发射器;
(2)根据所布放的基阵与发射器,调节信号发射设备形成测试声场;
(4)根据采集与存储的所有角度接收信号,利用离散傅里叶变换求取总声场傅里叶展开项系数Cm_n估计值;
基于结构的对称性,将基阵接收声场矢量进行傅里叶展开:
所述总声场傅里叶展开项系数Cm_n的估计值为:
(5)根据总声场傅里叶展开项系数估计值,将期望指向性输出与基阵实际输出分别展开;
将中心频率ω时各基元加权系数表示为:
将信号表达式限定到前M阶,则基阵实际输出的展开形式为:
期望指向性输出的展开形式为:
2.根据权利要求1所述的一种小尺寸矢量阵低频散射校正方法,其特征在于:所述步骤(1)布放基阵与发射器时以阵元间距d与波数k之间满足kd<<1设计小尺寸矢量阵,将小尺寸矢量阵固定在薄圆板上,并在圆板下方设计短圆柱腔以放置调理电路等电子设备。
3.根据权利要求2所述的一种小尺寸矢量阵低频散射校正方法,其特征在于:所述步骤(2)调节信号发射设备形成测试声场时布放声源与小尺寸矢量阵,二者之间距离满足声学远场条件,声源的主轴方向对准小尺寸矢量阵的正横方向;调节信号发生器产生CW脉冲信号,经功率放大器加载到声源上形成测试声场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811236822.4A CN109375197B (zh) | 2018-10-23 | 2018-10-23 | 一种小尺寸矢量阵低频散射校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811236822.4A CN109375197B (zh) | 2018-10-23 | 2018-10-23 | 一种小尺寸矢量阵低频散射校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109375197A CN109375197A (zh) | 2019-02-22 |
CN109375197B true CN109375197B (zh) | 2023-02-10 |
Family
ID=65401402
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811236822.4A Active CN109375197B (zh) | 2018-10-23 | 2018-10-23 | 一种小尺寸矢量阵低频散射校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109375197B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112578370B (zh) * | 2020-12-07 | 2023-07-21 | 江苏科技大学 | 一种远场目标强度的近场测试方法 |
CN113608188B (zh) * | 2021-06-21 | 2023-12-29 | 西安空间无线电技术研究所 | 一种基于天基tdm-mimo雷达空间碎片的测角方法及系统 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5357962A (en) * | 1992-01-27 | 1994-10-25 | Sri International | Ultrasonic imaging system and method wtih focusing correction |
AU4169497A (en) * | 1996-08-29 | 1998-04-14 | David T. Borup | Apparatus and method for imaging with wavefields using inverse scattering techniques |
JP4874254B2 (ja) * | 2004-09-10 | 2012-02-15 | ライオス テヒノロギー ゲーエムベーハー | 光fmcw後方散乱測定システムの校正 |
CN2814340Y (zh) * | 2005-04-30 | 2006-09-06 | 中国船舶重工集团公司第七一五研究所 | 一种矢量传感器 |
GB0619825D0 (en) * | 2006-10-06 | 2006-11-15 | Craven Peter G | Microphone array |
CN100580444C (zh) * | 2006-11-08 | 2010-01-13 | 哈尔滨工程大学 | 中高频下粘弹性材料声学参数测量方法 |
JP6013493B2 (ja) * | 2011-10-03 | 2016-10-25 | サーフ テクノロジー アクティーゼルスカブSurf Technology As | デュアルバンドパルス合成物による非線形イメージング |
CN103091666B (zh) * | 2011-11-07 | 2015-02-04 | 中国科学院电子学研究所 | 非理想定标器条件下的机载p波段极化sar定标方法 |
DE102013000684B3 (de) * | 2013-01-11 | 2014-01-23 | Klippel Gmbh | Anordnung und Verfahren zur holografischen Bestimmung des Direktschalles akustischer Quellen |
CN103983338B (zh) * | 2014-05-09 | 2016-08-24 | 哈尔滨工程大学 | 多极子矢量接收阵校准方法 |
WO2016106470A1 (zh) * | 2014-12-29 | 2016-07-07 | 哈尔滨工程大学 | 智能指向性全双工水声通信机及其通信方法 |
US20170023628A1 (en) * | 2015-07-20 | 2017-01-26 | Tektronix, Inc. | Time corrected time-domain reflectometer |
CN105158751A (zh) * | 2015-08-29 | 2015-12-16 | 许昌学院 | 一种声矢量阵快速方位估计方法 |
WO2017066198A1 (en) * | 2015-10-12 | 2017-04-20 | The Regents Of The University Of California | Spectroscopy imaging and analysis of live cells |
CN106093921B (zh) * | 2016-07-25 | 2019-04-26 | 中国电子科技集团公司第五十四研究所 | 基于稀疏分解理论的声矢量阵宽带测向方法 |
CN106546966B (zh) * | 2016-10-31 | 2019-03-29 | 西安电子科技大学 | 基于多项式拟合的杂波背景下雷达噪声功率估计方法 |
CN107037418B (zh) * | 2016-11-02 | 2019-10-18 | 哈尔滨工程大学 | 一种低频超指向性多极子矢量阵的有界空间校正方法 |
CN108226866A (zh) * | 2017-12-13 | 2018-06-29 | 贵州航天计量测试技术研究所 | 基于球谐函数展开的非共形面变换噪声源的定位方法 |
-
2018
- 2018-10-23 CN CN201811236822.4A patent/CN109375197B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109375197A (zh) | 2019-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7330396B2 (en) | Farfield analysis of noise sources | |
CN107037418B (zh) | 一种低频超指向性多极子矢量阵的有界空间校正方法 | |
CN106707250B (zh) | 基于互耦校正的雷达阵列自适应波束形成方法 | |
US7274622B1 (en) | Nonlinear techniques for pressure vector acoustic sensor array synthesis | |
Yang et al. | Functional delay and sum beamforming for three-dimensional acoustic source identification with solid spherical arrays | |
CN110109053B (zh) | 一种未知声速环境下快速doa估计方法 | |
CN109375197B (zh) | 一种小尺寸矢量阵低频散射校正方法 | |
CN110687538A (zh) | 一种基于近场聚焦的超波束形成方法 | |
CN107728112A (zh) | 在目标导向矢量严重失配情况下的稳健波束形成方法 | |
CN112285647A (zh) | 一种基于稀疏表示与重构的信号方位高分辨估计方法 | |
CN110736976B (zh) | 一种任意阵形的声纳波束形成器性能估计方法 | |
CN106855619A (zh) | 一种获取mimo成像雷达系统各方向的分辨率的方法 | |
CN103048641B (zh) | 具有联合约束优化形式的矢量阵稳健聚焦处理方法 | |
CN107202975B (zh) | 一种二维矢量阵阵元姿态误差校正方法 | |
CN101398482B (zh) | 声接收基阵被动宽带检测中的噪声场数值计算方法 | |
CN111983599A (zh) | 一种基于方位-俯仰字典的目标二维doa估计方法 | |
CN110398711A (zh) | 一种声呐共形阵基于阵列流形测量的方向图综合方法 | |
CN113624330B (zh) | 一种水下目标辐射噪声测量组合体积阵及测量方法 | |
CN111551943B (zh) | 一种基于damas2的稀疏布阵高分辨三维声成像方法及系统 | |
CN114609580A (zh) | 一种基于非圆信号的无孔互质阵列设计方法 | |
CN114563760A (zh) | 一种基于sca阵型的二阶超波束形成方法、设备及介质 | |
Hua et al. | Efficient two dimensional direction finding via auxiliary-variable manifold separation technique for arbitrary array structure | |
CN110824484A (zh) | 一种基于恒模算法的阵元位置估计方法 | |
CN114841024B (zh) | 一种基于矢量扩展互质阵的波束形成方法及装置 | |
CN110489780B (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 |