CN111352105A - 用于目标跟踪的测角方法 - Google Patents

用于目标跟踪的测角方法 Download PDF

Info

Publication number
CN111352105A
CN111352105A CN201910106239.XA CN201910106239A CN111352105A CN 111352105 A CN111352105 A CN 111352105A CN 201910106239 A CN201910106239 A CN 201910106239A CN 111352105 A CN111352105 A CN 111352105A
Authority
CN
China
Prior art keywords
target
range
sum
angle
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.)
Granted
Application number
CN201910106239.XA
Other languages
English (en)
Other versions
CN111352105B (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.)
Aerospace Information Research Institute 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 CN201910106239.XA priority Critical patent/CN111352105B/zh
Publication of CN111352105A publication Critical patent/CN111352105A/zh
Application granted granted Critical
Publication of CN111352105B publication Critical patent/CN111352105B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/66Radar-tracking systems; Analogous systems
    • G01S13/72Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
    • 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/28Details of pulse systems
    • G01S7/285Receivers
    • G01S7/292Extracting wanted echo-signals
    • G01S7/2923Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
    • G01S7/2927Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods by deriving and controlling a threshold value

Landscapes

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

Abstract

本公开提供一种用于目标跟踪的测角方法,包括:步骤S1:使用多个通道接收原始回波信号,并将原始回波信号合成和通道信号与差通道信号;步骤S2:将和差通道信号进行距离徙动校正,形成距离对齐的和差通道信号;步骤S3:利用该距离对齐的和差通道信号,构成潜在目标距离单元集合;步骤S4:对潜在目标距离单元集合进行角闪烁抑制,形成有效目标距离单元集合;步骤S5:利用该有效目标距离单元集合,使天线保持对目标的跟踪。本公开提供的用于目标跟踪的测角方法对和差通道信号进行距离徙动校正,消除了目标运动对目标检测的影响,增加了目标检测的效率。引入角闪烁抑制方法,减小了复杂目标的角闪烁引起的角度测量误差,提高目标的跟踪精度。

Description

