CN103558584B - 一种波达方向的检测方法及装置 - Google Patents

一种波达方向的检测方法及装置 Download PDF

Info

Publication number
CN103558584B
CN103558584B CN201310475827.3A CN201310475827A CN103558584B CN 103558584 B CN103558584 B CN 103558584B CN 201310475827 A CN201310475827 A CN 201310475827A CN 103558584 B CN103558584 B CN 103558584B
Authority
CN
China
Prior art keywords
arrival
signal
pulse compression
distance
antenna
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
CN201310475827.3A
Other languages
English (en)
Other versions
CN103558584A (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.)
Institute of Electronics of CAS
Original Assignee
Institute of Electronics of CAS
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 Institute of Electronics of CAS filed Critical Institute of Electronics of CAS
Priority to CN201310475827.3A priority Critical patent/CN103558584B/zh
Publication of CN103558584A publication Critical patent/CN103558584A/zh
Priority to EP14169310.1A priority patent/EP2860549B1/en
Application granted granted Critical
Publication of CN103558584B publication Critical patent/CN103558584B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904SAR modes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Direction-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/02Direction-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/14Systems for determining direction or deviation from predetermined direction
    • G01S3/46Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques

Abstract

本发明公开了一种波达方向的检测方法,对俯仰向各通道得到的回波数据做距离向脉冲压缩,得到脉冲压缩后的信号;对脉冲压缩后的信号在每个距离门上沿通道做傅立叶变换;对傅立叶变换后的结果进行插值,获得插值后频域幅值最大点处的频率值,根据所述频率值获得波达方向,本发明同时还公开了一种波达方向的检测装置。通过该方案能够减少波达方向检测时的运算量,并且克服目标海拔高度对波达方向检测的影响。

Description

一种波达方向的检测方法及装置
技术领域
本发明涉及合成孔径雷达(SAR,Synthetic Aperture Radar)技术领域,尤其涉及一种波达方向(DOA,Direction of Arrive)的检测方法及装置。
背景技术
SAR是一种高分辨率主动式微波成像雷达,能够使用较短的天线获得较高的距离向和方位向分辨率。数字波束合成(DBF,Digital Beam Forming)接收技术对于合成孔径雷达领域实现更高分辨率和更宽测绘带具有十分重要的应用价值,而在数字波束合成技术中波达方向的检测是影响数字波束合成效果的关键,因此,精确检测波达方向对于数字波束合成接收技术具有重要的意义。
目前,主要的波达方向检测技术都是基于统计的方法实现,主要有波束合成(beam-former)算法,capon算法和多重信号分类(MUSIC:Multiple SignalClassification)算法。一般来说,capon算法和MUSIC算法检测精度较高,但这三种算法都需要大量的数据做统计分析,运算量大,不利于实际工程实现。
另外,在实际的波达方向检测中,目标海拔高度的不同往往会带来指向误差,使得波达方向的检测出现偏差,从而使得SAR接收增益降低。
综上所述,如何在尽量减少运算量并克服目标海拔高度影响的情况下实现准确、快速的波达方向检测是目前亟待解决的问题。
发明内容
本发明实施例的主要目的在于提供一种波达方向的检测方法和装置,能够减少波达方向检测时的运算量,并克服目标海拔高度对波达方向检测的影响。
本发明的技术方案是这样实现的:
本发明实施例提供了一种波达方向的检测方法,所述方法包括:
对俯仰向各通道得到的回波数据做距离向脉冲压缩得到脉冲压缩后的信号;对脉冲压缩后的信号在每个距离门上沿通道做傅立叶变换;对傅立叶变换后的结果进行插值,获得插值后频域幅值最大点处的频率值,根据所述频率值获得波达方向。
上述方案中,所述脉冲压缩后的信号为
其中,τ是方位向慢时间,t是距离向快时间,σ是目标的后向散射系数,pr(i)是脉冲压缩后的点目标信号包络,t0为回波的双程延时,RT(τ)为发射天线至目标的瞬时斜距,α(t0)为波达方向的法向偏移角,dn为第n个接收天线相对于发射天线的位置值,dn的绝对值等于第n个接收天线与发射天线间的距离,λ为波长。
上述方案中,所述对脉冲压缩后得到的信号在每个距离门上沿通道作傅里叶变换包括:
在一个固定的距离门上,删除通道间的相同项,得到N维向量ν(n):
v ( n ) = exp [ j 2 π λ d n sin α ( t 0 ) ]
其中,t0为回波的双程延时,α(t0)为波达方向的法向偏移角,也就是波达方向,dn为第n个接收天线相对于发射天线的位置值,dn的绝对值等于第n个接收天线与发射天线间的距离,λ为波长;
将ν(n)变换到频域得到结果Γ(ξ):
Γ ( ξ ) = ∫ 0 L E exp { j 2 π sin [ α ( t 0 ) ] λ d n } exp [ - j 2 πξ d n ] d ( d n ) = = sin c { L E [ ξ - sin [ α ( t 0 ) ] λ ] } exp [ - jπξ L E ]
其中LE为方位向天线长度,ξ为频域参数值,t0为回波的双程延时,α(t0)为波达方向的法向偏移角,dn为第n个接收天线相对于发射天线的位置值,dn的绝对值等于第n个接收天线与发射天线间的距离,λ为波长。
上述方案中,所述插值后频域幅值最大点处的频率值ξmax为:其中,λ为波长,α(t0)为波达方向的法向偏移角。
具体的,所述Γ(ξ)频率范围为:其中d为相邻接收通道间距。
具体的,所述波达方向的法向偏移角α(t0)的绝对值小于其中,λ为波长,d为相邻接收通道间距。
上述方案中,所述根据所述频率值获得波达方向为获得波达方向的法向偏移角α(t0):α(t0)=arcsin(λξmax),其中,λ为波长,ξmax为插值后频域幅值最大点处的频率值。
上述方案中,所述方法还包括:根据波达方向的法向偏移角α(t0)获得各通道当前距离门的权值
其中,dn为第n个接收天线相对于发射天线的位置值,dn的绝对值等于第n个接收天线与发射天线间的距离,ξmax为插值后频域幅值最大点处的频率值;
根据所述权值进行数字波束DBF合成。
进一步的,所述根据所述权值进行DBF合成为:
以各通道权值对所述脉冲压缩后的信号进行加权求和,获得DBF合成后的信号G(τ,t):
其中,τ是方位向慢时间,t是距离向快时间。
本发明实施例还提供了一种波达方向的检测装置,所述装置包括:脉冲压缩模块、傅里叶变换模块、插值模块、以及波达方向获取模块,其中,
所述脉冲压缩模块,用于对俯仰向各通道得到的回波数据做距离向脉冲压缩,得到脉冲压缩后的信号;
所述傅里叶变换模块,用于对脉冲压缩后得到的信号在每个距离门上沿通道做傅里叶变换;
所述插值模块,用于对傅立叶变换后的结果进行插值;
所述波达方向获取模块,用于获得插值后频域幅值最大点处的频率值,并根据所述频率值获得波达方向。
进一步的,所述装置还包括:权值获取模块和波束合成模块,其中,
所述权值获取模块,用于根据波达方向获得各通道当前距离门的权值;
所述波束合成模块,用于以各通道权值对所述脉冲压缩后的信号进行加权求和,获得DBF合成后的信号。
本发明实施例所提供的波达方向的检测方法和装置,对俯仰向各通道得到的回波数据做距离向脉冲压缩,得到脉冲压缩后的信号;对脉冲压缩后的信号在每个距离门上沿通道做傅立叶变换;对傅立叶变换后的结果进行插值,获得插值后频域幅值最大点处的频率值,根据所述频率值能够获得波达方向;如此,能够减少波达方向检测时的运算量,并克服目标海拔高度对波达方向检测的影响。
附图说明
图1为本发明实施例提供的一种波达方向的检测方法的流程示意图;
图2为本发明实施例提供的一种波达方向的检测装置的结构示意图;
图3为本发明实施例提供的实验一效果对比图;
图4为本发明实施例提供的实验二效果对比图;
图5为本发明实施例提供的实验三效果对比图。
具体实施方式
本发明实施例中,对俯仰向各通道得到的回波数据做距离向脉冲压缩,得到脉冲压缩后的信号;对脉冲压缩后的信号在每个距离门上沿通道做傅立叶变换;对傅立叶变换后的结果进行插值,获得插值后频域幅值最大点处的频率值,根据所述频率值获得波达方向。
下面通过附图及具体实施例对本发明做进一步的详细说明。
本发明实施例实现一种波达方向的检测方法,如图1所示,该方法包括以下几个步骤:
步骤101:对俯仰向各通道得到的回波数据做距离向脉冲压缩,得到脉冲压缩后的信号;
这里,所述回波数据的生成具体包括:雷达发射脉冲信号后,由俯仰向均匀排列的接收阵列天线接收回波信号,回波信号进入每一个天线后经过低噪声放大器处理、混频、AD采样后生成回波数据。
具体的,所述对俯仰向各通道的回波数据做距离向脉冲压缩后得到的信号如下式所示:
其中,τ是方位向慢时间,t是距离向快时间,σ是目标的后向散射系数,pr(i)是脉冲压缩后的点目标信号包络,t0为回波的双程延时,RT(τ)为发射天线至目标的瞬时斜距,α(t0)为波达方向的法向偏移角,即为本发明实施例所要检测的波达方向,dn为第n个接收天线相对于发射天线的位置值,dn的绝对值等于第n个接收天线与发射天线之间的距离,λ为波长。
在实际的应用中,式(1)中的具体参数值与雷达脉冲的参数以及雷达与目标之间的几何关系有关。
具体的,这些数据的处理可以是由SAR系统内部的信号处理器实现或者经由SAR将数据通过宽带数据链等方式发送到地面接收站后在地面接收站中的信号处理器中进行处理。
步骤102:对脉冲压缩后的信号在每个距离门上沿通道做傅立叶变换;
在这一步骤中,所述对脉冲压缩后得到的信号在每个距离门上沿通道作傅里叶变换包括:在一个固定的距离门上,删除通道间的相同项,可以得到一个N维的向量ν(n),所述ν(n)具体为:
v ( n ) = exp [ j 2 π λ d n sin α ( t 0 ) ]
其中,t0为回波的双程延时,α(t0)为波达方向的法向偏移角,即为本发明实施例所要检测的波达方向,dn为第n个接收天线相对于发射天线的位置值,dn的绝对值等于第n个接收天线与发射天线间的距离,λ为波长。
将上述N维向量经傅里叶变换变换到频域,得到结果Γ(ξ)如下:
Γ ( ξ ) = ∫ 0 L E exp { j 2 π sin [ α ( t 0 ) ] λ d n } exp [ - j 2 πξ d n ] d ( d n ) = = sin c { L E [ ξ - sin [ α ( t 0 ) ] λ ] } exp [ - jπξ L E ]
其中,LE为方位向天线长度,ξ为频域参数值,t0为回波的双程延时,α(t0)为波达方向的法向偏移角,即为本发明实施例所要检测的波达方向,dn为第n个接收天线相对于发射天线的位置值,dn的绝对值等于第n个接收天线与发射天线间的距离,λ为波长。
下面通过一个具体示例对上述对脉冲压缩后的信号在每一个距离门上沿通道做傅里叶变换的具体实现方式进行介绍:
假如,距离向有10个接收通道,设某一接收窗得到的数据经脉冲压缩后的信号分别为S1(t)、S2(t)...S10(t),将10个通道的信号写成矩阵的形式如下:
S matri ( t ) = S 1 ( t ) S 2 ( t ) . . . S 10 ( t ) - - - ( 4 )
则,具体的,在每个距离门上沿通道对上述脉冲压缩后的信号做傅立叶变换时可以采用离散傅里叶变换的快速算法快速傅里叶变换(FFT)对矩阵Smatri(t)的每一列做10点的FFT。
步骤103:对傅立叶变换后的结果进行插值,获得插值后频域幅值最大点处的频率值,根据所述频率值获得波达方向;
在步骤102中通过式(3)将N个通道同一距离门上经过脉冲压缩后得到的信号做傅里叶变换后,易知Γ(ξ)的频域参数值ξ的取值范围为:
ξ ∈ [ - 1 2 d , 1 2 d ] - - - ( 5 )
其中d为相邻接收通道间距。
且从式(3)可以看出当时,Γ(ξ)取得最大值,其中,λ为波长,α(t0)为波达方向的法向偏移角,即为本发明实施例所要检测的波达方向。所以将傅里叶变换之后的频域结果插值后可以获得使频域幅值最大的点所对应的频域参数值,根据所述频域参数值可以获得使得接收增益最大的接收方向α(t0),即为波达方向,具体的,按照式(6)可以得出sin[α(t0)]=λξmax,通过对λξmax作反正弦变换即可获得波达方向α(t0):
α(t0)=arcsin(λξmax) (7)
其中,λ为波长,ξmax为插值后频域幅值最大点处的频率值,即:使Γ(ξ)取得最大值时所对应的频域参数值。
具体的,上述频域插值可以在傅里叶变换之前,通过时域补充零值样点的方法实现;或者,在傅里叶变换过程中,通过加密频域间隔实现。目前,频域插值算法比较多,在实际的应用中这一步骤可以采用各种频域插值算法实现,这里不作限制。
下面对本发明实施例所述方法能够获得的波达方向的范围进行分析。
本发明实施例所述方法主要是对回波数据进行脉冲压缩之后,再执行傅里叶变换,由于经过傅里叶变换后得到的信号Γ(ξ)在频域中是周期信号,因此,能够获得的波达方向将被限制在一个有限的范围之内,由式(5)可得波达方向的法向偏移角的正弦值范围为:
| sin [ α ( t 0 ) ] | ≤ λ 2 d - - - ( 8 )
其中,λ为波长,d为相邻接收通道间距。
所以波达方向的法向偏移角绝对值必须小于这个限制其实限制了测绘带的范围,而且该条件与测绘带内不出现栅瓣的条件一致。
可见,通过本发明实施例的实现方案可以将获得的波达方向限制在一定范围之内,这也将直接限制测绘带的范围,使得对波达方向的检测更加准确,这是因为:如果测绘带的范围太大,将使得波达方向的范围超出系统检测的能力范围,从而使得波达方向检测结果误差较大。
进一步的,本发明实施例还提供了依据上述求得的波达方向进行数字波束合成(DBF)的方法,该方法包括:
首先,根据步骤103中获得的波达方向,获得各通道当前距离门的权值如下:
其中,dn为第n个接收天线相对于发射天线的位置值,dn的绝对值等于第n个接收天线与发射天线间的距离,ξmax为插值后频域幅值最大点处的频率值。
之后,将各通道当前距离门的权值代入式(1)中所得到的信号中进行加权求和,得到DBF合成后的信号为:
其中,τ是方位向慢时间,t是距离向快时间。
通过上述方法,本发明实施例能够实现波达方向的检测、DBF权值的获取和数字波束的合成。
为了实现上述方法,本发明实施例还提供了一种波达方向的检测装置,如图2所示,所述装置包括:脉冲压缩模块21、傅里叶变换模块22、插值模块23以及波达方向获取模块24,其中,
所述脉冲压缩模块21,用于对俯仰向各通道得到的回波数据做距离向脉冲压缩,得到脉冲压缩后的信号;
所述傅里叶变换模块22,用于对脉冲压缩后的信号在每一个距离门上沿通道做傅里叶变换;
所述插值模块23,用于对傅立叶变换后的结果做插值;
所述波达方向获取模块24,用于获得插值后频域幅值最大点处的频率值,并根据所述频率值获得波达方向。
具体的,在上述装置中,所述脉冲压缩模块21对俯仰向各通道得到的回波数据做距离向脉冲压缩后得到的信号如式(1)所示。
具体的,所述傅里叶变换模块22对脉冲压缩后得到的信号在每一个距离门上沿通道做傅里叶变换,包括:
在一个固定的距离门上,略去通道间的相同项,获得一个N维的向量ν(n),所述ν(n)具体见式(2);
将上述N维向量ν(n)经傅里叶变换变换到频域,得到结果Γ(ξ),其中Γ(ξ)如(3)式所示;
具体的,所述插值模块23对傅立叶变换后的结果做插值,可以在傅里叶变换之前,通过时域补充零值样点的方法实现;或者,在傅里叶变换过程中,通过加密频域间隔实现。目前,频域插值算法比较多,在实际的应用中这一步骤可以采用各种频域插值算法实现,这里不作限制。
具体的,所述波达方向获取模块24,用于获得插值后频域幅值最大点处的频率值,并根据所述频率值获得波达方向,包括:
将傅里叶变换之后的频域结果插值后获得使频域幅值最大的点所对应点处的频率值,根据所述频率值获得使得接收增益最大的接收方向α(t0),也就是波达方向,其中α(t0)如式(7)所示。
为了在获得波达方向之后实现DBF合成,上述装置还可以进一步包括:权值获取模块25和波束合成模块26,其中
所述权值获取模块25,用于根据波达方向获得各通道当前距离门的权值;
所述波束合成模块26,用于以各通道当前距离门的权值对各通道回波数据进行加权求和获得DBF合成后的信号。
具体的,所述权值获取模块25根据波达方向获得各通道当前距离门的权值,包括:
根据获得的波达方向,获得各通道当前距离门的权值具体方法如(9)式所示;
具体的,所述波束合成模块26以各通道当前距离门的权值对各通道回波数据进行加权求和获得DBF合成后的信号,包括:
以各通道权值进行加权求和,得到DBF合成后的信号G(τ,t)如(10)式。
下面通过三种条件下的仿真实验结果来对本发明实施例的技术效果进行说明。
实验一
首先假设在测绘带内的每个距离门(以三个为例)上分别只有一个点目标,并且各目标均位于不同的海拔高度上,在这种情况下距离门斜距分别为:556Km、584Km、612Km,地距分别为:190Km、256Km、310Km,利用本发明实施例所述方法获得波达方向后得出新的权值,并且利用新的权值加权合成后的信号如图3(b)所示,不获得波达方向而直接DBF合成的结果如图3(a)。从图3可以看出,相较于不获得波达方向而直接DBF合成的结果,利用本发明实施例所述方法能够使每个距离门上的接收增益显著增强。
实验二
假设在测绘带内的每个距离门(以三个为例)上均有多个点目标,并且每个点目标均位于不同的海拔高度上,其中各点目标的海拔高度、斜距、地距各参数如表1所示:
目标 海拔 斜距(距离) 地距
P1 2000m 556Km 190Km
P2 4000m 584Km 256Km
P21 0m 584Km 256Km
P22 3000m 584Km 256Km
P3 6000m 612Km 310Km
P31 1000m 612Km 310Km
表1
利用本发明实施例所述方法获得波达方向后得出新的权值,并且利用新的权值加权合成后的信号如图4(b)所示,不通过检测获得波达方向而直接DBF合成的结果如图4(a)所示。从图4可以看出,本发明实施例在同一距离门上有多个目标时仍然可以有效的克服海拔差带来的增益损失。
实验三
加入分布于不同方位位置的点目标,在每个距离门上均分布有多个点目标,这些目标位于不同的方位位置和海拔高度,其中,在每个距离门上设置500个点目标,点目标的海拔在0-6000m之间均匀分布,点目标在方位向上的位置为:在天线波束照射的范围内均匀分布。仿真结果如图5所示。图5(a)为原始DBF的结果,图5(b)为利用本发明实施例检测出波达方向后进行DBF的结果。从图5可以看出,本发明实施例所述方案在考虑不同方位位置的目标时仍然具有很强的适用性。
通过上述实验及仿真结果可以看出,本发明实施例所述方案不仅能够准确、快速的获得波达方向,还能有效的克服海拔高度引起的波达方向检测误差,利用检测获得的波达方向更新权值之后进行DBF合成,可以克服海拔高度对DBF合成效果的影响。
以上所述,仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。

Claims (11)

1.一种波达方向的检测方法,其特征在于,所述方法包括:
对俯仰向各通道得到的回波数据做距离向脉冲压缩,得到脉冲压缩后的信号;对脉冲压缩后的信号在每个距离门上沿通道做傅立叶变换;对傅立叶变换后的结果进行插值,获得插值后频域幅值最大点处的频率值,根据所述频率值获得波达方向;
其中,所述对傅立叶变换后的结果进行插值在所述傅里叶变换之前,通过时域补充零值样点的方法实现;或者,在所述傅里叶变换过程中,通过加密频域间隔实现。
2.根据权利要求1所述的方法,其特征在于,所述脉冲压缩后的信号为
其中,τ是方位向慢时间,t是距离向快时间,σ是目标的后向散射系数,pr(□)是脉冲压缩后的点目标信号包络,t0为回波的双程延时,RT(τ)为发射天线至目标的瞬时斜距,α(t0)为波达方向的法向偏移角,dn为第n个接收天线相对于发射天线的位置值,dn的绝对值等于第n个接收天线与发射天线间的距离,λ为波长。
3.根据权利要求1所述的方法,其特征在于,所述对脉冲压缩后得到的信号在每个距离门上沿通道作傅里叶变换包括:
在一个固定的距离门上,删除通道间的相同项,得到N维向量ν(n):
v ( n ) = exp [ j 2 π λ d n s i n α ( t 0 ) ]
其中,t0为回波的双程延时,α(t0)为波达方向的法向偏移角,dn为第n个接收天线相对于发射天线的位置值,dn的绝对值等于第n个接收天线与发射天线间的距离,λ为波长;
将ν(n)变换到频域得到结果Γ(ξ):
Γ ( ξ ) = ∫ 0 L E exp { j 2 π sin [ α ( t 0 ) ] λ d n } exp [ - j 2 πξd n ] d ( d n ) = = sin c { L E [ ξ - sin [ α ( t 0 ) ] λ ] } exp [ - jπξL E ]
其中LE为方位向天线长度,ξ为频域参数值,t0为回波的双程延时,α(t0)为波达方向的法向偏移角,dn为第n个接收天线相对于发射天线的位置值,dn的绝对值等于第n个接收天线与发射天线间的距离,λ为波长。
4.根据权利要求1至3任一项所述的方法,其特征在于,所述插值后频域幅值最大点处的频率值ξmax为:其中,λ为波长,α(t0)为波达方向的法向偏移角。
5.根据权利要求3所述的方法,其特征在于,所述Γ(ξ)频率范围为:其中d为相邻接收通道间距。
6.根据权利要求3所述的方法,其特征在于,所述波达方向的法向偏移角α(t0)的绝对值小于其中,λ为波长,d为相邻接收通道间距。
7.根据权利要求1所述的方法,其特征在于,所述根据所述频率值获得波达方向为获得波达方向的法向偏移角α(t0):α(t0)=arcsin(λξmax),其中,λ为波长,ξmax为插值后频域幅值最大点处的频率值。
8.根据权利要求7所述的方法,其特征在于,所述方法还包括:根据波达方向的法向偏移角α(t0)获得各通道当前距离门的权值
其中,dn为第n个接收天线相对于发射天线的位置值,dn的绝对值等于第n个接收天线与发射天线间的距离,ξmax为插值后频域幅值最大点处的频率值;
根据所述权值进行数字波束DBF合成。
9.根据权利要求8所述的方法,其特征在于,所述根据所述权值进行DBF合成为:
以各通道权值对所述脉冲压缩后的信号进行加权求和,获得DBF合成后的信号G(τ,t):
其中,τ是方位向慢时间,t是距离向快时间。
10.一种波达方向的检测装置,其特征在于,所述装置包括:脉冲压缩模块、傅里叶变换模块、插值模块、以及波达方向获取模块,其中,
所述脉冲压缩模块,用于对俯仰向各通道得到的回波数据做距离向脉冲压缩,得到脉冲压缩后的信号;
所述傅里叶变换模块,用于对脉冲压缩后的信号在每个距离门上沿通道做傅里叶变换;
所述插值模块,用于对傅立叶变换后的结果进行插值;其中,所述对傅立叶变换后的结果进行插值在所述傅里叶变换模块进行所述傅里叶变换之前,通过时域补充零值样点的方法实现;或者,在所述傅里叶变换模块进行所述傅里叶变换过程中,通过加密频域间隔实现;
所述波达方向获取模块,用于获得插值后频域幅值最大点处的频率值,并根据所述频率值获得波达方向。
11.根据权利要求10所述的装置,其特征在于,所述装置还包括:权值获取模块和波束合成模块,其中,
所述权值获取模块,用于根据波达方向获得各通道当前距离门的权值;
所述波束合成模块,用于以各通道权值对所述脉冲压缩后的信号进行加权求和,获得DBF合成后的信号。
CN201310475827.3A 2013-10-12 2013-10-12 一种波达方向的检测方法及装置 Active CN103558584B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201310475827.3A CN103558584B (zh) 2013-10-12 2013-10-12 一种波达方向的检测方法及装置
EP14169310.1A EP2860549B1 (en) 2013-10-12 2014-05-21 Method and device for direction of arrival (doa) detection of sar signals

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310475827.3A CN103558584B (zh) 2013-10-12 2013-10-12 一种波达方向的检测方法及装置

Publications (2)

Publication Number Publication Date
CN103558584A CN103558584A (zh) 2014-02-05
CN103558584B true CN103558584B (zh) 2016-08-17

Family

ID=50012886

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310475827.3A Active CN103558584B (zh) 2013-10-12 2013-10-12 一种波达方向的检测方法及装置

Country Status (2)

Country Link
EP (1) EP2860549B1 (zh)
CN (1) CN103558584B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018104468A1 (en) * 2016-12-08 2018-06-14 Iee International Electronics & Engineering S.A. Direction of arrival estimation for automotive spread radar systems
CN109521426B (zh) * 2017-09-18 2022-01-07 比亚迪股份有限公司 基于汽车雷达获取目标的角度的方法及其装置
CN107741586B (zh) * 2017-09-29 2021-04-23 王辉 基于DBF-TOPS加权的星载Ka InSAR信号处理方法
CN109597021B (zh) * 2018-11-05 2021-02-05 北京航天控制仪器研究所 一种波达方向估计方法及装置
CN111814688B (zh) * 2020-07-09 2023-10-13 成都傅立叶电子科技有限公司 基于FFTc的数字波束形成角度捕获方法及装置、存储介质
CN113534179B (zh) * 2021-06-18 2023-08-01 中电海康集团有限公司 一种激光雷达测距方法及装置
CN114428225B (zh) * 2022-04-06 2022-06-14 中国电子科技集团公司第二十九研究所 一种基于准匹配滤波的多辐射源到达角度测量方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1388384A (zh) * 2001-05-26 2003-01-01 华为技术有限公司 一种提高接收信号波达方向估计准确度的方法
CN102013911A (zh) * 2010-12-02 2011-04-13 哈尔滨工程大学 一种基于门限检测的宽带信号波达方向估计方法
CN102012505A (zh) * 2010-10-15 2011-04-13 西安电子科技大学 雷达低仰角目标的波达方向估计方法
CN102621539A (zh) * 2012-04-24 2012-08-01 重庆大学 单次搜索宽带信号波达方向估计方法
CN103236885A (zh) * 2013-03-19 2013-08-07 南京大学 一种超宽带无频率色散的doa估计方法及其装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050195103A1 (en) * 2004-01-13 2005-09-08 Davis Dennis W. Phased arrays exploiting geometry phase and methods of creating such arrays

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1388384A (zh) * 2001-05-26 2003-01-01 华为技术有限公司 一种提高接收信号波达方向估计准确度的方法
CN102012505A (zh) * 2010-10-15 2011-04-13 西安电子科技大学 雷达低仰角目标的波达方向估计方法
CN102013911A (zh) * 2010-12-02 2011-04-13 哈尔滨工程大学 一种基于门限检测的宽带信号波达方向估计方法
CN102621539A (zh) * 2012-04-24 2012-08-01 重庆大学 单次搜索宽带信号波达方向估计方法
CN103236885A (zh) * 2013-03-19 2013-08-07 南京大学 一种超宽带无频率色散的doa估计方法及其装置

