CN108828585A - 一种空间平台目标无伺服跟踪测角方法 - Google Patents

一种空间平台目标无伺服跟踪测角方法 Download PDF

Info

Publication number
CN108828585A
CN108828585A CN201810571355.4A CN201810571355A CN108828585A CN 108828585 A CN108828585 A CN 108828585A CN 201810571355 A CN201810571355 A CN 201810571355A CN 108828585 A CN108828585 A CN 108828585A
Authority
CN
China
Prior art keywords
angle
target
value
vector
tracking
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
CN201810571355.4A
Other languages
English (en)
Other versions
CN108828585B (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201810571355.4A priority Critical patent/CN108828585B/zh
Publication of CN108828585A publication Critical patent/CN108828585A/zh
Application granted granted Critical
Publication of CN108828585B publication Critical patent/CN108828585B/zh
Expired - Fee Related 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/66Radar-tracking systems; Analogous systems
    • G01S13/68Radar-tracking systems; Analogous systems for angle tracking only

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)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明公开了一种空间平台目标无伺服跟踪测角方法,该方法首先利用角度搜索算法完成大视场角度初值估计,之后利用角度跟踪环进行目标角度跟踪,本发明能够减少目标角度跟踪计算中的计算量,提高精度。

Description

一种空间平台目标无伺服跟踪测角方法
技术领域
本发明属于微波雷达测量技术领域,具体涉及一种空间平台目标无伺服跟踪测角方法。
背景技术
对空间目标角度测量是天基空间目标探测与定位的关键。对一般的雷达,通常采用窄波束发射和接收信号,通过波束的机械扫描或电子扫描实现对探测空域的覆盖,因而雷达分时探测大视场不同位置目标的速率低。在星载条件下,雷达天线的机械扫描不仅是可靠性的薄弱环节,而且会给卫星的姿态稳定带来一定扰动。电扫相控阵天线实现复杂,高频段T/R组件效率低,体积、重量、功耗大,难于满足空间平台系统轻小型化与安装要求。
2011年韩月涛等人在《电子测量与仪器学报》第25卷第10期第842页至849页发表的“相位差矢量平均的干涉仪解模糊算法”一文中,提出了一种基于相位差矢量平均的逐次递推解模糊测角方法,该方法利用相位差矢量平均值采取逐次递推解模糊的方法进行测角。该方法降低了通道相位噪声的影响改善了测角精度,但该方法采用相位逐级递推解模糊处理,在低信噪比下解模糊错误会逐级传播。
2017年同钊等人在《空间电子技术》第6期第79页至88页发表的“基于相控阵天线的动态飞行器间高精度视线角测量方法”一文中,提出了一种基于和差波束的单脉冲比幅测角方法,该方法实现了对目标角度的动态测量。然而,该方法一方面需要增加额外的和差器以形成和差波束,多通道相控阵天线造成测角系统的体积、重量和功耗增加,不适于在对资源要求苛刻的空间平台应用;另一方面,该方法对接收信号的信噪比阈值要求高,限制了其应用范围。
发明内容
有鉴于此,本发明提供了一种空间平台目标无伺服跟踪测角方法,能够克服上述缺陷,减少目标角度跟踪计算中的计算量,提高精度。
实现本发明的技术方案如下:
一种空间平台目标无伺服跟踪测角方法,包括以下步骤:
步骤1、计算目标回波信号的平均相位矢量;
步骤2、构造本地参考信号的相位矢量;
步骤3、利用步骤1的平均相位矢量和步骤2的相位矢量构造导向矢量相关函数,在设定角度范围内搜索导向矢量相关函数的最大值,其对应的角度即为单次角度测量值;
步骤4、重复步骤3,基于多次角度测量值预测目标角度;
步骤5、再次计算目标回波信号的平均相位矢量;
步骤6、利用上一次计算得到的角度跟踪值的左偏角和右偏角分别构造本地参考信号的相位矢量;所述角度跟踪值的初始值为预测得到的目标角度;
步骤7、利用左偏角构造的相位矢量和步骤5得到的平均相位矢量进行矢量相关运算得到Pls(θ),利用右偏角构造的相位矢量和步骤5得到的平均相位矢量进行矢量相关运算得到Prs(θ),利用Pls(θ)和Prs(θ)计算得到角误差信号θpk
步骤8、对θpk进行环路滤波,滤波结果作为数字控制振荡器的更新控制字;数字控制振荡器的输出值为当前的角度跟踪值,数字控制振荡器的角度累加器的初始值为预测得到的目标角度。
进一步地,所述目标回波信号采用广角天线进行采集。
进一步地,步骤7中,利用Pls(θ)和Prs(θ)通过归一化角误差鉴别器计算得到角误差信号θpk
有益效果
1、本发明提供的一种空间平台目标无伺服跟踪测角方法,利用角度跟踪环进行目标角度跟踪,不但能够满足近距离和远距离目标的测角精度要求,并且计算量小,具有很好的远近目标适应能力。
2、本发明采用阵列广角天线进行目标回波信号的采集,无需在空间平台安装伺服系统控制天线机械扫描,避免了传统机械扫描天线跟踪测角对空间平台的扰动。
3、与相控阵雷达跟踪测角方法相比,本发明方法有效降低了体积、重量和功耗,更适于在对资源要求苛刻的空间平台应用。
4、本发明采用归一化角误差鉴别器对测角误差进行鉴别,消除了接收机热噪声功率起伏对鉴别器输出的影响,能够提供稳定的鉴别增益,提高了跟踪环路的稳定性。
附图说明
图1为本发明角跟踪环路结构框图。
图2为本发明近场目标回波模型示意图。
图3为本发明角跟踪环路一阶滤波器的数字实现。
图4为本地跟踪角NCO结构框图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种空间平台目标无伺服跟踪测角方法,该方法首先利用角度搜索算法完成大视场角度初值估计,之后利用角度跟踪环进行目标角度跟踪,下面参照附图和实施例,对本发明作进一步详细说明:
在本发明方法实施之前,根据目标可能出现的角度范围,设计选择广角(大视场)天线单元,并完成天线阵列的构型设计,天线阵列的构型无特定要求,该天线阵列可直接固定安装于空间平台,用于对空间目标回波或辐射无线信号的接收。
步骤1、计算远场或近场目标回波信号的平均相位矢量;
不失一般性考虑M维线阵,目标的入射角为θ,利用阵列接收信号载波相位差的N个测量值,构造N个相位矢量As(θ,n):
式(1)中,φm,n(θ)=πhmsinθ+Δφm,n表示线接收阵第m个测量通道第n个观测样本载波相位差测量值;hm=2dm/λ为第m个测量通道基线长度dm与载波信号半波长λ/2的比值;Δφm,n表示第m个测量通道第n个观测样本载波相位的测量误差。
计算N个相位矢量As(θ,n)的平均值
步骤2、构造本地参考信号的相位矢量;
对于远场目标:以宽步长在角度[-90°,90°]范围等间距划分2P个搜索格点,依据天线阵列流型,根据每一格点角度θi(i=1,2,…,2P)对应的参考信号载波相位差矢量[Φ1i),Φ2i),…,ΦMi)]T=[πh1sinθi,πh2sinθi,…,πhMsinθi]T(i=1,2,…,2P),构造出本地参考信号载波相位差的导向矢量A(θi):
对于近场目标:如图2所示为近场目标回波模型,以宽步长在角度范围[-90°,90°]等间距划分2P个搜索格点,依据天线阵列流型,根据每一格点角度θi(i=1,2,…,2P)计算目标S到各个通道的距离
式(4)中,为近场目标S到测角天线阵坐标原点O的距离估计值,天线阵坐标原点O为天线阵列的几何中心,为搜索角度θi时,目标S到天线第m通道的距离,λ为载波波长,hm为天线第m通道到坐标原点O距离半波长的倍数。
对应的基于距离辅助的参考信号载波相位差矢量为
式(5)中,为第m通道搜索角度θi时的相位差重构值。利用式(5)构造出本地参考信号载波相位差的导向矢量
步骤3、利用步骤1的平均相位矢量和步骤2的相位矢量构造导向矢量相关函数,在设定角度范围内搜索导向矢量相关函数的最大值,其对应的角度即为单次角度测量值;
定义并计算与A(θi)的相位矢量匹配相关函数C(θi):
遍历[-90°,90°]范围内2P个格点,计算C(θi)的值并搜索极大值点,极大值点对应的θi即为单次角度测量值
argmax(·)表示使括号内的表达式达到最大值时的变量的取值。
步骤4、重复步骤3,基于多次角度测量值预测目标角度;
重复步骤3累计2K次,得到2K个角度测量值并存储,为方便后文记为分别对前K个角度测量值和后K个角度测量值进行排序处理,分别得到排序后的测角值序列
sort(·)表示排序;
计算中L个点的平均值以及中L个点的平均值
计算差值的绝对值ε:
判断ε的大小,如果ε≤ξ(ξ为初始值预置门限),则角度测量初始值即为
如果ε>ξ,则认为中存在异常,继续重复步骤3~步骤4,直到满足步骤4的差值的绝对值ε≤ξ的条件。
步骤5、再次计算远场或近场目标回波信号的平均相位矢量;
在完成角度初始值的预测后,将预测值预置给角跟踪环路,作为跟踪环路本地角跟踪值的初始值,之后角跟踪环路通过本地角跟踪值,实现对目标入射角度的精确实时跟踪。附图1给出了角跟踪环路的基本组成,跟踪环路包括接收信号相位矢量累加器、角误差鉴别器、环路滤波器和本地角再生相位矢量重构单元、本地跟踪角NCO。
计算各个天线阵列接收通道之间的N′个时刻相位差的平均值;
式(13)中,φm,n′(θ)表示以天线第一通道为参考,第n′时刻第m通道的相位差测量值,m=1,…,M表示天线通道序号,n′=1,…,N′。
利用式(13)中的相位差积累结果构造接收相位的复指数矢量,得到接收信号的相位导向矢量平均值;
步骤6、利用上一次计算得到的角度跟踪值的左偏角和右偏角分别构造远场或近场本地参考信号的相位矢量;所述角度跟踪值的初始值为预测得到的目标角度;
角跟踪环提取当前的角度跟踪值作为入射角的估计值本地角相位导向矢量重构单元利用角度估计值作为角度中心值,根据天线阵的阵列流型,以第一通道为参考,得到本地相位差值。
式(15)中,为第m通道本地角估计值的无模糊相位重构值;为第m通道主值区间位于[-π,π)的相位重构值;为第m通道整周模糊值。hm=2dm/λ为基线长度与半波长的比值。
利用式(15)进行远场的本地相位导向矢量重构,得到
根据角跟踪结果,分别构造中心角左偏角右偏角对应的相位矢量。
式(17)中,为中心跟踪角相位矢量重构值;为左偏跟踪角相位矢量重构值;为右偏跟踪角相位矢量重构值。
对于近场的目标本地参考信号的相位矢量:以第一通道为参考,得到基于距离辅助的本地相位差值,利用式(4)进行本地相位复矢量重构,得到
根据角跟踪结果,分别构造中心角左偏角右偏角对应的相位矢量。
式(19)中,为中心跟踪角相位矢量重构值;为左偏跟踪角相位矢量重构值;为右偏跟踪角相位矢量重构值。
步骤7、利用左偏角构造的相位矢量和步骤5得到的平均相位矢量进行矢量相关运算得到Pls(θ),利用右偏角构造的相位矢量和步骤5得到的平均相位矢量进行矢量相关运算得到Prs(θ),利用Pls(θ)和Prs(θ)计算得到角误差信号θpk
构造合适的代价函数,分别计算出 与As′(θ)的代价函数Pps(θ)、Pls(θ)、Prs(θ)。常用的代价函数构造准则有:
相位矢量相关值最大准则
角误差鉴别器的输入为Pps(θ)、Pls(θ)、Prs(θ),输出为角误差信号θpk
步骤8、对θpk进行环路滤波,滤波结果作为数字控制振荡器的更新控制字;数字控制振荡器的输出值为当前的角度跟踪值,数字控制振荡器的角度累加器的初始值为预测得到的目标角度。
环路滤波器的作用是滤除角误差信号θpk中的高频谐波分量与相位噪声分量,滤波器的阶数和带宽决定了其对目标角度的动态响应与噪声性能。不失一般性,本发明采用一阶环路滤波器进行角度误差滤波,见附图3。
一阶数字角跟踪环滤波器在z域的离散传递函数F(z)为
式中,ωn为角跟踪环路滤波器的自然角频率ωn=2BL/(ξ+0.25/ξ),BL为环路带宽,ξ为环路阻尼系数。
根据式(22),得到对应角跟踪环路滤波器的时域表达式为
式(23)中,T为环路更新时间;K=KdKv为角跟踪环路总增益,Kd为角误差鉴别器增益,Kv=1为本地跟踪角再生NCO的增益;y(n)为环路滤波器的输出,该值作为本地跟踪角再生NCO的更新控制字。
本地跟踪角再生NCO模块由角控制字转化单元、角控制字寄存器和角度累加器组成,如附图4所示,工作过程如下:开始时,根据步骤4预测得到的角度初始值来预置角度累加器的初值;之后,将环路滤波器的输出量转化成更新控制字,角度累加器以更新控制字为步长作累加运算,得到角度跟踪值。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (3)

1.一种空间平台目标无伺服跟踪测角方法,其特征在于,包括以下步骤:
步骤1、计算目标回波信号的平均相位矢量;
步骤2、构造本地参考信号的相位矢量;
步骤3、利用步骤1的平均相位矢量和步骤2的相位矢量构造导向矢量相关函数,在设定角度范围内搜索导向矢量相关函数的最大值,其对应的角度即为单次角度测量值;
步骤4、重复步骤3,基于多次角度测量值预测目标角度;
步骤5、再次计算目标回波信号的平均相位矢量;
步骤6、利用上一次计算得到的角度跟踪值的左偏角和右偏角分别构造本地参考信号的相位矢量;所述角度跟踪值的初始值为预测得到的目标角度;
步骤7、利用左偏角构造的相位矢量和步骤5得到的平均相位矢量进行矢量相关运算得到Pls(θ),利用右偏角构造的相位矢量和步骤5得到的平均相位矢量进行矢量相关运算得到Prs(θ),利用Pls(θ)和Prs(θ)计算得到角误差信号θpk
步骤8、对θpk进行环路滤波,滤波结果作为数字控制振荡器的更新控制字;数字控制振荡器的输出值为当前的角度跟踪值,数字控制振荡器的角度累加器的初始值为预测得到的目标角度。
2.如权利要求1所述的一种空间平台目标无伺服跟踪测角方法,其特征在于,所述目标回波信号采用广角天线进行采集。
3.如权利要求1所述的一种空间平台目标无伺服跟踪测角方法,其特征在于,步骤7中,利用Pls(θ)和Prs(θ)通过归一化角误差鉴别器计算得到角误差信号θpk
CN201810571355.4A 2018-06-06 2018-06-06 一种空间平台目标无伺服跟踪测角方法 Expired - Fee Related CN108828585B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810571355.4A CN108828585B (zh) 2018-06-06 2018-06-06 一种空间平台目标无伺服跟踪测角方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810571355.4A CN108828585B (zh) 2018-06-06 2018-06-06 一种空间平台目标无伺服跟踪测角方法

Publications (2)

Publication Number Publication Date
CN108828585A true CN108828585A (zh) 2018-11-16
CN108828585B CN108828585B (zh) 2019-07-09

Family

ID=64144093

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810571355.4A Expired - Fee Related CN108828585B (zh) 2018-06-06 2018-06-06 一种空间平台目标无伺服跟踪测角方法

Country Status (1)

Country Link
CN (1) CN108828585B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109633625A (zh) * 2018-12-14 2019-04-16 北京遥测技术研究所 基于ofdm信号归一化一阶矩的单脉冲测角方法及跟踪应用

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102830386A (zh) * 2012-09-03 2012-12-19 西安建筑科技大学 一种强干扰背景下的任意阵型微弱信号源角度估计方法
CN103308909A (zh) * 2013-03-06 2013-09-18 西安电子科技大学 基于多参数交替搜索的米波雷达低仰角测高方法
CN103383452A (zh) * 2013-06-26 2013-11-06 西安电子科技大学 分布式阵列目标到达角估计方法
CN106019252A (zh) * 2016-05-18 2016-10-12 西安电子科技大学 一种基于Nested阵列的和差跟踪测角方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102830386A (zh) * 2012-09-03 2012-12-19 西安建筑科技大学 一种强干扰背景下的任意阵型微弱信号源角度估计方法
CN103308909A (zh) * 2013-03-06 2013-09-18 西安电子科技大学 基于多参数交替搜索的米波雷达低仰角测高方法
CN103383452A (zh) * 2013-06-26 2013-11-06 西安电子科技大学 分布式阵列目标到达角估计方法
CN106019252A (zh) * 2016-05-18 2016-10-12 西安电子科技大学 一种基于Nested阵列的和差跟踪测角方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHEN JIANJUN ET AL.: "Application of Particle Filter for Angle Measure", 《2010 IEEE INTERNATIONAL CONFERENCE ON INTELLIGENT COMPUTING AND INTELLIGENT SYSTEMS》 *
任晓松 等: "基于矢量最优估计的稳健测向方法", 《现代防御技术》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109633625A (zh) * 2018-12-14 2019-04-16 北京遥测技术研究所 基于ofdm信号归一化一阶矩的单脉冲测角方法及跟踪应用

Also Published As

Publication number Publication date
CN108828585B (zh) 2019-07-09

Similar Documents

Publication Publication Date Title
CN108871348B (zh) 一种利用天基可见光相机的低轨卫星自主定轨方法
CN106093892B (zh) 基于标校卫星同时开展雷达rcs标定与外测标定系统
CN102590839B (zh) 一种三天线定位定向的方法
CN101937077B (zh) 双基地前视、斜视合成孔径雷达多普勒中心频率的测定方法
CN101762812B (zh) 一种机载单站无源定位系统对目标辐射波长的测定方法
CN108181630B (zh) 一种北斗双天线旋转快速定向方法
CN103363988A (zh) 一种利用智能手机传感器实现地磁室内定位导航的方法
CN103792535B (zh) 一种利用sar卫星测量电离层tec值的方法
Hu et al. Modification of slant range model and imaging processing in GEO SAR
CN101666880B (zh) 一种构建双站线阵三维成像合成孔径雷达系统方法
CN108267736A (zh) 一种geo sar凝视成像模式方位模糊度确定方法
CN108828585A (zh) 一种空间平台目标无伺服跟踪测角方法
CN212845925U (zh) 一种船载sar系统
CN107248891B (zh) 一种用于移动通信天线指向监测的测向测姿装置
CN108594217A (zh) 一种空间目标俯仰和方位二维角度闭环跟踪系统
CN105137406B (zh) 天基雷达效能优化方法
Ding et al. Potential applications of the moon based synthetic aperture radar for earth observation
Tianyuan et al. Single observer passive location using phase rate of change with the extended Kalman particle filter
CN111025358A (zh) 基于导航卫星信号单短基线的定向方法
CN102262236B (zh) 一种运动载体方位标定装置及其标定方法
CN115128651B (zh) 一种基于卫星平台传感器系统的辐射源定位方法
Shang et al. A phase calibration method for single channel monopulse tracking systems based on curve fitting technique
Laine 2D model-based step-track procedure
Han et al. Design of Ship Tracking Space Situational Awareness Guidance System
Lanyi Total ionospheric electron content calibration using SERIES GPS satellite data

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Cui Wei

Inventor after: Li Yechen

Inventor after: Yao Jie

Inventor after: Han Yuetao

Inventor after: Wu Siliang

Inventor after: Liang Yibao

Inventor after: Shen Qing

Inventor before: Cui Wei

Inventor before: Yao Jie

Inventor before: Han Yuetao

Inventor before: Wu Siliang

Inventor before: Liang Yibao

Inventor before: Shen Qing

GR01 Patent grant
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: 20190709