CN104777478B - 一种相控阵雷达搜索捕获目标的方法 - Google Patents
一种相控阵雷达搜索捕获目标的方法 Download PDFInfo
- Publication number
- CN104777478B CN104777478B CN201510179314.7A CN201510179314A CN104777478B CN 104777478 B CN104777478 B CN 104777478B CN 201510179314 A CN201510179314 A CN 201510179314A CN 104777478 B CN104777478 B CN 104777478B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mover
- mtd
- angle
- 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
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
- G01S13/00—Systems 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/66—Radar-tracking systems; Analogous systems
- G01S13/68—Radar-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)
Abstract
本发明公开一种相控阵雷达搜索捕获目标的方法,属于相控阵跟踪雷达信号处理领域。本发明包括搜索过程和捕获过程,其中搜索过程通过波束引导信息按照两种波位编排状态对雷达信号进行搜索,如果接收信号的能量大于门限,判定符合转入捕获过程;捕获过程与普通相控阵雷达捕获方法相比,不仅进行了多次检测,以降低虚警概率,还增加了类单脉冲和差波束测角和α‑β滤波闭环的测角和跟踪滤波过程,利用跟踪滤波处理过程中新息过程中角度余弦残差为参量做判决,判断是否捕获成功,使得捕获过程有一个收敛的过程,在初始波束指向目标位置有较大误差时候,捕获到的雷达信号更稳定可靠,可使得整个相控阵自动跟踪***更加稳健。
Description
技术领域
本发明涉及相控阵跟踪雷达信号处理领域,特别涉及相控阵卫星跟踪***搜索/捕获的方法,以此过渡到稳定跟踪阶段。
背景技术
近年来,移动通信的不断发展,通信***对通信距离和范围的要求越来越高。移动卫星通信***成为一种机动通信的良好手段,可用于汽车,火车,舰船,飞机,导弹等各种移动载体上,因其覆盖范围广、对地域要求不敏感、通信容量大等优点,已经被广泛应用于军事和民用的多个领域。移动卫星通信***的关键设备是天线稳定跟踪***。
有源相控阵天线具有波束快速扫描的能力,采用电子移相方式控制波束方向,能够克服机械扫描天线波束指向转换的惯性及由此对***性能的限制。避免传统的机械雷达天线的机械伺服机构的运动速度较慢,跟不上载体的姿态和地理位置变化,使得在动态情况下天线的指向偏离卫星,导致通信质量下降或者造成通信中断的情况。
跟踪雷达是指那些能够自动连续的对目标进行跟踪,不断地对目标参数进行精确测量并能输出目标坐标位置参数的雷达,连续闭环式自动跟踪、目标坐标参数的高精度测量及高数据输出率是目标跟踪的主要特点。单脉冲雷达跟踪精度高,实现简单,数以典型的高精度跟踪雷达都选用单脉冲体制。
卫星通信相控阵天线自跟踪***,主要包括接收天线阵列、发射天线阵列、射频前端模块、自跟踪基带处理单元、波束引导控制单元、外部时钟和电源等部分。自跟踪接收***通过将天线阵列划分成4个子阵。各天线单元的信号合成四个子阵,进行4路A/D变换,4路子阵信号经模拟下变频变为中频信号并经ADC采样后,首先由预处理模块完成复信号产生、通道校正等处理工作。然后进入搜索捕获状态,根据载体平台提供的引导信息完成卫星信号的搜索捕获。捕获成功后***进入自跟踪状态,通过角度测量和跟踪滤波等处理使得天线波束自动对准卫星,实现接收阵列波束对卫星的自动跟踪。跟踪滤波过程中将实时检测目标是否跟丢,判定为跟丢后则重新转入搜索捕获状态。
由于运动平台设备获得的卫星位置预报信息和实际目标难免存在一定的角度误差,需要雷达在一定范围内进行搜索捕获过程,减小与目标位置之间的角度误差满足天线的跟踪环路跟踪状态条件。对于搜索捕获过程,一般的搜索捕获就是通过门限检测和判决,一种典型的技术是设置合适的门限,使得门限电压足够高,以防止大多数噪声尖峰超过门限,可是又要低的足以让弱信号通过,只是有的做了两次或多次判决确认,而本发明借助了闭环式跟踪滤波方法,运用了角度测量以及跟踪滤波,使得捕获过程有一个收敛的过程,更稳定,整个算法简易稳定,具有较高的搜索速度,可以运用到高速的移动平台的自跟踪***,保证自跟踪的准确性。
跟踪滤波技术是跟踪***的核心。国内外对该项技术也开展了深入研究。卡尔曼(Kalman)用状态方程来描述随机线性***,提出了卡尔曼滤波算法。该算法同时适用于平稳和非平稳随机过程。卡尔曼滤波在线性***中是最优的估计方法,但卡尔曼滤波需要实时求解增益矩阵,涉及到较复杂的矩阵求逆过程,为算法的工程实现增加了难度,针对非线性***,基于将非线性方程进行线性近似的思想,提出了扩展卡尔曼滤波等非线性滤波算法。但这类算法容易发散,工程应用受到限制。不敏卡尔曼滤波和粒子滤波两种较新的非线性滤波算法,但这些方法目前尚处于理论研究阶段,工程上难以实现。工程上通常采用α-β滤波算法,该算法是一种常增益的滤波方法,该方法避免了卡尔曼滤波中复杂的增益矩阵求解问题,易于实现,因而在工程上获得了较多应用。
发明内容
本发明提供了一种基于类单脉冲和差比幅测角的相控阵雷达搜索捕获方法,该方法与传统的搜索捕获方法相比,一般的搜索捕获通过门限检测和判决,只是有的做了两次或多次判决确认,而本方法借助了角度跟踪滤波,在搜索过程通过门限检测到信号后,进行类单脉冲和差比幅测角和α-β滤波的闭环跟踪过程使得捕获过程有一个收敛的过程,使得捕获更稳定可靠。
本发明具体采用如下技术方案:
一种相控阵雷达搜索捕获目标的方法,其流程如图1所示,具体包括以下步骤:
步骤1.相控阵雷达子阵划分;
将相控阵雷达天线阵列划分成四个子阵,分别记为子阵1、子阵2、子阵3及子阵4;
步骤2.设定能量检测门限Et;
记n1(n)、n2(n)、n3(n)、n4(n)分别为所述四个子阵均未接收到捕获目标的信号时相应的输出信号,则此时合成信号记为nΣ=n1(n)+n2(n)+n3(n)+n4(n);由此得噪声总能量:
求解所述噪声总能量的均值以及方差σΣ,得能量检测门限其中k1取值1.0到1.2、k2取值4到5;
步骤3.目标搜索过程;从载体平台提供的引导信息获得目标的初始位置后,通过控制相控阵雷达的波束位置指向该初始位置,并以该初始位置为中心进行波束搜索;波束搜索具体采用以下方式进行:
步骤3-1.构建两种不同范围的用于波束搜索的波位排列,分别记为排列1、排列2;
所述排列1如图4所示,包括5个波位:中心波位1,分别位于中心波位1的左上方、右上方、左下方及右下方的波位2、3、4及波位5;所述中心波位1与波位2、3、4、5相交叠,所述排列1中波位2、3、4、5的中心均与波位1的中心相距半个波束宽度;
所述排列2如图5所示,包括呈九宫格排列的九个波位:波位1位于九宫格中心,所述九宫格第一排由左至右分别为波位2、3及波位4,波位1的左、右分别为波位5、波位6,九宫格第三排由左至右分别为波位7、8、9;排列2中左右或上下相邻的两个波位中心之间均相距一个波束宽度;
步骤3-2.以目标的初始位置为中心波位1构建排列1所示的波束搜索排列方式,将相控阵雷达的波束指向调整至目标的初始位置,计算此时四个子阵合成信号的能量;
所述合成信号的能量具体由以下方法获得:
以子阵1为参考子阵,分别对每个子阵做波束形成,得到相应的输出信号y1(n)、y2(n)、y3(n)及y4(n),将所述四个输出信号合成得YΣ(n)=y1(n)+y2(n)+y3(n)+y4(n),由此求解出四个子阵合成信号的能量:
N表示时域上的总采样点数;
若所述能量不小于所述能量检测门限Et,则目标搜索完成,记此时相控阵雷达的波束指向为捕获用波束指向,执行步骤4;否则,将相控阵雷达的波束指向依次调整至波位2、3、4、5直至四个子阵合成信号的能量不小于所述能量检测门限Et时止;若将相控阵雷达的波束指向调整至波位5时,四个子阵合成信号的能量依旧小于所述能量检测门限Et,执行步骤3-3;
步骤3-3.以目标的初始位置为中心波位1构建排列2所示的波束搜索排列方式,将相控阵雷达的波束指向调整至目标的初始位置,计算此时四个子阵合成信号的能量;
若所述能量不小于所述能量检测门限Et,则目标搜索完成,记此时相控阵雷达的波束指向为捕获用波束指向,执行步骤4;否则,将相控阵雷达的波束指向依次调整至波位2、3、4、5、6、7、8、9直至四个子阵合成信号的能量不小于所述能量检测门限Et时止;若将相控阵雷达的波束指向调整至波位9时,四个子阵合成信号的能量依旧小于所述能量检测门限Et,则载体平台重新初始化获得最新的引导信息并向相控阵雷达提供目标的最新初始位置,并再次执行步骤3-1至步骤3-3;
步骤4.目标捕获过程;
步骤4-1.以步骤3所得的捕获用波束指向作为类单脉冲测角的初始值,利用和差波束比幅测角方法求解处目标信号的方位角和俯仰角
所述和差比幅测角方法采用的波束排列如图6所示,包括5个波位:中心波位5,分别位于中心波位5的正上方、正上方、左方及右方的波位1、2、3及波位4;所述中心波位5与波位1、2、3、4相交叠,所述排列中波位1、2、3、4的中心均与波位5的中心相距半个波束宽度;
步骤4-2.根据单脉冲和差波束比幅测角得到测量角度求得对应角度下的方向余弦:
以目标在三个方向的方向余弦及其变化的速度为状态向量:
其中为方向余弦变化的速度,应用α-β算法跟踪滤波得到入射角度估计值作为下一次测角的初始角度;
步骤4-3.在滤波过程中,受多种因素影响,观测数据中可能有部分数据与真实值之间存在很大误差,这些数据被称为野值,野值对跟踪性能会产生较大的影响,需要剔除;在进行α-β滤波时,将新息过程模值与某一门限W做判决,若满足:
观测数据有效,利用α-β滤波进行状态估计,并利用反三角函数求解出滤波过后的方位角和俯仰角作为下一次类单脉冲测角的初始值;否则,将当前观测值被判为野值,执行步骤4-4;
其中,门限值可选为W选为三个方向观测噪声标准差各自的3倍所组成的向量;
步骤4-4.若当前观测值被判为野值,则在状态估计时利用前一周期的状态向量估计值做线性外推,当连续3个观测时刻均检测到野值时判定为跟踪目标丢失,此时应当重新初始化α-β跟踪滤波器重新执行步骤3至步骤4所述的整个搜索和捕获过程;
步骤4-5.滤波过程中将角度余弦残差|αt|和|βt|与门限做判决,门限值通常取1/6至1/3波束宽度,当|αt|和|βt|均小于所述门限值时则认为捕获过程收敛,目标捕获完成可以转入跟踪过程,否则,将重复执行步骤4-1至步骤4-4直至捕获过程收敛,完成目标捕获。
本发明的有益效果是:传统方法在搜索过程通过门限检测到信号后进行两次或者多次检测的过程,尽管可以提高检测概率但依然无法保证每次进入自跟踪过程的波束指向满足自跟踪条件要求,如果波束指向信息有误会使得自跟踪结果误差大或者不能使用,本方法在搜索过程通过门限检测到信号后,捕获过程中进行了闭环跟踪过程,并且此时跟踪滤波器的采样率通常是正常跟踪状态的2~5倍可以使得捕获过程快速收敛,进入自跟踪过程波束指向信息满足其跟踪条件,保证跟踪结果准确,提高自跟踪精度。
附图说明
图1为本发明提供的相控阵雷达目标搜索捕获方法的流程图;
图2为相控阵雷达天线的坐标分布;
图3为相邻两波束示意图;
图4为搜索波位的排列分布1;
图5为搜索波位的排列分布2;
图6为和差比幅测角波束排列。
具体实施方式
本发明的搜索捕获方法包括搜索过程和捕获过程,***首先根据平台的波束引导信息进行目标位置搜索,搜索到有效信号后进入捕获过程。其中,搜索过程通过对两种波位编排状态进行波束扫描,对每个波位进行能量检测,以获得有效的检测信号;捕获过程,在搜索状态检测到目标后,阵列波束对检测到信号的方向进行重复接收和检测,排除虚警情况后,进行闭环的类单脉冲和差比幅测角以及α-β跟踪滤波,使得捕获过程快速收敛,进入自跟踪状态。
本发明具体采用如下技术方案:
一种相控阵雷达搜索捕获目标的方法,其流程如图1所示,具体包括以下步骤:
步骤1.相控阵雷达划分;
将相控阵雷达划分成四个子阵,分别记为子阵1、子阵2、子阵3及子阵4;
步骤2.设定能量检测门限Et;
记n1(n)、n2(n)、n3(n)、n4(n)分别为所述四个子阵均未接收到捕获目标的信号时相应的输出信号,则此时合成信号记为nΣ=n1(n)+n2(n)+n3(n)+n4(n);由此得噪声总能量:
求解所述噪声总能量的均值以及方差σΣ,得能量检测门限其中k1∈[1.0,1.2]、k2∈[4,5];
步骤3.目标搜索过程;从载体平台提供的引导信息获得目标的初始位置后,控制相控阵雷达的波束位置指向该初始位置,并以该初始位置为中心进行波束搜索;波束搜索具体采用以下方式进行:
步骤3-1.构建两种不同范围的用于波束搜索的波位排列,分别记为排列1、排列2;
所述排列1如图4所示,包括5个波位:中心波位1,分别位于中心波位1的左上方、右上方、左下方及右下方的波位2、3、4及波位5;所述中心波位1与波位2、3、4、5相交叠,所述排列1中波位2、3、4、5的中心均与波位1的中心相距半个波束宽度;
所述排列2如图5所示,包括呈九宫格排列的九个波位:波位1位于九宫格中心,所述九宫格第一排由左至右分别为波位2、3及波位4,波位1的左、右分别为波位5、波位6,九宫格第三排由左至右分别为波位7、8、9;排列2中左右或上下相邻的两个波位中心之间均相距一个波束宽度;
步骤3-2.以目标的初始位置为中心波位1构建排列1所示的波束搜索排列方式,将相控阵雷达的波束指向调整至目标的初始位置,计算此时四个子阵合成信号的能量;
所述合成信号的能量具体由以下方法获得:
以子阵1为参考子阵,分别对每个子阵做波束形成,得到相应的输出信号y1(n)、y2(n)、y3(n)及y4(n),将所述四个输出信号合成得YΣ(n)=y1(n)+y2(n)+y3(n)+y4(n),由此求解出四个子阵合成信号的能量:
其中N为时域上的总采样点数; (2)
若所述能量FΣ不小于所述能量检测门限Et,则目标搜索完成,记此时相控阵雷达的波束指向为捕获用波束指向,执行步骤4;否则,将相控阵雷达的波束指向依次调整至波位2、3、4、5直至四个子阵合成信号的能量FΣ不小于所述能量检测门限Et时止;若将相控阵雷达的波束指向调整至波位5时,四个子阵合成信号的能量FΣ依旧小于所述能量检测门限Et,执行步骤3-3;
步骤3-3.以目标的初始位置为中心波位1构建排列2所示的波束搜索排列方式,将相控阵雷达的波束指向调整至目标的初始位置,计算此时四个子阵合成信号的能量;
若所述能量FΣ不小于所述能量检测门限Et,则目标搜索完成,记此时相控阵雷达的波束指向为捕获用波束指向,执行步骤4;否则,将相控阵雷达的波束指向依次调整至波位2、3、4、5、6、7、8、9直至四个子阵合成信号的能量FΣ不小于所述能量检测门限Et时止;若将相控阵雷达的波束指向调整至波位9时,四个子阵合成信号的能量FΣ依旧小于所述能量检测门限Et,则载体平台重新初始化获得最新的引导信息并向相控阵雷达提供目标的最新初始位置,并再次执行步骤3-1至步骤3-3;
步骤4.目标捕获过程;
步骤4-1.以步骤3所得的捕获用波束指向作为类单脉冲测角的初始值,利用和差波束比幅测角方法求解处目标信号的方位角和俯仰角具体求解方法如下:
令雷达相控阵所在平面为x-o-y平面,捕获过程中相控阵雷达的入射信号初始指向与x轴及y轴的夹角分别为αx和αy;入射信号与z轴夹角为俯仰角θ,其取值范围为0°~90°,入射信号在x-o-y平面的投影与x轴夹角为方位角其取值范围为0°~360°,(αx,αy)与有如下关系:
4个天线子阵经移相后输出的信号分别记为y1(n)、y2(n)、y3(n)和y4(n),用第一通道的接收信号y1(n)分别与四个信号分别进行相关处理,就可得到子天线阵1与其余子阵之间的相关系数:R1,1、R1,2、R1,3和R1,4;这四个相关系数中包含有入射源的角度信息,具体的:
其中N时域上的总采样点数,将四个相关系数看成四个子阵输出的复幅度,可以得到和信号:
FΣi=R1,1+R1,2+R1,3+R1,4 (5)
以波束5为中心,角度余弦坐标系下波束5的指向为(α5,β5)=(α0,β0);沿α坐标方向左右各偏移该方向的半个波束宽度,顺序形成波束1和波束2,两波束指向分别为:
同样地,沿β坐标方向依次形成波束3和波束4,波束指向分别为
以上两式中,α3dB与β3dB分别表示角度余弦坐标系下α坐标方向和β坐标方向的半功率波束宽度。与入射信号方向越靠近的波束接收到的信号更强,据此可以判断目标信号偏离等信号轴的方向。
对5个波束方位进行接收信号合成,可以得到5个方位的和波束接收信号能量FΣ1~FΣ5。其中FΣi其中按照式(28)求得,下式求取α和β方向的误差电压:
分别记为目标方向偏离中心波束指向的角度余弦值,其中和为目标信号方向角度余弦,在一定范围内角误差信号uα与αt,uβ与βt近似成线性关系,即
其中斜率kα和kβ可以通过实验拟合得到,实际应用中,求解kα时假设β方向指向为目标方向,在α0附近左右偏离一个波束宽度范围里对应uα与αt的值并用直线逼近uα-αt的关系曲线,由实测数据拟合得到直线的斜率求解kα,β方向的角误差信号与角度余弦偏差的关系曲线斜率kβ与α方向的kα的求解方法相同,进而求解出信号方向的角度余弦
最终解出目标信号的方位角和俯仰角:
在αt≤α3dB/2,βt≤β3dB/2时,uα与αt,uβ与βt近似成线性关系,随着在αt>α3dB/2,βt>β3dB/2时,随着αt,βt的增大uα与αt,uβ与βt呈现出非线性关系,对非线性情况进行此时可采用分段拟合的处理方式或者限幅处理进行拟合求解目标信号的方位角和俯仰角。
步骤4-2.根据单脉冲和差波束比幅测角得到测量角度求得对应角度下的方向余弦为:
以目标在三个方向的方向余弦及其变化的速度为状态向量:
其中为方向余弦变化的速度,应用α-β算法跟踪滤波得到入射角度估计值作为下一次测角的初始角度;
α-β滤波器是一种针对匀速运动目标模型的常增益滤波器,其中的增益矩阵具有如下形式:
其中和分别为目标状态的位置和速度分量的常滤波增益,T为采样周期。先定义三个方向机动指标为:
其中σ1为过程噪声和为三个方向观测噪声的标准差。以和为例表达式如下式,和求解方法相同
状态转移矩阵为
观测矩阵为
下面为滤波具体步骤,离散化的目标运动状态方程为:
观测方程为z(n)=[xc,yc,zc],滤波步骤具体如下:
状态一步预测:
由观测信号z(n)计算新息过程:
对于α-β滤波算法,利用K(n),进行状态估计:
得到状态估计值后,提取向量中第一个,第三个,第五个参数的值,即提取滤波之后的角度余弦并通过反三角函数关系求解出作为下一次测角的输入值,进行闭环的测角和跟踪的过程。
步骤4-3.在滤波过程中,受多种因素影响,观测数据中可能有部分数据与真实值之间存在很大误差,这些数据被称为野值;野值对跟踪性能会产生较大的影响,需要剔除;在进行α-β滤波时,计算新息过程时应将其模值与某一门限W做判决,若满足:
则认为观测数据有效,否则判为野值;通常门限值可选为W选为三个方向观测噪声标准差各自的3倍所组成的向量;
步骤4-4.若当前观测值被判为野值,则在状态估计时利用前一周期的状态向量估计值做线性外推,当连续3个观测时刻均检测到野值时判定为跟踪目标丢失,此时应当重新初始化α-β跟踪滤波器重新执行步骤3至步骤4所述的整个搜索和捕获过程;
步骤4-5.滤波过程中将角度余弦残差|αt|和|βt|与门限做判决,门限值通常取1/6至1/3波束宽度,当|αt|和|βt|均小于所述门限值时则认为捕获过程收敛,目标捕获完成可以转入跟踪过程,否则,将重复执行步骤4-1至步骤4-4直至捕获过程收敛,完成目标捕获。
Claims (3)
1.一种相控阵雷达搜索捕获目标的方法,具体包括以下步骤:
步骤1.相控阵雷达子阵划分;
将相控阵雷达天线划分成四个子阵,分别记为子阵1、子阵2、子阵3及子阵4;
步骤2.设定能量检测门限Et;
记n1(n)、n2(n)、n3(n)、n4(n)分别为所述四个子阵均未接收到捕获目标的信号时相应的输出信号,则此时合成信号记为n∑=n1(n)+n2(n)+n3(n)+n4(n),由此得噪声总能量:
<mrow>
<msub>
<mi>E</mi>
<mi>&Sigma;</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msup>
<msub>
<mi>n</mi>
<mi>&Sigma;</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
求解所述噪声总能量的均值以及方差σ∑,得能量检测门限其中k1∈[1.0,1.2]、k2∈[4,5];
步骤3.目标搜索过程;从载体平台提供的引导信息获得目标的初始位置后,控制相控阵雷达天线的波束位置指向该初始位置,并以该初始位置为中心进行波束搜索;波束搜索具体采用以下方式进行:
步骤3-1.构建两种不同范围的用于波束搜索的波位排列,分别记为排列1、排列2;
所述排列1包括5个波位:中心波位1,分别位于中心波位1的左上方、右上方、左下方及右下方的波位2、3、4及波位5;所述中心波位1与波位2、3、4、5相交叠,所述排列1中波位2、3、4、5的中心均与波位1的中心相距半个波束宽度;
所述排列2包括呈九宫格排列的九个波位:波位1位于九宫格中心,所述九宫格第一排由左至右分别为波位2、3及波位4,波位1的左、右分别为波位5、波位6,九宫格第三排由左至右分别为波位7、8、9;排列2中左右或上下相邻的两个波位中心之间均相距一个波束宽度;
步骤3-2.以目标的初始位置为中心波位1构建排列1所示的波束搜索排列方式,将相控阵雷达的波束指向调整至目标的初始位置,计算此时四个子阵合成信号的能量F∑;
所述合成信号的能量具体由以下方法获得:
以子阵1为参考子阵,分别对每个子阵做波束形成,得到相应的输出信号y1(n)、y2(n)、y3(n)及y4(n),将所述四个输出信号合成得Y∑(n)=y1(n)+y2(n)+y3(n)+y4(n),由此求解出四个子阵合成信号的能量:
其中N为时域上的总采样点数;
若所述能量F∑不小于所述能量检测门限Et,则目标搜索完成,记此时相控阵雷达的波束指向为捕获用波束指向,执行步骤4;否则,将相控阵雷达的波束指向依次调整至波位2、3、4、5直至四个子阵合成信号的能量F∑不小于所述能量检测门限Et时止;若将相控阵雷达的波束指向调整至波位5时,四个子阵合成信号的能量F∑依旧小于所述能量检测门限Et,执行步骤3-3;
步骤3-3.以目标的初始位置为中心波位1构建排列2所示的波束搜索排列方式,将相控阵雷达的波束指向调整至目标的初始位置,计算此时四个子阵合成信号的能量;
若所述能量F∑不小于所述能量检测门限Et,则目标搜索完成,记此时相控阵雷达的波束指向为捕获用波束指向,执行步骤4;否则,将相控阵雷达的波束指向依次调整至波位2、3、4、5、6、7、8、9直至四个子阵合成信号的能量F∑不小于所述能量检测门限Et时止;若将相控阵雷达的波束指向调整至波位9时,四个子阵合成信号的能量F∑依旧小于所述能量检测门限Et,则载体平台重新初始化获得最新的引导信息并向相控阵雷达提供目标的最新初始位置,并再次执行步骤3-1至步骤3-3;
步骤4.目标捕获过程;
步骤4-1.以步骤3所得的捕获用波束指向作为类单脉冲测角的初始值,利用和差波束比幅测角方法求解目标信号的方位角和俯仰角
步骤4-2.根据和差波束比幅测角得到测量角度求得对应角度下的方向余弦为:
以目标在三个方向的方向余弦及其变化的速度为状态向量:
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
</mtd>
<mtd>
<mover>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
<mo>&CenterDot;</mo>
</mover>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mi>c</mi>
</msub>
</mtd>
<mtd>
<mover>
<msub>
<mi>y</mi>
<mi>c</mi>
</msub>
<mo>&CenterDot;</mo>
</mover>
</mtd>
<mtd>
<msub>
<mi>z</mi>
<mi>c</mi>
</msub>
</mtd>
<mtd>
<mover>
<msub>
<mi>z</mi>
<mi>c</mi>
</msub>
<mo>&CenterDot;</mo>
</mover>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
其中为方向余弦变化的速度,应用α-β算法跟踪滤波得到入射角度估计值作为下一次测角的初始角度;
步骤4-3.在进行α-β滤波时,计算新息过程α(n)时应将其模值与门限W做判决,若满足:
<mrow>
<mo>|</mo>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>=</mo>
<mo>|</mo>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mi>z</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>|</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo><</mo>
<mi>W</mi>
<mo>,</mo>
</mrow>
其中z(n)为α-β滤波过程中的观测信号;
H(n)为观测矩阵;
F(n-1)为状态转移矩阵;
<mrow>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
</mtd>
<mtd>
<mover>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
<mo>&CenterDot;</mo>
</mover>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mi>c</mi>
</msub>
</mtd>
<mtd>
<mover>
<msub>
<mi>y</mi>
<mi>c</mi>
</msub>
<mo>&CenterDot;</mo>
</mover>
</mtd>
<mtd>
<msub>
<mi>z</mi>
<mi>c</mi>
</msub>
</mtd>
<mtd>
<mover>
<msub>
<mi>z</mi>
<mi>c</mi>
</msub>
<mo>&CenterDot;</mo>
</mover>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>;</mo>
</mrow>
则认为观测数据有效,并进行状态估计,求解下一次测角的初始角度;否则新息过程α(n)判为野值;门限值选为W,W选为三个方向观测噪声标准差各自的3倍所组成的向量;
步骤4-4.若当前观测值被判为野值,则在状态估计时利用前一周期的状态向量估计值做线性外推,当连续3个观测时刻均检测到野值时判定为跟踪目标丢失,此时应当重新初始化α-β跟踪滤波器重新执行步骤3至步骤4所述的整个搜索和捕获过程;
步骤4-5.滤波过程中将角度余弦残差|αt|和|βt|与门限做判决,门限值通常取1/6至1/3波束宽度,当|αt|和|βt|均小于所述门限值时则认为捕获过程收敛,目标捕获完成可以转入跟踪过程,否则,将重复执行步骤4-1至步骤4-4直至捕获过程收敛,完成目标捕获。
2.根据权利要求1所述的相控阵雷达搜索捕获目标的方法,其特征在于,步骤4所述的方位角和俯仰角的具体求解方法如下:
其中,
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>&alpha;</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<msub>
<mi>u</mi>
<mi>&alpha;</mi>
</msub>
<mo>/</mo>
<msub>
<mi>k</mi>
<mi>&alpha;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>&beta;</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<msub>
<mi>u</mi>
<mi>&beta;</mi>
</msub>
<mo>/</mo>
<msub>
<mi>k</mi>
<mi>&beta;</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
2
所述和差比幅测角方法采用的波束排列包括5个波位:中心波位5,分别位于中心波位5的正上方、正下方、左方及右方的波位1、2、3及波位4;所述中心波位5与波位1、2、3、4相交叠,所述排列中波位1、2、3、4的中心均与波位5的中心相距半个波束宽度;
对5个波束方位进行接收信号合成,可以得到5个方位的和波束接收信号能量F∑1~F∑5,所述能量可按照式F∑的求解方法求得,求取误差电压:
<mrow>
<msub>
<mi>u</mi>
<mi>&alpha;</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&Sigma;</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>|</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>|</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&Sigma;</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&Sigma;</mi>
<mn>5</mn>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<mo>,</mo>
<msub>
<mi>u</mi>
<mi>&beta;</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&Sigma;</mi>
<mn>4</mn>
</mrow>
</msub>
<mo>|</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>|</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&Sigma;</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>&Sigma;</mi>
<mn>5</mn>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
令相控阵天线阵列所在平面为x-o-y平面,捕获过程中相控阵雷达的入射信号初始指向与x轴及y轴的夹角分别为αx和αy;入射信号与z轴夹角为俯仰角θ,其取值范围为0°~90°,入射信号在x-o-y平面的投影与x轴夹角为方位角其取值范围为0°~360°,(αx,αy)与有如下关系:
其中,斜率kα和kβ可以通过实验拟合得到:求解kα时假设β方向指向为目标方向,在α0附近左右偏离一个波束宽度范围里对应uα与αt的值并用直线逼近uα-αt的关系曲线,由实测数据拟合得到直线的斜率求解kα;β方向的角误差信号与角度余弦偏差的关系曲线斜率kβ与上述求解方法相同。
3.根据权利要求1所述的相控阵雷达搜索捕获目标的方法,其特征在于,所述α-β滤波方法具体如下:
离散化的目标运动状态方程为:
<mrow>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
</mtd>
<mtd>
<mover>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
<mo>&CenterDot;</mo>
</mover>
</mtd>
<mtd>
<msub>
<mi>y</mi>
<mi>c</mi>
</msub>
</mtd>
<mtd>
<mover>
<msub>
<mi>y</mi>
<mi>c</mi>
</msub>
<mo>&CenterDot;</mo>
</mover>
</mtd>
<mtd>
<msub>
<mi>z</mi>
<mi>c</mi>
</msub>
</mtd>
<mtd>
<mover>
<msub>
<mi>z</mi>
<mi>c</mi>
</msub>
<mo>&CenterDot;</mo>
</mover>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>;</mo>
</mrow>
观测信号为z(n)=[xc,yc,zc],滤波过程具体如下:
状态下一步预测:
<mrow>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>|</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中F(n-1)表示状态转移矩阵;
由观测信号z(n)计算新息过程:
<mrow>
<mover>
<mi>z</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>|</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>|</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mi>z</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>|</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中H(n)表示观测矩阵;
α-β滤波器是一种针对匀速运动目标模型的常增益滤波器,其增益矩阵为K(n),利用K(n),进行状态估计:
<mrow>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>|</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>|</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>&lsqb;</mo>
<mover>
<mi>z</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mover>
<mi>z</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>|</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
提取向量中第一、第三、第五个参数的值,即提取滤波之后的方向余弦并通过反三角函数关系求解出θ′,作为下一次测角的输入值,进行闭环的测角和跟踪的过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510179314.7A CN104777478B (zh) | 2015-04-16 | 2015-04-16 | 一种相控阵雷达搜索捕获目标的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510179314.7A CN104777478B (zh) | 2015-04-16 | 2015-04-16 | 一种相控阵雷达搜索捕获目标的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104777478A CN104777478A (zh) | 2015-07-15 |
CN104777478B true CN104777478B (zh) | 2017-10-17 |
Family
ID=53619057
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510179314.7A Active CN104777478B (zh) | 2015-04-16 | 2015-04-16 | 一种相控阵雷达搜索捕获目标的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104777478B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106932750A (zh) * | 2015-12-31 | 2017-07-07 | 北京橙鑫数据科技有限公司 | 目标设备方位确定方法及装置 |
CN105785354B (zh) * | 2016-03-11 | 2018-02-02 | 中国人民解放军国防科学技术大学 | 一种相控阵雷达探测高速运动目标的搜索参数优化方法 |
US10594358B2 (en) | 2017-04-21 | 2020-03-17 | Futurewei Technologies, Inc. | Leakage signal cancellation |
CN107329136B (zh) * | 2017-06-13 | 2020-12-29 | 电子科技大学 | 基于可变分析时刻的mimo雷达多目标自适应跟踪方法 |
CN108038483A (zh) * | 2017-10-13 | 2018-05-15 | 西安天和防务技术股份有限公司 | 一种主动式智能跟踪方法及***、存储介质及电子终端 |
CN107944337A (zh) * | 2017-10-13 | 2018-04-20 | 西安天和防务技术股份有限公司 | 一种低空目标智能跟踪方法及***、存储介质及电子终端 |
CN109782268A (zh) * | 2018-12-12 | 2019-05-21 | 武汉船舶通信研究所(中国船舶重工集团公司第七二二研究所) | 基于相控阵和差波束测角的目标跟踪捕获方法及装置 |
CN109949615B (zh) * | 2019-04-08 | 2021-03-12 | 中国长江电力股份有限公司 | 一种船舶低航速预警检测***及方法 |
CN110108301B (zh) * | 2019-05-14 | 2020-12-01 | 苏州大学 | 模值检测动基座鲁棒对准方法 |
CN110441748A (zh) * | 2019-06-24 | 2019-11-12 | 西安电子科技大学 | 一种基于幅度信息的α-β滤波方法 |
CN111077518B (zh) * | 2019-12-20 | 2020-11-06 | 哈尔滨工业大学 | 一种基于距离-多普勒量测的跟踪滤波方法及装置 |
CN111693985B (zh) * | 2020-06-24 | 2023-08-04 | 北京无线电测量研究所 | 雷达目标跟踪方法、雷达、计算机设备和存储介质 |
CN111708015B (zh) * | 2020-07-13 | 2022-06-21 | 哈尔滨工业大学 | 一种多径效应下的低空目标跟踪滤波方法 |
CN112235030B (zh) * | 2020-09-01 | 2023-05-05 | 航天科工空间工程发展有限公司 | 一种低轨通信卫星设计方法、装置、存储介质和设备 |
CN113447919B (zh) * | 2021-06-29 | 2022-09-02 | 重庆大学 | 扩展卡尔曼预测角度跟踪方法 |
CN113805169B (zh) * | 2021-08-11 | 2024-05-03 | 航天恒星科技有限公司 | 空间目标低功耗小卫星雷达搜索与跟踪方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100753844B1 (ko) * | 2006-09-29 | 2007-08-31 | 한국전자통신연구원 | 레이다 추적장치 및 그 방법 |
US20090315760A1 (en) * | 2007-06-01 | 2009-12-24 | Intelwaves Technologies Ltd. | Hybrid tracking control system and method for phased-array antennae |
CN101296018B (zh) * | 2008-05-26 | 2011-11-23 | 中国电子科技集团公司第五十四研究所 | 移动卫星通信相控阵天线波束形成跟踪方法 |
CN101887120B (zh) * | 2010-06-10 | 2012-09-05 | 哈尔滨工业大学 | 一种用于单波束机械扫描雷达测量目标方位角的方法 |
CN102522632B (zh) * | 2011-12-19 | 2013-11-27 | 东南大学 | 卫星移动通信相控阵天线的分子阵数字化波束跟踪方法 |
-
2015
- 2015-04-16 CN CN201510179314.7A patent/CN104777478B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104777478A (zh) | 2015-07-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104777478B (zh) | 一种相控阵雷达搜索捕获目标的方法 | |
CN102622520B (zh) | 一种机动目标跟踪的分布式多模型估计融合方法 | |
CN105572664A (zh) | 基于数据融合的组网导航雷达目标跟踪*** | |
CN103235297B (zh) | 基于宽带雷达观测的空间章动目标参数估计方法 | |
CN106980114A (zh) | 无源雷达目标跟踪方法 | |
CN105093198B (zh) | 一种分布式外辐射源雷达组网探测的航迹融合方法 | |
CN107064880A (zh) | 分布式多基雷达收发波束同步和波束控制的高精度方法 | |
CN105652255B (zh) | 雷达组网***的空间配准方法 | |
CN107300700B (zh) | 敏捷合成孔径雷达卫星聚束模式姿态机动需求计算方法 | |
CN101482605A (zh) | 加宽零陷的高动态卫星导航***干扰抑制方法 | |
CN110516193B (zh) | 笛卡尔坐标系下基于转换瑞利滤波器的机动目标跟踪方法 | |
CN110531315B (zh) | 一种基于信号强度变化率的卫星干扰源直接定位方法及定位装置 | |
CN106324591A (zh) | 一种基于相控阵体制雷达的目标多模式跟踪方法 | |
CN106772466A (zh) | 一种基于形状特征搜索的近地卫星目标自动捕获算法 | |
RU2699552C1 (ru) | Способ пассивной однопозиционной угломерно-доплеровской локации перемещающихся в пространстве радиоизлучающих объектов | |
CN109471100A (zh) | 一种sar多普勒调频率估计方法及*** | |
Chen et al. | Research on AIS and radar information fusion method based on distributed Kalman | |
CN112966366B (zh) | 遥感卫星对地姿态机动扫描的轨迹指向角计算方法及*** | |
CN106961837B (zh) | 多相位中心多波束合成孔径雷达方位向相位校正与信号处理方法 | |
CN105954730B (zh) | 一种sar回波快速时域生成方法 | |
CN108828585B (zh) | 一种空间平台目标无伺服跟踪测角方法 | |
CN104330772B (zh) | 基于多向寻优的全跟踪式ukf滤波算法的双站定位方法 | |
Wang et al. | Research on multi-maneuvering target tracking JPDA algorithm | |
CN106066471A (zh) | 一种移动目标的定位方法及*** | |
CN105242262A (zh) | 一种基于天线周期扫描时间差无源定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20200624 Address after: 610000 No. 1301-1304, 13th floor, No. 2 Building, Group B, 200 Tianquan Road, Chengdu High-tech Zone, Sichuan Province Patentee after: CHENGDU DIANKE HUIAN TECHNOLOGY Co.,Ltd. Address before: 611731 Chengdu province high tech Zone (West) West source Avenue, No. 2006 Patentee before: UNIVERSITY OF ELECTRONIC SCIENCE AND TECHNOLOGY |
|
TR01 | Transfer of patent right |