CN103293526A - 基于合成导向矢量的波束形成米波雷达测高方法 - Google Patents
基于合成导向矢量的波束形成米波雷达测高方法 Download PDFInfo
- Publication number
- CN103293526A CN103293526A CN2013101920660A CN201310192066A CN103293526A CN 103293526 A CN103293526 A CN 103293526A CN 2013101920660 A CN2013101920660 A CN 2013101920660A CN 201310192066 A CN201310192066 A CN 201310192066A CN 103293526 A CN103293526 A CN 103293526A
- Authority
- CN
- China
- Prior art keywords
- radar
- wave
- target
- reflection
- theta
- 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
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于合成导向矢量的波束形成米波雷达测高方法,主要解决现有技术测高时信号处理运算量大、计算复杂的问题,本发明从雷达回波中提取目标信号,确定目标仰角搜索范围;建立雷达波反射模型,用搜索仰角计算反射点与雷达中心的地面弧长和反射点与目标的地面弧长;计算目标和雷达中心的等效高度后计算反射角和反射波与直射波波程差;计算合成导向矢量;对回波数据波束形成后选取最大归一化功率值对应的角度作为目标仰角,并计算目标高度。本发明将波束形成和合成导向矢量相结合,并将目标反射模型进行简化处理,保证米波雷达的仰角估计精度,提高米波雷达测高的计算速度,简化测高计算过程,使米波雷达的实时信号处理更加容易。
Description
技术领域
本发明属于雷达信号处理技术领域,主要涉及米波雷达测高,具体地说是针对数字阵列米波雷达,提出一种基于合成导向矢量的波束形成米波雷达测高方法,在保证运算精度的同时降低运算量,用于雷达目标跟踪。
背景技术
米波雷达具有反隐身和抗反辐射导弹的能力,但因其波长较长,使其天线波束宽度较宽,角分辨率较低,在低仰角目标探测时地(海)面的多径反射现象更加严重。由于多径信号与直射波信号的相关性较强,导致米波雷达的测高问题一直是雷达界尚未很好解决的难题。
目前,国内外学者对于米波测高问题进行了大量的研究,取得了一些突破性的进展。陈伯孝等于2007年在《电子学报》和雷达年会上介绍了“基于波瓣***的米波雷达测高方法”,利用不同天线***波瓣的相位关系,确定目标所在仰角区间,对接收信号进行比幅处理提取归一化误差信号,最后根据归一化误差信号和仰角区间查表得到目标的高度。但该方法对阵地的平坦性要求较高,且测高精度只能达到距离的1%,难以满足一些精度较高的实际使用要求。杨雪亚于2011年在其博士论文中提出了基于多径信号模型的合成导向矢量最大似然测高算法,将精确多径反射模型引入到米波雷达测高中,提高了测高精度,在工程应用中效果较好。但此方法主要存在两个问题,一是在计算目标的反射角和波程差时应用了精确多径信号反射模型,导致计算时数值动态范围较大,单精度浮点数不能满足精度要求,而现有的DSP处理器(如ADSP-TS201)不能很好的支持双精度计算,计算所需时间大大加剧;二是最大似然方法需计算相关矩阵和投影矩阵,算法中需要矩阵相乘,计算量巨大,不利于提高测高计算速度,降低了雷达处理目标的能力。
波束形成算法是时域傅里叶谱估计的一种空域拓展形式,通过将各阵元输出进行加权求和,在一定时间内将天线阵列波束“导向”到一个方向上,对期望信号得到最大输出功率的导向位置给出了波达方向估计,是一种快速有效的算法。但这种方法会受到“瑞利限”的限制,即对位于一个波束宽度内的空间目标不可分辨。在低仰角米波雷达测高时,反射波和直射波夹角小于一个波束宽度,传统的波束形成算法将不可分辨。
米波雷达一般用作预警雷达,因此需要对目标高度进行快速精确的估计,上述几种方法均不能很好的解决米波雷达测高问题,制约着雷达目标测高的高效工作和工程化运行。
发明内容
本发明的目的在于克服上述已有技术的计算量大、测高速度低和测高精度低等不足,将合成导向矢量与波束形成相结合,解决了雷达测高时运算量大的问题,并且避免了普通波束形成在目标低仰角情况下不可分辨的问题,提供一种测高精度较高且运算量较小的基于合成导向矢量的波束形成米波雷达测高方法。
本发明是一种基于合成导向矢量的波束形成米波雷达测高方法,实现步骤包括:
步骤1:从雷达接收的回波中提取目标信号x(tk),并根据目标距离信息确定目标仰角搜索范围。
步骤2:建立雷达波反射模型,根据目标搜索仰角θd,计算反射点与雷达中心对应的地面弧长和反射点与目标对应的地面弧长,分别记为G1和G2,或者说是反射点分别与雷达中心和目标在地面投影点的弧长。
步骤3:分别计算目标和雷达中心在反射点切线上的等效高度,记为hr′和ht′,然后计算雷达波的反射角θi及反射波与直射波的波程差ΔR,目标等效高度为hr′,雷达中心等效高度为ht′。本发明通过模型的几何关系计算减小了数值范围,使单精度浮点数满足精度要求。
步骤4:利用直射波导向矢量a(θd)和反射波导向矢量a(θi)计算搜索角度的合成导向矢量as(θd)。
as(θd)=a(θd)+ρexp(-j2πΔR/λ)a(θi)
其中ρ为地面反射系数,λ为雷达波波长。
步骤5:对雷达回波进行波束形成处理,计算归一化阵列平均输出功率S(θd)。
步骤6:每个搜索仰角对应的归一化阵列平均输出功率计算完成后,选取最大归一化阵列平均输出功率对应的角度为目标仰角θ,并利用余弦定理计算目标的高度ht:
其中Rd为目标与雷达距离,hr为雷达中心高度,re为等效地球半径。ht即为所求的目标高度。
本发明利用合成导向矢量作为波束形成的权值,将直射波和反射波通过加权求和作为一个波束处理,因此避免了多波束夹角问题,解决了普通波束形成“瑞利限”对雷达低仰角测高时的限制。在计算反射角和波程差时采用近似多径反射模型,其满足计算精度同时降低数值范围,使此步骤在DSP中能够通过单精度计算实现,大大减少了运算时间,提高了雷达处理目标的能力。
本发明的实现还在于:其中步骤2中根据搜索仰角计算反射点与雷达中心对应的地面弧长G1和反射点与目标对应的地面弧长G2,按如下步骤进行:
(2a)建立雷达波反射模型,首先计算目标搜索角度θd对应的目标高度ht:
(2b)计算目标与雷达中心的地心夹角φ和地面投影点的弧长G:
G=φ·re
(2c)计算反射点与雷达中心的地面弧长G1和反射点与目标的地面弧长G2:
G2=G-G1
其中p和ξ均为中间变量,无具体物理含义。
本发明的实现还在于:其中步骤3中分别计算目标和雷达中心在反射点切线上的等效高度,记为hr′和ht′,然后计算雷达波的反射角θi及反射波与直射波的波程差ΔR,按如下步骤进行:
(3a)计算雷达中心等效高度hr′和目标等效高度ht′。
(3b)计算雷达波的反射角θi和波程差ΔR。
θi≈-asin(sin(θd)+2·hr′/Rd)
此步骤是本发明的关键性步骤之一,通过此步骤,利用近似模型减少数值范围,并得到构造合成导向矢量时所需的反射角和波程差,减少了运算量,提高了DSP运算速度。
本发明的实现还在于:其中步骤5中通过对雷达回波进行波束形成处理,计算归一化的阵列平均输出功率S(θd),按如下步骤进行:
(5a)计算第k个脉冲目标信号x(tk)对应搜索角度θd的波束形成值y(θd,k);
(5b)计算第k个脉冲目标信号x(tk)对应搜索角度θd的阵列输出功率p(θd,k);
p(θd,k)=y(θd,k)y(θd,k)H
(5c)每个脉冲目标信号对应的阵列输出功率计算完成后,计算所有K个脉冲的阵列平均输出功率P(θ);
(5d)对阵列输出功率消除导向矢量模值的影响,计算归一化阵列输出平均功率S(θd)。
此步骤是本发明的关键性步骤之一,通过此步骤,避免了最大似然算法中相关矩阵和投影矩阵的计算,减少了运算量,提高了目标角度估计的速度。
本发明结合了波束形成的运算简单和合成导向矢量对测高模型的良好匹配,能够在实现精确测高的同时大大减少运算量,提高雷达信号处理能力。
本发明与现有技术相比具有如下优点:
基于波瓣***的米波雷达测高方法计算简单,但测高精度受阵地影响较大,且测高精度较低;合成导向矢量最大似然算法测高精度较高,但运算量巨大,不利于雷达实时处理;普通波束形成算法测高精度较低,且当目标在低仰角时无法正常工作。
(1)本发明在目标角度估计时应用了波束形成算法,避免了最大似然算法中相关矩阵和投影矩阵计算时的巨大运算量,极大的减少了运算量,提高了测高速度。
(2)本发明使用直射波和反射波及波程差构造合成导向矢量,对米波雷达测高模型进行良好的匹配,将反射波和直射波当作一个波束处理,避免了目标低仰角时反射波和直射波夹角过小问题,解决了普通波束形成“瑞利限”对雷达低仰角测高时的限制,消除了多径效应对米波雷达测高的影响,提高了测角精度。
(3)本发明对米波雷达测高模型进行近似处理,保证运算精度的同时减少了计算雷达波反射角和波程差时的数值范围,满足单精度浮点数数值精度,适应目前DSP的计算结构,极大的减少了运算时间。
附图说明
图1是本发明的测高流程图;
图2是本发明中雷达接收目标信号模型图;
图3a是本发明计算反射角随直射角变化的曲线图;
图3b是本发明计算波程差随直射角变化的曲线图;
图4a是本发明中雷达波反射模型与精确雷达波反射模型的反射角相对误差曲线图;
图4b是本发明中雷达波反射模型与精确雷达波反射模型的波程差相对误差曲线图;
图5a是本发明针对某雷达实测数据处理的测角结果图;
图5b是本发明针对某雷达实测数据处理的测高结果图。
具体实施方式
下面结合附图详细说明本发明的内容和效果。
实施例1:
本发明是一种基于合成导向矢量的波束形成米波雷达测高方法,米波雷达在反隐身、反辐射导弹方面发挥着重要作用,而米波雷达测高问题是其中急需解决的难题之一,可用于雷达目标跟踪,参照图1,本发明的基于合成导向矢量的波束形成米波雷达测高方法包括以下步骤:
步骤1:从雷达接收的回波中提取目标信号x(tk),并根据目标距离信息确定目标仰角搜索范围。
从雷达接收的回波中提取目标信号x(tk)的过程通常是:对雷达接收的回波进行模数转换、幅相校正、IQ通道分离、脉冲压缩、杂波抑制和消除、动目标检测、动目标显示和恒虚警等处理后,得到目标信号和目标距离信息。
米波雷达一般用作预警雷达,探测距离一般为几百千米,目前飞机飞行高度一般在两万米以下。根据目标的距离信息确定测高时仰角搜索范围:目标距离大于300千米时,仰角搜索范围为0°~3°;目标距离为200千米~300千米时,仰角搜索范围为0°~6°;目标距离为100千米~200千米时,仰角搜索范围为0°~15°;目标距离小于100千米时,仰角搜索范围为0°~25°。
步骤2:根据雷达和目标的地理位置关系建立雷达波反射模型,使用目标搜索仰角θd计算反射点与雷达中心对应的地面弧长G1和反射点与目标对应的地面弧长G2,参见图2。
步骤3:分别计算目标和雷达中心在反射点切线上的等效高度,记为hr′和ht′,然后计算雷达波的反射角θi及反射波与直射波的波程差ΔR,参见图2。
步骤4:利用直射波导向矢量a(θd)和反射波导向矢量a(θi)计算搜索角度的合成导向矢量as(θd),通过加权求和将反射波和直射波作为一个波束处理,ρexp(-j2πΔR/λ)为权值,代表反射波与直射波由于地面反射和波程差引起的相位差。
as(θd)=a(θd)+ρexp(-j2πΔR/λ)a(θi)
其中ρ为地面反射系数,对于水平极化波,ρ一般取值为-0.95~-1之间,λ为雷达波波长。
步骤5:对雷达回波进行波束形成处理,利用合成导向矢量将雷达各阵元回波进行加权求和,将雷达波束“导向”到搜索仰角方向,计算归一化阵列平均输出功率S(θd),归一化阵列平均输出功率由每个回波脉冲的雷达输出功率进行平均后对合成导向矢量进行归一化处理得到,因此本发明测高计算运算量小,便于雷达实时信号处理。
步骤6:每个搜索仰角对应的归一化阵列平均输出功率计算完成后,选取最大归一化阵列平均输出功率对应的角度为目标仰角θ,并利用余弦定理计算目标的高度ht:
其中Rd为目标与雷达距离,hr为雷达中心高度,re为等效地球半径,ht即为所求目标高度,为进一步雷达信息处理提供了准确数据。
在现有技术中,基于波瓣***的米波雷达测高方法和普通波束形成算法均计算简单但精度较低,尤其在目标低仰角情况下,测高效果均不理想;基于导向矢量最大似然算法虽然测高效果较好但运算量较大,在现代雷达高实时性、目标高密集性的情况下不能满足战场需求。本发明利用合成导向矢量作为权值,对目标角度估计时采用波束形成算法,并对反射模型进行了近似处理,在对目标低仰角反射进行很好匹配的同时极大的减少了运算量,提高了测高速度。
实施例2:
基于合成导向矢量的波束形成米波雷达测高方法同实施例1,
步骤1:从雷达接收的回波中提取目标信号x(tk),并根据目标距离信息确定目标仰角搜索范围。
本发明中雷达接收目标信号的模型如图2所示。图2中一个远场的窄带信号入射到M个阵元组成的均匀线阵,选取雷达阵列中心A为参考阵元,hr为雷达中心高度,d为雷达阵元间隔,θa为天线的倾斜角,re为等效地球半径,Rd为目标到雷达中心的直射波波程,θd为直射波到雷达中心的搜索仰角,θi为反射波到雷达中心的反射角,Ri=Ri1+Ri2为目标到雷达中心的反射波波程,ΔR=Ri-Rd为雷达中心反射波和直射波的波程差,K为脉冲总数目,k为当前脉冲数。
雷达阵列接收的信号矢量为:
x(tk)=[a(θd)a(θi)][Sd(tk)Si(tk)]T+n(tk),k=1,…,K
其中a(θd)=[ad(1),ad(2),…,ad(M)]T是以雷达中心A为参考点的直射波导向矢量,a(θi)=[ai(1),ai(2),…,ai(M)]T为对应的反射波导向矢量,上标T表示转置。上式中
其中m=1,2,…,M,m表不第m个阵元,M表示阵元个数。
Sd(tk)和Si(tk)分别是雷达中心接收的直射波和反射波复包络;n(tk)是均值为零、方差为σ2的高斯白噪声。
对雷达中心的参考阵元来说,si(tk)=ρexp(-j2πΔR/λ)sd(tk),其中ρ为地面反射系数,因此有
x(tk)=[a(θd)+ρexp(-j2πΔR/λ)a(θi)]sd(tk)+n(tk),k=1,…,K
步骤2:建立雷达波反射模型,根据搜索仰角θd,计算反射点与雷达中心在地面投影点的弧长G1和反射点与目标在地面投影点的弧长G2。
图2中B点为反射波与地面交点,即反射点。ht为目标高度,φ为目标与雷达中心的地心夹角,G1为反射点与雷达中心地面投影点的弧长,G2为反射点与目标地面投影点的弧长。计算反射点与雷达中心的地面弧长G1和反射点与目标的地面弧长G2的步骤包括:
(2a)建立雷达波反射模型,应用余弦定理计算目标搜索角度θd对应的目标高度ht。
(2b)应用余弦定理计算目标与雷达中心的地心夹角φ和地面弧长G。
G=φ·re
(2c)应用fishback方程计算反射点与雷达中心在地面投影点的弧长G1和反射点与目标在地面投影点的弧长G2。
G2=G-G1
步骤3:利用得到的两弧长G1和G2再分别计算目标和雷达中心在反射点切线上的等效高度,记为hr′和ht′,然后计算雷达波的反射角θi及反射波与直射波的波程差ΔR。
(3a)计算雷达中心等效高度hr′和目标等效高度ht′。
由 可得:
等效雷达中心高度hr′和目标高度ht′可以应用正弦定理计算为
(3b)计算雷达波的反射角θi和波程差ΔR。
sin(θd)=(ht′-hr′)/Rd
sin(θi)=-(ht′+hr′)/Rs
在低仰角情况下,Rd≈Ri,且由于θd≈θi可得
θi≈-asin(sin(θd)+2·hr′/Rd)
步骤4:利用直射波导向矢量a(θd)和反射波导向矢量a(θi)计算搜索角度的合成导向矢量as(θd)。
直射波导向矢量为a(θd)=[ad(1),ad(2),…,ad(M)]T,其中
反射波导向矢量为a(θi)=[ai(1),ai(2),…,ai(M)]T,其中
则合成导向矢量为
as(θd)=a(θd)+ρexp(-j2πΔR/λ)a(θi)
步骤5:通过对雷达回波进行波束形成处理,计算归一化阵列平均输出功率S(θd)。
(5a)计算第k个脉冲目标信号x(tk)对应搜索角度θd的波束形成值y(θd,k)。
(5b)计算第k个脉冲目标信号x(tk)对应搜索角度θd的阵列输出功率p(θd,k)。
p(θd,k)=y(θd,k)y(θd,k)H
(5c)每个脉冲目标信号对应的阵列输出功率计算完成后,计算所有K个脉冲的阵列平均输出功率P(θ)。
(5d)由于合成导向矢量模值不为1,因此需要对阵列输出功率进行归一化处理,消除此影响,得到归一化阵列输出平均功率S(θd)。
步骤6:每个搜索仰角对应的归一化阵列平均输出功率计算完成后,根据波束形成原理,最大功率值指向的角度即为目标角度,因此选取最大归一化的阵列平均输出功率对应的角度为目标仰角θ。
计算出目标仰角后,利用余弦定理可计算目标高度ht。
其中Rd为目标与雷达距离,hr为雷达中心高度,re为等效地球半径。
参见图5b,图5b是采用本发明测高方法获得的目标高度随目标飞行的实时测高曲线图。
本发明使用简化的雷达反射模型,对米波雷达低仰角测高进行良好的匹配,避免了最大似然算法中由于多径效应导致的多维搜索,同时避免了“瑞利限”导致普通波束形成无法工作的限制,消除了多径效应对米波雷达测高的影响,提高了测角精度,为雷达测量***提供了精度保证。
实施例3:
基于合成导向矢量的波束形成米波雷达测高方法同实施例1-2,
本发明的效果可以通过以下仿真进一步说明。
仿真内容:雷达波反射模型分析。
仿真条件:
假设目标距离100千米,目标仰角范围0.25°~16°,角度步进0.05°,共计算316点,DSP采用ADSP-TS201,内核频率为600MHz,定义反射角相对误差为 波程差相对误差为 分别利用本发明和4/3地球精确曲面反射模型计算目标反射角和波程差,并计算两者的相对误差,统计DSP计算指令周期数和计算时间。
仿真结果:本发明对雷达波反射模型进行了近似处理,在受地球曲率影响较大的部分采用曲面反射模型,然后应用低仰角测高时的近似关系计算目标和雷达在平面反射模型的等效高度,并利用平面反射模型计算反射角和波程差,减少了数值动态范围,提高了运算速度。
图3a是采用本发明得到的目标反射角随入射角变化的曲线,其中横坐标为入射角,纵坐标为反射角。从图3a中可以看出,反射角在-0.25~-16度之间,反射角约等于入射角。
图3b是采用本发明得到的波程差随入射角变化的曲线,其中横坐标为入射角,纵坐标为波程差。从图3b中可以看出,波程差0.13~7.05米之间变化,与入射角成线性关系。
图4a是本发明与精确雷达波反射模型的反射角相对误差曲线图,其中横坐标为入射角,纵坐标为反射角相对误差。从图4a中可以看出,采用本发明计算反射角与4/3地球精确曲面反射模型的相对误差较小,均为0.1以下,满足计算精度要求。
图4b是本发明与精确雷达波反射模型的波程差相对误差曲线图,其中横坐标为入射角,纵坐标为波程差相对误差。从图4b中可以看出,采用本发明计算波程差时与4/3地球精确曲面反射模型的相对误差也为0.1以下,满足计算精度要求。
表1为本发明模型与精确反射模型的计算复杂度比较。
表1本发明模型和精确反射模型的计算复杂度比较
算法 | DSP指令周期数 | 运行时间(ms) |
本发明 | 1012148 | 1.69 |
精确反射模型 | 4397456 | 7.33 |
从表1中可以看出,本发明模型在计算反射角和波程差时的计算复杂度比精确反射模型降低接近3/4,在保证测高精度的基础上,提高了雷达测高的速度。
实施例4:
基于合成导向矢量的波束形成米波雷达测高方法同实施例1-2,
本发明的效果可以通过以下计算量分析进一步说明。
仿真内容:雷达测高计算复杂度比较。
仿真条件:假设雷达天线阵元数为M,回波脉冲数为K。
仿真结果:表2列出本发明和合成导向矢量最大似然算法计算复杂度对比。
表2本发明和合成导向矢量最大似然算法的计算复杂度对比
从表2中可以看出,本发明复杂度为O{723},而合成导向矢量最大似然算法的复杂度为O{16176},可以看出本发明比现有测高算法的计算量降低一个数量级。随着雷达阵元数目和回波脉冲数的增多,二者计算量差别将更明显。可见本发明能大大降低米波雷达低仰角测高处理的运算量,提高目标仰角估计的计算速度,为米波雷达测高提供了一种工程化利用的快速测高方法。
实施例5:
基于合成导向矢量的波束形成米波雷达测高方法同实施例1-2,
本发明的效果可以通过以下实测数据处理结果进一步说明。
仿真内容:某雷达实测数据处理结果。
仿真条件:雷达阵元个数为24,阵元间距d=λ/2,快拍数为13,DSP采用ADSP-TS201,内核频率为600MHz。分别利用本发明方法、普通波束形成方法和合成导向矢量最大似然算法对某条民航航迹进行高度估计。
仿真结果:测角和测高处理结果如图5a和5b所示,图中横轴表示目标的点迹号,纵轴分别表示测得的仰角和高度。图5a是本发明针对某雷达实测数据处理的测角结果图,从图5a中可以看出,当目标仰角低于2°时,普通波束形成算法由于反射波和直射波夹角小于一个波束宽度无法正常工作,本发明和合成导向矢量最大似然算法测角精度较高。图5b是本发明针对某雷达实测数据处理的测高结果图,从图5b中可以看出,当目标点迹号大于80时,此时由于目标仰角过小普通波束形成测高误差巨大,本发明和合成导向矢量最大似然算法仍能有效工作,测高精度较高。
综上,发明的基于合成导向矢量的波束形成米波雷达测高方法,主要解决现有技术测高时信号处理运算量大、计算复杂的问题。本发明从雷达回波中提取目标信号,确定目标仰角搜索范围;建立雷达波反射模型,用搜索仰角计算反射点与雷达中心的地面弧长和反射点与目标的地面弧长;计算目标和雷达中心的等效高度后计算反射角和反射波与直射波波程差;计算合成导向矢量;对回波数据波束形成后选取最大归一化功率值对应的角度作为目标仰角,并计算目标高度。本发明将波束形成和合成导向矢量相结合,并将目标反射模型进行简化处理,保证米波雷达的仰角估计精度,提高米波雷达对目标仰角估计的计算速度,也就是测高速度,简化测高计算过程,运算简单,使米波雷达的实时信号处理更加容易。
Claims (4)
1.一种基于合成导向矢量的波束形成米波雷达测高方法,包括以下步骤:
步骤1:从雷达接收的回波中提取目标信号x(tk),并确定目标仰角搜索范围;
步骤2:建立雷达波反射模型,根据目标搜索仰角θd,计算反射点与雷达中心的地面弧长和反射点与目标的地面弧长,分别记为G1和G2;
步骤3:分别计算目标和雷达中心在反射点切线上的等效高度,记为hr′和ht′,然后计算雷达波的反射角θi及反射波与直射波的波程差ΔR;
步骤4:利用直射波导向矢量a(θd)和反射波导向矢量a(θi)计算搜索角度的合成导向矢量as(θd);
as(θd)=a(θd)+ρexp(-j2πΔR/λ)a(θi)
其中ρ为地面反射系数,λ为雷达波波长;
步骤5:对雷达回波进行波束形成处理,计算归一化阵列平均输出功率S(θd);
步骤6:选取最大归一化阵列平均输出功率对应的角度为目标仰角θ,并利用余弦定理计算目标的高度ht:
其中Rd为目标与雷达距离,hr为雷达中心高度,re为等效地球半径。
2.根据权利要求1所述的米波雷达测高方法,其中步骤2中根据搜索仰角计算计算反射点与雷达中心的地面弧长G1和反射点与目标的地面弧长G2,按如下步骤进行:
(2a)建立雷达波反射模型,计算目标搜索角度θd对应的目标高度ht;
(2b)计算目标与雷达中心的地心夹角φ和地面投影点的弧长G;
G=φ·re
(2c)计算反射点与雷达中心的地面弧长G1和反射点与目标的地面弧长G2。
G2=G-G1
其中p和ξ均为中间变量,无具体物理含义。
3.根据权利要求2所述的米波雷达测高方法,其中步骤3中分别计算目标和雷达中心在反射点切线上的等效高度,记为hr′和ht′,然后计算雷达波的反射角θi及反射波与直射波的波程差ΔR,按如下步骤进行:
(3a)计算雷达中心等效高度hr′和目标等效高度ht′;
(3b)计算雷达波的反射角θi和波程差ΔR。
θi≈-asin(sin(θd)+2·hr′/Rd)
4.根据权利要求3所述的米波雷达测高方法,其中步骤5中对雷达回波进行波束形成处理,计算归一化阵列平均输出功率S(θd),按如下步骤进行:
(5a)计算第k个脉冲目标信号x(tk)对应搜索角度θd的波束形成值y(θd,k);
(5b)计算第k个脉冲目标信号x(tk)对应搜索角度θd的阵列输出功率p(θd,k);
p(θd,k)=y(θd,k)y(θd,k)H
(5c)计算K个脉冲的阵列平均输出功率P(θ);
(5d)计算归一化阵列输出平均功率S(θd)
式中H表示共轭转制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310192066.0A CN103293526B (zh) | 2013-05-13 | 2013-05-13 | 基于合成导向矢量的波束形成米波雷达测高方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310192066.0A CN103293526B (zh) | 2013-05-13 | 2013-05-13 | 基于合成导向矢量的波束形成米波雷达测高方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103293526A true CN103293526A (zh) | 2013-09-11 |
CN103293526B CN103293526B (zh) | 2015-06-03 |
Family
ID=49094719
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310192066.0A Expired - Fee Related CN103293526B (zh) | 2013-05-13 | 2013-05-13 | 基于合成导向矢量的波束形成米波雷达测高方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103293526B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105182325A (zh) * | 2015-08-13 | 2015-12-23 | 西安电子科技大学 | 基于秩1约束的米波mimo雷达低仰角目标测高方法 |
CN105785337A (zh) * | 2016-01-22 | 2016-07-20 | 西安电子科技大学 | 复杂地形下的米波雷达低仰角目标测高方法 |
CN105954746A (zh) * | 2016-04-29 | 2016-09-21 | 西安电子科技大学 | 基于广播自动相互监视信号的地形修正米波雷达测高方法 |
CN106646346A (zh) * | 2016-09-12 | 2017-05-10 | 西安电子科技大学 | 米波雷达有多径反射时对高仰角目标的仰角测量方法 |
CN108535698A (zh) * | 2018-04-04 | 2018-09-14 | 西安电子科技大学 | 基于波束空间的米波雷达低仰角估计方法 |
RU2720948C1 (ru) * | 2019-12-08 | 2020-05-15 | Закрытое акционерное общество Научно-исследовательский центр «РЕЗОНАНС» (ЗАО НИЦ «РЕЗОНАНС») | Способ измерения угла места воздушного объекта в метровом диапазоне электромагнитных волн с использованием электронной карты местности |
CN111814688A (zh) * | 2020-07-09 | 2020-10-23 | 成都傅立叶电子科技有限公司 | 基于FFTc的数字波束形成角度捕获方法及装置、存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5557282A (en) * | 1988-10-11 | 1996-09-17 | Itt Corporation | Height finding antenna apparatus and method of operation |
CN101995566A (zh) * | 2010-10-15 | 2011-03-30 | 西安电子科技大学 | 二维数字阵列雷达数字波束形成***及方法 |
CN102288944A (zh) * | 2011-05-12 | 2011-12-21 | 西安电子科技大学 | 基于地形匹配的数字阵列米波雷达超分辨测高方法 |
CN102830107A (zh) * | 2012-09-04 | 2012-12-19 | 南京信息工程大学 | 测量云中固、液态水含量的激光雷达探测方法与*** |
-
2013
- 2013-05-13 CN CN201310192066.0A patent/CN103293526B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5557282A (en) * | 1988-10-11 | 1996-09-17 | Itt Corporation | Height finding antenna apparatus and method of operation |
CN101995566A (zh) * | 2010-10-15 | 2011-03-30 | 西安电子科技大学 | 二维数字阵列雷达数字波束形成***及方法 |
CN102288944A (zh) * | 2011-05-12 | 2011-12-21 | 西安电子科技大学 | 基于地形匹配的数字阵列米波雷达超分辨测高方法 |
CN102830107A (zh) * | 2012-09-04 | 2012-12-19 | 南京信息工程大学 | 测量云中固、液态水含量的激光雷达探测方法与*** |
Non-Patent Citations (1)
Title |
---|
陈伯孝等: "约束MUSIC算法在米波雷达低仰角测高中的应用", 《电子学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105182325A (zh) * | 2015-08-13 | 2015-12-23 | 西安电子科技大学 | 基于秩1约束的米波mimo雷达低仰角目标测高方法 |
CN105785337A (zh) * | 2016-01-22 | 2016-07-20 | 西安电子科技大学 | 复杂地形下的米波雷达低仰角目标测高方法 |
CN105785337B (zh) * | 2016-01-22 | 2018-05-15 | 西安电子科技大学 | 复杂地形下的米波雷达低仰角目标测高方法 |
CN105954746A (zh) * | 2016-04-29 | 2016-09-21 | 西安电子科技大学 | 基于广播自动相互监视信号的地形修正米波雷达测高方法 |
CN106646346A (zh) * | 2016-09-12 | 2017-05-10 | 西安电子科技大学 | 米波雷达有多径反射时对高仰角目标的仰角测量方法 |
CN106646346B (zh) * | 2016-09-12 | 2019-02-19 | 西安电子科技大学 | 米波雷达有多径反射时对高仰角目标的仰角测量方法 |
CN108535698A (zh) * | 2018-04-04 | 2018-09-14 | 西安电子科技大学 | 基于波束空间的米波雷达低仰角估计方法 |
CN108535698B (zh) * | 2018-04-04 | 2021-07-20 | 西安电子科技大学 | 基于波束空间的米波雷达低仰角估计方法 |
RU2720948C1 (ru) * | 2019-12-08 | 2020-05-15 | Закрытое акционерное общество Научно-исследовательский центр «РЕЗОНАНС» (ЗАО НИЦ «РЕЗОНАНС») | Способ измерения угла места воздушного объекта в метровом диапазоне электромагнитных волн с использованием электронной карты местности |
CN111814688A (zh) * | 2020-07-09 | 2020-10-23 | 成都傅立叶电子科技有限公司 | 基于FFTc的数字波束形成角度捕获方法及装置、存储介质 |
CN111814688B (zh) * | 2020-07-09 | 2023-10-13 | 成都傅立叶电子科技有限公司 | 基于FFTc的数字波束形成角度捕获方法及装置、存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN103293526B (zh) | 2015-06-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103293526B (zh) | 基于合成导向矢量的波束形成米波雷达测高方法 | |
CN103308909B (zh) | 基于多参数交替搜索的米波雷达低仰角测高方法 | |
CN103018730B (zh) | 分布式子阵波达方向估计方法 | |
CN102288944B (zh) | 基于地形匹配的数字阵列米波雷达超分辨测高方法 | |
CN103760546B (zh) | 一种雷达用低空目标波达方向估计方法 | |
CN103353596A (zh) | 基于压缩感知的波束空间域米波雷达测高方法 | |
CN103364772B (zh) | 基于实数域广义多重信号分类算法的目标低仰角估计方法 | |
CN103885054B (zh) | 一种基于分布源反射模型的米波雷达低仰角测高方法 | |
CN107340495A (zh) | 一种基于阵列雷达的目标波达方向快速估计方法 | |
CN104155648A (zh) | 基于阵列数据重排的高频地波雷达单次快拍music测向方法 | |
CN102393525A (zh) | 子空间投影的导航干扰抑制与信号增强方法 | |
CN103197294A (zh) | 多频融合最大似然低空目标仰角估计方法 | |
CN110196414A (zh) | 一种基于补偿天线方向图误差的天线波束指向方法 | |
CN107607943A (zh) | 基于干涉相位辅助的延迟多普勒雷达高度表的测高方法 | |
CN102662158B (zh) | 一种对传感器天线阵列接收信号的快速处理方法 | |
CN102749067A (zh) | 自适应多级计算的精确反演海面高程遥测装置 | |
CN104898119A (zh) | 一种基于相关函数的动目标参数估计方法 | |
CN102608587B (zh) | 基于非线性最小二乘的空中机动目标检测方法 | |
CN106950529A (zh) | 声矢量近场源esprit和music参数估计方法 | |
CN109188019A (zh) | 基于多重信号分类算法的三维风速风向测量方法 | |
CN108535698A (zh) | 基于波束空间的米波雷达低仰角估计方法 | |
CN109597047A (zh) | 基于有监督深度神经网络的米波雷达doa估计方法 | |
CN104793177A (zh) | 基于最小二乘法的麦克风阵列测向方法 | |
CN104166129A (zh) | 一种实波束雷达迭代最小均方误差角超分辨方法 | |
CN108983169B (zh) | 一种基于数字高程模型的米波雷达地形修正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150603 Termination date: 20200513 |