CN109597021B - 一种波达方向估计方法及装置 - Google Patents
一种波达方向估计方法及装置 Download PDFInfo
- Publication number
- CN109597021B CN109597021B CN201811307881.6A CN201811307881A CN109597021B CN 109597021 B CN109597021 B CN 109597021B CN 201811307881 A CN201811307881 A CN 201811307881A CN 109597021 B CN109597021 B CN 109597021B
- Authority
- CN
- China
- Prior art keywords
- arrival
- wave
- fourier transform
- acoustic
- space
- 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
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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/80—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using ultrasonic, sonic or infrasonic waves
- G01S3/802—Systems for determining direction or deviation from predetermined direction
- G01S3/8027—By vectorial composition of signals received by plural, differently-oriented transducers
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
Description
技术领域
本发明属于信号处理技术领域,具体涉及一种基于空域和时域傅里叶变换的波达方向估计方法及装置。
背景技术
阵列信号处理是信号处理的一个重要分支,在雷达、声纳、通信、地震数据处理和医用成像等众多领域得到广泛的应用。波达方向(DOA)估计是阵列信号处理的一个主要研究内容,在雷达、声纳、通信等领域广泛应用。DOA估计是为了确定信号的方位,从接收数据中测出信号方向。阵列信号处理首先需要借助阵列获取多通道信号,阵列通常由多个传感器按照一定的拓扑结构组成,典型的阵列包括线阵、方阵或圆阵等。
可用于波达方向估计的方法有多种,传统的包括时延估计方法、波束形成估计方法以及极大似然估计方法等。这些传统方法由于一些缺点和局限性,逐渐被较成熟的基于多信号分类(MUSIC)估计方法和Sawada-ICA法所取代。
MUSIC是一种高分辨率的DOA算法,其不足之处在于:估计结果易受外来干扰噪声的影响,在强噪声干扰环境下甚至可能导致错误的估计结果。
Sawada-ICA法是基于时-频域独立分量分析(ICA)的DOA估计方法,其不足之处在于:没有考虑外加噪声分量,导致该方法的实用性较差,需要人为选取频响函数进行估计,效率比较低。
发明内容
本发明所解决的技术问题是:克服现有技术的不足,本发明提出一种波达方向估计方法及装置,可实现二维波达方向估计,测量装置简单,使用灵活、实用性好,可广泛应用在雷达、声纳、通信等领域。
本发明所采用的技术方案是:一种波达方向估计方法,包括步骤如下:
步骤一、以三列线阵两两正交构成传感器阵列,采集被测声源在X、Y、Z轴方向的声信号,并完成声信号的模数转化;X、Y、Z轴分别沿三列线阵方向;
步骤二、将步骤一中采集到的声信号按照X、Y、Z轴划分为u(x,t)、u(y,t)、u(z,t)三个二维时空t-s矩阵,t表示时间,s表示空间;
步骤三、对步骤二中经过划分的三个t-s矩阵分别进行二维傅里叶变换,得到对应的三个频率波数f-k矩阵,f表示频率,k表示波数;
步骤四、对步骤三得到的f-k矩阵进行峰值检测,获得三个f-k矩阵的峰值,峰值对应的波数值分别为kx、ky、kz;
所述步骤一中的传感器阵列由两两正交的线阵构成,根据空间乃奎斯特采样定律,线阵中传感器的间隔小于声波波长的一半,每条线阵长度大于两个声波波长。
所述步骤三中对二维时空t-s矩阵进行二维傅里叶变换指对声信号的空域和时域进行傅里叶变换,其中,对X、Y、Z轴数据的空域和时域傅里叶变换公式分别为:
式中,U(kx,f)、U(ky,f)、U(kz,f)分别表示经过空域和时域傅里叶变换得到X、Y、Z轴方向的f-k矩阵。
所述步骤四中,对f-k矩阵进行峰值检测指对步骤三中的f-k矩阵中阵列采集的声场信息进行频域和波数域滤波,获得来波在X、Y、Z轴方向的波数分量kx、ky、kz。
一种波达方向估计装置,包括:声波传感模块、AD采集模块、数据处理模块、波达方向计算模块;
声波传感模块,包括三列两两正交的传感器阵列,采集被测声源在X、Y、Z轴方向的声信号;
AD采集模块,将声波传感模块采集的模拟信号通过模数转换转变成数字信号;
数据处理模块,将采集到的声信号按照X、Y、Z轴划分为u(x,t)、u(y,t)、u(z,t)三个二维时空t-s矩阵;对经过划分的三个t-s矩阵分别进行二维傅里叶变换,得到对应的三个频率波数f-k矩阵;对f-k矩阵进行峰值检测,获得三个f-k矩阵的峰值和对应的波数值kx、ky、kz;其中,t表示时间,s表示空间,f表示频率,k表示波数;
所述声波传感模块中的传感器阵列由两两正交的线阵构成,根据空间乃奎斯特采样定律,线阵中传感器的间隔小于声波波长的一半,每条线阵长度大于两个声波波长。
所述数据处理模块对二维时空t-s矩阵进行二维傅里叶变换指对声信号的空域和时域进行傅里叶变换,其中,对X、Y、Z轴数据的空域和时域傅里叶变换公式分别为:
式中,U(kx,f)、U(ky,f)、U(kz,f)分别表示经过空域和时域傅里叶变换得到X、Y、Z轴方向的f-k矩阵。
本发明与现有技术相比的有益效果是:
(1)本发明提供了一种新的波达方向估计方法,采用三轴声阵列采集被测信号,可实现二维波达方向估计,计算简单有效,可用于在线实时跟踪,便于应用推广;
(2)本发明采用了基于二维傅里叶变换的阵列信号处理方法,通过频域和波数域滤波,有效滤除外来干扰噪声的影响,波达方向估计结果稳定,精度高,实用性好;
(3)本发明提供了一种新的波达方向估计装置,包括声波传感模块、AD采集模块、数据处理模块、波达方向计算模块,具有装置简单、使用灵活、实用性好的特点,在多个装置组阵使用时,可快速实现声源定位。
(4)本发明中的数据处理模块对采集数据进行二维傅里叶变换操作,该数据处理模块为集成快速傅里叶变换算法的定制硬件,实现数据处理的硬件加速,快速估计声源信号波达方向。
附图说明
图1是本发明提供的波达方向估计流程框图;
图2是本发明的实施示意图;
图3是本发明的波达方向估计装置的系统组成框图。
具体实施方式
下面结合附图对本发明进行进一步说明。
图1为本方法的流程框图,一种波达方向估计方法,包括步骤如下:
步骤一:布设传感器阵列,如图2所示,以三条线阵两两正交构成传感器阵列,根据空间乃奎斯特采样定律,线阵中传感器的间隔小于声波波长的一半,为了充分采集声场空间信息,每条线阵长度大于两个声波波长。传感器阵列同步采集被测声源在X、Y、Z三个方向的声信号,并通过AD采集模块完成声信号的模数转化;
步骤二:将步骤一中采集到的声信号数据按照X、Y、Z轴进行分组,分为u(x,t)、u(y,t)、u(z,t)三个二维时空t-s矩阵,即其中一维为时间采样,另一维为空间采样;t表示时间,s表示空间。x,y,z分别表示X、Y、Z轴上对应的采样点坐标值;
步骤三:对步骤二中经过划分的三个二维矩阵分别进行空域和时域二维傅里叶变换,得到对应的三个频率波数f-k矩阵,对X、Y、Z轴数据的空域和时域傅里叶变换公式分别为:
式中,U(kx,f)、U(ky,f)、U(kz,f)分别表示经过空域和时域傅里叶变换得到X、Y、Z轴方向的f-k矩阵。f表示频率,k表示波数。
步骤四:对步骤三得到的f-k矩阵进行峰值检测,获得三个f-k矩阵的峰值,而峰值对应的波数值分别为kx、ky、kz;
步骤五:将步骤四计算得到的标量波数分量kx、ky、kz合成矢量波矢波矢波矢的方向即为波达方向。进一步,将波达方向分为方位角θ和俯仰角方位角θ定义为波达方向在XY平面内投影与X轴的夹角,俯仰角定义为波达方向与XY平面的夹角,利用公式:
如图3所示,一种波达方向估计装置,包括:声波传感模块、AD采集模块、数据处理模块、波达方向计算模块;
声波传感模块,包括三列两两正交的传感器阵列,采集被测声源在X、Y、Z轴方向的声信号;
AD采集模块,将声波传感模块采集的模拟信号通过模数转换转变成数字信号;
数据处理模块,将采集到的声信号按照X、Y、Z轴划分为u(x,t)、u(y,t)、u(z,t)三个二维时空t-s矩阵;对经过划分的三个t-s矩阵分别进行二维傅里叶变换,得到对应的三个频率波数f-k矩阵;对f-k矩阵进行峰值检测,获得三个f-k矩阵的峰值和对应的波数值kx、ky、kz;其中,t表示时间,f表示频率,k表示波数;
以上所述,仅为本发明最佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (7)
1.一种波达方向估计方法,其特征在于,包括步骤如下:
步骤一、以三列线阵两两正交构成传感器阵列,采集被测声源在X、Y、Z轴方向的声信号,并完成声信号的模数转化;X、Y、Z轴分别沿三列线阵方向;
步骤二、将步骤一中采集到的声信号按照X、Y、Z轴划分为u(x,t)、u(y,t)、u(z,t)三个二维时空t-s矩阵,t表示时间,s表示空间;
步骤三、对步骤二中经过划分的三个t-s矩阵分别进行二维傅里叶变换,得到对应的三个频率波数f-k矩阵,f表示频率,k表示波数;
步骤四、对步骤三得到的f-k矩阵进行峰值检测,获得三个f-k矩阵的峰值,峰值对应的波数值分别为kx、ky、kz;
所述步骤一中的传感器阵列由两两正交的线阵构成,根据空间乃奎斯特采样定律,线阵中传感器的间隔小于声波波长的一半,每条线阵长度大于两个声波波长;
所述步骤三中对二维时空t-s矩阵进行二维傅里叶变换指对声信号的空域和时域进行傅里叶变换,其中,对X、Y、Z轴数据的空域和时域傅里叶变换公式分别为:
式中,U(kx,f)、U(ky,f)、U(kz,f)分别表示经过空域和时域傅里叶变换得到X、Y、Z轴方向的f-k矩阵。
2.根据权利要求1所述的一种波达方向估计方法,其特征在于:所述步骤四中,对f-k矩阵进行峰值检测指对步骤三中的f-k矩阵中阵列采集的声场信息进行频域和波数域滤波,获得来波在X、Y、Z轴方向的波数分量kx、ky、kz。
4.一种波达方向估计装置,其特征在于,包括:声波传感模块、AD采集模块、数据处理模块、波达方向计算模块;
声波传感模块,包括三列两两正交的传感器阵列,采集被测声源在X、Y、Z轴方向的声信号;
AD采集模块,将声波传感模块采集的模拟信号通过模数转换转变成数字信号;
数据处理模块,将采集到的声信号按照X、Y、Z轴划分为u(x,t)、u(y,t)、u(z,t)三个二维时空t-s矩阵;对经过划分的三个t-s矩阵分别进行二维傅里叶变换,得到对应的三个频率波数f-k矩阵;对f-k矩阵进行峰值检测,获得三个f-k矩阵的峰值和对应的波数值kx、ky、kz;其中,t表示时间,s表示空间,f表示频率,k表示波数;
5.根据权利要求4所述的一种波达方向估计装置,其特征在于:所述声波传感模块中的传感器阵列由两两正交的线阵构成,根据空间乃奎斯特采样定律,线阵中传感器的间隔小于声波波长的一半,每条线阵长度大于两个声波波长。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811307881.6A CN109597021B (zh) | 2018-11-05 | 2018-11-05 | 一种波达方向估计方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811307881.6A CN109597021B (zh) | 2018-11-05 | 2018-11-05 | 一种波达方向估计方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109597021A CN109597021A (zh) | 2019-04-09 |
CN109597021B true CN109597021B (zh) | 2021-02-05 |
Family
ID=65957482
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811307881.6A Active CN109597021B (zh) | 2018-11-05 | 2018-11-05 | 一种波达方向估计方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109597021B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112285641B (zh) * | 2020-09-16 | 2023-12-29 | 西安空间无线电技术研究所 | 一种基于ica的波达方向doa的估计方法及装置 |
CN113126068A (zh) * | 2021-04-22 | 2021-07-16 | 浙江大学 | 一种基于波数域滤波的水下航行器自噪声抑制方法 |
CN113561911B (zh) * | 2021-08-12 | 2023-03-17 | 森思泰克河北科技有限公司 | 车辆控制方法、装置、毫米波雷达及存储介质 |
CN114545324B (zh) * | 2022-04-24 | 2022-07-12 | 南京宇安防务科技有限公司 | 一种适用于非均匀阵列的快速测向方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101470187A (zh) * | 2007-12-26 | 2009-07-01 | 中国科学院声学研究所 | 一种用于线列阵的高精度测向方法 |
CN101644773A (zh) * | 2009-03-20 | 2010-02-10 | 中国科学院声学研究所 | 一种实时频域超分辨方位估计方法及装置 |
KR20140015893A (ko) * | 2012-07-26 | 2014-02-07 | 삼성테크윈 주식회사 | 음원 위치 추정 장치 및 방법 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2162065A1 (en) * | 1994-11-08 | 1996-05-09 | Jian-Cheng Zhang | Method of separating predetermined signal components from wireline acoustic array log data |
CN101051083B (zh) * | 2006-04-07 | 2010-11-03 | 中国科学院声学研究所 | 一种二维波达方向估计声纳信号处理方法 |
CN103558584B (zh) * | 2013-10-12 | 2016-08-17 | 中国科学院电子学研究所 | 一种波达方向的检测方法及装置 |
CN103744061A (zh) * | 2014-01-15 | 2014-04-23 | 西安电子科技大学 | 基于迭代最小二乘方法的mimo雷达doa估计方法 |
CN103777190B (zh) * | 2014-02-26 | 2016-03-30 | 南京信息工程大学 | 一种双基地mimo雷达高速高机动目标的角度估计方法 |
US9928213B2 (en) * | 2014-09-04 | 2018-03-27 | Qualcomm Incorporated | Event-driven spatio-temporal short-time fourier transform processing for asynchronous pulse-modulated sampled signals |
CN105204001A (zh) * | 2015-10-12 | 2015-12-30 | Tcl集团股份有限公司 | 一种声源定位的方法及系统 |
CN106950530B (zh) * | 2017-02-22 | 2019-08-06 | 西安电子科技大学 | 基于互协方差的声矢量传感器阵列解相干参数估计方法 |
CN107037398B (zh) * | 2017-06-23 | 2020-03-31 | 电子科技大学 | 一种二维music算法估计波达方向的并行计算方法 |
CN108132457A (zh) * | 2017-12-22 | 2018-06-08 | 景晖 | 一种确定位置的语音来波方向估计方法和装置 |
CN108387875A (zh) * | 2018-03-01 | 2018-08-10 | 南京理工大学 | 基于二维fft的时分mimo雷达相位误差矫正方法 |
CN108683619B (zh) * | 2018-05-18 | 2021-05-14 | 浙江大学城市学院 | 一种低复杂度的大规模mimo信道参数估计方法 |
-
2018
- 2018-11-05 CN CN201811307881.6A patent/CN109597021B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101470187A (zh) * | 2007-12-26 | 2009-07-01 | 中国科学院声学研究所 | 一种用于线列阵的高精度测向方法 |
CN101644773A (zh) * | 2009-03-20 | 2010-02-10 | 中国科学院声学研究所 | 一种实时频域超分辨方位估计方法及装置 |
KR20140015893A (ko) * | 2012-07-26 | 2014-02-07 | 삼성테크윈 주식회사 | 음원 위치 추정 장치 및 방법 |
Also Published As
Publication number | Publication date |
---|---|
CN109597021A (zh) | 2019-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109597021B (zh) | 一种波达方向估计方法及装置 | |
CN107644650B (zh) | 一种基于渐进串行正交化盲源分离算法的改进声源定位方法及其实现系统 | |
CN101957443B (zh) | 声源定位方法 | |
CN110927669B (zh) | 一种用于无线声传感器网络的cs多声源定位方法及系统 | |
CN110534126B (zh) | 一种基于固定波束形成的声源定位和语音增强方法及系统 | |
CN109581275B (zh) | 基于非圆信号和三维正交阵的二维水下doa估计方法和装置 | |
CN109669160B (zh) | 一种水下瞬态声信号的检测方法 | |
CN104407328A (zh) | 基于空间脉冲响应匹配的封闭空间声源定位方法及系统 | |
US7397427B1 (en) | Phase event detection and direction of arrival estimation | |
CN112285647B (zh) | 一种基于稀疏表示与重构的信号方位高分辨估计方法 | |
CN103792513A (zh) | 一种雷声定位系统及方法 | |
CN108398659B (zh) | 一种矩阵束与求根music结合的波达方向估计方法 | |
CN111812581B (zh) | 基于原子范数的球面阵列声源波达方向估计方法 | |
CN110515034B (zh) | 一种声信号方位角测量系统及方法 | |
CN109581274B (zh) | 基于夹角可调三维阵的非圆信号水下doa估计方法和装置 | |
CN109521392B (zh) | 基于非圆信号和l型线阵的水下一维doa估计方法和装置 | |
CN103048641B (zh) | 具有联合约束优化形式的矢量阵稳健聚焦处理方法 | |
CN114563756A (zh) | 双通道相关干涉仪测向样本线性插值方法 | |
CN111175727A (zh) | 一种基于条件波数谱密度的宽带信号方位估计的方法 | |
CN110907892B (zh) | 一种球麦克风阵列语音信号到达角估计方法 | |
CN110208736B (zh) | 基于四阶累量的非圆信号均匀阵列波达方向角估计方法 | |
CN111352075B (zh) | 一种基于深度学习的水下多声源定位方法及系统 | |
US20140269198A1 (en) | Beamforming Sensor Nodes And Associated Systems | |
CN116299182A (zh) | 一种声源三维定位的方法及装置 | |
CN107171748B (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 |