CN108957389A - 一种实数域多通道信号目标方位估计方法 - Google Patents
一种实数域多通道信号目标方位估计方法 Download PDFInfo
- Publication number
- CN108957389A CN108957389A CN201810510960.0A CN201810510960A CN108957389A CN 108957389 A CN108957389 A CN 108957389A CN 201810510960 A CN201810510960 A CN 201810510960A CN 108957389 A CN108957389 A CN 108957389A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- array
- azimuth
- channel signals
- multi channel
- 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.)
- Pending
Links
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/02—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 radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
- G01S3/143—Systems for determining direction or deviation from predetermined direction by vectorial combination of signals derived from differently oriented antennae
-
- 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
Abstract
本发明提供一种实数域多通道信号目标方位估计方法,包括:获得多通道信号的频谱用X(ω),将负频率对应的Xi(ω),i=1,2,…,M的幅度置零,结果用Y(ω)表示,计算多通道信号互谱RY(ω)=E[Y(ω)YH(ω)];通过在互谱上进行常规波束形成,计算方位谱P(θ),P(θ)=aH(θ)RY(ω)a(θ);搜索方位谱的最大值,对应θ值即为估计的目标方位。本发明直接在频域进行波束形成计算,通过避免使用Hilbert变换来降低计算成本,解决传统方案在实际应用中的局限性。
Description
技术领域
本发明涉及一种方位估计领域的方法,尤其涉及一种实数域多通道信号目标方位估计方法。
背景技术
方位估计是阵列信号处理的重要组成部分,在雷达和声纳探测领域具有广泛的应用。现有的方位估计算法往往要求输入信号必须是复数的多通道阵列数据,然而现实的传感器阵列输出多为实数,因此进行方位估计前将实数信号转化为复数信号。对实窄带信号来说,最广泛应用的算法是希尔伯特变换法,该算法能将实窄带信号转化为复的解析信号,解析信号在在频域上与原信号的频谱分布相同,只是没有负频率。实际中的物理信号被传感器接收和测量是都是实数,虽然可以用希尔波特变换将接收或测量的实数信号转化为复数信号,但是希尔伯特本身就是卷积的形式,运算量较大,在实际应用中具有较大的局限性。
发明内容
本发明的目的是为了直接在频域进行波束形成计算,通过避免使用Hilbert变换来降低计算成本,解决传统方案在实际应用中的局限性而提供一种实数域多通道信号目标方位估计方法。
本发明的目的是这样实现的:步骤如下:
步骤一:令表示传感器阵列输出的多通道时域信号,其中xi(t),i=1,…M表示传感器阵列第i个阵元输出的信号,M为阵列的阵元数;对x(t)进行N点FFT计算,获得多通道信号的频谱用X(ω):
其中:FFT[·]表示对中括号内的矩阵每一行进行快速傅里叶变换,Xi(ω),i=1,…M为第i阵元输出信号对应的频谱;
步骤二:将负频率对应的Xi(ω),i=1,2,…,M的幅度置零,结果用Y(ω)表示,即:
步骤三:计算多通道信号互谱:
RY(ω)=E[Y(ω)YH(ω)]
其中:Rz(ω)表示M×M维的互谱矩阵,E[·]表示对括号内的矩阵求期望,H代表矩阵共轭转置;
步骤四:通过在互谱上进行常规波束形成,计算方位谱P(θ):
P(θ)=aH(θ)RY(ω)a(θ)
其中:为阵列在θ方向的响应向量,k(θ)表示波数矢量,r代表阵列阵元的坐标矢量;
步骤五:搜索方位谱的最大值,对应θ值即为估计的目标方位。
本发明还包括这样一些结构特征:
1.步骤五还可以是:在[0°,360°)的水平角度内等按照一定的角度间隔计算P(θ),然后通过插值和峰值搜索来估计目标的方位。
与现有技术相比,本发明的有益效果是:本发明能够提供了一种新的实数域多通道信号目标方位估计方法,在进行波数形成前,不必进行Hilbert变换,也即不需要将多通道实数信号转变到复数域,直接在互谱域进行波数形成即可,解决了常规传感器阵列输出为实数而波束形成或方位估计算法需要复数的矛盾问题,另外发明的算法具有计算成本低的优点,易于在实际中进行推广应用。
附图说明
图1是本发明算法步骤的流程图;
图2是均匀线阵模型;
图3是均匀圆阵模型;
图4是在目标方位角为60°、信噪比5dB、工作频率800Hz情况下均匀线阵的方位谱;
图5是在目标方位角为60°、信噪比5dB、工作频率800Hz情况下均匀圆阵的方位谱。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述。
结合图1至图5,本发明不需要将阵列输出的多通道实数转变到复数域,直接在频域互谱域进行波束形成,进而对目标进行方位估计的方法;其具体步骤如下:
步骤一:令
表示传感器阵列输出的多通道时域信号。其中xi(t),i=1,…M表示传感器阵列第i个阵元输出的信号,M为阵列的阵元数。对x(t)进行N点FFT计算,获得多通道信号的频谱,用X(ω)表示,即:
其中FFT[·]表示对中括号内的矩阵每一行进行快速傅里叶变换,Xi(ω),i=1,…M为第i阵元输出信号对应的频谱。
步骤二:将负频率对应的Xi(ω),i=1,2,…,M的幅度置零,结果用Y(ω)表示,即:
步骤三:计算多通道信号互谱,即:
RY(ω)=E[Y(ω)YH(ω)] (4)
其中Rz(ω)表示M×M维的互谱矩阵,E[·]表示对括号内的矩阵求期望,H代表矩阵共轭转置。
步骤四:通过在互谱上进行常规波束形成,计算方位谱P(θ),即:
P(θ)=aH(θ)RY(ω)a(θ) (5)
其中为阵列在θ方向的响应向量,k(θ)表示波数矢量,r代表阵列阵元的坐标矢量。
步骤五:搜索方位谱的最大值,对应θ值即为估计的目标方位。可选地,也可在[0°,360°)的水平角度内等按照一定的角度间隔计算P(θ),然后通过插值和峰值搜索来估计目标的方位。
为验证算法有效性,结合具体参数给出本发明的两个具体实施例。
实施例一:阵元数M=12的均匀线阵,阵元间距d为半波长,如图2为均匀线阵模型,目标方位θ0=60°,信噪比为SNR=5dB,信号频率fc=0.8kHz,声速c=1500m/s,采样频率fs=5kHz。具体步骤如下:
步骤1:获取阵列输出的时域数据,接收模型为:
其中,s(t)=cos(2πfct)为目标信号,θ0位目标的期望方向,a(θ0)目标信号的方向矢量,n(t)代表加性高斯白噪声,si(t),i=1,2…,M为每个阵元接收到的目标信号,ni(t),i=1,2…,M为每个阵元接收到的噪声信号。
步骤2:对阵列输出的时域信号进行N点FFT变换,N=1024,即:
步骤3:对负频率对应的频谱进行置零处理,即:
步骤4:计算互谱矩阵
RY(ω)=E[Y(ω)YH(ω)]
步骤5:通过在互谱上进行常规波束形成,计算方位谱P(θ),即:
P(θ)=aH(θ)Rz(ω)a(θ) (5)
其中为方向矢量,如图4是目标方位角为60°,信噪比5dB,工作频率800Hz情况下均匀线阵的方位谱。
步骤6:寻找P(θ)的最大值,对应θ值即为目标的方位,另外也可通过预成角度的波数输出结合插值法计算目标的方位。
实施例二:阵元数M=12的均匀圆阵,阵元间距d为半波长,如图3为均匀线阵模型,目标方位θ0=60°,信噪比为SNR=10dB,信号频率fc=0.8kHz,声速c=1500m/s,采样频率fs=5kHz。具体步骤如下:
步骤1:获取阵列输出的时域数据,接收模型为:
其中,s(t)=cos(2πfct)为目标信号,θ0位目标的期望方向,a(θ0)目标信号的方向矢量,n(t)代表加性高斯白噪声,si(t),i=1,2…,M为每个阵元接收到的目标信号,ni(t),i=1,2…,M为每个阵元接收到的噪声信号。
步骤2:对阵列输出的时域信号进行N点FFT变换,N=1024,即:
步骤3:对负频率对应的频谱进行置零处理,即:
步骤4:计算互谱矩阵
RY(ω)=E[Y(ω)YH(ω)]
步骤5:通过在互谱上进行常规波束形成,计算方位谱P(θ),即:
P(θ)=aH(θ)Rz(ω)a(θ) (5)
其中为方向矢量,其中R为圆阵半径,如图5是目标方位角为60°,信噪比5dB,工作频率0.8kHz情况下均匀圆阵的方位谱。
步骤6:寻找P(θ)的最大值,对应θ值即为目标的方位,另外也可通过预成角度的波数输出结合插值法计算目标的方位。
综上,方位估计是阵列信号处理的重要组成部分,在雷达和声纳探测领域具有广泛的应用。现实的阵列输出和测量都是实数,无法直接用于方位估计算法,对于窄带信号常用的方法就是希尔伯特变换,但是运算量较大,实际使用有局限性。本发明要解决的技术问题是提供一种实数域的方位估计方法,直接在频域进行波束形成计算,通过避免使用Hilbert变换来降低计算成本,解决传统方案在实际应用中的局限性。
Claims (2)
1.一种实数域多通道信号目标方位估计方法,其特征在于:步骤如下:
步骤一:令表示传感器阵列输出的多通道时域信号,其中xi(t),i=1,…M表示传感器阵列第i个阵元输出的信号,M为阵列的阵元数;对x(t)进行N点FFT计算,获得多通道信号的频谱用X(ω):
其中:FFT[·]表示对中括号内的矩阵每一行进行快速傅里叶变换,Xi(ω),i=1,…M为第i阵元输出信号对应的频谱;
步骤二:将负频率对应的Xi(ω),i=1,2,…,M的幅度置零,结果用Y(ω)表示,即:
步骤三:计算多通道信号互谱:
RY(ω)=E[Y(ω)YH(ω)]
其中:Rz(ω)表示M×M维的互谱矩阵,E[·]表示对括号内的矩阵求期望,H代表矩阵共轭转置;
步骤四:通过在互谱上进行常规波束形成,计算方位谱P(θ):
P(θ)=aH(θ)RY(ω)a(θ)
其中:为阵列在θ方向的响应向量,k(θ)表示波数矢量,r代表阵列阵元的坐标矢量;
步骤五:搜索方位谱的最大值,对应θ值即为估计的目标方位。
2.根据权利要求1所述的一种实数域多通道信号目标方位估计方法,其特征在于:步骤五还可以是:在[0°,360°)的水平角度内等按照一定的角度间隔计算P(θ),然后通过插值和峰值搜索来估计目标的方位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810510960.0A CN108957389A (zh) | 2018-05-25 | 2018-05-25 | 一种实数域多通道信号目标方位估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810510960.0A CN108957389A (zh) | 2018-05-25 | 2018-05-25 | 一种实数域多通道信号目标方位估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108957389A true CN108957389A (zh) | 2018-12-07 |
Family
ID=64492123
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810510960.0A Pending CN108957389A (zh) | 2018-05-25 | 2018-05-25 | 一种实数域多通道信号目标方位估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108957389A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109856627A (zh) * | 2019-01-11 | 2019-06-07 | 中国船舶重工集团公司第七二四研究所 | 基于空间补零插值的三角阵相控阵2d-fft测向方法 |
CN109975743A (zh) * | 2019-04-12 | 2019-07-05 | 哈尔滨工程大学 | 一种线谱目标的互质阵互谱测向方法 |
CN112649787A (zh) * | 2020-11-30 | 2021-04-13 | 海鹰企业集团有限责任公司 | 一种基于低频圆环阵的目标方位估计方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644773A (zh) * | 2009-03-20 | 2010-02-10 | 中国科学院声学研究所 | 一种实时频域超分辨方位估计方法及装置 |
CN103364772A (zh) * | 2013-07-14 | 2013-10-23 | 西安电子科技大学 | 基于实数域广义多重信号分类算法的目标低仰角估计方法 |
-
2018
- 2018-05-25 CN CN201810510960.0A patent/CN108957389A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101644773A (zh) * | 2009-03-20 | 2010-02-10 | 中国科学院声学研究所 | 一种实时频域超分辨方位估计方法及装置 |
CN103364772A (zh) * | 2013-07-14 | 2013-10-23 | 西安电子科技大学 | 基于实数域广义多重信号分类算法的目标低仰角估计方法 |
Non-Patent Citations (2)
Title |
---|
林旺生: "基于矢量声场的水下被动探测与定位技术研究", 《中国博士学位论文全文数据库(电子期刊) 工程科技II辑》 * |
高云超: "希尔伯特-黄变换在水声信号处理中的应用研究", 《中国博士学位论文全文数据库(电子期刊) 信息科技辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109856627A (zh) * | 2019-01-11 | 2019-06-07 | 中国船舶重工集团公司第七二四研究所 | 基于空间补零插值的三角阵相控阵2d-fft测向方法 |
CN109975743A (zh) * | 2019-04-12 | 2019-07-05 | 哈尔滨工程大学 | 一种线谱目标的互质阵互谱测向方法 |
CN109975743B (zh) * | 2019-04-12 | 2020-11-03 | 哈尔滨工程大学 | 一种线谱目标的互质阵互谱测向方法 |
CN112649787A (zh) * | 2020-11-30 | 2021-04-13 | 海鹰企业集团有限责任公司 | 一种基于低频圆环阵的目标方位估计方法 |
CN112649787B (zh) * | 2020-11-30 | 2022-04-12 | 海鹰企业集团有限责任公司 | 一种基于低频圆环阵的目标方位估计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108375763B (zh) | 一种应用于多声源环境的分频定位方法 | |
CN110007148B (zh) | 一种基于离散频谱相位和幅值综合内插的单频信号频率估计方法 | |
CN109799495A (zh) | 一种用于高保真阵列处理的宽带时延估计方法 | |
Zhang et al. | Multitarget AOA estimation using wideband LFMCW signal and two receiver antennas | |
CN108957389A (zh) | 一种实数域多通道信号目标方位估计方法 | |
CN101644773A (zh) | 一种实时频域超分辨方位估计方法及装置 | |
CN104502904B (zh) | 一种鱼雷自导波束锐化方法 | |
CN107271955B (zh) | 一种宽带线性调频信号的时差和尺度差估计方法 | |
CN106646350A (zh) | 一种单只矢量水听器各通道幅度增益不一致时的修正方法 | |
CN105005038A (zh) | 一种改进的声矢量阵相干源doa估计算法 | |
CN105158751A (zh) | 一种声矢量阵快速方位估计方法 | |
CN113281727A (zh) | 一种基于水平线列阵的输出增强的波束形成方法及其系统 | |
CN110196407B (zh) | 一种基于频率预估的单矢量水听器信号来波方向估计方法 | |
CN109405954A (zh) | 一种基于矢量水听器阵列的uuv移动水声通信技术 | |
CN109932679B (zh) | 一种传感器列系统最大似然角度分辨率估计方法 | |
CN111175727B (zh) | 一种基于条件波数谱密度的宽带信号方位估计的方法 | |
CN108089146A (zh) | 一种对预估角误差鲁棒的高分辨宽带波达方向估计方法 | |
CN109541526A (zh) | 一种利用矩阵变换的圆环阵方位估计方法 | |
CN109814065B (zh) | 基于相位因子加权的波束形成方法 | |
Chen et al. | Insight into split beam cross-correlator detector with the prewhitening technique | |
Zhao et al. | Passive estimation of target velocity based on cross-spectrum histogram | |
CN111722178B (zh) | 一种基于指向性模型数值求解的远场窄带信号来波方向估计方法 | |
KR101801325B1 (ko) | 가상 채널 기반의 레이더 장치 및 그것을 이용한 가상 채널 생성 방법 | |
CN110231590B (zh) | 一种基于dft变换的阵列目标角度测量方法 | |
CN113075633A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181207 |