CN106646346B - 米波雷达有多径反射时对高仰角目标的仰角测量方法 - Google Patents

米波雷达有多径反射时对高仰角目标的仰角测量方法 Download PDF

Info

Publication number
CN106646346B
CN106646346B CN201610817425.0A CN201610817425A CN106646346B CN 106646346 B CN106646346 B CN 106646346B CN 201610817425 A CN201610817425 A CN 201610817425A CN 106646346 B CN106646346 B CN 106646346B
Authority
CN
China
Prior art keywords
elevation
measurement
steering vector
indicate
beam position
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
CN201610817425.0A
Other languages
English (en)
Other versions
CN106646346A (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN201610817425.0A priority Critical patent/CN106646346B/zh
Publication of CN106646346A publication Critical patent/CN106646346A/zh
Application granted granted Critical
Publication of CN106646346B publication Critical patent/CN106646346B/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
    • 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
    • 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开一种米波雷达有多径反射时对高仰角目标的仰角测量方法,该方法包括如下步骤:(1)确定仰角测量的左、右波束指向;(2)确定仰角测量的中心指向;(3)确定凹口区间;(4)计算特征值对应的特征矢量;(5)构成信号子空间;(6)确定加权值;(7)建立鉴角曲线;(8)确定接收信号内目标的仰角值。本发明可用于米波雷达存在多径反射信号的干扰时,对高仰角目标的仰角进行测量。本发明提高了米波雷达存在多径反射信号时,对高仰角目标的仰角测量精度。

Description

米波雷达有多径反射时对高仰角目标的仰角测量方法
技术领域
本发明属于通信技术领域,更进一步涉及雷达技术领域中的一种米波雷达有多径反射时对高仰角目标的仰角测量方法。本发明可用于改善米波雷达对高仰角目标测仰角过程中,由于存在地面或海面引起的多径反射回波所导致的目标仰角估计精度下降的问题。
背景技术
米波雷达由于波束宽,在对低仰角目标定位时,雷达波束打地,使得目标的直达波和地(海)面反射的多径反射波在天线波束主瓣内叠加,这两个信号强相关,这组强相关的信号同时被雷达天线接收,使得雷达不能正确的估计出目标的仰角,在对高仰角目标进行仰角测量时,由于天线波束旁瓣内存在多径反射信号,仰角测量精度低,在对仰角测量精度要求高的雷达系统中,多径反射信号的干扰仍然不可忽略,值得重视。
张光义院士在其发表的论著《相控阵雷达原理[M]》(国防工业出版社,1994北京)中公开了一种单脉冲测角方法。该单脉冲测角方法是利用单个脉冲的回波信号测出目标真实位置的测角方法。传统的单脉冲测角方法中的比幅单脉冲测角方法是利用在观测方向上形成的两个指向不同的接收波束,比较这两个接收波束的幅度差异,查表可知每个目标的估计仰角。传统的单脉冲测角方法忽略了米波雷达的天线方向图内的多径反射信号的干扰,因此用传统单脉冲测角方法在米波雷达在存在多径反射信号时对高仰角目标进行仰角测量存在测量精度低的问题。
发明内容
本发明针对上述现有技术存在的不足,提出一种米波雷达有多径反射时对高仰角目标的仰角测量方法。本发明提出一种计算左、右波束(左波束即下波束,右波束即上波束)指向方向图的加权值的方法,通过直接抑制米波雷达的多径反射信号的干扰,以解决现有技术在米波雷达存在多径反射信号时对高仰角目标的仰角测量精度低的问题。
本发明的基本思路是:确定仰角测量的左、右波束指向和中心指向,确定凹口区间,计算特征值对应的特征矢量,构成信号子空间,确定加权值,建立鉴角曲线,确定接收信号内目标的仰角值。
实现本发明目的的具体步骤如下:
(1)确定仰角测量的左、右波束指向:
(1a)在有多径反射时,米波雷达发射脉冲中任意选定两个相互交叠的波束;
(1b)将左边的波束作为仰角测量的左波束,右边的波束作为仰角测量的右波束,将仰角测量的左波束的指向作为仰角测量的左波束指向,仰角测量的右波束的指向作为仰角测量的右波束指向;
(2)确定仰角测量的中心指向:
将仰角测量的左波束与右波束交叠位置的指向作为仰角测量的中心指向;
(3)确定凹口区间:
按照下式,确定仰角测量的左、右波束指向方向图的凹口区间:
其中,D表示仰角测量的左、右波束指向方向图的凹口区间,θL表示仰角测量的左波束指向,θR表示仰角测量的右波束指向,Δ表示仰角测量的左、右波束指向方向图的凹口区间大小的二分之一;
(4)计算特征值对应的特征矢量:
(4a)按照导向矢量公式,计算仰角测量的左、右波束指向方向图的凹口区间上的导向矢量中的任意一个导向矢量;
(4b)按照自相关矩阵公式,计算仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵;
(4c)按照特征矢量公式,计算仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵的特征值及该特征值对应的特征矢量;
(5)构成信号子空间:
将仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵的所有特征值从大到小排序,选取排序后的特征值中的前m个特征值所对应的特征矢量,m的取值取决于雷达的工作环境,将所选取特征矢量组成信号子空间;
(6)确定加权值:
按照下式,计算仰角测量的左、右波束指向方向图的加权值:
w1=[I-G*(GH*G)-1*GH]*w3
w2=[I-G*(GH*G)-1*GH]*w4
其中,w1、w2分别表示仰角测量的左、右波束指向方向图的加权值,I表示单位矩阵,G表示信号子空间,*表示矩阵相乘操作,H表示共轭转置操作,(·)-1表示矩阵求逆操作,w3、w4分别表示采用传统比幅单脉冲测角方法确定的仰角测量的左、右波束指向方向图的加权值;
(7)建立鉴角曲线:
(7a)在有多径反射时,米波雷达发射脉冲中任意选定两个相互交叠的波束,将该交叠位置的指向作为仰角测量的中心指向,取一段以仰角测量的中心指向为中心的区间作为鉴角曲线上仰角测量区间;
(7b)按照导向矢量公式,计算鉴角曲线上仰角测量区间上的导向矢量中的任意一个导向矢量;
(7c)按照加权求和公式,计算鉴角曲线上仰角测量区间上的导向矢量中的任意一个导向矢量的左、右加权求和值;
(7d)按照下式,计算鉴角曲线上仰角测量区间上的导向矢量中的任意一个导向矢量的比值:
其中,fk表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的比值,pk1、pk2分别表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的左、右加权求和值;
(7e)连接鉴角曲线上仰角测量区间内上所有的导向矢量的比值,得到鉴角曲线;
(8)确定接收信号内目标的仰角值:
(8a)按照加权求和公式,计算在有多径反射时,米波雷达接收信号的左、右加权求和值;
(8b)按照下式,计算接收信号内目标的比值:
其中,g表示接收信号内目标的比值,q1、q2分别表示在有多径反射时,米波雷达接收信号的左、右加权求和值;
(8c)按照下式,计算鉴角曲线上仰角测量区间上的导向矢量中的任意一个导向矢量的比值与接收信号内目标的比值的差值:
ck=|fk-g|
其中,ck表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的比值与接收信号内目标的比值的差值,|·|表示取绝对值操作,fk表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的比值,g表示接收信号内目标的比值;
(8d)取鉴角曲线上仰角测量区间上所有的导向矢量的比值与接收信号内目标的比值的差值中最小的差值,将该差值对应的导向矢量的仰角值作为在有多径反射时,米波雷达接收信号内目标的仰角值。
本发明与现有的技术相比具有以下优点:
第一,由于本发明在仰角测量的左、右波束指向方向图上确定了凹口区间,使得米波雷达能够抑制仰角测量的左、右波束指向方向图的旁瓣内的多径反射信号的干扰,克服了现有技术忽略方向图内的多径反射信号的干扰,在对高仰角目标仰角测量时精度低的问题,使得本发明具有对高仰角目标仰角测量精度高的优点。
第二,由于本发明采用了由仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵的选取的特定的特征矢量构成的信号子空间,减少了计算量,使得本发明具有利于工程实现的优点。
附图说明
图1为本发明的流程图;
图2为本发明的米波雷达多径反射信号示意图;
图3为本发明的仿真实验中左、右波束指向方向图;
图4为本发明的仿真实验中鉴角曲线示意图;
图5为本发明的仿真实验中仰角测量精度图。
具体实施方式
下面结合附图对本发明作进一步的详细描述。
参照附图1,本发明的具体步骤如下:
步骤1,确定仰角测量的左、右波束指向。
在有多径反射时,米波雷达发射脉冲中任意选定两个相互交叠的波束。
将左边的波束作为仰角测量的左波束,右边的波束作为仰角测量的右波束,将仰角测量的左波束的指向作为仰角测量的左波束指向,仰角测量的右波束的指向作为仰角测量的右波束指向。
步骤2,确定仰角测量的中心指向。
将仰角测量的左波束与右波束交叠位置的指向作为仰角测量的中心指向。
步骤3,确定凹口区间。
按照下式,确定仰角测量的左、右波束指向方向图的凹口区间:
其中,D表示仰角测量的左、右波束指向方向图的凹口区间,θL表示仰角测量的左波束指向,θR表示仰角测量的右波束指向,Δ表示仰角测量的左、右波束指向方向图的凹口区间大小的二分之一。
步骤4,计算特征值对应的特征矢量。
按照导向矢量公式,计算仰角测量的左、右波束指向方向图的凹口区间上的导向矢量中的任意一个导向矢量:
其中,a(θi)表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量中的第i个导向矢量,θi表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量中的第i个导向矢量的仰角值,e表示指数符号,j表示虚数,π表示圆周率,d表示在有多径反射时,米波雷达天线的阵元间距,N表示在有多径反射时,米波雷达天线的阵元总数,sin表示正弦操作,λ表示在有多径反射时,米波雷达天线的发射信号的波长,T表示转置操作。
按照自相关矩阵公式,计算仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵:
R=a*aH
其中,R表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵,a表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量,*表示矩阵相乘操作,H表示共轭转置操作。
a的表达式如下:
a=[a(θ1),a(θ2)..........a(θn)]
其中,a(θ1),a(θ2),a(θn)分别表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量中的第1个,第2个,第n个导向矢量,θ1,θ2,θn分别表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量中的第1个,第2个,第n个导向矢量的仰角值。
按照特征矢量公式,计算仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵的特征值及该特征值对应的特征矢量:
R*Vi=λiVi
其中,R表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵,*表示矩阵相乘操作,Vi表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵的第i个特征值对应的特征矢量,λi表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵的第i个特征值。
步骤5,构成信号子空间。
将仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵的所有特征值从大到小排序,选取排序后的特征值中的前m个特征值所对应的特征矢量,m的取值取决于雷达的工作环境,将所选取特征矢量组成信号子空间。
步骤6,确定加权值。
按照下式,计算仰角测量的左、右波束指向方向图的加权值:
w1=[I-G*(GH*G)-1*GH]*w3
w2=[I-G*(GH*G)-1*GH]*w4
其中,w1、w2分别表示仰角测量的左、右波束指向方向图的加权值,I表示单位矩阵,G表示信号子空间,*表示矩阵相乘操作,H表示共轭转置操作,(·)-1表示矩阵求逆操作,w3、w4分别表示采用传统比幅单脉冲测角方法确定的仰角测量的左、右波束指向方向图的加权值。
G的表达式如下:
G=[V1,V2,..,Vm-1,Vm]
其中,G表示信号子空间,V1,V2,..,Vm-1,Vm分别表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵的选定的m个特征值所对应的m个特征矢量。
w3、w4的表达式分别如下:
w3、w4分别表示采用传统比幅单脉冲测角方法确定的仰角测量的左、右波束指向方向图的加权值,e表示指数符号,j表示虚数,π表示圆周率,d表示在有多径反射时,米波雷达天线的阵元间距,N表示在有多径反射时,米波雷达天线的阵元总数,sin表示正弦操作,θL表示仰角测量的左波束指向,θR表示仰角测量的右波束指向,λ表示在有多径反射时,米波雷达天线的发射信号的波长,T表示转置操作。
步骤7,建立鉴角曲线。
在有多径反射时,米波雷达发射脉冲中任意选定两个相互交叠的波束,将该交叠位置的指向作为仰角测量的中心指向,取一段以仰角测量的中心指向为中心的区间作为鉴角曲线上仰角测量区间。
按照导向矢量公式,计算鉴角曲线上仰角测量区间上的导向矢量中的任意一个导向矢量:
其中,b(θk)表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量中的第k个导向矢量,θk表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量中的第k个导向矢量的仰角值,e表示指数符号,j表示虚数,π表示圆周率,d表示在有多径反射时,米波雷达天线的阵元间距,N表示在有多径反射时,米波雷达天线的阵元总数,sin表示正弦操作,λ表示在有多径反射时,米波雷达天线的发射信号的波长,T表示转置操作。
按照加权求和公式,计算鉴角曲线上仰角测量区间上的导向矢量中的任意一个导向矢量的左、右加权求和值:
pk1=|w1 H*b(θk)|
pk2=|w2 H*b(θk)|
其中,pk1、pk2分别表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的左、右加权求和值,|·|表示取绝对值操作,w1、w2分别表示仰角测量的左、右波束指向方向图的加权值,H表示共轭转置操作,*表示矩阵相乘操作,b(θk)表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量,θk表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的仰角值。
按照下式,计算鉴角曲线上仰角测量区间上的导向矢量中的任意一个导向矢量的比值:
其中,fk表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的比值,pk1、pk2分别表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的左、右加权求和值。
连接鉴角曲线上仰角测量区间内上所有的导向矢量的比值,得到鉴角曲线。
步骤8,确定接收信号内目标的仰角值。
按照加权求和公式,计算在有多径反射时,米波雷达接收信号的左、右加权求和值:
q1=|w1 H*X|
q2=|w2 H*X|
其中,q1、q2分别表示在有多径反射时,米波雷达接收信号的左、右加权求和值,|·|表示取绝对值操作,w1、w2分别表示仰角测量的左、右波束指向方向图的加权值,H表示共轭转置操作,*表示矩阵相乘操作,X表示米波雷达接收的信号矩阵。
按照下式,计算接收信号内目标的比值:
其中,g表示接收信号内目标的比值,q1、q2分别表示在有多径反射时,米波雷达接收信号的左、右加权求和值。
按照下式,计算鉴角曲线上仰角测量区间上的导向矢量中的任意一个导向矢量的比值与接收信号内目标的比值的差值:
ck=|fk-g|
其中,ck表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的比值与接收信号内目标的比值的差值,|·|表示取绝对值操作,fk表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的比值,g表示接收信号内目标的比值。
取鉴角曲线上仰角测量区间上所有的导向矢量的比值与接收信号内目标的比值的差值中最小的差值,将该差值对应的导向矢量的仰角值作为在有多径反射时,米波雷达接收信号内目标的仰角值。
下面结合仿真图对本发明的效果做进一步的描述。
1.仿真条件:
本发明的仿真是在MATLAB R2010a的软件环境下进行的。#
2.仿真内容:
本发明的仿真实验所使用的天线为垂直放置的均匀线阵,阵元数为10,阵元间距为米1.6,天线中心高度为25米,波长为3.2米。多径反射系数为0.95。双波束指向中心为12°,波束1指向为9°,波束2指向为15°。凹口置零区间为-17°到-8°。鉴角曲线仰角测量区间为8.5°到16.5°。实际搜索目标的仰角区间为10°到14°,目标距离为100千米,目标信噪比为21dB,进行500次蒙特卡洛实验。
3.仿真效果分析:
图2为本发明的米波雷达多径反射信号示意图。图2中的hr表示米波雷达天线的架高,Rd表示目标与米波雷达的距离,θ1表示目标的仰角,θ2表示多径反射信号的反射角。
从图2中可以看出米波雷达既能接收到目标的信号,也能接收到多径反射信号。
图3(a)为本发明的仿真实验中左波束指向方向图。图3(a)是左波束指向方向图。图3(a)中的横坐标表示仰角范围,纵坐标表示幅值。图3(b)为本发明的仿真实验中右波束指向方向图。图3(b)是右波束指向方向图。图3(b)中的横坐标表示仰角范围,纵坐标表示幅值。
从图3(a)、图3(b)中可以看出左、右波束指向方向图均在各自波束指向处有最大的幅值,同时均在-17°到-8°仰角范围内形成凹口,达到抑制多径反射信号的目的。
图4为本发明的仿真实验得到的鉴角曲线,该鉴角曲线是利用本发明的仿真参数,通过计算仰角测量区间上的导向矢量中的任意一个导向矢量的左、右加权求和值,计算仰角测量区间上的导向矢量中的任意一个导向矢量的比值,连接仰角测量区间内上所有的导向矢量的比值,得到的鉴角曲线。图4中的横坐标表示仰角测量区间内上所有的导向矢量的仰角值,纵坐标表示仰角测量区间内上所有的导向矢量的比值。
从图4中可以看出在鉴角曲线上,纵坐标的每个比值只对应一个仰角值,满足单调性。
图5为本发明的仿真实验得到的仰角测量精度图。该仰角测量精度图是利用本发明的仿真参数,计算本发明方法下的目标仰角值和传统比幅单脉冲测角方法下的目标仰角值,计算各个方法下计算的目标仰角值与目标的真实仰角值的误差值,连接所有的误差值,绘制得到的仰角测量精度图。
图5中的横坐标表示目标的真实仰角,纵坐标表示仰角测量误差值。
从图5中可以看出,本发明的仰角测量误差值低于传统比幅单脉冲测角方法的仰角测量误差值,仰角的测量精度高。

Claims (6)

1.一种米波雷达有多径反射时对高仰角目标的仰角测量方法,包括如下步骤:
(1)确定仰角测量的左、右波束指向:
(1a)在有多径反射时,米波雷达发射脉冲中任意选定两个相互交叠的波束;
(1b)将左边的波束作为仰角测量的左波束,右边的波束作为仰角测量的右波束,将仰角测量的左波束的指向作为仰角测量的左波束指向,仰角测量的右波束的指向作为仰角测量的右波束指向;
(2)确定仰角测量的中心指向:
将仰角测量的左波束与右波束交叠位置的指向作为仰角测量的中心指向;
(3)确定凹口区间:
按照下式,确定仰角测量的左、右波束指向方向图的凹口区间:
其中,D表示仰角测量的左、右波束指向方向图的凹口区间,θL表示仰角测量的左波束指向,θR表示仰角测量的右波束指向,Δ表示仰角测量的左、右波束指向方向图的凹口区间大小的二分之一;
(4)计算特征值对应的特征矢量:
(4a)按照导向矢量公式,计算仰角测量的左、右波束指向方向图的凹口区间上的导向矢量中的任意一个导向矢量;
(4b)按照自相关矩阵公式,计算仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵;
(4c)按照特征矢量公式,计算仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵的特征值及该特征值对应的特征矢量;
(5)构成信号子空间:
将仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵的所有特征值从大到小排序,选取排序后的特征值中的前m个特征值所对应的特征矢量,m的取值取决于雷达的工作环境,将所选取特征矢量组成信号子空间;
(6)确定加权值:
按照下式,计算仰角测量的左、右波束指向方向图的加权值:
w1=[I-G*(GH*G)-1*GH]*w3
w2=[I-G*(GH*G)-1*GH]*w4
其中,w1、w2分别表示仰角测量的左、右波束指向方向图的加权值,I表示单位矩阵,G表示信号子空间,*表示矩阵相乘操作,H表示共轭转置操作,(·)-1表示矩阵求逆操作,w3、w4分别表示采用传统比幅单脉冲测角方法确定的仰角测量的左、右波束指向方向图的加权值;
(7)建立鉴角曲线:
(7a)在有多径反射时,米波雷达发射脉冲中任意选定两个相互交叠的波束,将该交叠位置的指向作为仰角测量的中心指向,取一段以仰角测量的中心指向为中心的区间作为鉴角曲线上仰角测量区间;
(7b)按照导向矢量公式,计算鉴角曲线上仰角测量区间上的导向矢量中的任意一个导向矢量;
(7c)按照加权求和公式,计算鉴角曲线上仰角测量区间上的导向矢量中的任意一个导向矢量的左、右加权求和值;
(7d)按照下式,计算鉴角曲线上仰角测量区间上的导向矢量中的任意一个导向矢量的比值:
其中,fk表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的比值,pk1、pk2分别表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的左、右加权求和值;
(7e)连接鉴角曲线上仰角测量区间内上所有的导向矢量的比值,得到鉴角曲线;
(8)确定接收信号内目标的仰角值:
(8a)按照加权求和公式,计算在有多径反射时,米波雷达接收信号的左、右加权求和值;
(8b)按照下式,计算接收信号内目标的比值:
其中,g表示接收信号内目标的比值,q1、q2分别表示在有多径反射时,米波雷达接收信号的左、右加权求和值;
(8c)按照下式,计算鉴角曲线上仰角测量区间上的导向矢量中的任意一个导向矢量的比值与接收信号内目标的比值的差值:
ck=|fk-g|
其中,ck表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的比值与接收信号内目标的比值的差值,|·|表示取绝对值操作,fk表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的比值,g表示接收信号内目标的比值;
(8d)取鉴角曲线上仰角测量区间上所有的导向矢量的比值与接收信号内目标的比值的差值中最小的差值,将该差值对应的导向矢量的仰角值作为在有多径反射时,米波雷达接收信号内目标的仰角值。
2.根据权利要求1所述的米波雷达有多径反射时对高仰角目标的仰角测量方法,其特征在于:步骤(4a)、步骤(7b)中所述的导向矢量公式如下:
其中,a(θi)表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量中的第i个导向矢量,θi表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量中的第i个导向矢量的仰角值,e表示指数符号,j表示虚数,π表示圆周率,d表示在有多径反射时,米波雷达天线的阵元间距,N表示在有多径反射时,米波雷达天线的阵元总数,sin表示正弦操作,λ表示在有多径反射时,米波雷达天线的发射信号的波长,T表示转置操作。
3.根据权利要求1所述的米波雷达有多径反射时对高仰角目标的仰角测量方法,其特征在于:步骤(4b)中所述的自相关矩阵公式如下:
R=a*aH
其中,R表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵,a表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量,*表示矩阵相乘操作,H表示共轭转置操作。
4.根据权利要求1所述的米波雷达有多径反射时对高仰角目标的仰角测量方法,其特征在于:步骤(4c)中所述的特征矢量公式如下:
R*Vi=λiVi
其中,R表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵,*表示矩阵相乘操作,Vi表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵的第i个特征值对应的特征矢量,λi表示仰角测量的左、右波束指向方向图的凹口区间上的导向矢量的自相关矩阵的第i个特征值。
5.根据权利要求1所述的米波雷达有多径反射时对高仰角目标的仰角测量方法,其特征在于:步骤(6)中所述左、右波束指向方向图的加权值是由下式得到的:
其中,w3、w4分别表示采用传统比幅单脉冲测角方法确定的仰角测量的左、右波束指向方向图的加权值,e表示以自然数为底的指数操作,j表示虚数符号,π表示圆周率,d表示在有多径反射时,米波雷达天线的阵元间距,N表示在有多径反射时,米波雷达天线的阵元总数,sin表示正弦操作,θL表示仰角测量的左波束指向,θR表示仰角测量的右波束指向,λ表示在有多径反射时,米波雷达天线的发射信号的波长,T表示转置操作。
6.根据权利要求1所述的米波雷达有多径反射时对高仰角目标的仰角测量方法,其特征在于:步骤(7c)、步骤(8a)中所述的加权求和公式如下:
pk1=|w1 H*b(θk)|
pk2=|w2 H*b(θk)|
其中,pk1、pk2分别表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的左、右加权求和值,|·|表示取绝对值操作,w1、w2分别表示仰角测量的左、右波束指向方向图的加权值,H表示共轭转置操作,*表示矩阵相乘操作,b(θk)表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量,θk表示鉴角曲线上仰角测量区间上的导向矢量中的第k个导向矢量的仰角值。
CN201610817425.0A 2016-09-12 2016-09-12 米波雷达有多径反射时对高仰角目标的仰角测量方法 Active CN106646346B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610817425.0A CN106646346B (zh) 2016-09-12 2016-09-12 米波雷达有多径反射时对高仰角目标的仰角测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610817425.0A CN106646346B (zh) 2016-09-12 2016-09-12 米波雷达有多径反射时对高仰角目标的仰角测量方法

Publications (2)

Publication Number Publication Date
CN106646346A CN106646346A (zh) 2017-05-10
CN106646346B true CN106646346B (zh) 2019-02-19

Family

ID=58852073

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610817425.0A Active CN106646346B (zh) 2016-09-12 2016-09-12 米波雷达有多径反射时对高仰角目标的仰角测量方法

Country Status (1)

Country Link
CN (1) CN106646346B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107991659B (zh) * 2017-11-17 2021-05-18 西安电子科技大学 基于字典学习的米波雷达低仰角目标测高方法
CN108535698B (zh) * 2018-04-04 2021-07-20 西安电子科技大学 基于波束空间的米波雷达低仰角估计方法
CN114664068B (zh) * 2022-03-31 2024-01-12 军事科学院军事医学研究院环境医学与作业医学研究所 一种体力劳动疲劳预警系统及方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103293526A (zh) * 2013-05-13 2013-09-11 西安电子科技大学 基于合成导向矢量的波束形成米波雷达测高方法
CN103728614A (zh) * 2014-01-15 2014-04-16 西安电子科技大学 基于机扫米波雷达的改进单脉冲测角方法
CN103744077A (zh) * 2014-01-15 2014-04-23 西安电子科技大学 机扫米波雷达在多目标情况下的测角方法
CN104796208A (zh) * 2015-04-03 2015-07-22 西安科技大学 正交化搜索的邻近强弱信号波达角估计方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103293526A (zh) * 2013-05-13 2013-09-11 西安电子科技大学 基于合成导向矢量的波束形成米波雷达测高方法
CN103728614A (zh) * 2014-01-15 2014-04-16 西安电子科技大学 基于机扫米波雷达的改进单脉冲测角方法
CN103744077A (zh) * 2014-01-15 2014-04-23 西安电子科技大学 机扫米波雷达在多目标情况下的测角方法
CN104796208A (zh) * 2015-04-03 2015-07-22 西安科技大学 正交化搜索的邻近强弱信号波达角估计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于压缩感知的米波雷达低空测角算法;王园园 等;《系统工程与电子技术》;20140430;第36卷(第4期);667-671
米波MIMO雷达低空目标波达方向估计新方法;刘源 等;《电子与信息学报》;20160331;第38卷(第3期);622-628
米波雷达测高多径模型研究;胡晓琴 等;《电波科学学报》;20080831;第23卷(第4期);651-657

Also Published As

Publication number Publication date
CN106646346A (zh) 2017-05-10

Similar Documents

Publication Publication Date Title
CN102156279B (zh) 基于mimo的双基地雷达地面动目标检测方法
CN107656255B (zh) 基于多径回波的超宽带雷达动目标二维定位方法
CN104237844B (zh) 基于相位补偿的分布式米波平面阵列雷达方位角测量方法
CN102288944A (zh) 基于地形匹配的数字阵列米波雷达超分辨测高方法
CN103235292B (zh) 平面相控阵调零保形校正的全维和差测角方法
CN103901410B (zh) 基于稀疏恢复的机载双基地mimo雷达杂波抑制方法
CN103364772B (zh) 基于实数域广义多重信号分类算法的目标低仰角估计方法
CN106646346B (zh) 米波雷达有多径反射时对高仰角目标的仰角测量方法
CN105182322B (zh) 基于反射信号相位差的被动式定位方法
CN106291540A (zh) 一种基于doa估计的多输入多输出探地雷达逆向投影目标成像方法
CN103293526B (zh) 基于合成导向矢量的波束形成米波雷达测高方法
CN103308909A (zh) 基于多参数交替搜索的米波雷达低仰角测高方法
CN103744077B (zh) 机扫米波雷达在多目标情况下的测角方法
CN103885054B (zh) 一种基于分布源反射模型的米波雷达低仰角测高方法
CN103869298B (zh) 一种分布式mimo天波超视距雷达海杂波仿真方法
CN106066473A (zh) 正交波形下mimo雷达目标的多波束联合测角和点迹融合方法
CN109188019A (zh) 基于多重信号分类算法的三维风速风向测量方法
CN103837867A (zh) 一种利用ais信息进行高频雷达天线通道校正的方法
CN103412286A (zh) 基于mimo雷达的发射极化优化doa估计方法
CN105738887A (zh) 基于多普勒通道划分的机载雷达杂波功率谱的优化方法
CN104502904A (zh) 一种鱼雷自导波束锐化方法
CN105182325A (zh) 基于秩1约束的米波mimo雷达低仰角目标测高方法
CN103245942A (zh) 一种基于mimo阵列的无失真扇扫成像方法
CN106338724B (zh) 机扫米波mimo三坐标雷达威力的获取方法
RU2711341C1 (ru) Способ двухмерного пеленгования

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