CN103809170A - 一种多普勒频偏修正方法及系统 - Google Patents
一种多普勒频偏修正方法及系统 Download PDFInfo
- Publication number
- CN103809170A CN103809170A CN201210451109.8A CN201210451109A CN103809170A CN 103809170 A CN103809170 A CN 103809170A CN 201210451109 A CN201210451109 A CN 201210451109A CN 103809170 A CN103809170 A CN 103809170A
- Authority
- CN
- China
- Prior art keywords
- compensation
- frequency
- array
- signal
- target
- 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
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
本发明涉及一种多普勒频偏修正方法,包括:对目标辐射噪声进行特性分析;构造接收模型,为拖曳线列阵构建直角坐标系;预设频点补偿总数;在某一初始时刻以及拖曳线列阵运动一定时间后的另一时刻,对信号做FFT变换以及频域波束形成,在当前频点上对频域波束形成的输出结果进行累加;计算空间移动补偿与时间延迟补偿,得到总的相位修正补偿,根据前一次补偿的结果与总的相位修正补偿得到扩展虚拟阵列的输出结果;获得时间方位历程图,利用该时间方位历程图对被测目标信号进行检测,判断当前的频点补偿次数是否已经达到了频点补偿总数,如果未达到,改变当前频点值的大小,重新采样,直至当前的频点补偿次数达到预设的频点补偿总数。
Description
技术领域
本发明涉及水声信号处理领域,特别涉及一种多普勒频偏修正方法及系统。
背景技术
随着低噪声舰艇和声隐身技术的发展,舰船等航行器的水下辐射噪声和目标强度进一步下降,传统的被动声呐对目标的探测更加困难,且声呐的检测能力随着安静型目标的出现和迅猛发展而受到严重挑战。
随着科技的发展,基于被动合成孔径技术的拖线阵声呐的研究成为水声信号处理领域的研究热点问题。拖曳线列阵声呐由于可以远离本舰自噪声、可以不受孔径限制,并能变深拖曳利用水文条件,被广泛用于水下目标特性的测量,且由小孔径阵列合成虚拟大孔径阵列的方法为解决低频信号大孔径声呐系统的实现提供了有效途径。被动合成孔径声呐正是通过拖曳线列阵运动合成得到比实际孔径大得多的合成孔径阵列,依靠短阵的机动突破阵列孔径的限制,获得更高的增益和更高的方位分辨力。依据水下舰船等目标特性分析得知螺旋桨的桨叶切割海水产生单频信号分量等线谱成分,为合成一个比物理孔径大得多的有效孔径提供了可能。
拖线阵声呐载体通常以匀速直线运动的方式对目标采取运动接收,由于实际水下复杂环境下介质和路径扰动引起的阵列移动偏差,及海浪、潮汐、风等因素的影响,声呐载体在进行孔径合成时,由于水听器阵列接收时与目标存在相对运动,产生多普勒运动频偏,会对目标的检测性能引起下降和处理增益带来损失,因此解决多普勒频偏是被动合成孔径声呐运动相位补偿的主要问题之一。考虑对弱目标信号检测时,频率未知,孔径合成时合理估计并补偿修正多普勒频偏也是一个难题,对于由于声呐载体与静止/运动目标的多普勒频偏带来的处理增益损失及检测性能的降低等问题,现有技术中并未有任何能够克服上述问题的技术方案。
拖线阵相对于目标以确定性规律运动是被动合成孔径声呐的基础,现有技术的被动合成孔径声呐算法一般用于理想情况下目标信号的检测和分辨,但实际应用时会受到目标的运动及水声信道等复杂因素的影响,而现有的被动合成孔径声呐算法在处理有关问题时忽略多普勒频偏的补偿,因而不具有可行性、通用性及实用性,不能在工程实际中得到合理应用。
发明内容
本发明的目的在于克服实际声呐载体采用匀速直线的运动方式进行孔径合成时,由于水听器接收阵与目标存在的相对运动产生的多普勒频偏,进而在孔径合成时对目标检测性能引起的下降及处理增益带来的损失等问题,从而提供一种对多普勒频率偏移进行补偿的方法及系统。
为了实现上述目的,本发明提供了一种多普勒频偏修正方法,用于包括多个阵元的拖曳线列阵,该方法包括:
步骤1)、对目标辐射噪声进行特性分析;
步骤2)、根据目标辐射噪声特性分析的结果,构造所述拖曳线列阵运动时对被测目标的信号进行声场空间采样的接收模型,然后为所述拖曳线列阵构建直角坐标系;
其中,在构建直角坐标系时,将所述拖曳线列阵中第一个阵元的初始位置作为直角坐标系的原点,将辐射目标作为点声源;
步骤3)、预设频点补偿总数;
步骤4)、某一时刻,所述拖曳线列阵声呐进行采样操作,然后对采样时各个阵元所接收到的信号做FFT变换以及频域波束形成,接着在当前频点上对各个阵元的频域波束形成的输出结果进行累加;
步骤5)、所述拖曳线列阵声呐在运动一定时间后,所述拖曳线列阵声呐进行采样操作,然后对采样时各个阵元所接收到的信号做FFT变换以及频域波束形成,接着在当前频点上对各个阵元的频域波束形成的输出结果进行累加;
步骤6)、根据步骤5)所得到的结果计算空间移动补偿与时间延迟补偿,将空间移动补偿与时间延迟补偿相加,得到所述拖曳线列阵声呐移动前后总的相位修正补偿,然后根据前一次补偿的结果与所述总的相位修正补偿得到扩展虚拟阵列的输出结果;
步骤7)、将步骤6)得到的补偿后的扩展虚拟阵列的输出结果进行相干累加,获得时间方位历程图,利用该时间方位历程图对被测目标信号进行检测,然后判断当前的频点补偿次数是否已经达到了步骤3)中预设的频点补偿总数,如果未达到,改变当前频点值的大小,然后重新执行步骤4)-步骤6),直至当前的频点补偿次数达到预设的频点补偿总数。
上述技术方案中,在所述的步骤4)与步骤5)中,所述FFT变换时所需的点数由补偿精度确定。
上述技术方案中,在所述的步骤7)中,对被测目标信号进行检测包括:
将不同频带内目标波束对应的检测方位进行记录,将各频带内存储的方位进行二次拟合,根据计算结果计算方位估计方差,将计算的方位估计方差与设定的检测方差门限进行比较,若小于门限,则检测到的水声目标信号结果属实,接收到的水下目标信号中确实含有水声目标信号,否则检测结果为虚警,接收到的信号中不含目标信号,继而完成对目标信号的检测。
上述技术方案中,在所述的步骤6)中,所述空间移动补偿根据移动的距离得到。
上述技术方案中,在所述的步骤6)中,所述的时间延迟补偿根据移动相应距离的时间得到。
本发明还提供了一种多普勒频偏修正系统,用于包括多个阵元的拖曳线列阵,该系统包括噪声特性分析模块、预处理模块、频点补偿总数预设模块、信号处理模块、延迟补偿模块、检测模块;其中,
所述的噪声特性分析模块用于对目标辐射噪声进行特性分析;
所述的预处理模块根据目标辐射噪声特性分析的结果,构造所述拖曳线列阵运动时对被测目标的信号进行声场空间采样的接收模型,然后为所述拖曳线列阵构建直角坐标系;其中,在构建直角坐标系时,将所述拖曳线列阵中第一个阵元的初始位置作为直角坐标系的原点,将辐射目标作为点声源;
所述的频点补偿总数预设模块用于预设频点补偿总数;
所述的信号处理模块在某一初始时刻以及拖曳线列阵运动一定时间后的另一时刻分别做如下操作:对拖曳线列阵声呐采样时各个阵元所接收到的信号做FFT变换以及频域波束形成,接着在当前频点上对各个阵元的频域波束形成的输出结果进行累加;
所述的延迟补偿模块根据所述的信号处理模块所得到的结果计算空间移动补偿与时间延迟补偿,将空间移动补偿与时间延迟补偿相加,得到所述拖曳线列阵声呐移动前后总的相位修正补偿,然后根据前一次补偿的结果与所述总的相位修正补偿得到扩展虚拟阵列的输出结果;
所述的检测模块将所述的延迟补偿模块得到的补偿后的扩展虚拟阵列的输出结果进行相干累加,获得时间方位历程图,利用该时间方位历程图对被测目标信号进行检测,然后判断当前的频点补偿次数是否已经达到了所述的频点补偿总数预设模块所预设的频点补偿总数,如果未达到,改变当前频点值的大小,然后重新启动信号处理模块,直至当前的频点补偿次数达到预设的频点补偿总数。
本发明的优点在于:
本发明的方法可以有效避免多普勒频偏对孔径合成带来的处理增益损失和检测性能下降,通过频率搜索、渐次补偿在波束域上进行相位修正,提高了方位估计精度,进而实现了对目标信号的有效检测。
附图说明
图1为一个实施例中利用拖曳线列阵声呐运动下对水下目标的辐射噪声进行检测的示意图;
图2为拖曳线列阵声呐构建直角坐标系的示意图;
图3为本发明的多普勒频偏修正方法的流程图。
具体实施方式
现结合附图对本发明作进一步的描述。
在对本发明的方法做详细说明之前,为了便于理解,首先对本发明中所涉及的一些概念做统一的说明。
时间方位历程图:时间方位历程图通过对连续的采样数据进行短时傅里叶变换而构成。将原始信号的采样序列进行分帧处理(如分为100帧,每帧2000个点),帧间相互重叠50%,然后分别对各帧信号进行波束形成处理,得到各帧对应目标方位,将各帧的时间联合即得到时间方位历程图。采用一张时间方位历程图可以读出目标所在方位,将其方位记录进行后续拟合处理及门限判别,可以实现对目标信号的有效检测。
拖曳线列阵声呐:也称“拖曳阵声呐”,它是将水听器镶嵌在电缆上形成线列阵,由拖曳电缆拖在舰艇尾后水中探测目标的声呐。
时间延迟补偿:反映了拖曳线列阵声呐运动某一时间后,时间变化对不同时间的两次接收时频域波束形成的影响。
空间移动补偿:反映了拖曳线列阵声呐运动某一时间后,空间位置变化对不同时间的两次接收时频域波束形成的影响。
以上是对本发明中所涉及的相关概念的描述,下面对本发明的方法加以说明。
如图3所示,本发明的多普勒频偏修正方法包括以下步骤:
步骤1)、对目标辐射噪声进行特性分析;
步骤2)、根据目标辐射噪声特性分析的结果,构造拖曳线列阵运动时对被测目标的信号进行声场空间采样的接收模型,细化分析各阵元所接收信号,构建直角坐标系;
其中,在构建直角坐标系时,将拖曳线列阵中第一个阵元的初始位置作为直角坐标系的原点,将辐射目标作为点声源。
步骤3)、预设频点补偿总数;
步骤4)、某一时刻,拖曳线列阵声呐进行采样操作,然后对采样时各个阵元所接收到的信号做FFT变换以及频域波束形成,接着在当前频点上对各个阵元的频域波束形成的输出结果进行累加,从而实现信号由阵元域处理转为波束域处理。
其中,所述FFT变换时所涉及的FFT点数与所要求的补偿精度相对应;所述当前频点为所要补偿的频率范围内的一个值,在第一次做频点补偿时,当前频点为所要补偿的频率范围的一个端值,例如,补偿在-0.5Hz~0.5Hz上进行,则当前频点为-0.5Hz,在后续的频点补偿操作中,当前频点发生变化,如由-0.5Hz变为-0.4Hz。
步骤5)、对拖曳线列阵声呐在运动一定时间后,拖曳线列阵声呐进行采样操作,然后对采样时各个阵元所接收到的信号做FFT变换以及频域波束形成,接着在当前频点上对各个阵元的频域波束形成的输出结果进行累加;
其中,所述FFT变换时所涉及的FFT点数与步骤4)中的FFT点数相同。
步骤6)、根据步骤5)所得到的结果计算空间移动补偿与时间延迟补偿,将空间移动补偿与时间延迟补偿相加,得到阵列移动前后总的相位修正补偿,根据前一次补偿的结果与所述总的相位修正补偿得到扩展虚拟阵列的输出结果;
其中,所述空间移动补偿根据移动的距离得到,一般预设为阵元重叠数的一半,或者根据实际移动距离作空间移动补偿;所述的时间延迟补偿根据移动相应距离的时间得到。
步骤7)、将步骤6)得到的补偿后的扩展虚拟阵列的输出结果进行相干累加,获得时间方位历程图,利用该时间方位历程图对被测目标信号进行检测,然后考察当前的频点补偿次数是否已经达到了步骤3)中预设的频点补偿总数,如果未达到,改变当前频点值的大小,然后重新执行步骤4)-步骤6),直至当前的频点补偿次数达到预设的频点补偿总数。
下面对上述步骤做进一步的描述。
在所述的步骤1)中,对目标辐射噪声进行特性分析,这一特性分析的实现不仅为水下螺旋桨切割形成线谱信号的被动合成孔径技术奠定基础,而且为后续信号检测处理中对拖船噪声抵消等技术提供了理论参考。
经研究发现:船体的机械噪声、螺旋桨噪声以及水动力噪声是舰船噪声的三种主要来源,其中机械噪声和螺旋桨噪声是舰船噪声的主要声源。机械噪声是航行或行业舰船上的各种机械振动通过船体向水中辐射而形成的噪声,由于各种机械运动形式的不同,其产生的水下辐射噪声性质也就不同,此噪声可以看成是强线谱和弱连续谱的迭加;螺旋桨噪声是由旋转着的螺旋桨所辐射的噪声,包括螺旋桨空化噪声和螺旋桨叶片振动所产生的噪声,往往是舰船辐射噪声高频段的主要部分;水动力噪声由不规则的、起伏的海流流过运动船只表面而形成,是水流动力作用于舰船的结果,在强度方面一般被机械噪声和螺旋桨噪声所掩盖。
水下舰船的辐射噪声谱特性与频率、航速、深度有关,低航速时,谱的低频端主要包括机械噪声和螺旋桨叶片速率谱线,随着频率增高,该谱线不规则地降低;航速较高时螺旋桨空化噪声的连续谱更为重要,掩盖了很多线谱。舰船辐射噪声信号的总的功率谱可以写成宽带连续谱与线谱的和:
G(t,f)=Gs(f)+GL(f)+2m(t)m(f)Gs(f)
式中,Gs(f)表示平稳连续功率谱,GL(f)表示线谱部分,2m(t)m(f)Gs(f)表示非平稳时变谱,m(t)和m(f)分别为调制函数和调制深度谱。对于舰船辐射噪声而言,在给定航速和深度下,谱的主要成分与临界频率有关,通常的舰船临界频率约在100~1000Hz之间,低于此频率时,谱的主要成分时船的机械和螺旋桨的线谱,高于此频率时,谱的主要成分则是螺旋桨空化的连续噪声谱。依据研究表明,舰船辐射噪声的频域特性是宽带连续谱与窄带线谱的迭加,窄带谱线主要集中在1kHz以下,而宽带连续谱则覆盖了多个倍频程。在实际检测时,被动合成孔径声呐感兴趣的信号是目标辐射噪声的线谱成分。因此在后续的步骤中,拖曳线列阵主要对窄带线谱信号进行检测。
在所述步骤2)中所构建的接收模型包括:由拖缆连接的若干个等间距的各向同性水听器形成一线列阵(在线列阵中的水听器也被称为阵元),该线列阵作匀速直线运动,以某一时刻第1号水听器接收到的被测目标信号为参考,其它水听器相对于被测目标均有不同的距离,即其它水听器接收到的被测目标信号相对于1号水听器接收到的被测目标信号有相应的传播时延,运动一定时间间隔后,各个水听器再次对被测目标信号进行接收处理,根据各水听器与运动目标的距离得到各水听器接收的信号。
在步骤4)与步骤5)中,需要对各个水听器所接收到的信号做FFT变换,在FFT变换时所需的点数与补偿精度有关。如补偿精度为0.1Hz时,FFT点数为实际采样率fs的10倍。当补偿精度发生变化,则相应的FFT变换也会对应改变。
所述的步骤7)具体包括以下步骤:
步骤7-1)、依据虚拟阵列补偿后波束累加获得的输出结果进行处理,结合各时间段信号处理获得时间方位历程图;
步骤7-2)、对整个时间方位历程图输出进行检测,得到目标波束对应方位,具体为:将不同频带内目标波束对应的检测方位进行记录,将各频带内存储的方位进行二次拟合,根据计算结果计算方位估计方差,将计算的方位估计方差与设定的检测方差门限进行比较,若小于门限,则检测到的水声目标信号结果属实,接收到水下目标信号确实含有水声目标信号,否则检测结果为虚警,接收到的信号中不含目标信号,继而完成对目标信号的检测;
步骤7-3)、判断当前的频点补偿次数是否已经达到了步骤3)中预设的频点补偿总数,如果未达到,改变当前频点值的大小,然后重新执行步骤4)-步骤6),直至当前的频点补偿次数达到预设的频点补偿总数。
下面结合一个具体的例子,对本发明方法做进一步说明。
图1为一个实施例中利用拖曳线列阵声呐运动下对水下目标的辐射噪声进行检测的示意图。从图中可以看出,所述拖曳线列阵声呐中有N个水听器,各水听器之间的距离为d,整个拖曳线列阵声呐以速度v做匀速直线运动。当舰船或其他类型的被测目标经过该水平拖曳线列阵声呐时,该声呐中的水听器能够对目标辐射噪声信号进行采样,完成对噪声源的检测。拖曳线列阵声呐在检测时,由于声呐中的水听器与被检测目标之间存在相对运动,使得声呐中的各个水听器所接收的信号与辐射声源存在频率不一致的情况,即:相对声源运动接收频率变高,背离声源运动接收频率变低。这一现象就是多普勒频偏。所述多普勒频偏会在孔径合成时造成目标检测性能的下降。
在下面的描述中就以该拖曳线列阵声呐为例,对如何实现多普勒频偏修正加以说明。
在图1中,R1与RN分别表示拖船与1号、N号水听器的距离。进一步细化分析,构建直角坐标系如图2示,单独考虑1号接收水听器,r0与ri表示移动前后拖船与1号水听器距离,辐射目标视为点声源,所在位置为(x0,y0),则:
即
式中Ts为点目标辐射信号间隔,t1为1号水听器接收信号间隔,由于拖线阵的水平移动引起r0与ri不同,进而导致接收信号采样间隔与辐射信号时间间隔不等,亦即产生多普勒运动频偏。阵列移动前后有阵元相互重叠,通过相位修正,可以得到虚拟阵元的扩展。拖船继续按速度v运动,可连续不断的将空间信息合成为虚拟阵元,当相位修正因子估计精确时,虚拟得到的合成孔径可以等效为实际阵元构成的物理孔径。
将t=0时刻n阵元接收到的信号记为xn(ti),对xn(ti)做M点FFT:
对应第k阵元为X0(k,K),K=1,...,M,做频域波束形成后将各阵元的输出结果进行累加,将阵元所接收到的信号由阵元域转换为波束域:
上述公式中fs表示采样频率。
将拖线列阵运动τ秒后,n阵元所接收到的信号记为xn(ti+τ),对其做M点FFT,对应第k阵元为X1(k,K),K=1,…,M,做频域波束形成后将各阵元的输出结果进行累加,有:
分别对时间延迟及空间位置移动进行相位估计,为保持相位一致,分别对空间位置移动及时间延迟变化予以考虑。对运动τ秒后S1进行空间移动补偿,已知运动距离为v·τ(v·τ可以为各阵元间距d的非整数倍),则空间移动补偿为:
已知运动时间为τ秒,则时间延迟补偿为:
联合考虑空间及时间上变化,即时间延迟及空间位置移动相位修正,将P1和P2相加,得到阵列移动前后总的相位修正补偿:
P0=P1+P2
由S0与P0加权累加得到的对应合成输出结果与角度有关,虚拟阵元长度为vτ,已知初始时N阵元基阵的长度为Nd,在一次合成后,N阵元基阵扩展后所得到的阵元总长度为Nd+vτ。若有多次合成,即再运动Jτ秒后,N阵元基阵扩展为长度为Nd+vJτ的线列阵。根据补偿后的波束输出累加获得的时间方位历程图对被测目标信号进行检测,按预设频点补偿次数逐次单频点变化,完成各频点的相位补偿,重复上述步骤,直至实现对目标信号的检测。
本发明还公开了一种多普勒频偏修正系统,用于包括多个阵元的拖曳线列阵,该系统包括噪声特性分析模块、预处理模块、频点补偿总数预设模块、信号处理模块、延迟补偿模块、检测模块;其中,
所述的噪声特性分析模块用于对目标辐射噪声进行特性分析;
所述的预处理模块根据目标辐射噪声特性分析的结果,构造所述拖曳线列阵运动时对被测目标的信号进行声场空间采样的接收模型,然后为所述拖曳线列阵构建直角坐标系;其中,在构建直角坐标系时,将所述拖曳线列阵中第一个阵元的初始位置作为直角坐标系的原点,将辐射目标作为点声源;
所述的频点补偿总数预设模块用于预设频点补偿总数;
所述的信号处理模块在某一初始时刻以及拖曳线列阵运动一定时间后的另一时刻分别做如下操作:对拖曳线列阵声呐采样时各个阵元所接收到的信号做FFT变换以及频域波束形成,接着在当前频点上对各个阵元的频域波束形成的输出结果进行累加;
所述的延迟补偿模块根据所述的信号处理模块所得到的结果计算空间移动补偿与时间延迟补偿,将空间移动补偿与时间延迟补偿相加,得到所述拖曳线列阵声呐移动前后总的相位修正补偿,然后根据前一次补偿的结果与所述总的相位修正补偿得到扩展虚拟阵列的输出结果;
所述的检测模块将所述的延迟补偿模块得到的补偿后的扩展虚拟阵列的输出结果进行相干累加,获得时间方位历程图,利用该时间方位历程图对被测目标信号进行检测,然后判断当前的频点补偿次数是否已经达到了所述的频点补偿总数预设模块所预设的频点补偿总数,如果未达到,改变当前频点值的大小,然后重新启动信号处理模块,直至当前的频点补偿次数达到预设的频点补偿总数。
本发明的方法及系统具有以下优点:
一方面,本发明的方法及系统不需要知道信号源的方位角、确切频率和准确的拖曳阵速度,通过相位估计修正能够得到虚拟阵列接收数据,采用常规波束形成等技术对虚拟阵列的输出进行相关累加可以得到大孔径阵列的波束输出,通过时间延迟及空间位置移动可获得相位修正补偿因子,被动合成孔径算法应用于信号检测不受孔径位置的约束限制,在波束域代替其在阵元域上进行相位修正,直接用于相位修正补偿;
另一方面,本发明的方法及系统充分考虑实际水下复杂环境下介质和路径扰动引起的阵列移动偏差,及海浪、潮汐、风等因素的影响,结合频率搜索,渐次补偿修正的方法解决多普勒频偏带来的影响,不要求阵元在相邻两次移动是位置完全重合;
最后,本发明的方法及系统在复杂水声环境下使被动合成孔径声呐对水下弱目标信号检测能力更优异,且检测效果较好,速度较快,即应用较短的基阵获得几倍到几十倍的实际物理孔径阵的增益和分辨力,综合提高对水下远距离、安静型弱目标信号的检测能力,具有一定的工程实用性。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种多普勒频偏修正方法,用于包括多个阵元的拖曳线列阵,该方法包括:
步骤1)、对目标辐射噪声进行特性分析;
步骤2)、根据目标辐射噪声特性分析的结果,构造所述拖曳线列阵运动时对被测目标的信号进行声场空间采样的接收模型,然后为所述拖曳线列阵构建直角坐标系;
其中,在构建直角坐标系时,将所述拖曳线列阵中第一个阵元的初始位置作为直角坐标系的原点,将辐射目标作为点声源;
步骤3)、预设频点补偿总数;
步骤4)、某一时刻,所述拖曳线列阵声呐进行采样操作,然后对采样时各个阵元所接收到的信号做FFT变换以及频域波束形成,接着在当前频点上对各个阵元的频域波束形成的输出结果进行累加;
步骤5)、所述拖曳线列阵声呐在运动一定时间后,所述拖曳线列阵声呐进行采样操作,然后对采样时各个阵元所接收到的信号做FFT变换以及频域波束形成,接着在当前频点上对各个阵元的频域波束形成的输出结果进行累加;
步骤6)、根据步骤5)所得到的结果计算空间移动补偿与时间延迟补偿,将空间移动补偿与时间延迟补偿相加,得到所述拖曳线列阵声呐移动前后总的相位修正补偿,然后根据前一次补偿的结果与所述总的相位修正补偿得到扩展虚拟阵列的输出结果;
步骤7)、将步骤6)得到的补偿后的扩展虚拟阵列的输出结果进行相干累加,获得时间方位历程图,利用该时间方位历程图对被测目标信号进行检测,然后判断当前的频点补偿次数是否已经达到了步骤3)中预设的频点补偿总数,如果未达到,改变当前频点值的大小,然后重新执行步骤4)-步骤6),直至当前的频点补偿次数达到预设的频点补偿总数。
2.根据权利要求1所述的多普勒频偏修正方法,其特征在于,在所述的步骤4)与步骤5)中,所述FFT变换时所需的点数由补偿精度确定。
3.根据权利要求1所述的多普勒频偏修正方法,其特征在于,在所述的步骤7)中,对被测目标信号进行检测包括:
将不同频带内目标波束对应的检测方位进行记录,将各频带内存储的方位进行二次拟合,根据计算结果计算方位估计方差,将计算的方位估计方差与设定的检测方差门限进行比较,若小于门限,则检测到的水声目标信号结果属实,接收到的水下目标信号中确实含有水声目标信号,否则检测结果为虚警,接收到的信号中不含目标信号,继而完成对目标信号的检测。
4.根据权利要求1所述的多普勒频偏修正方法,其特征在于,在所述的步骤6)中,所述空间移动补偿根据移动的距离得到。
5.根据权利要求1所述的多普勒频偏修正方法,其特征在于,在所述的步骤6)中,所述的时间延迟补偿根据移动相应距离的时间得到。
6.一种多普勒频偏修正系统,其特征在于,用于包括多个阵元的拖曳线列阵,该系统包括噪声特性分析模块、预处理模块、频点补偿总数预设模块、信号处理模块、延迟补偿模块、检测模块;其中,
所述的噪声特性分析模块用于对目标辐射噪声进行特性分析;
所述的预处理模块根据目标辐射噪声特性分析的结果,构造所述拖曳线列阵运动时对被测目标的信号进行声场空间采样的接收模型,然后为所述拖曳线列阵构建直角坐标系;其中,在构建直角坐标系时,将所述拖曳线列阵中第一个阵元的初始位置作为直角坐标系的原点,将辐射目标作为点声源;
所述的频点补偿总数预设模块用于预设频点补偿总数;
所述的信号处理模块在某一初始时刻以及拖曳线列阵运动一定时间后的另一时刻分别做如下操作:对拖曳线列阵声呐采样时各个阵元所接收到的信号做FFT变换以及频域波束形成,接着在当前频点上对各个阵元的频域波束形成的输出结果进行累加;
所述的延迟补偿模块根据所述的信号处理模块所得到的结果计算空间移动补偿与时间延迟补偿,将空间移动补偿与时间延迟补偿相加,得到所述拖曳线列阵声呐移动前后总的相位修正补偿,然后根据前一次补偿的结果与所述总的相位修正补偿得到扩展虚拟阵列的输出结果;
所述的检测模块将所述的延迟补偿模块得到的补偿后的扩展虚拟阵列的输出结果进行相干累加,获得时间方位历程图,利用该时间方位历程图对被测目标信号进行检测,然后判断当前的频点补偿次数是否已经达到了所述的频点补偿总数预设模块所预设的频点补偿总数,如果未达到,改变当前频点值的大小,然后重新启动信号处理模块,直至当前的频点补偿次数达到预设的频点补偿总数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210451109.8A CN103809170B (zh) | 2012-11-12 | 2012-11-12 | 一种多普勒频偏修正方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210451109.8A CN103809170B (zh) | 2012-11-12 | 2012-11-12 | 一种多普勒频偏修正方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103809170A true CN103809170A (zh) | 2014-05-21 |
CN103809170B CN103809170B (zh) | 2016-08-03 |
Family
ID=50706222
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210451109.8A Expired - Fee Related CN103809170B (zh) | 2012-11-12 | 2012-11-12 | 一种多普勒频偏修正方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103809170B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107462317A (zh) * | 2016-06-02 | 2017-12-12 | 中国科学院声学研究所 | 一种复杂场景下的水面目标辐射噪声模拟方法 |
CN110068331A (zh) * | 2018-01-24 | 2019-07-30 | 北京致感致联科技有限公司 | 水下导航定位设备及系统 |
CN111024217A (zh) * | 2019-12-27 | 2020-04-17 | 武昌船舶重工集团有限公司 | 一种船舶水下自噪声监测方法 |
CN111654338A (zh) * | 2020-05-07 | 2020-09-11 | Oppo广东移动通信有限公司 | 一种多普勒参数估计方法、设备及存储介质 |
CN113281727A (zh) * | 2021-06-02 | 2021-08-20 | 中国科学院声学研究所 | 一种基于水平线列阵的输出增强的波束形成方法及其系统 |
CN113589298A (zh) * | 2021-06-16 | 2021-11-02 | 中国船舶重工集团公司第七一五研究所 | 一种拖线阵湍流边界层噪声测试数据线谱净化方法 |
CN114460587A (zh) * | 2022-04-11 | 2022-05-10 | 杭州应用声学研究所(中国船舶重工集团公司第七一五研究所) | 一种主动声呐全景接触目标快速辨识方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5047993A (en) * | 1965-06-08 | 1991-09-10 | The United States Of America As Represented By The Secretary Of The Navy | Detection and classification for multi-beam sonar systems |
CN102565775B (zh) * | 2012-01-02 | 2013-06-05 | 浙江大学 | 基于fpga的实时并行时延相移波束形成方法 |
-
2012
- 2012-11-12 CN CN201210451109.8A patent/CN103809170B/zh not_active Expired - Fee Related
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107462317A (zh) * | 2016-06-02 | 2017-12-12 | 中国科学院声学研究所 | 一种复杂场景下的水面目标辐射噪声模拟方法 |
CN107462317B (zh) * | 2016-06-02 | 2019-07-09 | 中国科学院声学研究所 | 一种复杂场景下的水面目标辐射噪声模拟方法 |
CN110068331A (zh) * | 2018-01-24 | 2019-07-30 | 北京致感致联科技有限公司 | 水下导航定位设备及系统 |
CN111024217B (zh) * | 2019-12-27 | 2021-08-24 | 武昌船舶重工集团有限公司 | 一种船舶水下自噪声监测方法 |
CN111024217A (zh) * | 2019-12-27 | 2020-04-17 | 武昌船舶重工集团有限公司 | 一种船舶水下自噪声监测方法 |
CN111654338A (zh) * | 2020-05-07 | 2020-09-11 | Oppo广东移动通信有限公司 | 一种多普勒参数估计方法、设备及存储介质 |
CN111654338B (zh) * | 2020-05-07 | 2022-07-19 | Oppo广东移动通信有限公司 | 一种多普勒参数估计方法、设备及存储介质 |
CN113281727A (zh) * | 2021-06-02 | 2021-08-20 | 中国科学院声学研究所 | 一种基于水平线列阵的输出增强的波束形成方法及其系统 |
CN113281727B (zh) * | 2021-06-02 | 2021-12-07 | 中国科学院声学研究所 | 一种基于水平线列阵的输出增强的波束形成方法及其系统 |
CN113589298A (zh) * | 2021-06-16 | 2021-11-02 | 中国船舶重工集团公司第七一五研究所 | 一种拖线阵湍流边界层噪声测试数据线谱净化方法 |
CN113589298B (zh) * | 2021-06-16 | 2023-08-11 | 中国船舶重工集团公司第七一五研究所 | 一种拖线阵湍流边界层噪声测试数据线谱净化方法 |
CN114460587A (zh) * | 2022-04-11 | 2022-05-10 | 杭州应用声学研究所(中国船舶重工集团公司第七一五研究所) | 一种主动声呐全景接触目标快速辨识方法 |
CN114460587B (zh) * | 2022-04-11 | 2022-07-19 | 杭州应用声学研究所(中国船舶重工集团公司第七一五研究所) | 一种主动声呐全景接触目标快速辨识方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103809170B (zh) | 2016-08-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103529441B (zh) | 一种被动合成孔径目标信号检测和分辨方法及系统 | |
CN104101871A (zh) | 一种用于被动合成孔径的抑制窄带干扰方法及系统 | |
CN103809170B (zh) | 一种多普勒频偏修正方法及系统 | |
CN103675819B (zh) | 一种可用于被动合成孔径阵列偏航的目标检测方法及系统 | |
CN102981146B (zh) | 一种单矢量水听器被动定位方法 | |
CN103438987B (zh) | 基于超指向性小孔径圆柱阵的舰船辐射噪声源分辨方法 | |
Turgut et al. | Broadband source localization using horizontal-beam acoustic intensity striations | |
CN101813772B (zh) | 一种快速宽带频域扩展拖曳阵波束形成方法 | |
CN104678384B (zh) | 一种波束域的声压差互相关谱分析水下目标速度估计方法 | |
KR100195576B1 (ko) | 이동 물체 속도 측정 장치 | |
Wang et al. | A method of velocity estimation using composite hyperbolic frequency-modulated signals in active sonar | |
Morgado et al. | Design and experimental evaluation of an integrated USBL/INS system for AUVs | |
CN101915922A (zh) | 拖曳线列阵被动测距方法 | |
CN103076594A (zh) | 一种基于互相关的水声脉冲信号双阵元定位的方法 | |
RU2626380C1 (ru) | Система селекции движущихся целей с измерением дальности, радиальной скорости и направления движения | |
Kim et al. | Temporal domain processing for a synthetic aperture array | |
Urazghildiiev et al. | Maximum likelihood estimators and Cramér–Rao bound for estimating azimuth and elevation angles using compact arrays | |
Boltryk et al. | An ultrasonic transducer array for velocity measurement in underwater vehicles | |
CN104111349A (zh) | 一种多普勒测速器 | |
Lo | A matched-field processing approach to ranging surface vessels using a single hydrophone and measured replica fields | |
Yang et al. | Bayesian passive acoustic tracking of a cooperative moving source in shallow water | |
Song et al. | Track prediction of moving target for interference fringes without nearby points | |
Chen et al. | Passive range localization of the acoustic moving source using the demon spectrum | |
Zhu et al. | A coastal acoustic tomography experiment in the Sanmen Bay, China | |
Kang et al. | Measurement of Tonal Noise of Underwater Vehicle using Virtual Synthetic Array in Near-field |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160803 Termination date: 20181112 |
|
CF01 | Termination of patent right due to non-payment of annual fee |