CN104180891B - 一种基于声阵的海洋声传播损失实时测量方法 - Google Patents
一种基于声阵的海洋声传播损失实时测量方法 Download PDFInfo
- Publication number
- CN104180891B CN104180891B CN201410422855.3A CN201410422855A CN104180891B CN 104180891 B CN104180891 B CN 104180891B CN 201410422855 A CN201410422855 A CN 201410422855A CN 104180891 B CN104180891 B CN 104180891B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mfrac
- mtable
- 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
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明提出一种基于声阵的海洋声传播损失测量方法,包括利用水声接收阵系统实时实现传播损失测量,并且在阵元域的信噪比小于6dB情况下,仍可通过声接收阵的波束域输出结果,实时实现不同频段水下声传播损失的计算。与传统的基于浮标的静态传播损失测量技术相比,基于声阵的海洋声传播损失实时测量技术可以实时进行阵元域和波束域的多频段的水下声传播损失测量。并且在阵元域信噪比低于6dB而无法进行传播损失测量的情况下,利用声基阵的阵增益,通过波束域可以实现多频段的水下声传播损失测量。
Description
技术领域
本发明属于水声信号处理技术领域,具体涉及一种海洋声传播损失的实时测量方法。
背景技术
随着水下目标隐身技术的发展,作为探测设备的主/被动声纳,其低频、大孔径、大功率特征愈加明显。主/被动拖曳线列阵声纳已成为最主要的声纳探测设备,它具有远离本艇干扰、深度可变、低频大孔径、主动声源级高等优点,已成为探测安静型潜艇的最重要设备。为充分发挥声纳设备,尤其是主被动拖曳线列阵声纳在各种海区的使用效能,获取工作海区的海洋声学环境十分必要。现代海洋斗争的特点,已由冷战时期的深海对抗转为既兼顾深海同时又兼顾近岸复杂浅水域,对声纳的探测及目标参数估计提出了更高的性能要求。因此,对复杂多变海洋声学环境的充分了解,直接关系到声纳的性能预报、最佳效能发挥、设计参数优化。
水声信号的传播衰减直接关系声纳作用距离,传播衰减与工作频段、声速剖面、介质/界面吸收/反射等因素密切相关,测量并获取传播衰减曲线,可以为声纳性能的预报和设计提供依据,将实测曲线与模型估计曲线同步比对,可提高声纳性能预报的准确度。
发明内容
本发明所要解决的技术问题是:利用水声接收阵系统实时实现传播损失测量,在阵元域的信噪比小于6dB情况下,利用声阵的阵增益,通过声接收阵的波束域输出结果,实时实现不同频段水下声传播损失的计算。
为解决以上技术问题,本发明是通过以下技术方案实现的:
一种基于声阵的海洋声传播损失实时测量方法,其特征在于:采用水平拖曳线列阵实现,所述水平拖曳线列阵系统由16段声学段组成,各声学段由其内置的水听器阵列、前置与传输部分一起,完成声信号的声电转换、前置预处理、信号采集和传输功能;阵元数为184元,为实现波束域不同频段传播损失计算,采取布阵方式是:前64个声通道,阵元间距1.5m,中间24个声通道,阵元间距24m,最后96个声通道,阵元间距0.4m;声基阵总长758.4m;
用于单阵元传播损失实时测量,具体包括以下步骤:
步骤一:对实时接收的数据进行数据累计,累积8批数据得到x1(i,n)和x1_1(i,n),x1(i,n)和x1_1(i,n)分别为累积8批得到的当前测量数据和之前数据;
步骤二:对接收的数据进行傅立叶变化处理:
其中,16384为傅立叶变换的长度,处理数据有87.5%的重叠,即每次更新2048点,频率分辨力为0.7324Hz;
步骤三:分频段进行能量计算:
步骤四:对数据进行能量判断,以确定有无信号存在;在计算传播损失时,如果是脉冲信号,需进行信噪比判断:如果Est(fb)-Ent(fb)>Th,Th为门限,则输出传播损失;否则不输出;如果是连续波,则不需要进行能量判断;
步骤五:进行传播损失计算:
TL(fb)=SL(fb)-Est(fb)-HS-AG (5)
其中SL为声源级,事先标定;HS表示水听器灵敏度,单位为dB;AG表示放大倍数,单位为dB。
一种基于声阵的海洋声传播损失实时测量方法,其特征在于:
用于波束域传播损失实时测量,利用64元细线阵或96元粗线阵输出信号测量阵波束传播损失;接收数据表示为x1_64(i,n),i=0,……63;x1_96(i,n),i=0,……95,其中,i表示水听器通道号,x1_64(i,n)表示前64个水听器接收的阵数据,x1_96(i,n)表示后96个水听器接收的阵数据;n=0,……,2048为时间采样序列;具体包括以下步骤:
步骤一:对实时接收的数据进行滤波:
其中h1(n)为带通滤波器,滤波器通带为50-500Hz,过渡带为[10 50Hz]和[500600Hz],通带波纹小于等于0.5dB,阻带衰减大于等于40dB;h2(n)为带通滤波器,滤波器通带为500-2000Hz,过渡带为[400 500Hz]和[2000 2100Hz],通带波纹小于等于0.5dB,阻带衰减大于等于40dB;
步骤二:对接收的数据进行傅立叶变换处理:
步骤三:分别对两个阵的接收数据进行波束形成处理,在余弦坐标系中覆盖0°~180°水平空间:
其中0≤l≤128,k=17,...,85;
其中0≤l≤128,k=85,...,341;数据重叠50%,基元间距d1=1.5m,d2=0.4m,c为声速;
步骤四:波束选取;根据显控输入的人工选取方位或自动选取的方位lmax选择信号波束;自动选取方位过程如下:
公式(12)表示对100~500Hz频段信号计算波束图;
公式(13)表示对500~2000Hz频段信号计算波束图;利用公式(14)判断选取的方位;
步骤五:对选定的lmax波束数据进行逆傅立叶变换处理:
其中,2048为逆傅立叶变换的长度;
步骤六:对人工选取方位或自动选取的方位数据累积8批,得到的Bftime1_64(n)和Bftime1_96(n),分别为指定方位的累积8批得到的当前测量数据,Bftime1_1_64(n)和Bftime1_1_96(n)分别为指定方位的累积8批得到的之前噪声数据;
步骤七:对目标方位数据累积进行傅立叶变换处理:
其中,16384为傅立叶变换的长度,处理数据有87.5%的重叠,即每次更新2048点;频率分辨力为0.7324Hz;
步骤八:对波束形成并积累后的数据进行分频段能量计算:
步骤九:对接收信号进行能量判断;如果是脉冲信号需进行信噪比判断:如果Est_64(fb)-Ent_64(fb)>Th,Est_96(fb)-Ent_96(fb)>Th,Th为门限,则输出传播损失;否则,不输出;如果是连续波信号不需要进行能量判断;
步骤十:传播损失计算:
TL_64(fb)=SL(fb)-Est_64(fb)-HS-AG (25)
TL_96(fb)=SL(fb)-Est_96(fb)-HS-AG (26)
其中SL为配合舰船发射声源或爆炸声源的声源级,事先标定;HS表示水听器灵敏度,单位为dB;AG表示放大倍数,单位为dB。
本发明能带来以下有益效果:
1、基于声阵的海洋声传播损失实时测量技术可以实时进行阵元域和波束域的多频段的水下声传播损失测量。
2、基于声阵的海洋声传播损失实时测量技术在阵元域信噪比低于6dB而无法进行传播损失测量的情况下,利用声基阵的阵增益,通过波束域可以实现多频段的水下声传播损失测量。
附图说明
图1本发明的系统原理框图;
图2海上试验630-800Hz传播损失实时测量结果;
图3海上试验1000-1250Hz传播损失实时测量结果;
图4海上试验1600-2000Hz传播损失实时测量结果。
具体实施方式
下面结合具体实施例子和附图对本发明做进一步的描述:
本发明采用水平拖曳线列阵实现,由16段声学段组成,各声学段由其内置的水听器阵列、前置与传输部分一起,完成声信号的声电转换、前置预处理、信号采集和传输功能;阵元数为184元,为实现波束域不同频段传播损失计算,采取布阵方式是:前64个声通道,阵元间距1.5m,中间24个声通道,阵元间距24m,最后96个声通道,阵元间距0.4m;声基阵总长758.4m。
图1是本发明的系统原理框图,从图中可以看出本发明的实施过程:对水平拖曳线列阵的接收阵数据,首先进行放大、滤波、进行A/D转换,将接收的模拟信号进行数字采样。然后分别进行阵元域和波束域的传播损失计算,最后根据现场测量情况,操作者可以选择进行阵元域或波束域的传播损失测量结果输出。
1、单阵元传播损失实时测量方法
步骤一:对实时接收的数据进行数据累计,累积8批数据得到x1(i,n)和x1_1(i,n);
x1(i,n)和x1_1(i,n)分别为累积8批得到的当前测量数据和之前数据;
步骤二:对接收的数据进行傅立叶变化处理:
Nfft=16384为FFT长度,处理数据有87.5%的重叠,即每次更新2048点,频率分辨力为0.7324Hz;
步骤三:分频段进行能量计算:
步骤四:对数据进行能量判断,以确定有无信号存在;在计算传播损失时,如果是脉冲信号,需进行信噪比判断:如果Est(i,kn)-Ent(i,kn)>Th,Th为门限,则输出传播损失;否则不输出;如果是连续波,则不需要进行能量判断;
步骤五:进行传播损失计算:
TL(fb)=SL(fb)-Est(fb)-HS-AG (5)
其中SL为声源级,事先标定;HS表示水听器灵敏度,单位为dB;AG表示放大倍数,单位为dB。
2、波束域传播损失实时测量方法
利用64元细线阵或96元粗线阵输出信号测量阵波束传播损失;接收数据表示为x1_64(i,n),i=0,……64;x1_96(i,n)为水听器通道;n=0,……,2048为时间采样序列。
步骤一:对实时接收的数据进行滤波:
其中h1(n)为带通滤波器,滤波器通带为50-500Hz,过渡带为[10 50Hz]和[500600Hz],通带波纹小于等于0.5dB,阻带衰减大于等于40dB;h2(n)为带通滤波器,滤波器通带为500-2000Hz,过渡带为[400 500Hz]和[2000 2100Hz],通带波纹小于等于0.5dB,阻带衰减大于等于40dB;
步骤二:对接收的数据进行傅立叶变换处理:
步骤三:分别对两个阵的接收数据进行波束形成处理,在余弦坐标系中覆盖0°~180°水平空间:
其中0≤l≤128,k=17,...,85;
其中0≤l≤128,k=85,...,341;数据重叠50%,基元间距d1=1.5m,d2=0.4m,c为声速;
步骤四:波束选取;根据显控输入的人工选取方位或自动选取的方位lmax选择信号波束;自动选取方位过程如下:
公式(12)表示对100~500Hz频段信号计算波束图;
公式(13)表示对500~2000Hz频段信号计算波束图;利用公式(14)判断选取的方位;
步骤五:对选定的lmax波束数据进行逆傅立叶变换处理:
其中Nfft=2048为逆傅立叶变换长度;
步骤六:对自动或手动输入的目标方位数据累积8批,得到的Bftime1_64(n)和Bftime1_96(n),分别为指定方位的累积8批得到的当前测量数据,Bftime1_1_64(n)和Bftime1_1_96(n)分别为指定方位的累积8批得到的之前噪声数据;
步骤七:对目标方位数据累积进行傅立叶变换处理:
Nfft=16384为傅立叶变换长度,处理数据有87.5%的重叠,即每次更新2048点;频率分辨力为0.7324Hz;
步骤八:对波束形成并积累后的数据进行分频段能量计算:
步骤九:对接收信号进行能量判断;如果是脉冲信号需进行信噪比判断:如果Est_64(fb)-Ent_64(fb)>Th,Est_96(fb)-Ent_96(fb)>Th,Th为门限,则输出传播损失;否则,不输出;如果是连续波信号不需要进行能量判断;
步骤十:传播损失计算:
TL_64(fb)=SL(fb)-Est_64(fb)-HS-AG (25)
TL_96(fb)=SL(fb)-Est_96(fb)-HS-AG (26)
其中SL为配合舰船发射声源或爆炸声源的声源级,事先标定;HS表示水听器灵敏度,单位为dB;AG表示放大倍数,单位为dB。
图2是海上试验630-800Hz传播损失实时测量结果。
图3是海上试验1000-1250Hz传播损失实时测量结果。
图4是海上试验1600-2000Hz传播损失实时测量结果。
水平拖曳线列阵的阵型由于海洋环境因素的影响,会造成阵型失真,对波束域的传播损失测量会造成一定影响,但实际应用中,当拖船拖速满足测量要求时,可以忽略阵型失真的影响。而其他的阵型,如圆柱阵、舷侧阵等,则不会出现阵型失真的影响。
3、显示/数据记录系统
显示输出包括:单通道和波束域传播损失输出,包括18个频段。
数据记录系统:数据记录仪记录原始的阵元数据、场景数据和传播损失计算结果,同时,可以利用数据回放功能实现离线的数据处理。
本发明不局限于上述实施方式,不论其实施方式作任何变化,凡是采用本发明所提供的思路,都是本发明的一种变形,均应认为在发明的保护范围之内。
Claims (3)
1.一种基于声阵的海洋声传播损失实时测量方法,其特征在于:
用于单阵元传播损失实时测量,具体包括以下步骤:
步骤一:对实时接收的数据进行数据累计,累积8批数据得到x1(i,n)和x1_1(i,n),x1(i,n)和x1_1(i,n)分别为累积8批得到的当前测量数据和之前数据;
步骤二:对接收的数据进行傅立叶变化处理:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mrow>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
</mrow>
<mn>16384</mn>
</munder>
<mo>{</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mn>1</mn>
<mo>_</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mrow>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
</mrow>
<mn>16384</mn>
</munder>
<mo>{</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>_</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,16384为傅立叶变换的长度,处理数据有87.5%的重叠,即每次更新2048点,频率分辨力为0.7324Hz;
步骤三:分频段进行能量计算:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>3</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mn>2</mn>
</munderover>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>N</mi>
<mi>k</mi>
</msub>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<msub>
<mi>k</mi>
<mi>H</mi>
</msub>
</munderover>
<mo>|</mo>
<msub>
<mi>Y</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mi>b</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>18</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>n</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>3</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mn>2</mn>
</munderover>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>N</mi>
<mi>k</mi>
</msub>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<msub>
<mi>k</mi>
<mi>H</mi>
</msub>
</munderover>
<mo>|</mo>
<msub>
<mi>Y</mi>
<mrow>
<mn>1</mn>
<mo>_</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mi>b</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>18</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤四:对数据进行能量判断,以确定有无信号存在;在计算传播损失时,如果是脉冲信号,需进行信噪比判断:如果Est(fb)-Ent(fb)>Th,Th为门限,则输出传播损失;否则不输出;如果是连续波,则不需要进行能量判断;
步骤五:进行传播损失计算:
TL(fb)=SL(fb)-Est(fb)-HS-AG (5)
其中SL为声源级,事先标定;HS表示水听器灵敏度,单位为dB;AG表示放大倍数,单位为dB。
2.一种基于声阵的海洋声传播损失实时测量方法,其特征在于:
用于波束域传播损失实时测量,利用64元细线阵或96元粗线阵输出信号测量阵波束传播损失;接收数据表示为x1_64(i,n),i=0,……63;x1_96(i,n),i=0,……95,其中,i表示水听器通道号,x1_64(i,n)表示前64个水听器接收的阵数据,x1_96(i,n)表示后96个水听器接收的阵数据;n=0,……,2048为时间采样序列;具体包括以下步骤:
步骤一:对实时接收的数据进行滤波:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>h</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>&CircleTimes;</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<mn>63</mn>
<mo>;</mo>
<mn>0</mn>
<mo>&le;</mo>
<mi>n</mi>
<mo>&le;</mo>
<mn>2047</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>h</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>&CircleTimes;</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<mi>i</mi>
<mo>&le;</mo>
<mn>95</mn>
<mo>;</mo>
<mn>0</mn>
<mo>&le;</mo>
<mi>n</mi>
<mo>&le;</mo>
<mn>2047</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中h1(n)为带通滤波器,滤波器通带为50-500Hz,过渡带为[10 50Hz]和[500 600Hz],通带波纹小于等于0.5dB,阻带衰减大于等于40dB;h2(n)为带通滤波器,滤波器通带为500-2000Hz,过渡带为[400 500Hz]和[2000 2100Hz],通带波纹小于等于0.5dB,阻带衰减大于等于40dB;
步骤二:对接收的数据进行傅立叶变换处理:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mn>1</mn>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2048</mn>
</mfrac>
<munder>
<mrow>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
</mrow>
<mn>2048</mn>
</munder>
<mo>{</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mo>...</mo>
<mo>...</mo>
<mo>,</mo>
<mn>63</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
1
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mn>1</mn>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2048</mn>
</mfrac>
<munder>
<mrow>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
</mrow>
<mn>2048</mn>
</munder>
<mo>{</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>2</mn>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mo>...</mo>
<mo>...</mo>
<mo>,</mo>
<mn>95</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤三:分别对两个阵的接收数据进行波束形成处理,在余弦坐标系中覆盖0°~180°水平空间:
<mrow>
<msub>
<mi>Bf</mi>
<mrow>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>64</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mn>63</mn>
</munderover>
<msub>
<mi>Y</mi>
<mrow>
<mn>1</mn>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mo>&times;</mo>
<msub>
<mi>f</mi>
<mi>k</mi>
</msub>
<mo>&times;</mo>
<mi>d</mi>
<mn>1</mn>
<mo>&times;</mo>
<mi>i</mi>
<mo>&times;</mo>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mi>l</mi>
<mn>64</mn>
</mfrac>
</mrow>
<mo>)</mo>
<mo>/</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中0≤l≤128,k=17,...,85;
<mrow>
<msub>
<mi>Bf</mi>
<mrow>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>96</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mn>95</mn>
</munderover>
<msub>
<mi>Y</mi>
<mrow>
<mn>1</mn>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mo>&times;</mo>
<msub>
<mi>f</mi>
<mi>k</mi>
</msub>
<mo>&times;</mo>
<mi>d</mi>
<mn>2</mn>
<mo>&times;</mo>
<mi>i</mi>
<mo>&times;</mo>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mi>l</mi>
<mn>64</mn>
</mfrac>
</mrow>
<mo>)</mo>
<mo>/</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
其中0≤l≤128,k=85,...,341;数据重叠50%,基元间距d1=1.5m,d2=0.4m,c为声速;
步骤四:波束选取;根据显控输入的人工选取方位或自动选取的方位lmax选择信号波束;自动选取方位过程如下:
<mrow>
<msub>
<mi>BP</mi>
<mrow>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>9</mn>
</mrow>
<mn>85</mn>
</munderover>
<mo>|</mo>
<msub>
<mi>Bf</mi>
<mrow>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
公式(12)表示对100~500Hz频段信号计算波束图;
<mrow>
<msub>
<mi>BP</mi>
<mrow>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>85</mn>
</mrow>
<mn>341</mn>
</munderover>
<mo>|</mo>
<msub>
<mi>Bf</mi>
<mrow>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
公式(13)表示对500~2000Hz频段信号计算波束图;利用公式(14)判断选取的方位;
步骤五:对选定的lmax波束数据进行逆傅立叶变换处理:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Bftime</mi>
<mrow>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>2048</mn>
<mo>*</mo>
<munder>
<mrow>
<mi>I</mi>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
</mrow>
<mn>2048</mn>
</munder>
<mo>{</mo>
<msub>
<mi>Bf</mi>
<mrow>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>9</mn>
<mo>,</mo>
<mn>10</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>85</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Bftime</mi>
<mrow>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>2048</mn>
<mo>*</mo>
<munder>
<mrow>
<mi>I</mi>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
</mrow>
<mn>2048</mn>
</munder>
<mo>{</mo>
<msub>
<mi>Bf</mi>
<mrow>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>85</mn>
<mo>,</mo>
<mn>86</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>341</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,2048为逆傅立叶变换的长度;
步骤六:对人工选取方位或自动选取的方位数据累积8批,得到的Bftime1_64(n)和Bftime1_96(n),分别为指定方位的累积8批得到的当前测量数据,Bftime1_1_64(n)和Bftime1_1_96(n)分别为指定方位的累积8批得到的之前噪声数据;
步骤七:对目标方位数据累积进行傅立叶变换处理:
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>16384</mn>
</mfrac>
<munder>
<mrow>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
</mrow>
<mn>16384</mn>
</munder>
<mo>{</mo>
<mi>B</mi>
<mi>f</mi>
<mi>t</mi>
<mi>i</mi>
<mi>m</mi>
<mi>e</mi>
<msub>
<mn>1</mn>
<mrow>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>Y</mi>
<mo>_</mo>
<msub>
<mn>1</mn>
<mrow>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>16384</mn>
</mfrac>
<munder>
<mrow>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
</mrow>
<mn>16384</mn>
</munder>
<mo>{</mo>
<mi>B</mi>
<mi>f</mi>
<mi>t</mi>
<mi>i</mi>
<mi>m</mi>
<mi>e</mi>
<mn>1</mn>
<mo>_</mo>
<msub>
<mn>1</mn>
<mrow>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>Y</mi>
<mrow>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>16384</mn>
</mfrac>
<munder>
<mrow>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
</mrow>
<mn>16384</mn>
</munder>
<mo>{</mo>
<mi>B</mi>
<mi>f</mi>
<mi>t</mi>
<mi>i</mi>
<mi>m</mi>
<mi>e</mi>
<msub>
<mn>1</mn>
<mrow>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>Y</mi>
<mo>_</mo>
<msub>
<mn>1</mn>
<mrow>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>16384</mn>
</mfrac>
<munder>
<mrow>
<mi>F</mi>
<mi>F</mi>
<mi>T</mi>
</mrow>
<mn>16384</mn>
</munder>
<mo>{</mo>
<mi>B</mi>
<mi>f</mi>
<mi>t</mi>
<mi>i</mi>
<mi>m</mi>
<mi>e</mi>
<mn>1</mn>
<mo>_</mo>
<msub>
<mn>1</mn>
<mrow>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>20</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,16384为傅立叶变换的长度,处理数据有87.5%的重叠,即每次更新2048点;频率分辨力为0.7324Hz;
步骤八:对波束形成并积累后的数据进行分频段能量计算:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>N</mi>
<mi>k</mi>
</msub>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<msub>
<mi>k</mi>
<mi>H</mi>
</msub>
</munderover>
<mo>|</mo>
<msub>
<mi>Y</mi>
<mrow>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mi>b</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>8</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>21</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>n</mi>
<mi>t</mi>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>N</mi>
<mi>k</mi>
</msub>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<msub>
<mi>k</mi>
<mi>H</mi>
</msub>
</munderover>
<mo>|</mo>
<mi>Y</mi>
<mo>_</mo>
<msub>
<mn>1</mn>
<mrow>
<mo>_</mo>
<mn>64</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mi>b</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>8</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>22</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mi>t</mi>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>N</mi>
<mi>k</mi>
</msub>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<msub>
<mi>k</mi>
<mi>H</mi>
</msub>
</munderover>
<mo>|</mo>
<msub>
<mi>Y</mi>
<mrow>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mi>b</mi>
<mo>=</mo>
<mn>9</mn>
<mo>,</mo>
<mn>10</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>14</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>23</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>n</mi>
<mi>t</mi>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>N</mi>
<mi>k</mi>
</msub>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>L</mi>
</msub>
</mrow>
<msub>
<mi>k</mi>
<mi>H</mi>
</msub>
</munderover>
<mo>|</mo>
<mi>Y</mi>
<mo>_</mo>
<msub>
<mn>1</mn>
<mrow>
<mo>_</mo>
<mn>96</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>f</mi>
<mi>b</mi>
<mo>=</mo>
<mn>9</mn>
<mo>,</mo>
<mn>10</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>14</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>24</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤九:对接收信号进行能量判断;如果是脉冲信号需进行信噪比判断:如果Est_64(fb)-Ent_64(fb)>Th,Est_96(fb)-Ent_96(fb)>Th,Th为门限,则输出传播损失;否则,不输出;如果是连续波信号不需要进行能量判断;
步骤十:传播损失计算:
TL_64(fb)=SL(fb)-Est_64(fb)-HS-AG (25)
TL_96(fb)=SL(fb)-Est_96(fb)-HS-AG (26)
其中SL为配合舰船发射声源或爆炸声源的声源级,事先标定;HS表示水听器灵敏度,单位为dB;AG表示放大倍数,单位为dB。
3.根据权利要求1或2所述的一种基于声阵的海洋声传播损失实时测量方法,其特征在于:采用水平拖曳线列阵实现,所述水平拖曳线列阵由16段声学段组成,各声学段由其内置的水听器阵列、前置与传输部分一起,完成声信号的声电转换、前置预处理、信号采集和传输功能;阵元数为184元,为实现波束域不同频段传播损失计算,采取布阵方式是:前64个声通道,阵元间距1.5m,中间24个声通道,阵元间距24m,最后96个声通道,阵元间距0.4m;声基阵总长758.4m。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410422855.3A CN104180891B (zh) | 2014-08-26 | 2014-08-26 | 一种基于声阵的海洋声传播损失实时测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410422855.3A CN104180891B (zh) | 2014-08-26 | 2014-08-26 | 一种基于声阵的海洋声传播损失实时测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104180891A CN104180891A (zh) | 2014-12-03 |
CN104180891B true CN104180891B (zh) | 2017-09-05 |
Family
ID=51962120
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410422855.3A Active CN104180891B (zh) | 2014-08-26 | 2014-08-26 | 一种基于声阵的海洋声传播损失实时测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104180891B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104568113B (zh) * | 2014-12-31 | 2017-06-27 | 浙江大学 | 一种基于模型的海洋声传播调查爆炸波自动截取方法 |
CN105301114B (zh) * | 2015-10-20 | 2018-01-12 | 浙江大学 | 一种基于多通道空时逆滤波技术的声学覆盖层插入损失测量方法 |
CN108008445B (zh) * | 2017-10-11 | 2019-10-18 | 中国船舶重工集团公司第七一五研究所 | 一种改进的浅海传播损失预报模型的计算方法 |
CN109000779B (zh) * | 2018-07-12 | 2019-12-20 | 青岛中海潮科技有限公司 | 一种海洋声传播可靠性模型快速匹配系统 |
CN110109087B (zh) * | 2019-05-07 | 2021-02-02 | 中国科学院声学研究所 | 一种声纳不规则探测范围显示方法及系统 |
CN111024208B (zh) * | 2019-11-26 | 2022-03-25 | 中国船舶重工集团有限公司第七一0研究所 | 一种垂直阵声压梯度波束形成与信号检测方法 |
-
2014
- 2014-08-26 CN CN201410422855.3A patent/CN104180891B/zh active Active
Non-Patent Citations (5)
Title |
---|
A modified transfer matrix method for prediction of transmission loss of multilayer acoustic materials;C.-M. Lee 等;《Journal of Sound and Vibration》;20090925;第326卷(第1-2期);第290-301页 * |
Comparison and implementation of the various numerical methods used for calculating transmission loss in silencer systems;S. Bilawchuk 等;《Applied Acoustics》;20030930;第64卷(第9期);第903-916页 * |
一种用能量谱去噪声计算声传播损失的方法;刘贞文等;《厦门大学学报(自然科学版)》;20090531;第48卷(第3期);第378-381页 * |
大陆坡海域二维声传播研究;秦继兴等;《声学学报》;20140331;第39卷(第2期);第145-153页 * |
空气中声源激发的浅海水下声场传播实验研究;王光旭等;《声学学报》;20111130;第36卷(第6期);第588-595页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104180891A (zh) | 2014-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104180891B (zh) | 一种基于声阵的海洋声传播损失实时测量方法 | |
CN101852854B (zh) | 一种水下多波束测探系统及其探测方法 | |
Dalen et al. | The influence of wind‐induced bubbles on echo integration surveys | |
CN103076594B (zh) | 一种基于互相关的水声脉冲信号双阵元定位的方法 | |
CN103616693B (zh) | 探渔声纳及声纳回波信号处理方法 | |
CN106291564B (zh) | 一种海底冷泉水体回声反射探测系统和方法 | |
CN105301114A (zh) | 一种基于多通道空时逆滤波技术的声学覆盖层插入损失测量方法 | |
CN112269164A (zh) | 深海可靠声路径下基于干涉结构匹配处理弱目标定位方法 | |
Fialkowski et al. | Methods for identifying and controlling sonar clutter | |
CN104656073B (zh) | 三维成像声纳波束形成方法及在多核处理器上的实现方法 | |
CN108107436A (zh) | 一种基于可靠声路径的水下目标主动分类与定位方法 | |
CN114167423A (zh) | 基于深度回归网络的雷达海浪参数测量方法 | |
CN104678384A (zh) | 一种波束域的声压差互相关谱分析水下目标速度估计方法 | |
RU137126U1 (ru) | Гидроакустический комплекс надводного корабля | |
CN110763327A (zh) | 一种平面障板条件下单矢量水听器信号检测方法 | |
CN109061654A (zh) | 一种深海环境下单圆环阵主动三维定位方法 | |
CN109100711A (zh) | 一种深海环境下单基地主动声纳低运算量三维定位方法 | |
CN101231342A (zh) | 拖曳式水平渔探仪 | |
CN112255607A (zh) | 一种海杂波的抑制方法 | |
CN101762824B (zh) | 基于单向水声测距的海洋地震拖缆位置测量方法 | |
Lunkov et al. | Using discrete low-frequency components of shipping noise for gassy sediment characterization in shallow water | |
CN103926586B (zh) | 一种使用发射子阵的mimo阵列测深方法 | |
CN106680823B (zh) | 利用抹香鲸叫声脉冲的目标距离和速度探测方法与装置 | |
Vahidpour et al. | An automated approach to passive sonar classification using binary image features | |
CN100473959C (zh) | 海洋声场相速度的水听器测试方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |