CN114114240B - 超稀疏阵列在栅瓣影响下的三维目标跟踪方法和装置 - Google Patents
超稀疏阵列在栅瓣影响下的三维目标跟踪方法和装置 Download PDFInfo
- Publication number
- CN114114240B CN114114240B CN202111295109.9A CN202111295109A CN114114240B CN 114114240 B CN114114240 B CN 114114240B CN 202111295109 A CN202111295109 A CN 202111295109A CN 114114240 B CN114114240 B CN 114114240B
- Authority
- CN
- China
- Prior art keywords
- radar
- kalman filter
- equation
- measurement
- target
- 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
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000005259 measurement Methods 0.000 claims abstract description 83
- 239000011159 matrix material Substances 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 230000026676 system process Effects 0.000 claims description 2
- 230000003190 augmentative effect Effects 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 12
- 238000012545 processing Methods 0.000 abstract description 10
- 230000005856 abnormality Effects 0.000 description 7
- 238000001514 detection method Methods 0.000 description 7
- 230000003416 augmentation Effects 0.000 description 6
- 238000005070 sampling Methods 0.000 description 6
- 238000013461 design Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000005855 radiation Effects 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
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/72—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
- G01S13/723—Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
- G01S13/726—Multiple target tracking
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
- G01S7/352—Receivers
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
技术领域
本发明属于雷达数据处理技术领域,具体涉及超稀疏阵列在栅瓣影响下的三维目标跟踪方法。
背景技术
目标跟踪是雷达***的主要功能,而栅瓣是一种对目标航迹估计的准确性影响比较大的因子。栅瓣的产生与雷达天线自身的设计密不可分。当对雷达***的分辨率有较高要求时,往往需要增大天线阵列孔径,在阵列信号处理中等效于减小波束宽度,从而提高方位估计的精度。最直观的方法是按照处理信号的频带的上限来设定阵元间距,通过直接增加阵元数目的方式来增大阵列孔径,然而这种大面积的密集阵元布置方式成本过于高昂,往往难以实现。为了降低天线阵列的制造成本,降低其复杂度,只能采用稀疏布置天线单元的设计,通过减少阵元数及通道数的方式实现用较少的阵元数目达到较大的阵列孔径,用较低的成本实现期望的阵列性能。然而较大的阵元间距会导致栅瓣效应的出现,栅瓣也会产生对应的来波,使得跟踪目标时出现方位模糊,从而无法准确判定目标的波达方向。当存在多目标时,强目标的栅瓣会对信号特征较弱的目标的检测与估计造成严重干扰,甚至可能导致弱目标无法被检测到。而在三维目标跟踪中,方位角和俯仰角都会受到栅瓣多值的影响,出现测角多值的现象。
雷达进行角度测量通常采用的方法是相位干涉仪测向定位技术,通过比较两个天线的相位来获得方向。单基线相位干涉原理如图1所示。
来波方向与法线方向的夹角为φ,则平面波前到达天线单元1和天线单元2的时间就有先有后,体现在固定频率信号上就存在相位差。为了提高波达角的估计精度,通常采用增大基线长度的方法。但是当基线长度大于半波长时,由于栅瓣的存在,不能分辨来波的真正方向,会产生多个可能结果,这就是角度模糊现象。此时相位差是一个以2π为模糊的观测值,角度观测值与实际值的对于关系为其中k为多个未知的模糊值,它的个数与栅瓣的个数有关。这种情况的出现意味着在同一个采样点会同时出现k+1个观测值,其中一个为实际值,其他为模糊值。考虑到栅瓣往往是成对出现,而且对称分布于主瓣两侧,k的取值往往也是成对的正负数。
设阵元间距为d,波长为λ,φ0为主瓣对应角度,当的时候,这些φ对应的方向上也会存在与主瓣幅度相同的栅瓣。为了避免栅瓣出现,必须要求。假设目标方位与阵列所成夹角为φs,根据自然指向性函数可以预测栅瓣出现位置和宽度,第m个栅瓣出现的方位为
φm=arccos(cosφs±mλ/d),m=1,2,…,0°<θm<180°
φm=arccos(cos(2π-φs)±mλ/d),m=1,2,…,180°<θm<360°
或者
φm=arcsin(sinφs±mλ/d),m=1,2,…,-90°<θm<90°
φm=arcsin(sin(2π-φs)±mλ/d),m=1,2,…,90°<θm<270°
当m=0时为主极大位置,测量区域范围内的其他整数m出现的位置都为栅瓣位置,在测量区域内可能存在多个栅瓣。当阵元间距为波长一半,即d=λ/2,则没有栅瓣存在。当d=2λ,而且扫描范围为0。到180。,对应的m=-2,-1,1,2,共计4个栅瓣,对应的角度为180。,120。,60。,0。。当d=3λ,对应的栅瓣则会有6个,以此类推。
当目标位置变化的维度是三维的时候,会有两个角度用于表示目标位置,一个为方位角,另一个为俯仰角,而栅瓣会对方位角和俯仰角都产生测角多值的影响。
目前应对针对栅瓣引发的测角模糊现象的主要处理思路是对雷达天线的阵元布置方式进行优化,从而避免栅瓣的出现,这是从直接消除栅瓣这一思路上解决栅瓣影响问题。通过天线阵元优化布置的方式实现阵列栅瓣影响抑制,包括以下几类方法。
一类常用的方法是利用组合阵法或优化稀布阵技术的手段避免稀疏天线阵引起的栅瓣影响。对于此类方法的研究有包括:1)采取优化天线阵面排布,阵元形成非周期排列的方式来降低稀疏阵列引起的栅瓣影响。采用子阵重叠的方式布阵是非周期单元布阵的一种方式,这种方法是使子阵方向图在电扫范围内尽可能接***顶阵因子方向图,而在扫描范围外无能量辐射,目的就是再次分散栅瓣的能量,如随机错位的子阵防止栅瓣产生,或利用子阵的随机排布来打破阵面布局的周期性,或采用单元不等间距布阵形式来打破子阵的规则性、采用不同阵元间距的接收阵列和发射阵列避免栅瓣的产生。2)通过对稀疏阵列进行优化设计,使得利用最少的有源阵元数将栅瓣谱级抑制到满足要求的范围内。打乱阵面的周期结构以形成非周期布阵,并将阵中辐射单元设计成高效辐射单元,二者之间有机结合,共同抑制栅瓣出现。这类方法需要对阵列的形状进行设计,其对抗栅瓣影响的性能主要由阵列形状决定。
另一类方法是通过虚拟内插阵元法,在信号处理中降低阵列稀疏程度来抵抗阵稀疏引起的栅瓣影响,通过合理的设计收发阵元的布局,可以虚拟出更多数量的有效虚拟阵元数,通过这种方式可以达到更多的有效阵元数,利用阵列的优化配置和虚拟孔径的方法来实现对栅瓣的抑制。其它的处理方法还有将时域互相关法和空间域处理方法相结合的适用于***子阵的处理方法、阵列声强器法以及对步进信号的栅瓣影响抑制算法、在栅瓣位置处采用空间滤波器抑制电磁波的功率传输等等。
上述提到的方法都是从根源上解决问题,通过消除栅瓣的方式避免问题,然而在有些情况下栅瓣的出现是无法避免的,例如机动式的分布式雷达,这种雷达***可看作是一个随时间改变阵元排布的巨型天线阵列,其形状很难维持在不出现栅瓣的状态。此外阵元布置优化是针对特定波长进行对应的设置,而雷达波束的波长并不是固定不变的,当波长变短,而原有的阵列没有变化时,也可能会导致栅瓣的出现。当栅瓣已经出现时,要考虑消除栅瓣引发的测角模糊,识别出主瓣对应的测角,使得雷达能够在栅瓣影响下依然能够实现对目标的持续跟踪,并保证较高的估计精度。
发明内容
针对上述不足,提出了一种基于卡尔曼滤波器组的目标跟踪方法,利用卡尔曼滤波器的不同表现区分主瓣和栅瓣对应的测角,从而提高跟踪精度。
超稀疏阵列雷达在栅瓣模糊下的三维目标跟踪方法,其特征在于,包括:
分别建立超稀疏阵列雷达在栅瓣模糊下的雷达状态方程和雷达量测方程;
根据所述雷达状态方程和所述雷达量测方程,建立用于辨别栅瓣和估计目标状态的卡尔曼滤波器;
初始化所述卡尔曼滤波器;
获取所述超稀疏阵列雷达在各个时刻的雷达目标量测信息;
根据初始设置值、雷达状态方程和雷达量测方程,进行目标跟踪,直至目标跟踪结束。
超稀疏阵列雷达在栅瓣模糊下的三维目标跟踪装置,其特征在于,所述装置包括:
方程构建模块,用于分别建立超稀疏阵列雷达在栅瓣模糊下的雷达状态方程和雷达量测方程;
卡尔曼滤波器构建模块,用于根据所述雷达状态方程和所述雷达量测方程,建立用于辨别栅瓣和估计目标状态的卡尔曼滤波器;
初始化模块,用于初始化所述卡尔曼滤波器;
量测信息获取模块,用于获取所述超稀疏阵列雷达在各个时刻的雷达目标量测信息;
目标跟踪模块,根据初始设置值与雷达目标量测信息,进行目标跟踪,直至目标跟踪结束。
本发明创造的优点:
基于卡尔曼滤波器组的测角模糊分辨方法,实现方式是利用多个卡尔曼滤波器对栅瓣和主瓣共同产生的多个测角一一对应分别进行参数估计,根据不同的卡尔曼滤波器得出的偏差值区分主瓣和栅瓣,从而使得后续的信号处理直接使用主瓣对应的信息,避免了栅瓣造成的影响。依据这种方法,可以使得栅瓣的影响在信号处理过程进行之前就被排除,从而消除了测角模糊现象,并获得较为准确的目标方位估计信息,实现较为精确的三维目标跟踪。
附图说明
图1是单基线相位干涉原理图;
图2是***结构示意图;
图3是仿真结果示意图;
具体实施方式
下面结合附图对本发明作进一步描述,应当理解,此处所描述的内容仅用于说明和解释本发明,并不用于限定本发明。
本申请的应用场景在于当栅瓣的出现无法避免时,利用卡尔曼滤波器组区分主瓣和栅瓣对应的信号,从而在信号处理过程开始前消除栅瓣造成的影响,实现三维目标跟踪,保证对目标方位的估计结果具有较高的精度。解决现有技术中缺乏当栅瓣出现后的应对方法的问题,为分布式超稀疏阵列的三维目标跟踪提供保障。
根据栅瓣形成机理,利用卡尔曼滤波从测角模糊现象中分辨出哪个测角为主瓣对应的观测值,需要对测角模糊现象进行建模,从而与卡尔曼滤波器相结合。测角模糊现象可看作是在同一个采样点存在多个不同的测量角,但在实际测量中,由于雷达天线扫描范围的限制,往往在同一个采样点存在的测量角的数量不满足模糊值成对出现的情况,具有一定的随机分布特性。
在雷达天线实际工作过程中,由于雷达扫描范围与主瓣和栅瓣位置的关系,可能会出现以下三种情况:1、在一个采样点上得到的测量值包含对应主瓣的测量值,而栅瓣对应的测量值是成对出现的;2、在一个采样点上得到的测量值包含对应主瓣的测量值,而栅瓣对应的测量值不成对出现的;3、在一个采样点上得到的测量值不包含对应主瓣的测量值。对于卡尔曼滤波器组而言,只能分辨哪个是主瓣,从而根据主瓣信息计算目标方位,而无法从栅瓣信息中得出目标的方位,因此只能针对前两种情况完成主瓣的辨识工作。对于第三种情况,考虑到进行目标跟踪时需要保证主瓣始终能够处于扫描范围之内,连主瓣信息都无法得到的时候,则目标跟踪可以被认为是无法进行,因此第三种情况在本申请算法中不考虑。
对于雷达扫描范围变化导致测角数量变化,可以针对不同的扫描范围,构建不同的卡尔曼滤波器组。考虑到栅瓣的数量以及分布是能够根据相关公式进行推算的,因此在雷达构造方式或者多雷达组合的分布方式能够确定的前提下可以直接推算出扫描范围内包含多少个栅瓣,从而确定卡尔曼滤波器的数量,完成后续的算法设计工作。并且在卡尔曼滤波器估计过程中,雷达扫描范围应当维持在初始设定状态,否则扫描范围的变化会导致卡尔曼滤波器组和实际情况无法对应,导致估计工作无法完成。
在异常检测里一般使用的卡尔曼滤波器,本质上是一类以仿真算法为基础的异常检测方法,此类方法实施流程存在两个部分:一方面是出现残差;另一方面是残差决策。残差的出现代表着以装置模型建立的一个或一组函数为基础去表现异常事件量的多少;残差的决策代表着以已确定的残差值为基础,建立合适的阈值及决策标准去探究异常来源,这是残差决策的基本流程。以卡尔曼滤波器为基础的异常检测流程包括以下内容:首先是比较卡尔曼滤波器的预测值和实际测定值的偏差,得到残差信号,也就是残差加权平方和(WSSR),其数量代表了异常状态的数量,以此判断装置有没有出现故障。其次是以残差的异常为基础进行评价,如果装置中不存在异常,装置残差值在理论上会较为平稳的处于某个小范围之间;但是如果存在异常,那么残差即将出现大幅度的误差,使用残差值能够对判断异常状况是否存在,判断出异常的类型、范围等有关数据,完成对异常的检测。
本申请中是针对测角多值的异常检测问题,对应的方法是以若干个输入信息存在区别的卡尔曼滤波器为基础构建的卡尔曼滤波器组去处理。卡尔曼滤波器组中涵盖有n个卡尔曼滤波器,每个卡尔曼滤波器的增广状态变量都去除了与模糊值k相关的部分,其中n为雷达接收到的测角组(每个测角组包含一个方位角与一个俯仰角)的总数,每个卡尔曼滤波器接收1个测角组的检测值。任意测角组对应的是主瓣还是栅瓣均通过1个与之对应的卡尔曼滤波器进行判断。卡尔曼滤波器组通过残差的测算之后,会得到n个残差。当第i个测角为主瓣对应测角时,第i个卡尔曼滤波器得到的信息不包含模糊值的干扰,其估计结果更为准确,残差值约等于0。但是其他的卡尔曼滤波器对应的测角则是栅瓣对应的测角,其估计过程会受到模糊值的干扰,导致最终的估计结果与测量值存在一定的偏差,由于预测偏差造成残差值的增大。按照n个卡尔曼滤波器的残差值数量级不同的比较,能够识别出其中的最小值,最小值和另外的残差值存在数量级的不同,据此可确定主瓣对应的测角。
在卡尔曼滤波器组进行测角模糊处理时,阈值的取值是主瓣判断的关键。设定WSSR对应的阈值在较小的取值上可以提高主瓣分辨成功的几率,但是由于噪声的干扰,会增加误判的几率,将主瓣也判别为栅瓣,而阈值设定在较大的取值上则会导致相反的结果,因此阈值的设定对混合卡尔曼滤波器组的主瓣分辨功能的表现有很大的影响。一种比较合理的阈值设定方式是根据雷达在不同监测区域上的量测精度的不同,建立对应的一系列阈值,从而提升测角模糊处理的准确性。
其于上述原理,本申请提供一种基于分布式稀疏阵列在栅瓣测角模糊影响下的目标跟踪方法,其实现步骤包括:
步骤S00:建立存在栅瓣影响情况下的雷达状态方程。
当栅瓣存在时,由于栅瓣的影响,雷达会产生多个测角值,不同的测角值可对应不同的模糊数,由于模糊数的出现并没有对被监测物体的状态参数产生影响,因此建立的状态方程不设置模糊数,对应的被监测物体的状态方程可以表示为下式所示
其中x、y和z为被监测目标的方位信息,和/>为被监测目标在x轴、y轴和z轴的速度分量,/>状态方程中包含***过程噪音W,其协方差矩阵为
Q=E(WWT)
步骤S01:建立具有栅瓣的分布式阵列雷达的量测方程。
雷达能够测量到的参数包括目标距离,目标速度以及由于栅瓣产生的多个测角。由于一开始无法确定哪个是主瓣哪个是栅瓣,因此多个不同瓣的测角可以进行编号处理,其中第i个瓣对应的输出方程为
其中r是被监测物体和雷达之间的距离,是被监测物体的速度,/>是被监测物体的方位角,θ是被监测物体的俯仰角,kfw,i和kfy,i分别为方位角和俯仰角对应的模糊值,输出方程中包含测量噪音W,其协方差矩阵为
R=E(WWT)
步骤S02:根据步骤S00状态方程,以及步骤S01所述的量测方程,建立用于辨别栅瓣和估计目标状态的卡尔曼滤波器模型。针对每个测角构建一个卡尔曼滤波器,然后组成卡尔曼滤波器组。
对应第i个测角的卡尔曼滤波器的状态空间模型是在步骤S00所述的状态方程和步骤S01所述的量测方程的基础上构建而成。与第i个瓣所对应的测角的卡尔曼滤波器的结构为
其中xkal,i、ykal,i和zkal,i为第i个卡尔曼滤波器估计的被监测目标的方位信息, 和/>为第i个卡尔曼滤波器估计的被监测目标在x轴、y轴和z轴的速度分量,rkal,i是卡尔曼滤波器估计的被监测物体和雷达之间的距离,/>是卡尔曼滤波器估计的被监测物体的速度,/>是卡尔曼滤波器估计的被监测物体的方位角,θkal,i是卡尔曼滤波器估计的被监测物体的俯仰角。WSSRi为第i个卡尔曼滤波器的估计结果与测量结果之间的残差加权平方和。根据不同卡尔曼滤波器计算的残差加权平方和的差异,判断哪个滤波器对应的测角是主瓣测角,从而避免栅瓣影响。
从上式中可以看出,卡尔曼滤波器的增广状态变量中并没有包含两个模糊值,这是因为对三维目标进行跟踪时,如果用卡尔曼滤波器直接估计栅瓣造成的模糊值,由于方位角和俯仰角都存在模糊值,这就使得卡尔曼滤波器需要估计五个增广状态变量,包括x、y、z三个轴对应的速度以及两个模糊值,而卡尔曼滤波器可用的可测输出却只有方位角、俯仰角、速度和距离这四个参数,如果增广状态变量包含两个模糊值,就会使得需要估计的参数多于卡尔曼滤波器的输入参数,进而导致卡尔曼滤波器的估计工作无法进行。
为了能够在不引入增广状态变量的前提下进行主瓣和栅瓣的区分,采用的方法是用其他的数值来间接体现模糊值所造成的影响。
该滤波器中增益矩阵K的计算过程如下所示。由于卡尔曼滤波器的增益矩阵无法直接通过非线性输出方程计算得到,因此需要对非线性方程进行小偏差线性化计算,从而得到
其中r是被监测物体和雷达之间的距离,是被监测物体的速度,/>是被监测物体的方位角,θ是被监测物体的俯仰角,V为测量噪声。H是输出方程
考虑到卡尔曼滤波器的鲁棒性,卡尔曼滤波器应满足对应的二次型指标,卡尔曼滤波器的增益需要使得二次型指标最小化。根据极值原理,可以导出最优的增益矩阵公式为
K=PHTR-1
其中P为时间更新矩阵,可以通过线性***对应的黎卡提方程进行计算
其中Q为***噪声W的协方差矩阵,Q=E(WWT),R为测量噪声V的协方差矩阵,R=E(VVT)。
由于构成卡尔曼滤波器的***为非线性***,因此增益矩阵K需要根据***运行状态的变化进行更新,具体的更新方程为
P=FPFT+Q
K=PHT(HPH+R)-1
P=P-KHP
步骤S03:分布式阵列雷达对目标进行观测,获得目标各个时刻的观测值。
步骤S04:根据步骤S02所述的卡尔曼滤波器和步骤S03所述的目标雷达观测值,对目标跟踪过程进行初始设置。
这一步骤的主要目的是进行卡尔曼滤波器组的初始化设定。首先根据接收到的观测值中的测角数量确定卡尔曼滤波器的增广状态变量的数量;然后在卡尔曼滤波器的估计工作开始之前赋予卡尔曼滤波器一组目标信息的初始值,并根据初始值计算卡尔曼滤波器内增益矩阵K的数值。
步骤S05:根据步骤S04的初始设置值和步骤S02所构建的卡尔曼滤波器组,进行由k时刻到k+1时刻的目标跟踪,直至目标跟踪结束。
由于初始值不可能与目标的实际状态完全吻合,会存在一定偏差,因此卡尔曼滤波器需要在后续的估计过程中计算新的估计值,使其不断逼近观测值,从而消除这种偏差。然而只有主瓣所对应的测角由于没有受到模糊值的干扰,主瓣对应的卡尔曼滤波器的估计偏差能够达到最小,它所对应的卡尔曼滤波器计算得到的加权残差平方和WSSR在数量级上是最小的,其他卡尔曼滤波器计算的WSSR都会比主瓣对应的计算结果在数量级上大一些,从而分辨出这些测角中哪个是主瓣对应的测角。
本申请提出的针对分布式稀疏阵列存在栅瓣时的目标跟踪方法,能够区分主瓣和栅瓣所对应的测角,并能够获得较为准确的目标信息。
实际应用中,如图2所示,某分布式基地雷达***可以看作是一个大型的稀疏天线阵,由于栅瓣的影响,雷达***接收到的目标信号包含多个不同的测角,卡尔曼滤波器根据接收到的信息对目标方位以及模糊数进行估计,当滤波器的估计信息与观测值具有较高吻合度时,滤波器所估计的状态参数能够比较准确地反应物体当前的状态。
以仿真实验验证该方法。仿真中被跟踪目标以固定速率运行,而雷达持续对目标进行跟踪。仿真中对和/>施加了均方差为0.05的过程噪声干扰,对距离r施加了均方差为1的测量噪声干扰,对方位角/>和俯仰角θ施加了均方差为0.05的测量噪声干扰,对速度V施加了均方差为0.05的测量噪声干扰。卡尔曼滤波器组需要在噪声干扰下处理测角模糊,分辨主瓣和栅瓣对应的测角,对目标实现有效的定位以及跟踪。
从仿真结果对应的图3可以看出,基于主瓣信息的卡尔曼滤波器计算得到的WSSR明显小于基于栅瓣信息得到的WSSR,根据WSSR的差异确定哪个是主瓣,从而避免栅瓣的影响。由于没有模糊值的干扰,主瓣对应的卡尔曼滤波器的估计结果更为准确,偏差较小,因此对应的WSSR更接近于0,而栅瓣对应的卡尔曼滤波器计算的WSSR数值则会比较大,通过这种差异,可以从多个测角中判断出哪个是主瓣对应的测角。
最后应说明的是:以上所述仅为本发明的解释,并不用于限制本发明,尽管对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.超稀疏阵列雷达在栅瓣模糊下的三维目标跟踪方法,其特征在于,包括:
分别建立超稀疏阵列雷达在栅瓣模糊下的雷达状态方程和雷达量测方程;
根据所述雷达状态方程和所述雷达量测方程,建立用于辨别栅瓣和估计目标状态的卡尔曼滤波器;
初始化所述卡尔曼滤波器;
获取所述超稀疏阵列雷达在各个时刻的雷达目标量测信息;
根据初始设置值、雷达状态方程和雷达量测方程,进行目标跟踪,直至目标跟踪结束;
所述建立超稀疏阵列雷达在栅瓣模糊下的雷达状态方程,
所述雷达状态方程满足下述关系式:
其中x、y和z为被监测目标的方位信息,和/>为被监测目标在x轴、y轴和z轴的速度分量,/>状态方程中包含***过程噪音W;
所述雷达量测方程满足下述关系式:第i个瓣对应的输出方程为
其中r是被监测物体和雷达之间的距离,是被监测物体的速度,/>是被监测物体的方位角,θ是被监测物体的俯仰角,kfw,i和kfy,i分别为方位角和俯仰角对应的模糊值,输出方程中包含测量噪音V;
所述卡尔曼滤波器满足下述关系式:建立一个卡尔曼滤波器组,每个滤波器对应一个测角,滤波器组中构建的与第i个瓣的测角对应的非线性卡尔曼滤波器的结构为
其中xkal,i、ykal,i和zkal,i为第i个卡尔曼滤波器估计的被监测目标的方位信息, 和/>为第i个卡尔曼滤波器估计的被监测目标在x轴、y轴和z轴的速度分量,rkal,i是卡尔曼滤波器估计的被监测物体和雷达之间的距离,/>是卡尔曼滤波器估计的被监测物体的速度,/>是卡尔曼滤波器估计的被监测物体的方位角,θkal,i是卡尔曼滤波器估计的被监测物体的俯仰角,WSSRi为第i个卡尔曼滤波器的估计结果与测量结果之间的残差加权平方和,根据不同卡尔曼滤波器计算的残差加权平方和的差异,判断哪个滤波器对应的测角是主瓣测角,从而避免栅瓣影响,K为增益矩阵。
2.根据权利要求1所述的方法,其特征在于,该滤波器中增益矩阵K的计算过程如下,对非线性输出方程进行小偏差线性化计算,从而得到
其中r是被监测物体和雷达之间的距离,是被监测物体的速度,/>是被监测物体的方位角,θ是被监测物体的俯仰角,V为测量噪声,H是对应的线性化参数矩阵,
其中H为输出方程的偏导数矩阵,增益矩阵公式为
K=PHTR-1
其中P为时间更新矩阵,通过线性***对应的黎卡提方程进行计算
其中Q为***噪声W的协方差矩阵,Q=E(WWT),R为测量噪声V的协方差矩阵,R=E(VVT)。
3.根据权利要求2所述的方法,其特征在于,增益矩阵K需要根据***运行状态的变化进行更新,更新方程为P=FPFT+Q
K=PHT(HPH+R)-1
P=P-KHP。
4.根据权利要求1所述的方法,其特征在于,卡尔曼滤波器组的初始化设定时,首先根据接收到的观测值中的测角数量确定卡尔曼滤波器的增广状态变量的数量;然后在卡尔曼滤波器的估计工作开始之前赋予卡尔曼滤波器一组目标信息的初始值,并根据初始值计算卡尔曼滤波器内增益矩阵K的数值。
5.根据权利要求4所述的方法,其特征在于,主瓣所对应的测角没有受到模糊值的干扰,所对应的卡尔曼滤波器计算得到的残差加权平方和WSSR在数量级上最小,其他卡尔曼滤波器计算的WSSR会比主瓣对应的计算结果在数量级上大,从而分辨出主瓣对应的测角。
6.一种采用权利要求1的方法的超稀疏阵列雷达在栅瓣模糊下的三维目标跟踪装置,其特征在于,所述装置包括:
方程构建模块,用于分别建立超稀疏阵列雷达在栅瓣模糊下的雷达状态方程和雷达量测方程;
卡尔曼滤波器构建模块,用于根据所述雷达状态方程和所述雷达量测方程,建立用于辨别栅瓣和估计目标状态的卡尔曼滤波器;
初始化模块,用于初始化所述卡尔曼滤波器;
量测信息获取模块,用于获取所述超稀疏阵列雷达在各个时刻的雷达目标量测信息;
目标跟踪模块,根据初始设置值、雷达状态方程和雷达量测方程,进行目标跟踪,直至目标跟踪结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111295109.9A CN114114240B (zh) | 2021-11-03 | 2021-11-03 | 超稀疏阵列在栅瓣影响下的三维目标跟踪方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111295109.9A CN114114240B (zh) | 2021-11-03 | 2021-11-03 | 超稀疏阵列在栅瓣影响下的三维目标跟踪方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114114240A CN114114240A (zh) | 2022-03-01 |
CN114114240B true CN114114240B (zh) | 2024-02-27 |
Family
ID=80380644
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111295109.9A Active CN114114240B (zh) | 2021-11-03 | 2021-11-03 | 超稀疏阵列在栅瓣影响下的三维目标跟踪方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114114240B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112014835B (zh) * | 2020-09-01 | 2023-05-26 | 中国电子科技集团公司信息科学研究院 | 分布式稀疏阵列雷达在栅瓣模糊下的目标跟踪方法和装置 |
CN114779203B (zh) * | 2022-06-20 | 2022-09-09 | 中国电子科技集团公司第五十四研究所 | 广域随机稀布阵空域扫描能量信息匹配的目标定位方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002181926A (ja) * | 2000-12-15 | 2002-06-26 | Furuno Electric Co Ltd | 目標追尾方法およびその装置 |
CN101813764A (zh) * | 2010-03-25 | 2010-08-25 | 电子科技大学 | 一种针对均匀超稀疏阵列天线波束指向模糊的抑制方法 |
CN107045125A (zh) * | 2017-03-17 | 2017-08-15 | 电子科技大学 | 一种基于预测值量测转换的交互多模型雷达目标跟踪方法 |
WO2019037996A1 (de) * | 2017-08-21 | 2019-02-28 | Astyx Gmbh | Abbildendes radarsystem mit einem empfangsarray zur winkelbestimmung von objekten in zwei dimensionen durch eine gespreizte anordnung der empfangsantennen einer dimension |
EP3588128A1 (en) * | 2018-06-26 | 2020-01-01 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Method for detection and height and azimuth estimation of objects in a scene by radar processing using sparse reconstruction with coherent and incoherent arrays |
CN110988811A (zh) * | 2019-12-20 | 2020-04-10 | 中国电子科技集团公司第二十研究所 | 一种应用于稀疏布阵宽带相控阵的栅瓣解模糊方法 |
CN111596290A (zh) * | 2020-06-01 | 2020-08-28 | 中国电子科技集团公司信息科学研究院 | 一种基于最大相关熵扩展卡尔曼滤波的雷达目标跟踪方法 |
CN111679266A (zh) * | 2020-06-18 | 2020-09-18 | 成都纳雷科技有限公司 | 汽车毫米波雷达稀疏阵列栅瓣虚假目标识别方法及*** |
CN112014835A (zh) * | 2020-09-01 | 2020-12-01 | 中国电子科技集团公司信息科学研究院 | 分布式稀疏阵列雷达在栅瓣模糊下的目标跟踪方法和装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100945466B1 (ko) * | 2008-05-07 | 2010-03-08 | 서강대학교산학협력단 | 초음파 영상 시스템의 비균일 가중 주기 희박 어레이의최적 설계 방법 |
US9562968B2 (en) * | 2012-10-22 | 2017-02-07 | Saab-Sensis Corporation | Sensor system and method for determining target location using sparsity-based processing |
-
2021
- 2021-11-03 CN CN202111295109.9A patent/CN114114240B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002181926A (ja) * | 2000-12-15 | 2002-06-26 | Furuno Electric Co Ltd | 目標追尾方法およびその装置 |
CN101813764A (zh) * | 2010-03-25 | 2010-08-25 | 电子科技大学 | 一种针对均匀超稀疏阵列天线波束指向模糊的抑制方法 |
CN107045125A (zh) * | 2017-03-17 | 2017-08-15 | 电子科技大学 | 一种基于预测值量测转换的交互多模型雷达目标跟踪方法 |
WO2019037996A1 (de) * | 2017-08-21 | 2019-02-28 | Astyx Gmbh | Abbildendes radarsystem mit einem empfangsarray zur winkelbestimmung von objekten in zwei dimensionen durch eine gespreizte anordnung der empfangsantennen einer dimension |
EP3588128A1 (en) * | 2018-06-26 | 2020-01-01 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Method for detection and height and azimuth estimation of objects in a scene by radar processing using sparse reconstruction with coherent and incoherent arrays |
CN110988811A (zh) * | 2019-12-20 | 2020-04-10 | 中国电子科技集团公司第二十研究所 | 一种应用于稀疏布阵宽带相控阵的栅瓣解模糊方法 |
CN111596290A (zh) * | 2020-06-01 | 2020-08-28 | 中国电子科技集团公司信息科学研究院 | 一种基于最大相关熵扩展卡尔曼滤波的雷达目标跟踪方法 |
CN111679266A (zh) * | 2020-06-18 | 2020-09-18 | 成都纳雷科技有限公司 | 汽车毫米波雷达稀疏阵列栅瓣虚假目标识别方法及*** |
CN112014835A (zh) * | 2020-09-01 | 2020-12-01 | 中国电子科技集团公司信息科学研究院 | 分布式稀疏阵列雷达在栅瓣模糊下的目标跟踪方法和装置 |
Non-Patent Citations (5)
Title |
---|
Multi-Phase-Center Sidelobe Suppression Method for Circular GBSAR Based on Sparse Spectrum;Yanping Wang et al.;IEEE Access;第8卷;全文 * |
一种随机稀布MIMO雷达的旁瓣抑制方法;陈金立;朱艳萍;李家强;;科学技术与工程(18);全文 * |
分布式卫星编队的误差校正方法及运动目标检测性能分析;刘颖;廖桂生;张涛麟;;电子与信息学报(10);全文 * |
基于最大相关熵的雷达扩展卡尔曼滤波算法研究;李春霞等;中国传媒大学学报(自然科学版);第27卷(第3期);全文 * |
角度模糊对天基稀疏孔径GMTI雷达的影响与抑制方法;陆必应;梁甸农;;电子学报(03);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114114240A (zh) | 2022-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114114240B (zh) | 超稀疏阵列在栅瓣影响下的三维目标跟踪方法和装置 | |
CN107290709B (zh) | 基于范德蒙分解的互质阵列波达方向估计方法 | |
CN112014835B (zh) | 分布式稀疏阵列雷达在栅瓣模糊下的目标跟踪方法和装置 | |
CN104777468B (zh) | 利用雷达天线和差通道抑制雷达旁瓣截获的方法 | |
CN108663653B (zh) | 基于l形电磁矢量传感器阵列的波达方向估计方法 | |
CN108872926B (zh) | 一种基于凸优化的幅相误差校正及doa估计方法 | |
CN106054153A (zh) | 一种基于分数阶变换的海杂波区目标检测与自适应杂波抑制方法 | |
KR100902559B1 (ko) | 레이더 간섭계 및 그것을 이용한 표적위치 추정방법 | |
CN106980105B (zh) | 电磁矢量传感器阵列空间旋转解相干测向方法 | |
CN117110980A (zh) | 基于fpga的自适应单脉冲测向方法 | |
CN109783960B (zh) | 一种基于网格部分细化的波达方向估计方法 | |
Zhang et al. | Detection of sea-surface small targets masked by range sidelobes of large objects | |
RU2615491C1 (ru) | Способ одновременного измерения двух угловых координат цели в обзорной амплитудной моноимпульсной радиолокационной системе с антенной решеткой и цифровой обработкой сигнала | |
Li et al. | Identifying unambiguous frequency patterns for two-target localization using frequency diverse array | |
RU2303274C1 (ru) | Способ радиопеленгования и радиопеленгатор для его осуществления | |
CN114740434A (zh) | 一种等距分布子阵列***及抗伴飞式干扰的方法 | |
CN115575906A (zh) | 非均匀背景下干扰智能抑制的融合检测方法 | |
CN115906372A (zh) | 一种基于稀疏阵幅相误差校正的无网格参数估计方法 | |
CN114609580A (zh) | 一种基于非圆信号的无孔互质阵列设计方法 | |
CN114460531A (zh) | 一种均匀线阵music空间谱估计方法 | |
CN118228189A (zh) | 一种水下主被动联合自导的多级融合方法 | |
CN111983555B (zh) | 一种基于神经网络的干涉仪角度解算方法 | |
CN117706467B (zh) | 一种干涉仪布阵方法 | |
Di Claudio et al. | Robust direction estimation of UWB sources in Laguerre-Gauss beamspaces | |
CN114460536A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |