CN107340495A - 一种基于阵列雷达的目标波达方向快速估计方法 - Google Patents
一种基于阵列雷达的目标波达方向快速估计方法 Download PDFInfo
- Publication number
- CN107340495A CN107340495A CN201710510050.8A CN201710510050A CN107340495A CN 107340495 A CN107340495 A CN 107340495A CN 201710510050 A CN201710510050 A CN 201710510050A CN 107340495 A CN107340495 A CN 107340495A
- Authority
- CN
- China
- Prior art keywords
- mrow
- angle measurement
- msub
- signal data
- digital signal
- 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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于阵列雷达的目标波达方向快速估计方法,其主要思路为:确定阵列雷达,所述阵列雷达包含N个阵元;阵列雷达接收雷达模拟信号,并对该雷达模拟信号进行A/D采样和脉冲压缩处理,得到脉冲压缩回波数字信号数据矢量;设定脉冲压缩回波数字信号数据矢量的来波方向θ,并得到θ方向上的合成导向矢量;根据θ方向上的合成导向矢量,得到N个阵元的DBF权矢量;对脉冲压缩回波数字信号数据矢量进行DBF粗测,得到DBF粗测结果,所述DBF粗测结果为脉冲压缩回波数字信号数据矢量的来波方向θ的仰角位置;根据DBF粗测结果对脉冲压缩回波数字信号数据矢量进行SVML精测,得到基于阵列雷达的目标波达方向。
Description
技术领域
本发明属于米波雷达信号处理技术领域,特别涉及一种基于阵列雷达的目标波达方向快速估计方法,适用于实际工程应用。
背景技术
米波雷达测角测高是信号处理领域的重要分支,雷达发射电磁波并接收来自其威力覆盖范围内的目标回波信号,然后从接收的回波信号中提取目标的位置和其他信息,以用于探测、定位以及目标识别;米波雷达因为其工作频段和波长较长的特点,具有良好的反隐身和反辐射导弹能力,但同时因为米波波束较宽,存在波束“打地”以及多径效应等问题,严重影响了米波雷达测角精度。
早期目标的波达方向估计主要采用数字波束形成(DBF)算法,而DBF算法的角度分辨率和天线的孔径成反比,为了达到较好的角度分辨率需要增大雷达天线的阵列孔径,这一点是难以实现的;对于米波雷达而言,在低仰角区域的直达波和多径反射波位于同一波束宽度内,不可分辨,会出现角度模糊;也就是说,采用DBF测角算法简单,运算速度快,能够测出目标方向,但是测角精度不够。
为了达到较高的测角精度,出现了一类子空间拟合算法:最大似然(ML)算法,最大似然(ML)算法在目标信噪比较低、快拍数较少的条件下,相比于其他算法能够达到更好的测角性能;但是ML算法需要进行空域的多维搜索,运算量较大,为了减少运算量,学者们对其进行了改进,最终形成了合成导向矢量(SVML)算法,SVML算法对直达角和反射角都进行搜索,虽然较好的解决了多径反射波和直达波在一个波束宽度瓣内DBF算法无法进行目标波达方向估计的问题,但却存在两个不足:第一,需要对直达角和反射角分别进行二维估计,计算量偏大,无法确保信号处理的实时性;第二,大量角度的搜索,增大了运算量,很难满足测角的实时性。
发明内容
针对上述现有技术存在的不足,本发明的目的在于提出一种基于阵列雷达的目标波达方向快速估计方法,该种基于阵列雷达的目标波达方向快速估计方法是一种基于CUDA架构的DBF与SVML结合测向算法,能够根据目标回波信息和地形准确快速估计目标波达方向,进而能够提高米波雷达测角测高的速率和精度。
本发明的主要思路:本发明结合DBF算法与SVML算法进行测角,主要是因为米波雷达波束较宽,在低仰角区域的直达波和多径反射波位于同一波束宽度内,不可分辨,虽然采用DBF算法会出现角度模糊,但是DBF测角算法比较简单,运算速度快,而SVML算法的波瓣比DBF算法的波瓣更窄,角度分辨率更高;所以将DBF算法与SVML算法结合进行测角能够优势互补,DBF算法主要负责粗测,SVML算法负责精测,这样就缩小了搜索范围,加快了搜索进程,可以同时计算多个角度的似然值,一方面可以提高测角速率,另一方面,可以进一步减小搜索步进间隔以提高测角精度。
为达到上述技术目的,本发明采用如下技术方案予以实现。
一种基于阵列雷达的目标波达方向快速估计方法,包括以下步骤:
步骤1,确定阵列雷达,所述阵列雷达包含N个阵元;阵列雷达接收雷达模拟信号,并对该雷达模拟信号进行A/D采样和脉冲压缩处理,得到脉冲压缩回波数字信号数据矢量;
设定脉冲压缩回波数字信号数据矢量的来波方向θ,并得到θ方向上的合成导向矢量
步骤2,根据θ方向上的合成导向矢量,得到N个阵元的DBF权矢量;
步骤3,对脉冲压缩回波数字信号数据矢量进行DBF粗测,得到DBF粗测结果,所述DBF粗测结果为脉冲压缩回波数字信号数据矢量的来波方向θ的仰角位置;
步骤4,根据DBF粗测结果对脉冲压缩回波数字信号数据矢量进行SVML精测,得到基于阵列雷达的目标波达方向。
本发明的有益效果:
第一,本发明方法实现了将DBF算法和SVML算法相结合,使得SVML方法搜索范围变小,减小了运算量,加快了运算速率。
第二,本发明在提高目标角度测量精度方面更有优势:一方面,增加DBF算法粗测,即首先采用DBF算法测出目标波达方向的大致范围,提高了测角速率;然后采用SVML算法进行精测,从而减少角度搜索范围,提高测角速率,解决了低仰角测角模糊问题;另一方面,采用并行处理的实现方法进行角度测量,不仅加快了运算速率,提高信号处理速率,同时能够利用并行这一特点在同样时间内以更小步进对粗测范围进行搜索,提高测角测高的准确度。
附图说明
下面结合附图说明和具体实施方式对本发明作进一步详细说明。
图1为本发明的一种基于阵列雷达的目标波达方向快速估计方法流程图;
图2(a)为使用本发明方法得到的DBF粗测结果图;
图2(b)为使用本发明方法得到的SVML精测结果图;
图2(c)为图2(b)的似然值转换成为单位dB的SVML谱图。
具体实施方式
参照图1,为本发明的一种基于阵列雷达的目标波达方向快速估计方法流程图;其中所述基于阵列雷达的目标波达方向快速估计方法,包括以下步骤:
步骤1,确定阵列雷达,所述阵列雷达包含N个阵元;阵列雷达接收雷达模拟信号,并对该雷达模拟信号进行A/D采样和脉冲压缩处理,得到脉冲压缩回波数字信号数据矢量。
步骤1的子步骤为:
1a)确定阵列雷达,阵列雷达为米博雷达;所述阵列雷达包含N个阵元;阵列雷达的检测范围内存在目标,且阵列雷达发射电磁波并接收雷达模拟信号,然后对接收到的雷达模拟信号进行A/D采样和脉冲压缩处理,得到脉冲压缩回波数字信号数据矢量;将所述脉冲压缩回波数字信号数据矢量记为X(t),X(t)=[x1(t),x2(t),...,xn(t),...,xN(t)],xn(t)表示第n个阵元的脉冲压缩回波数字信号数据,t表示时间变量;X(t)为N×P维,N表示阵列雷达包含的阵元总个数,P表示阵列雷达发射的脉冲总个数,且N、P分别为大于0的正整数。
1b)假设脉冲压缩回波数字信号数据矢量X(t)的来波方向为θ,则θ方向上的合成导向矢量为其表达式为:
其中,κ表示阵列雷达发射的电磁波波数,κ=2π/λ,λ表示阵列雷达发射的电磁波波长,d表示阵列雷达中相邻阵元的间距,N为阵列雷达包含的阵元总个数,上标T表示转置操作,e表示指数函数,j表示虚数单位,sin表示正弦函数。
步骤2,为了降低阵列的副瓣电平,需要对θ方向上的合成导向矢量进行加窗处理,得到加窗处理后的合成导向矢量,记为N个阵元的DBF权矢量W,其表达式为:
其中,W=[W1,W2,...,Wn,...,WN]T,n∈{1,2,…,N},Wn表示第n个阵元的DBF权矢量,上标T表示转置操作,Wwin表示长度为N的窗函数,上标T表示转置操作。
步骤3,对脉冲压缩回波数字信号数据矢量X(t)进行DBF粗测,得到DBF粗测结果,所述DBF粗测结果为脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的仰角位置。
对脉冲压缩回波数字信号数据矢量X(t)进行DBF粗测,其粗测过程为:
3.1初始化:假设脉冲压缩回波数字信号数据矢量X(t)的来波方向θ估计需要搜索的空间仰角范围为令为了使脉冲压缩回波数字信号数据矢量X(t)的来波方向θ每隔设定的步进△θ遍历搜索且对脉冲压缩回波数字信号数据矢量X(t)的来波方向θ估计需要搜索的空间仰角范围进行划分,得到ang_num_cu个测角,分别为
其中,表示取的最小整数,ang_num_cu为内包含的测角总个数;表示第个测角,的初始值为1,第1个测角为第ang_num_cu个测角为第个测角为
3.2计算第个测角的输出模值对脉冲压缩回波数字信号数据矢量X(t)进行加权求和,得到第个测角的输出模值其表达式为:
其中,表示第个测角的DBF权矢量,表示方向上的合成导向矢量,上标H表示共轭转置操作,Wwin表示长度为N的窗函数。
3.3令分别取1至ang_num_cu,并行执行3.2,进而同时得到第1个测角的输出模值y1(t)至第ang_num_cu个测角的输出模值yang_num_cu(t),记为ang_num_cu个测角的输出模值 其中,第1个测角的输出模值y1(t)至第ang_num_cu个测角的输出模值yang_num_cu(t)的计算过程为并行计算。
3.4得到DBF粗测结果:在ang_num_cu个测角的输出模值中搜索最大输出模值对应的测角,记为 为脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的仰角位置;也就是说,当N个阵元的DBF权矢量对脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的仰角位置方向上的信号实现同相叠加,即时,输出模值最大;表示方向上的合成导向矢量,W'表示方向上的DBF权矢量。
所述脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的仰角位置为DBF粗测结果。
步骤4,对根据DBF粗测结果对脉冲压缩回波数字信号数据矢量进行SVML精测,得到基于阵列雷达的目标波达方向。
(4a)确定脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的精测范围,其过程为:如果则脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的精测范围为如果则脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的精测范围为σ表示设定的角度,一般取σ≤2°。
为了便于说明,令脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的精测范围为[α1,α2],如果则如果α1=0,在脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的精测范围内令△η为步长,△η<△θ;令ang_num为[α1,α2]内包含的测角总个数, 表示取的最小整数。
初始化:令θt表示第t个测角,t∈{1,2,…,ang_num},t的初始值为1,第1个测角θ1为α1,第ang_num个测角θang_num为α2,第t个测角θt为α1+(t-1)△η,α2=α1+(ang_num-1)△η,ang_num表示脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的精测范围[α1,α2]内包含的测角总个数,ang_num为大于0的正整数。
(4b)计算第t个测角θt的直达波和第t个测角θt反射波的波程差△rt,以及第t个测角θt的反射波的角度θtr,其表达式分别为:
△rt=rt1+rt2-rd
其中,第t个测角θt的反射路径包括两段,分别为第t个测角θt经反射点反射前的入射段和第t个测角θt经反射点反射后的反射段,并分别记为第t个测角θt反射路径的入射段和第t个测角θt反射路径的反射段,rt1表示第t个测角θt反射路径的入射段,rt2表示第t个测角θt的反射路径的反射段,rd表示目标和阵列雷达的直线距离,上标-1表示求逆操作,cos表示求余弦操作。
然后计算第t个测角θt的合成导向矢量a(θt),其表达式为:
其中,ad(θt)表示第t个测角θt的直达波导向矢量,ar(θtr)表示第t个测角θt的反射波导向矢量,ρ表示设定的反射系数,常取-0.95;λ表示阵列雷达发射的电磁波波长,上标j表示虚数单位,e表示指数函数,下标r表示反射。
(4c)利用第t个测角θt的合成导向矢量a(θt),构造第t个测角θt的正交投影矩阵P(θt):
P(θt)=a(θt)[aH(θt)a(θt)](-1)aH(θt)
其中,a(θt)表示第t个测角θt的合成导向矢量,上标H表示共轭转置操作。
(4d)令t分别取1至ang_num,并行执行(4b)和(4c),进而同时第1个测角θ1的正交投影矩阵P(θ1)至第ang_num个测角θang_num的正交投影矩阵P(θang_num),记为ang_num个测角的正交投影矩阵P(θ),P(θ)=[P(θ1),P(θ2),...,P(θt),...,P(θang_num)];其中,第1个测角θ1的合成导向矢量a(θ1)至第ang_num个测角θang_num的合成导向矢量a(θang_num)的计算过程为并行计算。
(4e)根据ang_num个测角的正交投影矩阵P(θ),构造关于脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的最大似然函数Q(θ),其表达式为:
Q(θ)=tr{P(θ)R}
其中,tr{-}表示求迹操作,R表示脉冲压缩回波数字信号数据矢量X(t)的自相关矩阵,θ表示脉冲压缩回波数字信号数据矢量X(t)的来波方向。
(4f)对关于脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的最大似然函数Q(θ)进行谱峰搜索,计算得到脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的估计值
其中,表示求解关于脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的最大似然函数Q(θ)中最大谱峰所对应的来波方向操作;所述脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的估计值为基于阵列雷达的目标波达方向。
下面结合仿真实验对本发明效果做进一步的验证说明。
1.仿真条件
本实验的操作***为Windows 7 SP1,显卡是计算能力为2.0的NVIDIA TeslaC2050下进行的,该显卡支持双精度浮点运算。通过利用VS2013+CUDA 7.5作为仿真平台,实现基于阵列雷达的目标波达方向估计;本实验中目标真实角度为1.4°。
2.仿真内容
首先,进行DBF粗测,得到图2(a)的DBF粗测结果图;然后,利用DBF粗测结果,进行SVML精测,得到图2(b)的SVML谱图;将图2(b)的似然值转换为dB为单位,得到图2(c)的结果图。
3.附图说明
图2(a)为使用本发明方法得到的DBF粗测结果图;其中,横坐标表示直达波入射角度,纵坐标表示归一化幅度值,从图2(a)中可看出,当目标处在低仰角区时,DBF波瓣较宽,无法准确测出目标仰角。
图2(b)是在DBF粗测结果上进行SVML精测,得到的SVML谱图,图2(b)为使用本发明方法得到的SVML精测结果图;其中,横坐标表示直达波入射角度,纵坐标表示SVML似然函数值,由图2(b)中可看出,在低仰角区时,SVML波瓣较窄,依然保持很好的测角精度,可以较为准确的得到目标仰角。
图2(c)为图2(b)的似然值转换成为单位dB的SVML谱图,其中,横坐标表示直达波入射角度,纵坐标为以dB为单位的SVML似然函数值。
通过图2(c)与图2(a)进行对比,可以直观的看出低仰角区,SVML算法相对DBF算法而言测角精度更高。
4.表格说明
表1为本发明实例中不同测高算法在不同平台上测角时间统计,在本发明实例中,设定低仰角区目标真实仰角为1.4°。
表1为本发明实例中不同测高算法在CPU与GPU所运行时间对比(单位:ms)
从表格数据可以看出,GPU比CPU的数据运算效率更高。要说明的是,表格中DBF粗测一栏括号内数字表示总时间0.242ms内有0.096ms用于GPU和CPU间数据传输,也就是说,GPU有很大一部分时间是用在数据传输上。所以,在GPU运算过程中,应当尽量减数据传输时间。
表1为本发明实例中不同测高算法在CPU与GPU所运行时间对比(单位:ms),其中,最后一列为GPU相对于CPU的加速比,从表格中可以看出,SVML算法加速比达到292.444倍,总的加速比达到39.748倍,证明了本方法的有效性。
Claims (5)
1.一种基于阵列雷达的目标波达方向快速估计方法,其特征在于,包括以下步骤:
步骤1,确定阵列雷达,所述阵列雷达包含N个阵元;阵列雷达接收雷达模拟信号,并对该雷达模拟信号进行A/D采样和脉冲压缩处理,得到脉冲压缩回波数字信号数据矢量;
设定脉冲压缩回波数字信号数据矢量的来波方向θ,并得到θ方向上的合成导向矢量;
步骤2,根据θ方向上的合成导向矢量,得到N个阵元的DBF权矢量;
步骤3,对脉冲压缩回波数字信号数据矢量进行DBF粗测,得到DBF粗测结果,所述DBF粗测结果为脉冲压缩回波数字信号数据矢量的来波方向θ的仰角位置;
步骤4,根据DBF粗测结果对脉冲压缩回波数字信号数据矢量进行SVML精测,得到基于阵列雷达的目标波达方向。
2.如权利要求1所述的一种基于阵列雷达的目标波达方向快速估计方法,其特征在于,在步骤1中,将所述脉冲压缩回波数字信号数据矢量记为X(t),X(t)=[x1(t),x2(t),...,xn(t),...,xN(t)],xn(t)表示第n个阵元的脉冲压缩回波数字信号数据,t表示时间变量;X(t)为N×P维,N表示阵列雷达包含的阵元总个数,P表示阵列雷达发射的脉冲总个数,且N、P分别为大于0的正整数;
所述θ方向上的合成导向矢量为其表达式为:
<mrow>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>,</mo>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>&kappa;</mi>
<mi>d</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>&kappa;</mi>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
其中,κ表示阵列雷达发射的电磁波波数,κ=2π/λ,λ表示阵列雷达发射的电磁波波长,d表示阵列雷达中相邻阵元的间距,N为阵列雷达包含的阵元总个数,上标T表示转置操作,e表示指数函数,j表示虚数单位,sin表示正弦函数。
3.如权利要求2所述的一种基于阵列雷达的目标波达方向快速估计方法,其特征在于,在步骤2中,将所述得到N个阵元的DBF权矢量记为W,其表达式为:
<mrow>
<mi>W</mi>
<mo>=</mo>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<msub>
<mi>W</mi>
<mrow>
<mi>w</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
其中,W=[W1,W2,...,Wn,...,WN]T,n∈{1,2,…,N},Wn表示第n个阵元的DBF权矢量,上标T表示转置操作,Wwin表示长度为N的窗函数,上标T表示转置操作。
4.如权利要求3所述的一种基于阵列雷达的目标波达方向快速估计方法,其特征在于,步骤3的子步骤为:
3.1初始化:假设脉冲压缩回波数字信号数据矢量X(t)的来波方向θ估计需要搜索的空间仰角范围为 令对脉冲压缩回波数字信号数据矢量X(t)的来波方向θ估计需要搜索的空间仰角范围进行划分,得到ang_num_cu个测角,分别为
其中, 表示取大于或者等于的最小整数,ang_num_cu为内包含的测角总个数;表示第个测角,的初始值为1,第1个测角为第ang_num_cu个测角为第个测角为
3.2计算第个测角的输出模值对脉冲压缩回波数字信号数据矢量X(t)进行加权求和,得到第个测角的输出模值其表达式为:
<mrow>
<msub>
<mi>y</mi>
<mo>&part;</mo>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<msub>
<mi>W</mi>
<mo>&part;</mo>
</msub>
<mi>H</mi>
</msup>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,表示第个测角的DBF权矢量,表示方向上的合成导向矢量,上标H表示共轭转置操作,Wwin表示长度为N的窗函数;
3.3令分别取1至ang_num_cu,并行执行3.2,进而同时得到第1个测角的输出模值y1(t)至第ang_num_cu个测角的输出模值yang_num_cu(t),记为ang_num_cu个测角的输出模值
3.4在ang_num_cu个测角的输出模值中搜索最大输出模值对应的测角,记为 为脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的仰角位置;所述脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的仰角位置为DBF粗测结果。
5.如权利要求4所述的一种基于阵列雷达的目标波达方向快速估计方法,其特征在于,步骤4的子步骤为:
(4a)确定脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的精测范围,其过程为:如果则脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的精测范围为如果则脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的精测范围为σ表示设定的角度;
令脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的精测范围为[α1,α2],如果则如果α1=0,在脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的精测范围内令△η为步长,△η<△θ;令ang_num为[α1,α2]内包含的测角总个数, 表示取大于或者等于的最小整数;
初始化:令θt表示第t个测角,t∈{1,2,…,ang_num},t的初始值为1,第1个测角θ1为α1,第ang_num个测角θang_num为α2,第t个测角θt为α1+(t-1)△η,α2=α1+(ang_num-1)△η,ang_num表示脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的精测范围[α1,α2]内包含的测角总个数,ang_num为大于0的正整数;
(4b)计算第t个测角θt的直达波和第t个测角θt反射波的波程差△rt,以及第t个测角θt的反射波的角度θtr,其表达式分别为:
△rt=rt1+rt2-rd
<mrow>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>t</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msup>
<mi>cos</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>(</mo>
<mfrac>
<mrow>
<msubsup>
<mi>r</mi>
<mi>d</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>t</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>r</mi>
<mrow>
<mi>t</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>r</mi>
<mi>d</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>t</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中,第t个测角θt的反射路径包括两段,分别为第t个测角θt经反射点反射前的入射段和第t个测角θt经反射点反射后的反射段,并分别记为第t个测角θt反射路径的入射段和第t个测角θt反射路径的反射段,rt1表示第t个测角θt反射路径的入射段,rt2表示第t个测角θt的反射路径的反射段,rd表示目标和阵列雷达的直线距离,上标-1表示求逆操作,cos表示求余弦操作;然后计算第t个测角θt的合成导向矢量a(θt),其表达式为:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>a</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>a</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&rho;a</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>t</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mi>&rho;e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;&Delta;r</mi>
<mi>t</mi>
</msub>
</mrow>
<mi>&lambda;</mi>
</mfrac>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,ad(θt)表示第t个测角θt的直达波导向矢量,ar(θtr)表示第t个测角θt的反射波导向矢量,ρ表示设定的反射系数,常取-0.95;λ表示阵列雷达发射的电磁波波长,上标j表示虚数单位,e表示指数函数,下标r表示反射;
(4c)利用第t个测角θt的合成导向矢量a(θt),构造第t个测角θt的正交投影矩阵P(θt):
P(θt)=a(θt)[aH(θt)a(θt)](-1)aH(θt)
其中,a(θt)表示第t个测角θt的合成导向矢量,上标H表示共轭转置操作;
(4d)令t分别取1至ang_num,并行执行(4b)和(4c),进而同时得到第1个测角θ1的正交投影矩阵P(θ1)至第ang_num个测角θang_num的正交投影矩阵P(θang_num),记为ang_num个测角的正交投影矩阵P(θ),P(θ)=[P(θ1),P(θ2),...,P(θt),...,P(θang_num)];
(4e)根据ang_num个测角的正交投影矩阵P(θ),构造关于脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的最大似然函数Q(θ),其表达式为:
Q(θ)=tr{P(θ)R}
其中,tr{-}表示求迹操作,R表示脉冲压缩回波数字信号数据矢量X(t)的自相关矩阵,θ表示脉冲压缩回波数字信号数据矢量X(t)的来波方向;
(4f)对关于脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的最大似然函数Q(θ)进行谱峰搜索,计算得到脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的估计值
其中,表示求解关于脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的最大似然函数Q(θ)中最大谱峰所对应的来波方向操作;所述脉冲压缩回波数字信号数据矢量X(t)的来波方向θ的估计值为基于阵列雷达的目标波达方向。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710510050.8A CN107340495A (zh) | 2017-06-28 | 2017-06-28 | 一种基于阵列雷达的目标波达方向快速估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710510050.8A CN107340495A (zh) | 2017-06-28 | 2017-06-28 | 一种基于阵列雷达的目标波达方向快速估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107340495A true CN107340495A (zh) | 2017-11-10 |
Family
ID=60217963
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710510050.8A Pending CN107340495A (zh) | 2017-06-28 | 2017-06-28 | 一种基于阵列雷达的目标波达方向快速估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107340495A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108710103A (zh) * | 2018-06-12 | 2018-10-26 | 西安电子科技大学 | 基于稀疏阵列的强弱多目标超分辨测向与信源数估计方法 |
WO2019095912A1 (zh) * | 2017-11-16 | 2019-05-23 | 华南理工大学 | 基于可调夹角均匀线阵的水下波达方向估计方法及装置 |
CN109901131A (zh) * | 2019-03-24 | 2019-06-18 | 西安电子科技大学 | 基于斜投影的多径利用相干波束形成方法 |
CN110058193A (zh) * | 2019-04-22 | 2019-07-26 | 西安空间无线电技术研究所 | 一种基于单接收通道的数字多波束测角方法及*** |
CN110940973A (zh) * | 2019-12-10 | 2020-03-31 | 成都纳雷科技有限公司 | 一种用于雷达目标检测的角度测量方法及装置 |
CN113325363A (zh) * | 2020-02-28 | 2021-08-31 | 加特兰微电子科技(上海)有限公司 | 确定波达方向的方法、装置及相关设备 |
CN113359196A (zh) * | 2021-05-26 | 2021-09-07 | 上海交通大学 | 基于子空间法和dbf的多目标生命体征探测方法 |
CN114428236A (zh) * | 2022-04-06 | 2022-05-03 | 南京楚航科技有限公司 | 一种车载毫米波雷达角度置信度估计方法 |
CN114817675A (zh) * | 2022-06-20 | 2022-07-29 | 中国电子科技集团公司第五十四研究所 | 利用广域随机稀布阵波束旁瓣特征的目标快速搜索方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103018730A (zh) * | 2012-11-27 | 2013-04-03 | 西安电子科技大学 | 分布式子阵波达方向估计方法 |
CN103353587A (zh) * | 2013-06-13 | 2013-10-16 | 西安电子科技大学 | 基于改进的自适应数字波束形成的干涉测角方法 |
US20140159945A1 (en) * | 2012-12-11 | 2014-06-12 | National Chiao Tung University | Method and Device for Estimating Direction of Arrival |
CN106093869A (zh) * | 2016-05-27 | 2016-11-09 | 上海无线电设备研究所 | 基于地形参数表的周期性合成导向矢量最大似然算法 |
-
2017
- 2017-06-28 CN CN201710510050.8A patent/CN107340495A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103018730A (zh) * | 2012-11-27 | 2013-04-03 | 西安电子科技大学 | 分布式子阵波达方向估计方法 |
US20140159945A1 (en) * | 2012-12-11 | 2014-06-12 | National Chiao Tung University | Method and Device for Estimating Direction of Arrival |
CN103353587A (zh) * | 2013-06-13 | 2013-10-16 | 西安电子科技大学 | 基于改进的自适应数字波束形成的干涉测角方法 |
CN106093869A (zh) * | 2016-05-27 | 2016-11-09 | 上海无线电设备研究所 | 基于地形参数表的周期性合成导向矢量最大似然算法 |
Non-Patent Citations (2)
Title |
---|
SHUAI LIU等: "An Improved Polarization and DOA Estimation Algorithm", 《2016 IEEE》 * |
汪坤: "某米波雷达测高方法的改进及工程实现", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019095912A1 (zh) * | 2017-11-16 | 2019-05-23 | 华南理工大学 | 基于可调夹角均匀线阵的水下波达方向估计方法及装置 |
CN108710103A (zh) * | 2018-06-12 | 2018-10-26 | 西安电子科技大学 | 基于稀疏阵列的强弱多目标超分辨测向与信源数估计方法 |
CN109901131A (zh) * | 2019-03-24 | 2019-06-18 | 西安电子科技大学 | 基于斜投影的多径利用相干波束形成方法 |
CN109901131B (zh) * | 2019-03-24 | 2022-12-02 | 西安电子科技大学 | 基于斜投影的多径利用相干波束形成方法 |
CN110058193A (zh) * | 2019-04-22 | 2019-07-26 | 西安空间无线电技术研究所 | 一种基于单接收通道的数字多波束测角方法及*** |
CN110940973B (zh) * | 2019-12-10 | 2022-11-04 | 成都纳雷科技有限公司 | 一种用于雷达目标检测的角度测量方法及装置 |
CN110940973A (zh) * | 2019-12-10 | 2020-03-31 | 成都纳雷科技有限公司 | 一种用于雷达目标检测的角度测量方法及装置 |
CN113325363A (zh) * | 2020-02-28 | 2021-08-31 | 加特兰微电子科技(上海)有限公司 | 确定波达方向的方法、装置及相关设备 |
CN113359196A (zh) * | 2021-05-26 | 2021-09-07 | 上海交通大学 | 基于子空间法和dbf的多目标生命体征探测方法 |
CN113359196B (zh) * | 2021-05-26 | 2023-01-20 | 上海交通大学 | 基于子空间法和dbf的多目标生命体征探测方法 |
CN114428236B (zh) * | 2022-04-06 | 2022-07-01 | 南京楚航科技有限公司 | 一种车载毫米波雷达角度置信度估计方法 |
CN114428236A (zh) * | 2022-04-06 | 2022-05-03 | 南京楚航科技有限公司 | 一种车载毫米波雷达角度置信度估计方法 |
CN114817675A (zh) * | 2022-06-20 | 2022-07-29 | 中国电子科技集团公司第五十四研究所 | 利用广域随机稀布阵波束旁瓣特征的目标快速搜索方法 |
CN114817675B (zh) * | 2022-06-20 | 2022-09-09 | 中国电子科技集团公司第五十四研究所 | 利用广域随机稀布阵波束旁瓣特征的目标快速搜索方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107340495A (zh) | 一种基于阵列雷达的目标波达方向快速估计方法 | |
US7315488B2 (en) | Methods and systems for passive range and depth localization | |
CN103308909B (zh) | 基于多参数交替搜索的米波雷达低仰角测高方法 | |
CN104007440B (zh) | 一种加速分解后向投影聚束合成孔径雷达成像方法 | |
CN104898119B (zh) | 一种基于相关函数的动目标参数估计方法 | |
US7961147B1 (en) | Long baseline phase interferometer ambiguity resolution using frequency differences | |
US5614912A (en) | Radar processing method and apparatus | |
CN103364772B (zh) | 基于实数域广义多重信号分类算法的目标低仰角估计方法 | |
US6594201B2 (en) | System and method for localizing targets using multiple arrays | |
CN106646373B (zh) | 基于浅海多途到达角和到达时延的近距离声源联合定位方法 | |
CN103760546B (zh) | 一种雷达用低空目标波达方向估计方法 | |
CN103885054B (zh) | 一种基于分布源反射模型的米波雷达低仰角测高方法 | |
CN107607943A (zh) | 基于干涉相位辅助的延迟多普勒雷达高度表的测高方法 | |
CN102393525A (zh) | 子空间投影的导航干扰抑制与信号增强方法 | |
CN103293526A (zh) | 基于合成导向矢量的波束形成米波雷达测高方法 | |
CN109188019A (zh) | 基于多重信号分类算法的三维风速风向测量方法 | |
CN103616685B (zh) | 基于图像特征的isar图像几何定标方法 | |
CN106199600A (zh) | 基于多普勒估计的方位多通道合成孔径雷达成像方法 | |
CN104793177B (zh) | 基于最小二乘法的麦克风阵列测向方法 | |
CN108919220A (zh) | 基于嵌入式gpu的弹载sar前侧视成像方法 | |
CN108089147A (zh) | 改进的短波单位定位方法 | |
CN105738887A (zh) | 基于多普勒通道划分的机载雷达杂波功率谱的优化方法 | |
CN109581320A (zh) | 一种毫米波雷达后端中频信号的谱峰快速搜索方法 | |
CN113064147A (zh) | 一种低信噪比下新型匹配场被动定位方法 | |
CN106872936B (zh) | 近场源l型声矢量传感器阵列解模糊多参数估计方法 |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20171110 |
|
WD01 | Invention patent application deemed withdrawn after publication |