CN112083430A - 一种适用于轨道角动量三维成像声呐的旁瓣抑制方法 - Google Patents
一种适用于轨道角动量三维成像声呐的旁瓣抑制方法 Download PDFInfo
- Publication number
- CN112083430A CN112083430A CN202010944466.2A CN202010944466A CN112083430A CN 112083430 A CN112083430 A CN 112083430A CN 202010944466 A CN202010944466 A CN 202010944466A CN 112083430 A CN112083430 A CN 112083430A
- Authority
- CN
- China
- Prior art keywords
- vortex
- wave
- phi
- angular momentum
- orbital angular
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000003384 imaging method Methods 0.000 title claims abstract description 22
- 230000001629 suppression Effects 0.000 title claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 9
- 239000000126 substance Substances 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 4
- 238000002372 labelling Methods 0.000 claims description 2
- 230000002401 inhibitory effect Effects 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
-
- 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/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种适用于轨道角动量三维成像声呐的旁瓣抑制方法,其技术特点是:对不同阶的轨道角动量涡旋声波的回波信号进行加权求和,获得涡旋波的波束;根据涡旋波的波束,构建旁瓣高的方位的点的集合S;对旁瓣高的方位的点的集合S进行处理,获得优化后的加权矢量;使用优化后的加权矢量对接收水听器接收的信号进行第一次波束形成;对步骤4得到的数据上进行加权求和获得第二次波束形成,遍历不同的空间方位后得到最终波束图。本发明通过优化波束形成的加权矢量,使得第一次波束形成在第二次波束旁瓣高的方位形成凹陷,降低轨道角动量三维成像声呐旁瓣级,从而达到抑制最终波束的旁瓣的目的。
Description
技术领域
本发明属于海洋声学装备技术领域,涉及一种旁瓣抑制方法,尤其是一种适用于轨道角动量三维成像声呐的旁瓣抑制方法。
背景技术
随着海洋技术的发展,人们对涡旋声波的研究也越来越重视。涡旋声波具有螺旋形的波前相位,可以对空间信息进行调制,提高声波的信息传递和获取能力。因此采用轨道角动量涡旋声波作为发射信号的相控阵三维成像声呐,不仅可以提高成像精度,还能够大大减少所需接收面阵的水听器数量。但是,接收面阵的稀疏性质,会带来接收面阵波束的旁瓣级的提高,例如,采用64个阵元的接收面阵的轨道角动量三维成像声呐,旁瓣级只有-9.8dB。因此,如何有效地抑制轨道角动量三维成像声呐的旁瓣级是目前迫切需要解决的问题。
发明内容
本发明的目的在于克服现有技术的不足,提供一种设计合理且能有效降低轨道角动量三维成像声呐旁瓣级的适用于轨道角动量三维成像声呐的旁瓣抑制方法。
本发明解决其技术问题是采取以下技术方案实现的:
一种适用于轨道角动量三维成像声呐的旁瓣抑制方法,包括以下步骤:
步骤1:对不同阶的轨道角动量涡旋声波的回波信号进行加权求和,获得涡旋波的波束;
步骤2:根据涡旋波的波束,构建旁瓣高的方位的点的集合S;
步骤3:采用凹槽噪声算法、优化最低旁瓣法或旁瓣约束高增益优化法,对旁瓣高的方位的点的集合S进行处理,获得优化后的加权矢量;
步骤4:使用优化后的加权矢量对接收水听器接收的信号进行第一次波束形成;
步骤5:对步骤4得到的数据上进行加权求和获得第二次波束形成,遍历不同的空间方位后得到最终波束图。
而且,所述步骤1的具体实现方法为:
在正方形的接收面阵上随机布放若干个接收水听器,以正方形的接收面阵的中心为圆心,布放一个半径为a的发射圆环阵,在发射圆环阵上均匀布放N个发射换能器,以正X轴上的发射换能器为第一个发射基元,对发射换能器进行标号:1,2,…,N;第n个发射换能器发射的信号为:
其中,φn=2πnl/N代表发射调相,l为发射涡旋波的阶数,ω为发射的频率;
利用发射圆环阵列依次发射不同阶数l的涡旋声波,对不同阶l的涡旋声波的回波信号进行加权求和,得到涡旋波的波束:
而且,所述步骤2的具体实现方法为:设置一个阈值∈,对涡旋波的波束归一化幅值超过该阈值的集合,对该集合离散化得到旁瓣高的方位的点的集合S:
S={(φ1,θ1),···,(φJ,θJ)}
这里(φi,θi)为旁瓣高于阈值∈的方位的俯仰角和方位角坐标。
而且,所述步骤3的凹槽噪声算法为:在集合S的方位放置虚拟信号干扰源,得到基阵数据的协方差矩阵为:
则优化后的加权矢量为:
而且,所述步骤4的第一次波束形成的结果为:
而且,所述步骤5中第二次波束形成的结果为:
本发明的优点和积极效果是:
本发明设计合理,其通过优化第一次波束形成的加权矢量,使得第一次波束形成在第二次波束旁瓣高的方位形成凹陷,降低轨道角动量三维成像声呐旁瓣级,从而达到抑制最终波束的旁瓣的目的,不仅可以提高成像精度,还能够大大减少所需的接收面阵的水听器的数量。
附图说明
图1为本发明采用的收发阵列示意图。
图2为采用常规算法仿真得到的波束图。
图3为按照本发明仿真得到的波束图。
具体实施方式
以下结合附图对本发明实施例做进一步详述:
一种适用于轨道角动量三维成像声呐的旁瓣抑制方法,包括以下步骤:
在本实施例中,考虑一个以50倍波长为边长的正方形的接收面阵,在正方形的接收面阵上随机布放64个接收水听器,以正方形的接收面阵的中心为圆心,布放一个半径为a的发射圆环阵,在发射圆环阵上均匀布放N个发射换能器,如图1所示为该收发阵列的示意图,以正X轴上的发射换能器为第一个发射基元,对发射换能器进行标号:1,2,…,N;第n个发射换能器发射的信号为:
其中,φn=2πnl/N代表发射调相,l为发射涡旋波的阶数,ω为发射的频率;
利用发射圆环阵列依次发射不同阶数l的涡旋声波,对不同阶l的涡旋声波的回波信号进行加权求和,得到涡旋波的波束:
步骤2:根据涡旋波的波束,构建旁瓣高的方位的点的集合S;
在本步骤中,设置一个阈值∈,选出步骤1得到的涡旋波的波束幅值超过该阈值的集合,并对该集合离散化得到旁瓣高的方位的点的集合S:
S={(φ1,θ1),···,(φJ,θJ)}。
步骤3:采用凹槽噪声算法、优化最低旁瓣法、旁瓣约束高增益优化法等算法任意一种方法,对旁瓣高的方位的点的集合S进行处理,获得优化后的加权矢量。
本实施例以凹槽噪声法为例对集合S中的方位进行处理。其在旁瓣区人为放置若干虚拟干扰源,可以获得旁瓣抑制的波束形成加权矢量。考虑在集合S的方位放置虚拟信号干扰源,则可以得到基阵数据的协方差矩阵为:
进而获得优化后的加权矢量为:
步骤4:用步骤3得到的优化后的加权矢量对接收水听器接收的信号进行第一次波束形成,第一次波束形成表示如下:
本步骤利用优化后的加权矢量,使得第一次的波束形成在方位点S处的波束产生凹陷,从而达到抑制最终波束的旁瓣的目的。
步骤5:在步骤4得到的数据上进行第二次波束形成,得到遍历不同的空间方位(φ,θ)后得到最终的波束图s2(φ,θ)。
遍历不同的空间方位得到最终的波束图。
如图2所示为采用常规算法仿真得到的波束图,图3为本方法仿真得到的波束图。仿真采用32个发射换能器发射涡旋信号,并由64个接收水听器组成的随机稀疏面阵进行回波接收,对准方向横向角和纵向角都为10°。图2为优化前的波束图,由于只采用了64个水听器接收回波信号,稀疏率仅为0.6%,但是也导致了旁瓣级较高为-9.68dB。图3为优化后的波束图,旁瓣级为-12.75dB。可以看到,采用了凹槽噪声算法,降低了轨道角动量三维成像声呐的波束的旁瓣级。
需要强调的是,本发明所述的实施例是说明性的,而不是限定性的,因此本发明并不限于具体实施方式中所述的实施例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。
Claims (6)
1.一种适用于轨道角动量三维成像声呐的旁瓣抑制方法,其特征在于:包括以下步骤:
步骤1:对不同阶的轨道角动量涡旋声波的回波信号进行加权求和,获得涡旋波的波束;
步骤2:根据涡旋波的波束,构建旁瓣高的方位的点的集合S;
步骤3:采用凹槽噪声算法、优化最低旁瓣法或旁瓣约束高增益优化法,对旁瓣高的方位的点的集合S进行处理,获得优化后的加权矢量;
步骤4:使用优化后的加权矢量对接收水听器接收的信号进行第一次波束形成;
步骤5:对步骤4得到的数据上进行加权求和获得第二次波束形成,遍历不同的空间方位后得到最终波束图。
2.根据权利要求1所述的适用于轨道角动量三维成像声呐的旁瓣抑制方法,其特征在于:所述步骤1的具体实现方法为:在正方形的接收面阵上随机布放若干个接收水听器,以正方形的接收面阵的中心为圆心,布放一个半径为a的发射圆环阵,在发射圆环阵上均匀布放N个发射换能器,以正X轴上的发射换能器为第一个发射基元,对发射换能器进行标号:1,2,…,N;第n个发射换能器发射的信号为:
其中,φn=2πnl/N代表发射调相,l为发射涡旋波的阶数,ω为发射的频率;
利用发射圆环阵列依次发射不同阶数l的涡旋声波,对不同阶l的涡旋声波的回波信号进行加权求和,得到涡旋波的波束:
3.根据权利要求1所述的适用于轨道角动量三维成像声呐的旁瓣抑制方法,其特征在于:所述步骤2的具体实现方法为:设置一个阈值∈,对涡旋波的波束归一化幅值超过该阈值的集合,对该集合离散化得到旁瓣高的方位的点的集合S:
S={(φ1,θ1),···,(φJ,θJ)}
这里(φi,θi)为旁瓣高于阈值∈的方位的俯仰角和方位角坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010944466.2A CN112083430B (zh) | 2020-09-10 | 2020-09-10 | 一种适用于轨道角动量三维成像声呐的旁瓣抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010944466.2A CN112083430B (zh) | 2020-09-10 | 2020-09-10 | 一种适用于轨道角动量三维成像声呐的旁瓣抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112083430A true CN112083430A (zh) | 2020-12-15 |
CN112083430B CN112083430B (zh) | 2023-11-07 |
Family
ID=73731752
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010944466.2A Active CN112083430B (zh) | 2020-09-10 | 2020-09-10 | 一种适用于轨道角动量三维成像声呐的旁瓣抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112083430B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112764041A (zh) * | 2021-01-19 | 2021-05-07 | 鹏城实验室 | 成像系统和方法 |
CN112904347A (zh) * | 2021-01-19 | 2021-06-04 | 鹏城实验室 | 成像系统和方法 |
CN112904345A (zh) * | 2021-01-19 | 2021-06-04 | 鹏城实验室 | 定位系统及定位方法 |
CN113009480A (zh) * | 2021-02-20 | 2021-06-22 | 鹏城实验室 | 一种成像方法及系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100013696A1 (en) * | 2008-07-18 | 2010-01-21 | Raytheon Company | System and method of orbital angular momentum (oam) diverse signal processing using classical beams |
CN103984103A (zh) * | 2014-02-18 | 2014-08-13 | 上海大学 | 一种利用计算全息图产生涡旋光束的方法 |
CN105842702A (zh) * | 2016-05-13 | 2016-08-10 | 南京航空航天大学 | 多波束成像声呐的旁瓣抑制方法、阵列稀疏方法 |
CN106329108A (zh) * | 2016-10-31 | 2017-01-11 | 宁夏大学 | 一种双环结构多模态oam电磁涡旋波阵列天线 |
CN106886020A (zh) * | 2017-03-06 | 2017-06-23 | 中国人民解放军国防科学技术大学 | 一种单天线接收条件下的电磁涡旋成像方法 |
CN108594221A (zh) * | 2018-04-27 | 2018-09-28 | 中国人民解放军国防科技大学 | 基于同心圆环阵列的涡旋电磁波产生与优化方法 |
CN108919199A (zh) * | 2018-07-12 | 2018-11-30 | 南京航空航天大学 | 多波束成像声纳稀疏阵列的旁瓣抑制方法、阵列稀疏方法 |
CN110824433A (zh) * | 2019-10-09 | 2020-02-21 | 清华大学 | 电磁波量子态轨道角动量雷达探测和方法 |
CN110988868A (zh) * | 2019-11-19 | 2020-04-10 | 南京理工大学 | 一种等效分数阶模式涡旋电磁波产生及成像方法 |
-
2020
- 2020-09-10 CN CN202010944466.2A patent/CN112083430B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100013696A1 (en) * | 2008-07-18 | 2010-01-21 | Raytheon Company | System and method of orbital angular momentum (oam) diverse signal processing using classical beams |
CN103984103A (zh) * | 2014-02-18 | 2014-08-13 | 上海大学 | 一种利用计算全息图产生涡旋光束的方法 |
CN105842702A (zh) * | 2016-05-13 | 2016-08-10 | 南京航空航天大学 | 多波束成像声呐的旁瓣抑制方法、阵列稀疏方法 |
CN106329108A (zh) * | 2016-10-31 | 2017-01-11 | 宁夏大学 | 一种双环结构多模态oam电磁涡旋波阵列天线 |
CN106886020A (zh) * | 2017-03-06 | 2017-06-23 | 中国人民解放军国防科学技术大学 | 一种单天线接收条件下的电磁涡旋成像方法 |
CN108594221A (zh) * | 2018-04-27 | 2018-09-28 | 中国人民解放军国防科技大学 | 基于同心圆环阵列的涡旋电磁波产生与优化方法 |
CN108919199A (zh) * | 2018-07-12 | 2018-11-30 | 南京航空航天大学 | 多波束成像声纳稀疏阵列的旁瓣抑制方法、阵列稀疏方法 |
CN110824433A (zh) * | 2019-10-09 | 2020-02-21 | 清华大学 | 电磁波量子态轨道角动量雷达探测和方法 |
CN110988868A (zh) * | 2019-11-19 | 2020-04-10 | 南京理工大学 | 一种等效分数阶模式涡旋电磁波产生及成像方法 |
Non-Patent Citations (2)
Title |
---|
XIAOLI YIN ET AL.: "Sidelobe suppression in the generation of beam carrying orbital angular momentum", 《2014 13TH INTERNATIONAL CONFERENCE ON OPTICAL COMMUNICATIONS AND NETWORKS (ICOCN)》, pages 1 - 4 * |
李志焕 等: "涡旋光束旁瓣抑制的研究进展", 《激光杂志》, vol. 34, no. 2, pages 1 - 4 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112764041A (zh) * | 2021-01-19 | 2021-05-07 | 鹏城实验室 | 成像系统和方法 |
CN112904347A (zh) * | 2021-01-19 | 2021-06-04 | 鹏城实验室 | 成像系统和方法 |
CN112904345A (zh) * | 2021-01-19 | 2021-06-04 | 鹏城实验室 | 定位系统及定位方法 |
WO2022156000A1 (zh) * | 2021-01-19 | 2022-07-28 | 鹏城实验室 | 定位系统及定位方法 |
CN112764041B (zh) * | 2021-01-19 | 2023-05-26 | 鹏城实验室 | 成像系统和方法 |
CN112904347B (zh) * | 2021-01-19 | 2023-05-26 | 鹏城实验室 | 成像系统和方法 |
CN113009480A (zh) * | 2021-02-20 | 2021-06-22 | 鹏城实验室 | 一种成像方法及系统 |
CN113009480B (zh) * | 2021-02-20 | 2022-04-19 | 鹏城实验室 | 一种成像方法及系统 |
WO2022174500A1 (zh) * | 2021-02-20 | 2022-08-25 | 鹏城实验室 | 一种成像方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN112083430B (zh) | 2023-11-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112083430B (zh) | 一种适用于轨道角动量三维成像声呐的旁瓣抑制方法 | |
CA2204298C (en) | Circularly symmetric, zero redundancy, planar array having broad frequency range applications | |
CN108828603B (zh) | 一种基于十字型的三维成像声纳阵列的稀疏优化方法 | |
CN107657070B (zh) | 基于遗传算法频率分集共形阵列波束的解耦合方法及天线 | |
CN112083431B (zh) | 基于声学轨道角动量的低频远程三维成像方法 | |
CN112083432B (zh) | 基于声学轨道角动量的超精细三维成像方法 | |
CN112162266A (zh) | 一种基于凸优化理论的共形阵二维波束优化方法 | |
Stergiopoulos | Implementation of adaptive and synthetic-aperture processing schemes in integrated active-passive sonar systems | |
CN110988868A (zh) | 一种等效分数阶模式涡旋电磁波产生及成像方法 | |
Chi et al. | High-resolution real-time underwater 3-D acoustical imaging through designing ultralarge ultrasparse ultra-wideband 2-D arrays | |
CN110736976B (zh) | 一种任意阵形的声纳波束形成器性能估计方法 | |
CN112698337B (zh) | 一种宽带三维成像声呐稀疏布阵方法 | |
CN109669172B (zh) | 基于主瓣内强干扰抑制的弱目标方位估计方法 | |
JP2021005861A (ja) | フェーズドアレイアンテナを用いた高速空間探索 | |
CN111381212A (zh) | 基于子阵划分的虚拟超短基线定位方法 | |
CN113126030B (zh) | 基于宽带声场干涉结构的深海直达声区目标深度估计方法 | |
CN111025272B (zh) | 具备隧道效应抑制能力的平面声学基阵超宽覆盖波束发射方法 | |
CN110244288A (zh) | 一种基于聚焦归位原理的水声阵列信号处理方法 | |
Elechi et al. | Analysis of a SONAR detecting system using multi-beamforming algorithm | |
CN113009518B (zh) | 一种卫星导航信号的多波束抗干扰方法 | |
Pillai et al. | Product beamforming and nested array in tandem for enhanced sonar performance | |
RU2515133C1 (ru) | Сферическая гидроакустическая антенна | |
CN115825883B (zh) | 星载gnss-s雷达实时直达信号干扰抑制方法 | |
RU2736097C1 (ru) | Способ обнаружения локального объекта на фоне распределенной помехи | |
CN112433218B (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 |