用于目标跟踪的测角方法
技术领域
本公开涉及雷达目标跟踪与目标跟踪技术领域,尤其涉及一种用于目标跟踪的测角方法。
背景技术
单脉冲技术是一种经典的测角方法,基于单脉冲技术的单脉冲雷达在目标跟踪领域具有广泛的应用。理论上,单脉冲雷达仅使用一个脉冲内的回波就可以实现对目标的角度测量。常规的单脉冲雷达具有多个通道,包括和通道与差通道,对和通道与差通道进行比幅处理,得到目标的角度。
雷达发射信号的带宽增加后,能够提高距离向的分辨率,对目标反射回来的宽带信号进行处理,可以获取目标的高分辨距离向。从而实现对复杂目标不同部位的检测与识别,获取更多的目标信息。
然而,在实现本公开的过程中,本申请发明人发现,为了提高测角精度,经常将一段时间内采集到的脉冲进行联合处理,称为脉冲积累技术。使用宽带雷达对运动目标进行跟踪时,由于距离分辨率高,即使目标的运动速度很慢,在较短的脉冲积累时间内,目标运动的距离也会超过一个距离分辨单元的长度,即发生距离徙动。并且,对复杂目标进行跟踪时,散射点之间的相互干扰会产生角闪烁现象,引起测角误差,严重情况下会导致目标丢失。
公开内容
(一)要解决的技术问题
基于上述技术问题,本公开提供一种用于目标跟踪的测角方法,以缓解现有技术中目标跟踪方法容易发生发成距离徙动以及产生角闪烁现象,引起测角误差,甚至导致目标丢失的技术问题。
(二)技术方案
本公开提供一种用于目标跟踪的测角方法,包括:步骤S1:使用多个通道接收宽带信号经过目标散射后得到的原始回波信号,并将所述原始回波信号合成和通道信号与差通道信号;步骤S2:将所述和通道信号与所述差通道信号进行距离徙动校正,形成距离对齐的和差通道信号;步骤S3:利用该距离对齐的和差通道信号,确定可能存在散射点的距离单元,构成潜在目标距离单元集合;步骤S4:对所述潜在目标距离单元集合进行角闪烁抑制,形成有效目标距离单元集合;步骤S5:利用该有效目标距离单元集合,确定目标偏离天线的角度,控制天线伺服转动相应的角度,使天线保持对目标的跟踪。
在本公开的一些实施例中,所述步骤S3包括:步骤S31:利用距离对齐的和差通道信号合成高分辨距离像,并进行非相干积累,得到非相干积累的距离功率像序列;以及步骤S32:对所述非相干积累的距离功率像序列进行阈值检测,确定可能存在散射点的距离单元,构成潜在目标距离单元集合。
在本公开的一些实施例中,所述步骤S32包括:步骤S32a:对所述非相干积累的距离功率像序列进行排序,计算其中值,并将该中值乘以门限因子η作为阈值;以及步骤S32b:检测功率值大于阈值的距离单元,视为可能存在散射点的距离单元,记录所有可能存在目标的距离单元序号l,得到潜在目标距离单元集合L0,以及潜在目标距离单元对应的距离值R(L0)。
在本公开的一些实施例中,所述步骤S5包括:步骤S51:对所述有效目标距离单元集合进行和差比幅测角,确定每个散射点偏离天线的角度;步骤S52:根据散射点偏离天线的角度计算出散射点的坐标,确定目标的几何中心;以及步骤S53:根据目标的几何中心,确定目标的几何中心偏离天线的角度,控制天线伺服转动相应的角度,使天线波束保持对目标的实时照射。
在本公开的一些实施例中,其中,所述原始回波信号满足下式:
Figure BDA0001966237100000021
其中,k=1,2,3,4表示通道的序号;m=1,2,3,...,M表示脉冲编号;n=1,2,3,...,N表示每一个脉冲中的第n个采样点,每个脉冲总共采N个采样点,
Figure BDA0001966237100000031
表示第n个采样点的频率,B为信号带宽;i为目标上散射点的序号,散射强度为σi,在第m个脉冲时刻,距雷达的距离为Rim,Fki为第i个散射点回波信号被第k个通道接收时的方向图系数;fc为信号中心频率,c=3e8m/s为电磁波的传播速度。
在本公开的一些实施例中,所述步骤S1中:将多个通道的所述原始回波信号进行数字和差处理,合成和通道信号、方位差通道信号,以及俯仰差通道信号,并分别满足下式:
Ssum(m,n)=S1(m,n)+S2(m,n)+S3(m,n)+S4(m,n)
Sdif_az(m,n)=S1(m,n)+S4(m,n)-S2(m,n)-S3(m,n)
Sdif_el(m,n)=S1(m,n)+S2(m,n)-S3(m,n)-S4(m,n)
其中Ssum(m,n),Sdif_az(m,n),Sdif_el(m,n)分别表示合成的和通道信号、方位差通道信号以及俯仰差通道信号。
在本公开的一些实施例中,所述步骤S31中,将距离对齐的和差通道信号在快时间向进行处理,合成高分辨距离像,该高分辨距离像满足下式:
Figure BDA0001966237100000032
Figure BDA0001966237100000033
Figure BDA0001966237100000034
其中Ysum(m,l),Ydif_az(m,l),Ydif_el(m,l)分别表示利用和通道信号、方位差通道信号以及俯仰差通道信号合成的高分辨距离像,As(m,l,i),Aaz(m,l,i),Ael(m,l,i)分别表示各个距离像单元中第i个散射点的幅度。
在本公开的一些实施例中,所述步骤S31中,对所述和通道高分辨距离像在慢时间向进行非相干积累,得到非相干积累的距离功率像序列,该非相干积累的距离功率像序列满足下式:
Figure BDA0001966237100000041
在本公开的一些实施例中,所述步骤S4包括:步骤S41:根据所述潜在目标距离单元集合提取和通道信号,方位差通道信号,俯仰差通道信号的高分辨距离像上相应距离单元内的信号;步骤S42:计算所述潜在目标距离单元集合内每个距离单元对应的方位向和差比与俯仰向和差比,并统计每个距离单元上每一列和差比值的虚部最大值与实部方差;以及步骤S43:将所述虚部最大值与预先设置的阈值ηim进行比较,将所述实部方差值与预先设置的阈值ηvar进行比较,将所述潜在目标距离单元集合中虚部最大值大于阈值的距离单元,与实部方差大于阈值的距离单元从潜在目标距离单元集合中剔除,构成有效目标距离单元集合L0eff
在本公开的一些实施例中,所述步骤S51中,根据所述有效目标距离单元集合内的距离单元序号,提取每个距离单元对应的方位向和差比与俯仰向和差比,利用下式进行和差比幅测角,确定有效距离单元上散射点偏离天线的角度,并使每个有效距离单元得到一个角度序列:
Figure BDA0001966237100000042
Figure BDA0001966237100000043
其中θbaz,kmaz分别表示方位向天线差波束的3dB波束宽度与单脉冲斜率,θbel,kmel分别表示俯仰向天线差波束的3dB波束宽度与单脉冲斜率,real表示取实部运算,l是所述有效目标距离单元集合L0eff内的距离单元序号。
在本公开的一些实施例中,所述步骤S52包括:步骤S52a:对所述有效目标距离单元集合内每个距离单元的角度序列进行线性拟合,得到每个距离单元上散射点在第M个脉冲时刻偏离天线的角度;步骤S52b:根据所述有效目标距离单元集合内每个距离单元上散射点的距离与角度,计算每个有效散射点对应的直角坐标;以及步骤S52c:对所述有效散射点的直角坐标计算平均值,确定目标的几何中心直角坐标。
在本公开的一些实施例中,其中:所述步骤S52a中,对于每一个有效目标距离单元l0,根据所述偏离天线的角度序列θaz(m,l0),由线性拟合关系得出一个第M个脉冲时刻的偏离角度θaz(l0);所述步骤S52b中,利用下式计算每个有效散射点对应的直角坐标:
x(l0)=R(l0)cos(θa_az+Δθaz(l0))cos(θa_el+Δθel(l0))
y(l0)=R(l0)cos(θa_az+Δθaz(l0))sin(θa_el+Δθel(l0))
z(l0)=R(l0)sin(θa_az+Δθaz(l0))
其中,θa_az表示天线中心的方位角,θa_el表示天线中心的俯仰角;
所述步骤S52c中,利用下式确定目标的几何中心直角坐标:
Figure BDA0001966237100000051
Figure BDA0001966237100000052
Figure BDA0001966237100000053
在本公开的一些实施例中,所述步骤S53中,利用下式,确定目标的几何中心偏离天线的角度:
Figure BDA0001966237100000054
Figure BDA0001966237100000055
Δθc_az=θc_aza_az
Figure BDA0001966237100000056
(三)有益效果
从上述技术方案可以看出,本公开提供的用于目标跟踪的测角方法具有以下有益效果的其中之一或其中一部分:
(1)对回波信号实施距离徙动校正,再进行非相干积累检测目标,避免了目标运动对目标检测的影响,增加了检测到散射点的数目,提供更多目标信息;
(2)对存在角闪烁现象的散射点进行剔除,对角闪烁现象进行抑制,减小目标角度的测量误差,提升目标跟踪精度。
附图说明
图1为本公开实施例提供的用于目标跟踪的高精度测角方法的步骤示意图。
图2为本发明实施例中步骤S3的具体步骤示意图。
图3为本发明实施例中步骤S32的具体步骤示意图。
图4为本发明实施例中步骤S5的具体步骤示意图。
图5为本发明实施例中步骤S4的具体步骤示意图。
图6为本发明实施例中步骤S52的具体步骤示意图。
图7为本发明实施例中仿真使用的多散射点目标模型。
图8为本发明实施例的步骤S2中进行距离徙动校正后的和通道高分辨距离像。
图9为本发明实施例的步骤S2中未进行距离徙动校正的和通道高分辨距离像。
图10为本发明实施例的步骤S2中距离徙动校正后对高分辨距离像进行非相干积累得到的距离功率像。
图11为本发明实施例的步骤S2中未进行距离徙动校正对高分辨距离像进行非相干积累得到的距离功率像。
图12为本发明实施例的步骤S5中所述的方位角跟踪曲线。
图13为本发明实施例的步骤S5中所述的俯仰角跟踪曲线。
具体实施方式
本公开提供的用于目标跟踪的测角方法对和差通道信号进行距离徙动校正,消除了目标运动对目标检测的影响,增加了目标检测的效率,同时引入角闪烁抑制方法,减小了复杂目标的角闪烁引起的角度测量误差,提高目标的跟踪精度。
为使本公开的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本公开进一步详细说明。
本公开实施例提供一种用于目标跟踪的测角方法,如图1至图6所示所示,包括:
步骤S1:使用多个通道接收宽带信号经过目标散射后得到的原始回波信号,并将原始回波信号合成和通道信号与差通道信号;
步骤S2:将和通道信号与所述差通道信号进行距离徙动校正,形成距离对齐的和差通道信号;
步骤S3:利用该距离对齐的和差通道信号,确定可能存在散射点的距离单元,构成潜在目标距离单元集合;
步骤S4:对潜在目标距离单元集合进行角闪烁抑制,形成有效目标距离单元集合;
步骤S5:利用该有效目标距离单元集合,确定目标偏离天线的角度,控制天线伺服转动相应的角度,使天线保持对目标的跟踪。
在本公开的一些实施例中,如图2所示,步骤S3包括:步骤S31:利用距离对齐的和差通道信号合成高分辨距离像,并进行非相干积累,得到非相干积累的距离功率像序列;以及步骤S32:对非相干积累的距离功率像序列进行阈值检测,确定可能存在散射点的距离单元,构成潜在目标距离单元集合。
在本公开的一些实施例中,如图3所示,步骤S32包括:步骤S32a:对非相干积累的距离功率像序列进行排序,计算其中值,并将该中值乘以门限因子η作为阈值;以及步骤S32b:检测功率值大于阈值的距离单元,视为可能存在散射点的距离单元,记录所有可能存在目标的距离单元序号l,得到潜在目标距离单元集合L0,以及潜在目标距离单元对应的距离值R(L0)。
在本公开的一些实施例中,如图4所示,步骤S5包括:步骤S51:对有效目标距离单元集合进行和差比幅测角,确定每个散射点偏离天线的角度;步骤S52:根据散射点偏离天线的角度计算出散射点的坐标,确定目标的几何中心;以及步骤S53:根据目标的几何中心,确定目标的几何中心偏离天线的角度,控制天线伺服转动相应的角度,使天线波束保持对目标的实时照射。
在本公开的一些实施例中,其中,原始回波信号满足下式:
Figure BDA0001966237100000071
其中,k=1,2,3,4表示通道的序号;m=1,2,3,...,M表示脉冲编号;n=1,2,3,...,N表示每一个脉冲中的第n个采样点,每个脉冲总共采N个采样点,
Figure BDA0001966237100000081
表示第n个采样点的频率,B为信号带宽;i为目标上散射点的序号,散射强度为σi,在第m个脉冲时刻,距雷达的距离为Rim,Fki为第i个散射点回波信号被第k个通道接收时的方向图系数;fc为信号中心频率,c=3e8m/s为电磁波的传播速度。
在本公开的一些实施例中,步骤S1中:将多个通道的所述原始回波信号进行数字和差处理,合成和通道信号、方位差通道信号,以及俯仰差通道信号,并分别满足下式:
Ssum(m,n)=S1(m,n)+S2(m,n)+S3(m,n)+S4(m,n)
Sdif_az(m,n)=S1(m,n)+S4(m,n)-S2(m,n)-S3(m,n)
Sdif_el(m,n)=S1(m,n)+S2(m,n)-S3(m,n)-S4(m,n)
其中Ssum(m,n),Sdif_az(m,n),Sdif_el(m,n)分别表示合成的和通道信号、方位差通道信号以及俯仰差通道信号。
在本公开的一些实施例中,步骤S31中,将距离对齐的和差通道信号在快时间向进行处理,合成高分辨距离像,该高分辨距离像满足下式:
Figure BDA0001966237100000082
Figure BDA0001966237100000083
Figure BDA0001966237100000084
其中Ysum(m,l),Ydif_az(m,l),Ydif_el(m,l)分别表示利用和通道信号、方位差通道信号以及俯仰差通道信号合成的高分辨距离像,As(m,l,i),Aaz(m,l,i),Ael(m,l,i)分别表示各个距离像单元中第i个散射点的幅度。
在本公开的一些实施例中,步骤S31中,对和通道高分辨距离像在慢时间向进行非相干积累,得到非相干积累的距离功率像序列,该非相干积累的距离功率像序列满足下式:
Figure BDA0001966237100000091
在本公开的一些实施例中,如图5所示,步骤S4包括:步骤S41:根据潜在目标距离单元集合提取和通道信号,方位差通道信号,俯仰差通道信号的高分辨距离像上相应距离单元内的信号;步骤S42:计算潜在目标距离单元集合内每个距离单元对应的方位向和差比与俯仰向和差比,并统计每个距离单元上每一列和差比值的虚部最大值与实部方差;以及步骤S43:将虚部最大值与预先设置的阈值ηim进行比较,将实部方差值与预先设置的阈值ηvar进行比较,将潜在目标距离单元集合中虚部最大值大于阈值的距离单元,与实部方差大于阈值的距离单元从潜在目标距离单元集合中剔除,构成有效目标距离单元集合LOeff
在本公开的一些实施例中,步骤S51中,根据有效目标距离单元集合内的距离单元序号,提取每个距离单元对应的方位向和差比与俯仰向和差比,利用下式进行和差比幅测角,确定有效距离单元上散射点偏离天线的角度,并使每个有效距离单元得到一个角度序列:
Figure BDA0001966237100000092
Figure BDA0001966237100000093
其中θbaz,kmaz分别表示方位向天线差波束的3dB波束宽度与单脉冲斜率,θbel,kmel分别表示俯仰向天线差波束的3dB波束宽度与单脉冲斜率,real表示取实部运算,l是所述有效目标距离单元集合L0eff内的距离单元序号。
在本公开的一些实施例中,如图6所示,步骤S52包括:步骤S52a:对有效目标距离单元集合内每个距离单元的角度序列进行线性拟合,得到每个距离单元上散射点在第M个脉冲时刻偏离天线的角度;步骤S52b:根据有效目标距离单元集合内每个距离单元上散射点的距离与角度,计算每个有效散射点对应的直角坐标;以及步骤S52c:对有效散射点的直角坐标计算平均值,确定目标的几何中心直角坐标。
在本公开的一些实施例中,步骤S52a中,对于每一个有效目标距离单元l0,根据偏离天线的角度序列θaz(m,l0),由线性拟合关系得出一个第M个脉冲时刻的偏离角度θaz(l0);
在本公开的一些实施例中,步骤S52b中,利用下式计算每个有效散射点对应的直角坐标:
x(l0)=R(l0)coS(θa_az+Δθaz(l0))cos(θa_el+Δθel(l0))
y(l0)=R(l0)coS(θa_az+Δθaz(l0))sin(θa_el+Δθel(l0))
z(l0)=R(l0)sin(θa_az+Δθaz(l0))
其中,θa_az表示天线中心的方位角,θa_el表示天线中心的俯仰角。
在本公开的一些实施例中,步骤S52c中,利用下式确定目标的几何中心直角坐标:
Figure BDA0001966237100000101
Figure BDA0001966237100000102
Figure BDA0001966237100000103
在本公开的一些实施例中,步骤S53中,利用下式,确定目标的几何中心偏离天线的角度:
Figure BDA0001966237100000104
Figure BDA0001966237100000105
Δθc_az=θc_aza_az
Figure BDA0001966237100000106
依据以上描述,本领域技术人员应当对本公开实施例提供的用于目标跟踪的测角方法有了清楚的认识。
以下以一具体实施例验证本公开提供的用于目标跟踪的测角方法的有效性:
在本发明实施例中,仿真对一个具有多个散射点的复杂目标的运动进行角度测量与跟踪,所使用的目标模型为一个由11个散射点构成的飞机轮廓,如图7所示,模型的尺寸为0.15×0.12m。雷达位于坐标原点,目标的初始坐标为(-10,30,1),运动速度为vx=2m/s,vy=0.5m/s,vz=1m/s。
本实施例提供的用于目标跟踪的测角方法的具体步骤为:
步骤S1:使用多个通道接收宽带信号经过目标散射后得到的原始回波信号,并将原始回波信号合成和通道信号与差通道信号;
步骤S2:将和通道信号与所述差通道信号进行距离徙动校正,形成距离对齐的和差通道信号,实际应用中,实现距离徙动校正的方法并非一种,在本实施例中,这一校正是通过Keystone变换实现的,其具体实现是对和差通道信号的每一列进行插值,其中需要用到尺度变换满足下式:
Figure BDA0001966237100000111
Figure BDA0001966237100000112
步骤S3:利用该距离对齐的和差通道信号,确定可能存在散射点的距离单元,构成潜在目标距离单元集合,具体包括:
步骤S31:利用距离对齐的和差通道信号合成高分辨距离像(在本实施例中,高分辨距离像是通过对距离对齐的和差通道信号在快时间向进行逆傅里叶变换实现的,如图8至图9所示,图8至图9为由本实施例中一个脉冲积累周期内的信号经处理得到的高分辨距离像,为了直观地对比距离徙动校正的效果,在本实施例中同时给出了未进行距离徙动校正时的和通道高分辨距离像,其中图8为进行距离徙动校正后的和通道高分辨距离像,图9为未进行距离徙动校正的和通道高分辨距离像),对和通道高分辨距离像取模得到高分辨距离像的幅度,对幅度计算平方得到功率,将功率信号在慢时间向进行累加,得到非相干积累的距离功率像,得到非相干积累的距离功率像序列(图10至图11为本实施例中一个脉冲积累周期内得到的非相干积累距离功率像,为了直观地对比距离徙动校正的效果,在本实施例中同时给出了未进行距离徙动校正时的和通道高分辨距离像,图10为距离徙动校正后对高分辨距离像进行非相干积累得到的距离功率像,图11为未进行距离徙动校正对高分辨距离像进行非相干积累得到的距离功率像);以及
步骤S32:对非相干积累的距离功率像序列进行阈值检测,确定可能存在散射点的距离单元,构成潜在目标距离单元集合,具体包括:
步骤S32a:对非相干积累的距离功率像序列进行排序,计算其中值,并将该中值乘以门限因子η作为阈值,本实施例中η=20;以及
步骤S32b:检测功率值大于阈值的距离单元,视为可能存在散射点的距离单元,记录所有可能存在目标的距离单元序号l,得到潜在目标距离单元集合L0,以及潜在目标距离单元对应的距离值R(L0);
步骤S4:对潜在目标距离单元集合进行角闪烁抑制,形成有效目标距离单元集合,具体包括:
步骤S41:根据潜在目标距离单元集合提取和通道信号,方位差通道信号,俯仰差通道信号的高分辨距离像上相应距离单元内的信号;
步骤S42:计算潜在目标距离单元集合内每个距离单元对应的方位向和差比与俯仰向和差比,并统计每个距离单元上每一列和差比值的虚部最大值与实部方差;以及
步骤S43:将虚部最大值与预先设置的阈值ηim进行比较,将实部方差值与预先设置的阈值ηvar进行比较,将潜在目标距离单元集合中虚部最大值大于阈值的距离单元,与实部方差大于阈值的距离单元从潜在目标距离单元集合中剔除,构成有效目标距离单元集合L0eff,在本实施例中,
Figure BDA0001966237100000121
α=0.2,其中α的值可以根据目标尺寸的先验信息进行确定,ηvar=0.5;
步骤S5:利用该有效目标距离单元集合,确定目标偏离天线的角度,控制天线伺服转动相应的角度,使天线保持对目标的跟踪,具体包括:
步骤S51:对有效目标距离单元集合进行和差比幅测角,确定每个散射点偏离天线的角度,具体为:根据有效目标距离单元集合内的距离单元序号,提取每个距离单元对应的方位向和差比与俯仰向和差比,进行和差比幅测角,确定有效距离单元上散射点偏离天线的角度,并使每个有效距离单元得到一个角度序列,在本实施例中,θbaz=7.236mrad,kbaz=1.522,θbel=7.236mrad,θkel=1.522;
步骤S52:根据散射点偏离天线的角度计算出散射点的坐标,确定目标的几何中心,具体包括:
步骤S52a:对有效目标距离单元集合内每个距离单元的角度序列进行线性拟合,得到每个距离单元上散射点在第M个脉冲时刻偏离天线的角度;
步骤S52b:根据有效目标距离单元集合内每个距离单元上散射点的距离与角度,计算每个有效散射点对应的直角坐标;以及
步骤S52c:对有效散射点的直角坐标计算平均值,确定目标的几何中心直角坐标;以及
步骤S53:根据目标的几何中心,确定目标的几何中心偏离天线的角度,控制天线伺服转动相应的角度,使天线波束保持对目标的实时照射,在本实施例中,对目标的角度跟踪曲线如图12至图13所示,其中,图12为方位角跟踪曲线,图13为俯仰角跟踪曲线。
综上所述,本公开提供的用于目标跟踪的测角方法对和差通道信号进行距离徙动校正的同时引入角闪烁抑制方法,消除了目标运动对目标检测的影响,增加了目标检测的效率,并且减小了复杂目标的角闪烁引起的角度测量误差,提高目标的跟踪精度。
还需要说明的是,实施例中提到的方向用语,例如“上”、“下”、“前”、“后”、“左”、“右”等,仅是参考附图的方向,并非用来限制本公开的保护范围。贯穿附图,相同的元素由相同或相近的附图标记来表示。在可能导致对本公开的理解造成混淆时,将省略常规结构或构造。
并且图中各部件的形状和尺寸不反映真实大小和比例,而仅示意本公开实施例的内容。另外,在权利要求中,不应将位于括号之间的任何参考符号构造成对权利要求的限制。
类似地,应当理解,为了精简本公开并帮助理解各个公开方面中的一个或多个,在上面对本公开的示例性实施例的描述中,本公开的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而,并不应将该公开的方法解释成反映如下意图:即所要求保护的本公开要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说,如前面的权利要求书所反映的那样,公开方面在于少于前面公开的单个实施例的所有特征。因此,遵循具体实施方式的权利要求书由此明确地并入该具体实施方式,其中每个权利要求本身都作为本公开的单独实施例。
以上所述的具体实施例,对本公开的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本公开的具体实施例而已,并不用于限制本公开,凡在本公开的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本公开的保护范围之内。

Claims (13)

1.一种用于目标跟踪的测角方法,包括:
步骤S1:使用多个通道接收宽带信号经过目标散射后得到的原始回波信号,并将所述原始回波信号合成和通道信号与差通道信号;
步骤S2:将所述和通道信号与所述差通道信号进行距离徙动校正,形成距离对齐的和差通道信号;
步骤S3:利用该距离对齐的和差通道信号,确定可能存在散射点的距离单元,构成潜在目标距离单元集合;
步骤S4:对所述潜在目标距离单元集合进行角闪烁抑制,形成有效目标距离单元集合;
步骤S5:利用该有效目标距离单元集合,确定目标偏离天线的角度,控制天线伺服转动相应的角度,使天线保持对目标的跟踪。
2.根据权利要求1所述的用于目标跟踪的测角方法,所述步骤S3包括:
步骤S31:利用距离对齐的和差通道信号合成高分辨距离像,并进行非相干积累,得到非相干积累的距离功率像序列;以及
步骤S32:对所述非相干积累的距离功率像序列进行阈值检测,确定可能存在散射点的距离单元,构成潜在目标距离单元集合。
3.根据权利要求2所述的用于目标跟踪的测角方法,所述步骤S32包括:
步骤S32a:对所述非相干积累的距离功率像序列进行排序,计算其中值,并将该中值乘以门限因子η作为阈值;以及
步骤S32b:检测功率值大于阈值的距离单元,视为可能存在散射点的距离单元,记录所有可能存在目标的距离单元序号l,得到潜在目标距离单元集合L0,以及潜在目标距离单元对应的距离值R(L0)。
4.根据权利要求2所述的用于目标跟踪的测角方法,所述步骤S5包括:
步骤S51:对所述有效目标距离单元集合进行和差比幅测角,确定每个散射点偏离天线的角度;
步骤S52:根据散射点偏离天线的角度计算出散射点的坐标,确定目标的几何中心;以及
步骤S53:根据目标的几何中心,确定目标的几何中心偏离天线的角度,控制天线伺服转动相应的角度,使天线波束保持对目标的实时照射。
5.根据权利要求4所述的用于目标跟踪的测角方法,其中,所述原始回波信号满足下式:
Figure FDA0001966237090000021
其中,k=1,2,3,4表示通道的序号;m=1,2,3,...,M表示脉冲编号;n=1,2,3,...,N表示每一个脉冲中的第n个采样点,每个脉冲总共采N个采样点,
Figure FDA0001966237090000022
表示第n个采样点的频率,B为信号带宽;i为目标上散射点的序号,散射强度为σi,在第m个脉冲时刻,距雷达的距离为Rim,Fki为第i个散射点回波信号被第k个通道接收时的方向图系数;fc为信号中心频率,c=3e8m/s为电磁波的传播速度。
6.根据权利要求5所述的用于目标跟踪的测角方法,所述步骤S1中:将多个通道的所述原始回波信号进行数字和差处理,合成和通道信号、方位差通道信号,以及俯仰差通道信号,并分别满足下式:
Ssum(m,n)=S1(m,n)+S2(m,n)+S3(m,n)+S4(m,n)
Sdif_az(m,n)=S1(m,n)+S4(m,n)-S2(m,n)-S3(m,n)
Sdif_el(m,n)=S1(m,n)+S2(m,n)-S3(m,n)-S4(m,n)
其中Ssum(m,n),Sdif_az(m,n),Sdif_el(m,n)分别表示合成的和通道信号、方位差通道信号以及俯仰差通道信号。
7.根据权利要求5所述的用于目标跟踪的测角方法,所述步骤S31中,将距离对齐的和差通道信号在快时间向进行处理,合成高分辨距离像,该高分辨距离像满足下式:
Figure FDA0001966237090000031
Figure FDA0001966237090000032
Figure FDA0001966237090000033
其中Ysum(m,l),Ydif_az(m,l),Ydif_el(m,l)分别表示利用和通道信号、方位差通道信号以及俯仰差通道信号合成的高分辨距离像,As(m,l,i),Aaz(m,l,i),Ael(m,l,i)分别表示各个距离像单元中第i个散射点的幅度。
8.根据权利要求7所述的用于目标跟踪的测角方法,所述步骤S31中,对所述和通道高分辨距离像在慢时间向进行非相干积累,得到非相干积累的距离功率像序列,该非相干积累的距离功率像序列满足下式:
Figure FDA0001966237090000034
9.根据权利要求6所述的用于目标跟踪的测角方法,所述步骤S4包括:
步骤S41:根据所述潜在目标距离单元集合提取和通道信号,方位差通道信号,俯仰差通道信号的高分辨距离像上相应距离单元内的信号;
步骤S42:计算所述潜在目标距离单元集合内每个距离单元对应的方位向和差比与俯仰向和差比,并统计每个距离单元上每一列和差比值的虚部最大值与实部方差;以及
步骤S43:将所述虚部最大值与预先设置的阈值ηim进行比较,将所述实部方差值与预先设置的阈值ηvar进行比较,将所述潜在目标距离单元集合中虚部最大值大于阈值的距离单元,与实部方差大于阈值的距离单元从潜在目标距离单元集合中剔除,构成有效目标距离单元集合L0eff
10.根据权利要求9所述的用于目标跟踪的测角方法,所述步骤S51中,根据所述有效目标距离单元集合内的距离单元序号,提取每个距离单元对应的方位向和差比与俯仰向和差比,利用下式进行和差比幅测角,确定有效距离单元上散射点偏离天线的角度,并使每个有效距离单元得到一个角度序列:
Figure FDA0001966237090000041
Figure FDA0001966237090000042
其中θbaz,kmaz分别表示方位向天线差波束的3dB波束宽度与单脉冲斜率,θbel,kmel分别表示俯仰向天线差波束的3dB波束宽度与单脉冲斜率,real表示取实部运算,l是所述有效目标距离单元集合L0eff内的距离单元序号。
11.根据权利要求10所述的用于目标跟踪的测角方法,所述步骤S52包括:
步骤S52a:对所述有效目标距离单元集合内每个距离单元的角度序列进行线性拟合,得到每个距离单元上散射点在第M个脉冲时刻偏离天线的角度;
步骤S52b:根据所述有效目标距离单元集合内每个距离单元上散射点的距离与角度,计算每个有效散射点对应的直角坐标;以及
步骤S52c:对所述有效散射点的直角坐标计算平均值,确定目标的几何中心直角坐标。
12.根据权利要求11所述的用于目标跟踪的测角方法,其中:
所述步骤S52a中,对于每一个有效目标距离单元l0,根据所述偏离天线的角度序列θaz(m,l0),由线性拟合关系得出一个第M个脉冲时刻的偏离角度θaz(l0);
所述步骤S52b中,利用下式计算每个有效散射点对应的直角坐标:
x(l0)=R(l0)cos(θa_az+Δθaz(l0))cos(θa_el+Δθel(l0))
y(l0)=R(l0)cos(θa_az+Δθaz(l0))sin(θa_el+Δθel(l0))
z(l0)=R(l0)sin(θa_az+Δθaz(l0))
其中,θa_az表示天线中心的方位角,θa_el表示天线中心的俯仰角;
所述步骤S52c中,利用下式确定目标的几何中心直角坐标:
Figure FDA0001966237090000051
Figure FDA0001966237090000052
Figure FDA0001966237090000053
13.根据权利要求12所述的用于目标跟踪的测角方法,所述步骤S53中,利用下式,确定目标的几何中心偏离天线的角度:
Figure FDA0001966237090000054
Figure FDA0001966237090000055
Δθc_az=θc_aza_az
Figure FDA0001966237090000056
CN201910106239.XA 2019-02-01 2019-02-01 用于目标跟踪的测角方法 Active CN111352105B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910106239.XA CN111352105B (zh) 2019-02-01 2019-02-01 用于目标跟踪的测角方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910106239.XA CN111352105B (zh) 2019-02-01 2019-02-01 用于目标跟踪的测角方法

