CN106772326B - 一种多子阵合成孔径声纳相位误差分析方法 - Google Patents

一种多子阵合成孔径声纳相位误差分析方法 Download PDF

Info

Publication number
CN106772326B
CN106772326B CN201611102892.1A CN201611102892A CN106772326B CN 106772326 B CN106772326 B CN 106772326B CN 201611102892 A CN201611102892 A CN 201611102892A CN 106772326 B CN106772326 B CN 106772326B
Authority
CN
China
Prior art keywords
array element
subsystem
phase
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.)
Active
Application number
CN201611102892.1A
Other languages
English (en)
Other versions
CN106772326A (zh
Inventor
张学波
肖军
陈晓辉
夏雪
彭成
杨博
陈迎春
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chinese People's Liberation Army 91388
Original Assignee
Chinese People's Liberation Army 91388
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chinese People's Liberation Army 91388 filed Critical Chinese People's Liberation Army 91388
Priority to CN201611102892.1A priority Critical patent/CN106772326B/zh
Publication of CN106772326A publication Critical patent/CN106772326A/zh
Application granted granted Critical
Publication of CN106772326B publication Critical patent/CN106772326B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

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 (2)

1.一种多子阵合成孔径声纳相位误差分析方法，其特征在于，包含以下步骤：
S1、针对每个接收阵元和发射阵元组成的子系统，计算双程斜距历程；
S2、针对每个接收阵元和发射阵元组成的子系统，基于数值分析方法计算相位驻留点，该相位驻留点计算公式如下：
其中v表示声纳载体的运动速度，c表示声波在水中的传播速度；i表示第i个接收阵元和发射阵元所组成子系统的索引，di表示第i个子系统的收、发阵元间距，表示对应于第i个子系统的相位驻留点；r表示目标斜距向的坐标；fτ和ft分别表示对应于距离向快时间τ与方位向慢时间t的瞬时频率；
S3、针对每个接收阵元和发射阵元所组成的子系统，基于数值分析方法的相位驻留点，计算二维频域系统函数，其相位计算公式如下：
其中fc表示宽频带信号的中心频率；P(fτ)表示系统所发射宽频带信号频谱的相位；表示信号所历经的双程斜距历程；
S4、针对每个接收阵元和发射阵元所组成的子系统，计算待评估模型的二维频域系统函数；
S5、针对每个接收阵元和发射阵元所组成的子系统，根据多子阵合成孔径声纳的方位空不变和距离空变特性，设置目标空间坐标；
S6、针对每个接收阵元和发射阵元所组成的子系统，在二维频域内计算相位误差，其二维频域相位误差计算公式如下：
其中表示针对第i个子系统所设置的第n个空间理想点目标的距离向坐标；表示待评估成像模型的二维频域系统函数。
2.根据权利要求1所述的一种多子阵合成孔径声纳相位误差分析方法，其特征在于，所述的步骤S5的目标空间坐标的设置，其依据是多子阵合成孔径声纳的方位空不变和距离空变特性，从而在二维平面设置一系列方位坐标固定、距离向坐标空变的理性点目标。
CN201611102892.1A 2016-12-05 2016-12-05 一种多子阵合成孔径声纳相位误差分析方法 Active CN106772326B (zh)

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 CN106772326A (zh) 2017-05-31
CN106772326B true 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)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108344987B (zh) * 2018-01-04 2022-06-14 中国人民解放军91388部队 多子阵合成孔径声纳频域函数基于数值计算的方法
CN108037497B (zh) * 2018-01-04 2022-05-17 中国人民解放军91388部队 多子阵合成孔径声纳数据的收发合置转换方法
CN109116364B (zh) * 2018-08-03 2023-01-03 中国人民解放军91388部队 基于时域顺序排列的多子阵sas回波信号等效单阵转换方法
CN116520335B (zh) * 2023-06-29 2023-09-05 海底鹰深海科技股份有限公司 一种多接收阵元合成孔径声呐波数域成像方法
CN116990791B (zh) * 2023-09-25 2023-11-28 海底鹰深海科技股份有限公司 一种多接收阵合成孔径声纳回波仿真方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
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
CN103713289B (zh) * 2013-12-15 2016-03-23 浙江大学 基于分布式Phased-MIMO联合处理的目标探测方法

Also Published As

Publication number Publication date
CN106772326A (zh) 2017-05-31

Similar Documents

Publication Publication Date Title
CN106772326B (zh) 一种多子阵合成孔径声纳相位误差分析方法
CN106842210B (zh) 一种新的多子阵合成孔径声纳快速成像算法
CN107976660B (zh) 弹载多通道雷达超低空目标分析与多径回波建模方法
CN103713285B (zh) 基于信息融合的分布式米波阵列雷达测角方法
CN107688178A (zh) 一种基于77GHz毫米波雷达的锯齿波测距测速方法
CN103616671A (zh) 一种相控阵雷达数字仿真系统及其仿真方法
CN110412571A (zh) 基于电磁涡旋波的合成孔径雷达三维成像方法
CN110286370B (zh) 航行状态下的船载地波雷达运动补偿方法
CN105974390B (zh) 基于多普勒信息的机扫米波雷达质量中心测角方法
CN108490443A (zh) 基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法
JP2014182010A (ja) レーダ装置
CN103792512A (zh) 一种基于水平直线阵的近场三维被动定位方法
CN106842162A (zh) 一种有源雷达对消隐身系统及方法
CN104991252A (zh) 双站圆迹sar快速时域成像方法
CN115856898B (zh) 一种适应于全海深的面阵三维成像声呐点位归算方法
CN106291529B (zh) 一种双基地雷达目标定位装置及其定位方法
JP5238531B2 (ja) レーダ装置、海洋レーダ観測装置およびドップラ周波数データ算出方法
CN113671443A (zh) 基于掠射角声线修正的水声传感器网络深海目标定位方法
CN107688176B (zh) 一种隐身目标探测装置及方法
CN103487793B (zh) 一种基于简正波理论的宽带混响波形仿真方法
CN102621557B (zh) 一种动态目标二维捕获的方法
CN108427111B (zh) 一种雷达测距方法及装置
CN103487796A (zh) 一种利用水声信道统计不变特征实现被动测距的方法
CN101598788B (zh) 合成孔径声纳信号的快速仿真方法
CN103969642B (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