CN110390073A - 一种矢量传感的多通道空间合成方位滤波方法 - Google Patents
一种矢量传感的多通道空间合成方位滤波方法 Download PDFInfo
- Publication number
- CN110390073A CN110390073A CN201910762326.0A CN201910762326A CN110390073A CN 110390073 A CN110390073 A CN 110390073A CN 201910762326 A CN201910762326 A CN 201910762326A CN 110390073 A CN110390073 A CN 110390073A
- Authority
- CN
- China
- Prior art keywords
- vector
- signal
- angle
- channel
- horizontal direction
- 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
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H3/00—Measuring characteristics of vibrations by using a detector in a fluid
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation 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),各通道的接收信号表示为:
其中s(t)为目标信号,np(t)、nx(t)、ny(t)、nz(t)表示各通道的噪声,为独立同分布的;ρ为环境密度,c为声速,θ为水平方向角,α为俯仰角,表示为向量形式r:
r=s(t)I+n (2)
式中I为振速的方向矢量,n为噪声方向矢量:
第二步:考虑二维平面,初始化波达方向的水平方向角θ的角度范围[0,2π];
第三步:设定波达方向的角度搜索范围间隔为Δθ,对应将给出水平方向角θ的角度范围[0,2π]的分割为多个搜索空间角度范围的集合{[θ1,θ2],[θ2,θ3],...,[θM-1,θM]},θ1=0,θM=2π,Δθ=2π/M为一个确定的值,实际中M的取值可根据要求的精度与算法复杂度进行折衷选取,一般Δθ=1°;
第四步:对多空间方位的每一个搜索范围间隔Δθ进行矢量合成,得到合成后的信号:
其中r(t,Δθ)表示合成信号,N表示对每一个搜索空间角度范围如[θ1,θ2]的等分离散化,i表示对搜索空间角度范围[θ1,θ2]的第i个等分离散化,θi=i(θ2-θ1)/N;
第五步:将合成r(t,Δθ)与声压信号p(t)做互相关,记录对应的相关度值
其中μr和μp分别为合成信号与声压信号的均值,重复第三步到第五步直到搜索完成全部的水平方向角θ的角度范围;
第六步:对第五步得到的相关度值进行排序,选择最大相关度值并输出其所对应的合成信号。
本发明的有益效果在于利用矢量传感接收信号的空间方位特性,对平面振速通道的离散空间合成,利用噪声的空间性及信号的方位性,通过对合成信号与声压通道的相关性实现方位滤波,滤除空间其他方位的干扰信息,解决了传统标量声传感技术下滤波降噪的难题。
附图说明
图1为本发明矢量被动空间感知示意图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
针对新型矢量传感器的快速发展,利用其空间特性,本发明提供一种单矢量传感的空间合成方位滤波方法。
本发明解决其技术问题所采用的技术方案的具体步骤如下:
第一步:利用矢量传感器采集4通道信号,矢量传感通常由声压传感器和三维振速传感器构成,测量空间同一点处的声压P及正交的三维振速分量。其中三维振速分量vx、vy、vz是对单位功率水声波阵面测量入射波场声振速v(t)的笛卡尔分量:
声压p(t)及正交的三维振速分量vx(t)、vy(t)、vz(t),各通道的接收信号表示为:
其中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为噪声方向矢量:
第二步:考虑二维平面,初始化波达方向的水平方向角θ的角度范围[0,2π];
第三步:设定波达方向的角度搜索范围间隔为Δθ,对应将给出水平方向角θ的角度范围[0,2π]的分割为多个搜索空间角度范围的集合{[θ1,θ2],[θ2,θ3],...,[θM-1,θM]},θ1=0,θM=2π,Δθ=2π/M为一个确定的值,实际中M的取值可根据要求的精度与算法复杂度进行折衷选取,一般Δθ=1°;
第四步:对多空间方位的每一个搜索范围间隔Δθ进行矢量合成,得到合成后的信号:
其中r(t,Δθ)表示合成信号,N表示对每一个搜索空间角度范围如[θ1,θ2]的等分离散化,i表示对搜索空间角度范围[θ1,θ2]的第i个等分离散化,θi=i(θ2-θ1)/N;
第五步:将合成r(t,Δθ)与声压信号p(t)做互相关,记录对应的相关度值
其中μr和μp分别为合成信号与声压信号的均值,重复第三步到第五步直到搜索完成全部的水平方向角θ的角度范围;
第六步:对第五步得到的相关度值进行排序,选择最大相关度值并输出其所对应的合成信号。
Claims (1)
1.一种矢量传感的多通道空间合成方位滤波方法,其特征在于包括下述步骤:
第一步:利用矢量传感器采集4通道信号,声压p(t)及正交的三维振速分量vx(t)、vy(t)、vz(t),各通道的接收信号表示为:
其中s(t)为目标信号,np(t)、nx(t)、ny(t)、nz(t)表示各通道的噪声,为独立同分布的;ρ为环境密度,c为声速,θ为水平方向角,α为俯仰角,表示为向量形式r:
r=s(t)I+n (2)
式中I为振速的方向矢量,n为噪声方向矢量:
第二步:考虑二维平面,初始化波达方向的水平方向角θ的角度范围[0,2π];
第三步:设定波达方向的角度搜索范围间隔为Δθ,对应将给出水平方向角θ的角度范围[0,2π]的分割为多个搜索空间角度范围的集合{[θ1,θ2],[θ2,θ3],...,[θM-1,θM]},θ1=0,θM=2π,Δθ=2π/M为一个确定的值;
第四步:对多空间方位的每一个搜索范围间隔Δθ进行矢量合成,得到合成后的信号:
其中r(t,Δθ)表示合成信号,N表示对每一个搜索空间角度范围如[θ1,θ2]的等分离散化,i表示对搜索空间角度范围[θ1,θ2]的第i个等分离散化,θi=i(θ2-θ1)/N;
第五步:将合成r(t,Δθ)与声压信号p(t)做互相关,记录对应的相关度值
其中μr和μp分别为合成信号与声压信号的均值,重复第三步到第五步直到搜索完成全部的水平方向角θ的角度范围;
第六步:对第五步得到的相关度值进行排序,选择最大相关度值并输出其所对应的合成信号。
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 true CN110390073A (zh) | 2019-10-29 |
CN110390073B 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) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112071297A (zh) * | 2020-09-07 | 2020-12-11 | 西北工业大学 | 一种矢量声的自适应滤波方法 |
CN113253196A (zh) * | 2021-06-01 | 2021-08-13 | 中国电子科技集团公司第三十六研究所 | 一种多信号测向方法、装置和电子设备 |
Citations (10)
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 | 中国电子科技集团公司第三研究所 | 多目标方位估计方法及水声垂直矢量阵列*** |
-
2019
- 2019-08-19 CN CN201910762326.0A patent/CN110390073B/zh active Active
Patent Citations (10)
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)
Title |
---|
邹吉武等: "海洋矢量噪声场试验研究", 《哈尔滨工程大学学报》 * |
高鑫等: "基于单矢量传感器的强干扰噪声抑制方法", 《传感器与微***》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112071297A (zh) * | 2020-09-07 | 2020-12-11 | 西北工业大学 | 一种矢量声的自适应滤波方法 |
CN112071297B (zh) * | 2020-09-07 | 2023-11-10 | 西北工业大学 | 一种矢量声的自适应滤波方法 |
CN113253196A (zh) * | 2021-06-01 | 2021-08-13 | 中国电子科技集团公司第三十六研究所 | 一种多信号测向方法、装置和电子设备 |
CN113253196B (zh) * | 2021-06-01 | 2024-05-21 | 中国电子科技集团公司第三十六研究所 | 一种多信号测向方法、装置和电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN110390073B (zh) | 2023-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102829770B (zh) | Gps浮标测波方法和测波*** | |
CN112987004B (zh) | 一种浅海环境下基于水平阵列的水面水下目标分类方法 | |
CN203714144U (zh) | 一种基于声学和gps智能定位的浮标装置 | |
CN106772246A (zh) | 基于声阵列的无人机实时检测和定位***及方法 | |
CN103076594B (zh) | 一种基于互相关的水声脉冲信号双阵元定位的方法 | |
CN112098938B (zh) | 一种基于六元锥矢量阵的水声目标降维匹配声场定位方法 | |
CN109283511A (zh) | 一种宽覆盖多波束接收基阵校准方法 | |
CN110390073A (zh) | 一种矢量传感的多通道空间合成方位滤波方法 | |
CN111273265B (zh) | 基于demon谱的单矢量水听器自主测距装置及方法 | |
CN103048642A (zh) | 基于频域最小二乘法的水声脉冲信号匹配场定位方法 | |
CN104569484B (zh) | 一种多输入多输出阵列式超声测风***及测量方法 | |
CN101915922A (zh) | 拖曳线列阵被动测距方法 | |
CN112269163B (zh) | 一种基于坐底单三维矢量水听器的水中声源方位深度协同跟踪方法 | |
CN102043145A (zh) | 基于声矢量传感器均匀直线阵的快速宽带频域波束形成方法 | |
CN108089155A (zh) | 一种深海环境下单水听器声源被动定位方法 | |
CN105629220A (zh) | 一种基于单水听器的深海水声被动测距方法 | |
CN108957462A (zh) | 一种基于平坦海底的多波束水体数据处理方法 | |
CN103076604A (zh) | 一种基于频散特征的低频水声脉冲信号距离的测量方法 | |
CN109061654A (zh) | 一种深海环境下单圆环阵主动三维定位方法 | |
CN101762824B (zh) | 基于单向水声测距的海洋地震拖缆位置测量方法 | |
CN109100799A (zh) | 一种起伏海水表面的电缆深度定位方法及处理终端 | |
CN102830424B (zh) | 一种检波器组合参数计算方法 | |
CN108414984A (zh) | 一种基于二阶干涉的水下目标定位方法 | |
CN106154220A (zh) | L型简化声矢量传感器阵列多参数联合估计四元数方法 | |
CN113153232B (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 |