Publications (2)

Publication Number Publication Date
CN111352105A true CN111352105A (zh) 2020-06-30
CN111352105B CN111352105B (zh) 2021-03-12

Family

ID=71195523

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910106239.XA Active CN111352105B (zh) 2019-02-01 2019-02-01 用于目标跟踪的测角方法

Country Status (1)

Country Link
CN (1) CN111352105B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111948643A (zh) * 2020-08-13 2020-11-17 中国科学院空天信息创新研究院 基于多波束非线性鉴角的目标跟踪方法
CN112415504A (zh) * 2020-11-02 2021-02-26 中国科学院空天信息创新研究院 一种雷达目标跟踪方法及装置
CN112965028A (zh) * 2021-02-10 2021-06-15 西南电子技术研究所(中国电子科技集团公司第十研究所) 多波束相控阵差和比角度估计方法
CN114465638A (zh) * 2022-01-19 2022-05-10 上海交通大学 单通道单脉冲跟踪方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102749620A (zh) * 2012-06-13 2012-10-24 北京航空航天大学 一种弹载/机载雷达单脉冲前视成像处理方法
CN105607055A (zh) * 2015-12-24 2016-05-25 西安电子科技大学 一种基于天线方向图的机载雷达单脉冲前视成像方法
CN106526583A (zh) * 2016-10-21 2017-03-22 北京无线电测量研究所 一种基于天线方向图信息的地面运动目标定位方法
CN107515398A (zh) * 2017-08-23 2017-12-26 哈尔滨工业大学 基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法
JP2018205174A (ja) * 2017-06-06 2018-12-27 株式会社東芝 レーダ装置及びそのレーダ信号処理方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102749620A (zh) * 2012-06-13 2012-10-24 北京航空航天大学 一种弹载/机载雷达单脉冲前视成像处理方法
CN105607055A (zh) * 2015-12-24 2016-05-25 西安电子科技大学 一种基于天线方向图的机载雷达单脉冲前视成像方法
CN106526583A (zh) * 2016-10-21 2017-03-22 北京无线电测量研究所 一种基于天线方向图信息的地面运动目标定位方法
JP2018205174A (ja) * 2017-06-06 2018-12-27 株式会社東芝 レーダ装置及びそのレーダ信号処理方法
CN107515398A (zh) * 2017-08-23 2017-12-26 哈尔滨工业大学 基于宽带回波的逆合成孔径雷达成像目标连续跟踪方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HOSSEIN DARVISHI等: "Adaptive hybrid method for low-angle target tracking in multipath", 《IET RADAR,SONAR&NAVIGATION》 *
LIANGTAN WAN等: "Joint Range-Doppler-Angle Estimation for Intelligent Tracking of Moving Aerial Targets", 《IEEE INTERNET OF THINGS JOURNAL》 *
李保国: "基于目标角闪烁特性的末制导雷达跟踪与识别技术", 《中国博士学位论文全文数据库信息科技辑》 *
逯程等: "基于位置信息的高分辨雷达角闪烁抑制方法", 《海军航空工程学院学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111948643A (zh) * 2020-08-13 2020-11-17 中国科学院空天信息创新研究院 基于多波束非线性鉴角的目标跟踪方法
CN111948643B (zh) * 2020-08-13 2021-03-26 中国科学院空天信息创新研究院 基于多波束非线性鉴角的目标跟踪方法
CN112415504A (zh) * 2020-11-02 2021-02-26 中国科学院空天信息创新研究院 一种雷达目标跟踪方法及装置
CN112965028A (zh) * 2021-02-10 2021-06-15 西南电子技术研究所(中国电子科技集团公司第十研究所) 多波束相控阵差和比角度估计方法
CN112965028B (zh) * 2021-02-10 2023-08-29 西南电子技术研究所(中国电子科技集团公司第十研究所) 多波束相控阵差和比角度估计方法
CN114465638A (zh) * 2022-01-19 2022-05-10 上海交通大学 单通道单脉冲跟踪方法及系统

Also Published As

Publication number Publication date
CN111352105B (zh) 2021-03-12

Similar Documents

Publication Publication Date Title
CN111352105B (zh) 用于目标跟踪的测角方法
CN111352107B (zh) 基于多通道数字和差的单脉冲跟踪与成像方法
CN111537966B (zh) 一种适用于毫米波车载雷达领域的阵列天线误差校正方法
CN109471101B (zh) 一种弹载sar图像校正方法
CN101498788B (zh) 一种逆合成孔径雷达的目标转角估计和横向定标方法
CN105445701B (zh) Ddma‑mimo雷达目标的单脉冲角度估计方法
CN106093932A (zh) 一种扫描波束高分辨率的雷达散射计
CN109655819B (zh) 一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法
CN112433210A (zh) 一种双站前视探地雷达快速时域成像方法
CN113671485B (zh) 基于admm的米波面阵雷达二维doa估计方法
CN113156444B (zh) 一种基于运动补偿的多波束声呐高精度成像方法
CN112147593B (zh) 一种高速密集爆炸物破片目标的四维参数估计方法
CN112415504A (zh) 一种雷达目标跟踪方法及装置
CN112698267A (zh) 一种干扰源测试定位方法
CN112433208A (zh) 一种双站圆周探地雷达快速时域成像方法及系统
CN116430344A (zh) 一种基于改进相位特征的海面慢速小目标检测方法
CN107037434B (zh) 基于和差波束单脉冲测角的成像方法
CN111948643B (zh) 基于多波束非线性鉴角的目标跟踪方法
CN110618403B (zh) 一种基于双波束雷达的着陆飞行器参数测量方法
Gao et al. Static background removal in vehicular radar: Filtering in azimuth-elevation-doppler domain
CN115390070A (zh) 一种机动平台单通道短孔径sar的海面目标聚焦定位方法
CN114994676A (zh) 一种一站固定式双站低频超宽带sar运动目标成像方法
CN114325700A (zh) 一种星载多通道sar动目标成像方法
CN108983192B (zh) 基于gps辐射源的雷达运动目标参数估计方法
CN111965637A (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
TA01 Transfer of patent application right

Effective date of registration: 20200730

Address after: 100190, No. 19 West Fourth Ring Road, Beijing, Haidian District

Applicant after: Research Institute of aerospace information innovation, Chinese Academy of Sciences

Address before: 100190, No. 19 West Fourth Ring Road, Beijing, Haidian District

Applicant before: Institute of Electronics, Chinese Academy of Sciences

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant