CN110390073B - 一种矢量传感的多通道空间合成方位滤波方法 - Google Patents

一种矢量传感的多通道空间合成方位滤波方法 Download PDF

Info

Publication number
CN110390073B
CN110390073B CN201910762326.0A CN201910762326A CN110390073B CN 110390073 B CN110390073 B CN 110390073B CN 201910762326 A CN201910762326 A CN 201910762326A CN 110390073 B CN110390073 B CN 110390073B
Authority
CN
China
Prior art keywords
signal
vector
angle
channel
theta
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
CN201910762326.0A
Other languages
English (en)
Other versions
CN110390073A (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201910762326.0A priority Critical patent/CN110390073B/zh
Publication of CN110390073A publication Critical patent/CN110390073A/zh
Application granted granted Critical
Publication of CN110390073B publication Critical patent/CN110390073B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H3/00Measuring characteristics of vibrations by using a detector in a fluid
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12Simultaneous equations, e.g. systems of linear equations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Computing Systems (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明提供了一种矢量传感的多通道空间合成方位滤波方法,利用矢量传感器采集4通道信号,设定波达方向的角度搜索范围间隔,对应将给出水平方向角的角度范围的分割为多个搜索空间角度范围的集合,对多空间方位的每一个搜索范围间隔Δθ进行矢量合成,将合成信号与声压信号p(t)做互相关,对相关度值进行排序,选择最大相关度值并输出其所对应的合成信号。本发明利用矢量传感接收信号的空间方位特性,对平面振速通道的离散空间合成,利用噪声的空间性及信号的方位性,通过对合成信号与声压通道的相关性实现方位滤波,滤除空间其他方位的干扰信息,解决了传统标量声传感技术下滤波降噪的难题。

Description

一种矢量传感的多通道空间合成方位滤波方法
技术领域
本发明涉及信号处理领域,尤其是一种新型水下矢量信号处理方法。
背景技术
21世纪以来,世界各国在政治、经济、军事方面围绕海洋领域的竞争愈演愈烈,世界各国均提出了相应的海洋发展战略,对海洋资源的保护、开发、利用已成为世界共同关注的焦点。当前,我国的国家核心利益主要体现在经济发展和安全权益两个方面,经济建设是改革开放以来的中心任务,而维护海洋权益是新时期实现海洋强国的根本保证。因此,研究先进的微弱信号检测方法对远距离水中目标探测和识别具有巨大的研究价值和现实意义。
矢量声信号是由矢量水听器对声场空间一点处的声压和质点振速的若干正交分量直接测量所得,与传统水听器测得的声压信号相比,矢量水听器(UAV)一般包含一个声场的声压(标量)通道以及三个声场振速的三维正交矢量通道,能够反映更加完整的声场信息,可以支持更丰富的信号处理方法及应用需求。
针对这种新型的矢量传感器,利用其独有的空间方位特性开展信号处理的方法是一个值得深入的研究方向,由于环境噪声具空间性,其本质是产生自不同方位、不同距离的合成,因此利用矢量合成来实现方位滤波方法可以得到强的空间指向性,进而提升滤波降噪性能,能够为空间方位相关的信号处理方法及应用奠定基础。
发明内容
为了克服现有技术的不足,针对常规方法依赖先验信息、在远距离低信噪比下效果变差的问题,本发明提供一种矢量传感的多通道空间合成方位滤波方法。
本发明解决其技术问题所采用的技术方案的具体步骤如下:
第一步:利用矢量传感器采集4通道信号,声压p(t)及正交的三维振速分量vx(t)、vy(t)、vz(t),各通道的接收信号表示为:
Figure BDA0002170723740000011
其中s(t)为目标信号,np(t)、nx(t)、ny(t)、nz(t)表示各通道的噪声,为独立同分布的;ρ为环境密度,c为声速,θ为水平方向角,α为俯仰角,表示为向量形式r:
r=s(t)I+n (2)
式中I为振速的方向矢量,n为噪声方向矢量:
Figure BDA0002170723740000021
第二步:考虑二维平面,初始化波达方向的水平方向角θ的角度范围[0,2π];
第三步:设定波达方向的角度搜索范围间隔为Δθ,对应将给出水平方向角θ的角度范围[0,2π]的分割为多个搜索空间角度范围的集合{[θ12],[θ23],...,[θM-1M]},θ1=0,θM=2π,Δθ=2π/M为一个确定的值,实际中M的取值可根据要求的精度与算法复杂度进行折衷选取,一般Δθ=1°;
第四步:对多空间方位的每一个搜索范围间隔Δθ进行矢量合成,得到合成后的信号:
Figure BDA0002170723740000022
其中r(t,Δθ)表示合成信号,N表示对每一个搜索空间角度范围如[θ12]的等分离散化,i表示对搜索空间角度范围[θ12]的第i个等分离散化,θi=i(θ21)/N;
第五步:将合成r(t,Δθ)与声压信号p(t)做互相关,记录对应的相关度值
Figure BDA0002170723740000023
Figure BDA0002170723740000024
其中μr和μp分别为合成信号与声压信号的均值,重复第三步到第五步直到搜索完成全部的水平方向角θ的角度范围;
第六步:对第五步得到的相关度值进行排序,选择最大相关度值并输出其所对应的合成信号。
本发明的有益效果在于利用矢量传感接收信号的空间方位特性,对平面振速通道的离散空间合成,利用噪声的空间性及信号的方位性,通过对合成信号与声压通道的相关性实现方位滤波,滤除空间其他方位的干扰信息,解决了传统标量声传感技术下滤波降噪的难题。
附图说明
图1为本发明矢量被动空间感知示意图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
针对新型矢量传感器的快速发展,利用其空间特性,本发明提供一种单矢量传感的空间合成方位滤波方法。
本发明解决其技术问题所采用的技术方案的具体步骤如下:
第一步:利用矢量传感器采集4通道信号,矢量传感通常由声压传感器和三维振速传感器构成,测量空间同一点处的声压P及正交的三维振速分量。其中三维振速分量vx、vy、vz是对单位功率水声波阵面测量入射波场声振速v(t)的笛卡尔分量:
Figure BDA0002170723740000031
声压p(t)及正交的三维振速分量vx(t)、vy(t)、vz(t),各通道的接收信号表示为:
Figure BDA0002170723740000032
其中s(t)为目标信号,np(t)、nx(t)、ny(t)、nz(t)表示各通道的噪声,为独立同分布的;ρ为环境密度,c为声速,θ为水平方向角,x轴正方向为0°,α为俯仰角,以水平面(xoy平面)为0°,t表示时间,表示为向量形式r:
r=s(t)I+n (2)
式中I为振速的方向矢量,n为噪声方向矢量:
Figure BDA0002170723740000033
第二步:考虑二维平面,初始化波达方向的水平方向角θ的角度范围[0,2π];
第三步:设定波达方向的角度搜索范围间隔为Δθ,对应将给出水平方向角θ的角度范围[0,2π]的分割为多个搜索空间角度范围的集合{[θ12],[θ23],...,[θM-1M]},θ1=0,θM=2π,Δθ=2π/M为一个确定的值,实际中M的取值可根据要求的精度与算法复杂度进行折衷选取,一般Δθ=1°;
第四步:对多空间方位的每一个搜索范围间隔Δθ进行矢量合成,得到合成后的信号:
Figure BDA0002170723740000041
其中r(t,Δθ)表示合成信号,N表示对每一个搜索空间角度范围如[θ12]的等分离散化,i表示对搜索空间角度范围[θ12]的第i个等分离散化,θi=i(θ21)/N;
第五步:将合成r(t,Δθ)与声压信号p(t)做互相关,记录对应的相关度值
Figure BDA0002170723740000042
Figure BDA0002170723740000043
其中μr和μp分别为合成信号与声压信号的均值,重复第三步到第五步直到搜索完成全部的水平方向角θ的角度范围;
第六步:对第五步得到的相关度值进行排序,选择最大相关度值并输出其所对应的合成信号。

Claims (1)

1.一种矢量传感的多通道空间合成方位滤波方法,其特征在于包括下述步骤:
第一步:利用矢量传感器采集4通道信号,声压p(t)及正交的三维振速分量vx(t)、vy(t)、vz(t),各通道的接收信号表示为:
Figure FDA0002170723730000011
其中s(t)为目标信号,np(t)、nx(t)、ny(t)、nz(t)表示各通道的噪声,为独立同分布的;ρ为环境密度,c为声速,θ为水平方向角,α为俯仰角,表示为向量形式r:
r=s(t)I+n (2)
式中I为振速的方向矢量,n为噪声方向矢量:
Figure FDA0002170723730000012
第二步:考虑二维平面,初始化波达方向的水平方向角θ的角度范围[0,2π];
第三步:设定波达方向的角度搜索范围间隔为Δθ,对应将给出水平方向角θ的角度范围[0,2π]的分割为多个搜索空间角度范围的集合{[θ12],[θ23],...,[θM-1M]},θ1=0,θM=2π,Δθ=2π/M为一个确定的值;
第四步:对多空间方位的每一个搜索范围间隔Δθ进行矢量合成,得到合成后的信号:
Figure FDA0002170723730000013
其中r(t,Δθ)表示合成信号,N表示对每一个搜索空间角度范围如[θ12]的等分离散化,i表示对搜索空间角度范围[θ12]的第i个等分离散化,θi=i(θ21)/N;
第五步:将合成r(t,Δθ)与声压信号p(t)做互相关,记录对应的相关度值
Figure FDA0002170723730000014
Figure FDA0002170723730000015
其中μr和μp分别为合成信号与声压信号的均值,重复第三步到第五步直到搜索完成全部的水平方向角θ的角度范围;
第六步:对第五步得到的相关度值进行排序,选择最大相关度值并输出其所对应的合成信号。
CN201910762326.0A 2019-08-19 2019-08-19 一种矢量传感的多通道空间合成方位滤波方法 Active CN110390073B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910762326.0A CN110390073B (zh) 2019-08-19 2019-08-19 一种矢量传感的多通道空间合成方位滤波方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910762326.0A CN110390073B (zh) 2019-08-19 2019-08-19 一种矢量传感的多通道空间合成方位滤波方法

Publications (2)

Publication Number Publication Date
CN110390073A CN110390073A (zh) 2019-10-29
CN110390073B true CN110390073B (zh) 2023-03-24

Family

ID=68288904

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910762326.0A Active CN110390073B (zh) 2019-08-19 2019-08-19 一种矢量传感的多通道空间合成方位滤波方法

Country Status (1)

Country Link
CN (1) CN110390073B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112071297B (zh) * 2020-09-07 2023-11-10 西北工业大学 一种矢量声的自适应滤波方法
CN113253196B (zh) * 2021-06-01 2024-05-21 中国电子科技集团公司第三十六研究所 一种多信号测向方法、装置和电子设备

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2010126808A (ru) * 2010-06-30 2012-01-10 Тихоокеанский океанологический институт им. В.И. Ильичева Дальневосточного отделения Российской академии наук (ТОИ ДВО РАН) (RU) Способ определения пеленга на шумящий объект
RU2474836C1 (ru) * 2011-07-20 2013-02-10 Учреждение Российской академии наук Институт проблем морских технологий Дальневосточного отделения РАН (ИПМТ ДВО РАН) Гидроакустический комплекс для измерения азимутального угла на источник звука в мелком море
CN103399312A (zh) * 2013-07-15 2013-11-20 哈尔滨工程大学 基于压缩感知技术的时空联合滤波高分辨方位估计方法
JP2014069038A (ja) * 2012-10-02 2014-04-21 Hitachi Medical Corp 多チャンネル画像合成方法及び磁気共鳴イメージング装置
CN104698431A (zh) * 2015-03-17 2015-06-10 河海大学 模糊分量空间角的估计及多通道sar方位解模糊方法
CN105207964A (zh) * 2015-08-17 2015-12-30 哈尔滨工程大学 一种基于单矢量传感器的水声自适应判决反馈均衡方法
CN106199600A (zh) * 2016-06-30 2016-12-07 西安电子科技大学 基于多普勒估计的方位多通道合成孔径雷达成像方法
CN106443672A (zh) * 2016-08-30 2017-02-22 西安电子科技大学 一种方位多通道sar信号自适应重构方法
RU2653585C1 (ru) * 2017-09-08 2018-05-11 Федеральное государственное бюджетное учреждение науки Институт проблем морских технологий Дальневосточного отделения Российской академии наук (ИМПТ ДВО РАН) Способ обнаружения шумящих, движущихся в море объектов
CN109283492A (zh) * 2018-10-29 2019-01-29 中国电子科技集团公司第三研究所 多目标方位估计方法及水声垂直矢量阵列系统

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2010126808A (ru) * 2010-06-30 2012-01-10 Тихоокеанский океанологический институт им. В.И. Ильичева Дальневосточного отделения Российской академии наук (ТОИ ДВО РАН) (RU) Способ определения пеленга на шумящий объект
RU2474836C1 (ru) * 2011-07-20 2013-02-10 Учреждение Российской академии наук Институт проблем морских технологий Дальневосточного отделения РАН (ИПМТ ДВО РАН) Гидроакустический комплекс для измерения азимутального угла на источник звука в мелком море
JP2014069038A (ja) * 2012-10-02 2014-04-21 Hitachi Medical Corp 多チャンネル画像合成方法及び磁気共鳴イメージング装置
CN103399312A (zh) * 2013-07-15 2013-11-20 哈尔滨工程大学 基于压缩感知技术的时空联合滤波高分辨方位估计方法
CN104698431A (zh) * 2015-03-17 2015-06-10 河海大学 模糊分量空间角的估计及多通道sar方位解模糊方法
CN105207964A (zh) * 2015-08-17 2015-12-30 哈尔滨工程大学 一种基于单矢量传感器的水声自适应判决反馈均衡方法
CN106199600A (zh) * 2016-06-30 2016-12-07 西安电子科技大学 基于多普勒估计的方位多通道合成孔径雷达成像方法
CN106443672A (zh) * 2016-08-30 2017-02-22 西安电子科技大学 一种方位多通道sar信号自适应重构方法
RU2653585C1 (ru) * 2017-09-08 2018-05-11 Федеральное государственное бюджетное учреждение науки Институт проблем морских технологий Дальневосточного отделения Российской академии наук (ИМПТ ДВО РАН) Способ обнаружения шумящих, движущихся в море объектов
CN109283492A (zh) * 2018-10-29 2019-01-29 中国电子科技集团公司第三研究所 多目标方位估计方法及水声垂直矢量阵列系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于单矢量传感器的强干扰噪声抑制方法;高鑫等;《传感器与微系统》;20100920(第09期);全文 *
海洋矢量噪声场试验研究;邹吉武等;《哈尔滨工程大学学报》;20110125(第01期);全文 *

Also Published As

Publication number Publication date
CN110390073A (zh) 2019-10-29

Similar Documents

Publication Publication Date Title
CN110390073B (zh) 一种矢量传感的多通道空间合成方位滤波方法
CN101339242B (zh) 声源定位的无线测量方法
CN112987004B (zh) 一种浅海环境下基于水平阵列的水面水下目标分类方法
CN110824429B (zh) 深海环境下利用非同步垂直阵的宽带声源被动定位方法
CN109669160B (zh) 一种水下瞬态声信号的检测方法
CN111273265B (zh) 基于demon谱的单矢量水听器自主测距装置及方法
CN103048642A (zh) 基于频域最小二乘法的水声脉冲信号匹配场定位方法
CN110133580B (zh) 一种随机共振增强的声矢量信号定向方法
CN112130112B (zh) 一种基于声矢量阵联合信息处理信源数估计方法
CN101915922A (zh) 拖曳线列阵被动测距方法
CN114280541B (zh) 一种基于深海分布式垂直线列阵的目标被动定位方法
CN109319074B (zh) 一种多正交信号发射无人潜器声引导回收系统
CN111537982A (zh) 一种畸变拖曳阵线谱特征增强方法及系统
CN104678384A (zh) 一种波束域的声压差互相关谱分析水下目标速度估计方法
CN111580079B (zh) 一种基于单矢量水听器的水下运动目标最近接近距离估计方法
CN111679245B (zh) 一种基于均匀圆阵的子阵重合的分裂波束定向方法
CN105353340B (zh) 一种双层圆柱阵水下被动目标检测方法
CN112098938A (zh) 一种基于六元锥矢量阵的水声目标降维匹配声场定位方法
CN110501669B (zh) 一种中心对称声矢量圆阵快速空间谱压缩超分辨方位估计方法
CN110095755B (zh) 一种声源定位方法
CN110502823B (zh) 一种矢量传感的多源混叠信号分离方法
CN108414984A (zh) 一种基于二阶干涉的水下目标定位方法
CN111504253A (zh) 一种波浪滑翔机的水下三维声学立体阵确定方法
CN112162235B (zh) 一种光滑分段随机共振增强的声矢量信号定向方法
CN105138743A (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