CN106772326A - 一种多子阵合成孔径声纳相位误差分析方法 - Google Patents
一种多子阵合成孔径声纳相位误差分析方法 Download PDFInfo
- Publication number
- CN106772326A CN106772326A CN201611102892.1A CN201611102892A CN106772326A CN 106772326 A CN106772326 A CN 106772326A CN 201611102892 A CN201611102892 A CN 201611102892A CN 106772326 A CN106772326 A CN 106772326A
- Authority
- CN
- China
- Prior art keywords
- array element
- subsystem
- analysis method
- phase error
- synthetic aperture
- 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
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
-
- 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
- G01S15/8902—Side-looking sonar
- G01S15/8904—Side-looking sonar using synthetic aperture techniques
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
本发明提供了一种多子阵合成孔径声纳相位误差分析方法,包含以下步骤:针对每个接收阵元和发射阵元组成的子系统,计算双程斜距历程;针对每个接收阵元和发射阵元组成的子系统,基于数值分析方法计算相位驻留点;针对每个接收阵元和发射阵元所组成的子系统,基于数值分析方法的相位驻留点,计算二维频域系统函数;计算待评估模型的二维频域系统函数;针对每个接收阵元和发射阵元所组成的子系统,根据多子阵合成孔径声纳的方位空不变和距离空变特性,设置目标空间坐标;针对每个接收阵元和发射阵元所组成的子系统,在二维频域内计算相位误差,克服了传统成像模型或成像算法评估周期长、数值指标不够直观等问题。
Description
技术领域
本发明属于信号处理领域,特别涉及一种多子阵合成孔径声纳相位误差分析方法。
背景技术
合成孔径声纳技术依靠小尺寸基元的匀速直线运动虚拟合成一副长度随探测距离成正比的大孔径基阵,通过空间不同位置采样信号的相干叠加得到与距离、频率无关的方位向高分辨率。而多子阵技术的采用更是解决了方位向高分辨和距离向测距之间的矛盾,使SAS技术开始真正走向工程实用化的道路。
成像算法是合成孔径声纳技术的核心技术,目前,国内外已有许多关于合成孔径声纳成像的几何模型,根据这些模型,就可以得到各种不同的成像算法。针对这些几何模型或者成像算法的优劣,国内外一般均依赖于名义分辨率、峰值旁瓣比和积分旁瓣比三个参数进行评价,这往往要历经目标回波仿真、目标成像处理和评估指标计算等步骤,性能评估周期较长,过程繁琐且相对于图形来说还不直观。另外,多子阵合成孔径声纳具有距离向空变特性,也就是说其二维频域相位误差具有距离向瞬时频率、方位向瞬时频率和距离三个自由度。目前,关于多子阵合成孔径声纳二维频域相位误差的分析方法,国内外还未见系统的报道。因此,开展包含距离向瞬时频率、方位向瞬时频率和距离三自由度的二维频域相位误差分析,对于各种成像几何模型或者成像算法性能优劣的分析以及成像算法后期研究的开展具有重要的实际意义。
发明内容
本发明的目的在于克服传统成像模型或成像算法评估周期长、数值指标不够直观等问题,提供一种多子阵合成孔径声纳相位误差分析方法,进而可根据各种模型的误差分析结果,有针对性地开展成像算法方面的研究工作。
为解决上述技术问题,本发明所采用的技术方案为:一种多子阵合成孔径声纳相位误差分析方法,包含以下步骤:
S1、针对每个接收阵元和发射阵元组成的子系统,计算双程斜距历程;
S2、针对每个接收阵元和发射阵元组成的子系统,基于数值分析方法计算相位驻留点;
S3、针对每个接收阵元和发射阵元所组成的子系统,基于数值分析方法的相位驻留点,计算二维频域系统函数;
S4、针对每个接收阵元和发射阵元所组成的子系统,计算待评估模型的二维频域系统函数;
S5、针对每个接收阵元和发射阵元所组成的子系统,根据多子阵合成孔径声纳的方位空不变和距离空变特性,设置目标空间坐标;
S6、针对每个接收阵元和发射阵元所组成的子系统,在二维频域内计算相位误差。
作为本发明的一个优选的技术方案,所述的步骤S2的相位驻留点计算公式如下:
其中v表示声纳载体的运动速度,c表示声波在水中的传播速度;i表示第i个接收阵元和发射阵元所组成子系统的索引,di表示第i个子系统的收、发阵元间距,表示对应于第i个子系统的相位驻留点;r表示目标斜距向的坐标;fτ和ft分别表示对应于距离向快时间τ与方位向慢时间t的瞬时频率。
作为本发明的一个优选的技术方案,所述的步骤S3的二维频域系统函数,其计算公式如下:
其中fc表示宽频带信号的中心频率;P(fτ)表示系统所发射宽频带信号频谱的相位;表示信号所经历的双程斜距历程。
作为本发明的一个优选的技术方案,所述的步骤S5的目标空间坐标的设置,其依据是多子阵合成孔径声纳的方位空不变和距离空变特性,从而在二维平面设置一系列方位坐标固定、距离向坐标空变的理性点目标。
作为本发明的一个优选的技术方案,所述的步骤S6的二维频域相位误差计算公式如下:
其中表示针对第i个子系统所设置的第n个空间理想点目标的距离向坐标;表示待评估成像模型的二维频域系统函数。
本申请实施例中提供的一个或多个技术方案,至少具有如下技术效果或优点:
克服了传统成像模型或成像算法图像聚焦性能评估周期长、数值指标不够直观等问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的二维频域相位误差计算流程。
图2为多子阵合成孔径声纳系统二维成像几何。
图3为收、发阵元最近时的二维频域相位误差。
其中(a)为斜距25米时的相位误差,(b)为斜距80米时的相位误差,(c)为斜距100米时的相位误差,(d)为斜距140米时的相位误差。
图4为收、发阵元最远时的二维频域相位误差。
其中(a)为斜距25米时的相位误差,(b)为斜距80米时的相位误差,(c)为斜距100米时的相位误差,(d)为斜距140米时的相位误差。
具体实施方式
为了更好的理解上述技术方案,下面将结合说明书附图以及具体的实施方式对上述技术方案进行详细的说明。
如图1所示,本实施例所述的一种多子阵合成孔径声纳相位误差分析方法,包含以下步骤:
S1、针对每个接收阵元和发射阵元组成的子系统,计算双程斜距历程;
S2、针对每个接收阵元和发射阵元组成的子系统,基于数值分析方法计算相位驻留点;
S3、针对每个接收阵元和发射阵元所组成的子系统,基于数值分析方法所求得的相位驻留点,计算二维频域系统函数;
S4、针对每个接收阵元和发射阵元所组成的子系统,计算待评估模型的二维频域系统函数;
S5、针对每个接收阵元和发射阵元所组成的子系统,根据多子阵合成孔径声纳的方位空不变和距离空变特性,设置目标空间坐标;
S6、针对每个接收阵元和发射阵元所组成的子系统,在二维频域内计算相位误差。
其中,在本实施例中,所述的步骤S2的相位驻留点计算公式如下:
其中v表示声纳载体的运动速度,c表示声波在水中的传播速度;i表示第i个接收阵元和发射阵元所组成子系统的索引,di表示第i个子系统的收、发阵元间距,表示对应于第i个子系统的相位驻留点;r表示目标斜距向的坐标;fτ和ft分别表示对应于距离向快时间τ与方位向慢时间t的瞬时频率。
作为本发明的一个优选的技术方案,所述的步骤S3的二维频域系统函数,其计算公式如下:
其中fc表示宽频带信号的中心频率;P(fτ)表示系统所发射宽频带信号频谱的相位;表示信号历经的双程斜距历程。
其中,在本实施例中,所述的步骤S5的目标空间坐标的设置,其依据是多子阵合成孔径声纳的方位空不变和距离空变特性,从而设置一系列方位坐标固定、距离向坐标空变的空间理性点目标。
其中,在本实施例中,所述的步骤S6的二维频域相位误差计算公式如下:
其中表示针对第i个子系统所设置的第n个空间理想点目标的距离向坐标;表示待评估成像模型的二维频域系统函数。
更具体的,如图2所示为多子阵合成孔径声纳系统二维成像几何,声纳平台在以速度v前进的过程中,发射阵元同时向正侧方向以固定的脉冲重复频率发射与位置无关的宽频带信号;经过t时间后,发射阵在方位向的位置为vt,此时发射阵元与空间坐标为(0,r)理想点目标之间的距离RT(t;r)为:
多接收阵的采用提高了传统合成孔径声纳系统的测绘速率,如果忽略声波传播时间τ*内平台的运动,那将会导致方位向聚焦性能的下降以及距离向目标走动现象的发生;当波束较窄时,τ*可用2r/c近似。不失一般性,在此仅考虑与发射阵相隔di的第i个接收阵元,那么点目标到接收阵的实际传播路径长度为:
信号所经历的双程斜距历程为:
发射信号的频谱为P(fτ),对于第i个接收阵元和发射阵元所组成的第i个子系统,根据相位驻留原理可得
即
采用数值分析方法中的二分法求解上式,便可得到对应于第i个子系统的相位驻留点依据所求解的相位驻留点计算二维频域系统函数为:
以待评估的Loffeld模型为例,针对第i个子系统,基于Loffeld模型的二维频域系统函数为:
根据多子阵合成孔径声纳的方位空不变和距离空变特性,在二维空间设置方位固定、距离空变的理想点目标,计算对应于第i个子系统距离空变的二维频域相位误差:
下面给出一组典型的合成孔径声纳系统参数:
依据所给出的典型参数,按照上述步骤,针对收、发阵元最近时的第1个子系统,不同斜距情况时的相位误差如图3所示;针对收、发阵元最远时的第32个子系统,不同斜距情况时的相位误差如图4所示,第三个关于距离的自由度可以分别联合图3或者图4中不同距离处的相位误差图得到。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (5)
1.一种多子阵合成孔径声纳相位误差分析方法,其特征在于,包含以下步骤:
S1、针对每个接收阵元和发射阵元组成的子系统,计算双程斜距历程;
S2、针对每个接收阵元和发射阵元组成的子系统,基于数值分析方法计算相位驻留点;
S3、针对每个接收阵元和发射阵元所组成的子系统,基于数值分析方法的相位驻留点,计算二维频域系统函数;
S4、针对每个接收阵元和发射阵元所组成的子系统,计算待评估模型的二维频域系统函数;
S5、针对每个接收阵元和发射阵元所组成的子系统,根据多子阵合成孔径声纳的方位空不变和距离空变特性,设置目标空间坐标;
S6、针对每个接收阵元和发射阵元所组成的子系统,在二维频域内计算相位误差。
2.根据权利要求1所述的一种多子阵合成孔径声纳相位误差分析方法,其特征在于,所述的步骤S2的相位驻留点计算公式如下:
其中v表示声纳载体的运动速度,c表示声波在水中的传播速度;i表示第i个接收阵元和发射阵元所组成子系统的索引,di表示第i个子系统的收、发阵元间距,表示对应于第i个子系统的相位驻留点;r表示目标斜距向的坐标;fτ和ft分别表示对应于距离向快时间τ与方位向慢时间t的瞬时频率。
3.根据权利要求1所述的一种多子阵合成孔径声纳相位误差分析方法,其特征在于,所述的步骤S3的二维频域系统函数,其相位计算公式如下:
其中fc表示宽频带信号的中心频率;P(fτ)表示系统所发射宽频带信号频谱的相位;表示信号所历经的双程斜距历程。
4.根据权利要求1所述的一种多子阵合成孔径声纳相位误差分析方法,其特征在于,所述的步骤S5的目标空间坐标的设置,其依据是多子阵合成孔径声纳的方位空不变和距离空变特性,从而在二维平面设置一系列方位坐标固定、距离向坐标空变的理性点目标。
5.根据权利要求1所述的一种多子阵合成孔径声纳相位误差分析方法,其特征在于,所述的步骤S6的二维频域相位误差计算公式如下:
其中表示针对第i个子系统所设置的第n个空间理想点目标的距离向坐标;表示待评估成像模型的二维频域系统函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611102892.1A CN106772326B (zh) | 2016-12-05 | 2016-12-05 | 一种多子阵合成孔径声纳相位误差分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611102892.1A CN106772326B (zh) | 2016-12-05 | 2016-12-05 | 一种多子阵合成孔径声纳相位误差分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106772326A true CN106772326A (zh) | 2017-05-31 |
CN106772326B CN106772326B (zh) | 2019-04-30 |
Family
ID=58884575
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611102892.1A Active CN106772326B (zh) | 2016-12-05 | 2016-12-05 | 一种多子阵合成孔径声纳相位误差分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106772326B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108037497A (zh) * | 2018-01-04 | 2018-05-15 | 中国人民解放军91388部队 | 多子阵合成孔径声纳数据的收发合置转换方法 |
CN108344987A (zh) * | 2018-01-04 | 2018-07-31 | 中国人民解放军91388部队 | 基于数值计算方法的多子阵合成孔径声纳频域系统函数 |
CN109116364A (zh) * | 2018-08-03 | 2019-01-01 | 中国人民解放军91388部队 | 基于时域顺序排列的多子阵sas回波信号等效单阵转换方法 |
CN116520335A (zh) * | 2023-06-29 | 2023-08-01 | 海底鹰深海科技股份有限公司 | 一种多接收阵元合成孔径声呐波数域成像方法 |
CN116990791A (zh) * | 2023-09-25 | 2023-11-03 | 海底鹰深海科技股份有限公司 | 一种多接收阵合成孔径声纳回波仿真方法 |
CN117828246A (zh) * | 2024-01-10 | 2024-04-05 | 海底鹰深海科技股份有限公司 | 多子阵合成孔径声纳的成像方法和成像系统以及计算装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6037892A (en) * | 1998-05-28 | 2000-03-14 | Multispec Corporation | Method for automatic focusing of radar or sonar imaging systems using high-order measurements |
CN103713289A (zh) * | 2013-12-15 | 2014-04-09 | 浙江大学 | 基于分布式Phased-MIMO联合处理的目标探测方法 |
-
2016
- 2016-12-05 CN CN201611102892.1A patent/CN106772326B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6037892A (en) * | 1998-05-28 | 2000-03-14 | Multispec Corporation | Method for automatic focusing of radar or sonar imaging systems using high-order measurements |
CN103713289A (zh) * | 2013-12-15 | 2014-04-09 | 浙江大学 | 基于分布式Phased-MIMO联合处理的目标探测方法 |
Non-Patent Citations (4)
Title |
---|
HAILIANG YANG ET AL: "A robust multiple-receiver Range-Doppler algorithm for synthetic aperture sonar imagery", 《IEEE》 * |
张学波 等: "合成孔径声呐多接收阵数据融合CS 成像算法", 《哈尔滨工程大学学报》 * |
张学波 等: "四阶模型的多接收阵合成孔径声呐距离-多普勒成像算法", 《电子与信息学报》 * |
张学波 等: "基于双基模型的多接收阵合成孔径声纳CS 成像算法", 《高技术通讯》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108037497A (zh) * | 2018-01-04 | 2018-05-15 | 中国人民解放军91388部队 | 多子阵合成孔径声纳数据的收发合置转换方法 |
CN108344987A (zh) * | 2018-01-04 | 2018-07-31 | 中国人民解放军91388部队 | 基于数值计算方法的多子阵合成孔径声纳频域系统函数 |
CN108037497B (zh) * | 2018-01-04 | 2022-05-17 | 中国人民解放军91388部队 | 多子阵合成孔径声纳数据的收发合置转换方法 |
CN108344987B (zh) * | 2018-01-04 | 2022-06-14 | 中国人民解放军91388部队 | 多子阵合成孔径声纳频域函数基于数值计算的方法 |
CN109116364A (zh) * | 2018-08-03 | 2019-01-01 | 中国人民解放军91388部队 | 基于时域顺序排列的多子阵sas回波信号等效单阵转换方法 |
CN109116364B (zh) * | 2018-08-03 | 2023-01-03 | 中国人民解放军91388部队 | 基于时域顺序排列的多子阵sas回波信号等效单阵转换方法 |
CN116520335A (zh) * | 2023-06-29 | 2023-08-01 | 海底鹰深海科技股份有限公司 | 一种多接收阵元合成孔径声呐波数域成像方法 |
CN116520335B (zh) * | 2023-06-29 | 2023-09-05 | 海底鹰深海科技股份有限公司 | 一种多接收阵元合成孔径声呐波数域成像方法 |
CN116990791A (zh) * | 2023-09-25 | 2023-11-03 | 海底鹰深海科技股份有限公司 | 一种多接收阵合成孔径声纳回波仿真方法 |
CN116990791B (zh) * | 2023-09-25 | 2023-11-28 | 海底鹰深海科技股份有限公司 | 一种多接收阵合成孔径声纳回波仿真方法 |
CN117828246A (zh) * | 2024-01-10 | 2024-04-05 | 海底鹰深海科技股份有限公司 | 多子阵合成孔径声纳的成像方法和成像系统以及计算装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106772326B (zh) | 2019-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106772326A (zh) | 一种多子阵合成孔径声纳相位误差分析方法 | |
CN106842210B (zh) | 一种新的多子阵合成孔径声纳快速成像算法 | |
CN101526614B (zh) | 基于子孔径和等效散射体的sar回波快速模拟方法 | |
CN107976660B (zh) | 弹载多通道雷达超低空目标分析与多径回波建模方法 | |
CN109444901B (zh) | 一种异构环境下多子阵sas子孔径成像方法 | |
CN102749621B (zh) | 一种双基地合成孔径雷达频域成像方法 | |
CN103197317A (zh) | 基于fpga的sar成像方法 | |
CN110836981A (zh) | 一种分层水流高分辨径向声学多普勒频率测量方法 | |
CN115856898B (zh) | 一种适应于全海深的面阵三维成像声呐点位归算方法 | |
CN112540348A (zh) | 一种基于空间尺度上的声线修正算法在长基线水声定位系统上的应用 | |
CN107064911A (zh) | 一种多子阵合成孔径声纳运动误差补偿方法 | |
CN110412587A (zh) | 一种基于解卷积的下视合成孔径三维成像方法及系统 | |
CN104656073A (zh) | 三维成像声纳波束形成方法及在多核处理器上的实现方法 | |
CN104820221B (zh) | 多基合成孔径雷达的目标三维定位方法 | |
CN112433208A (zh) | 一种双站圆周探地雷达快速时域成像方法及系统 | |
CN103487793B (zh) | 一种基于简正波理论的宽带混响波形仿真方法 | |
CN102735207B (zh) | 一种基于距离和载波相差的二维角度解算方法 | |
CN107505604A (zh) | 一种考虑收、发阵元指向性的回波仿真方法 | |
CN114184256A (zh) | 一种多目标背景下的水位测量方法 | |
CN106291529A (zh) | 一种双基地雷达目标定位装置及其定位方法 | |
CN108427111B (zh) | 一种雷达测距方法及装置 | |
CN107942307A (zh) | 调频连续波着陆雷达波束的地面入射角估计方法 | |
CN108761417A (zh) | 基于知识辅助最大似然的机载雷达杂波抑制方法 | |
CN108107415A (zh) | 基于机会约束的集中式mimo雷达多波束功率分配方法 | |
CN101598788B (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 |