Also Published As

Publication number Publication date
EP2860549A1 (en) 2015-04-15
CN103558584A (zh) 2014-02-05
EP2860549B1 (en) 2017-07-05

Similar Documents

Publication Publication Date Title
CN103558584B (zh) 一种波达方向的检测方法及装置
US9207313B2 (en) MIMO angle estimation with simultaneous mainlobe jammer cancellation
CN106021637B (zh) 互质阵列中基于迭代稀疏重构的doa估计方法
US8669901B2 (en) Method for determining azimuth and elevation angles of arrival of coherent sources
CN103176168B (zh) 一种机载非正侧视阵雷达近程杂波对消方法
CN108008386B (zh) 一种基于单快拍music算法的距离向处理方法
CN110488255A (zh) 一种相控阵雷达单脉冲高分辨测角系统及方法
CN104155648A (zh) 基于阵列数据重排的高频地波雷达单次快拍music测向方法
CN105137409B (zh) 基于幅相约束的目标信号稳健空时自适应处理方法
CN106855622A (zh) 一种子阵级相控阵雷达的测角方法
CN102087354A (zh) 无源雷达分组ls-clean微弱目标检测方法
CN102544755A (zh) 一种基于强散射点的均匀线阵校准方法
US20130176173A1 (en) Grating lobe mitigation in presence of simultaneous receive beams
CN108957406A (zh) 一种基于神经网络的雷达主瓣干扰抑制方法
CN110426670A (zh) 基于tls-cs的外辐射源雷达超分辨doa估计方法
CN107607915B (zh) 基于固定地物回波的有源相控阵雷达接收通道校正方法
CN113238211B (zh) 一种干扰条件下参数化自适应阵列信号检测方法与系统
CN107290732A (zh) 一种量子大爆炸的单基地mimo雷达测向方法
CN107576947A (zh) 基于时间平滑的l型阵对相干信源二维波达方向估计方法
CN110244273A (zh) 一种基于均匀分布式阵列的目标角度估计方法
CN111537997A (zh) 一种基于mimo和压缩感知技术的三维雷达成像方法
CN105334435A (zh) 一种基于任意阵形的自适应局部放电超声监测方法
CN109188373B (zh) 基于子阵阻塞矩阵预处理的抗主瓣干扰方法
CN100585429C (zh) 一种基于非直线天线阵列的无源通道校正方法
Yong et al. Echo DOA based high-resolution target location

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant