基于径向加速度的RA-Singer-EKF机动目标跟踪算法制造技术

本发明专利技术公开了一种基于径向加速度的RA-Singer-EKF机动目标跟踪算法，属于雷达机动目标跟踪领域。本发明专利技术的方法可快速准确地提供机动目标的径向加速度和径向速度信息，有效提高雷达对机动目标的跟踪性能。本发明专利技术的方法包括以下步骤：(一)对雷达接收信号进行采样，利用匹配追踪(OMP)方法得到目标径向加速度和径向速度；(二)在数据处理阶段将径向加速度和径向速度进行坐标转换并引入量测方程和状态方程；(三)采用Singer模型和扩展卡尔曼滤波(EKF)算法实现机动目标跟踪。本发明专利技术能够精确实时地反映目标的机动情况，提高了目标跟踪精度，改善了速度和加速度估计精度，工程实现容易，具有较强的工程应用价值和推广前景。

【技术实现步骤摘要】

本专利技术隶属于雷达机动目标跟踪领域，适用于高/中脉冲重复频率雷达(如机载脉冲多普勒雷达等)对机动目标的精确跟踪。
技术介绍
随着科学技术的发展，现代军用飞机和导弹的机动能力大大增强，对雷达的检测和跟踪性能提出了新的挑战，对高机动目标跟踪也成为雷达目标跟踪领域的难点和重点。现有的机动目标跟踪技术是在雷达测量信息(位置、多普勒速度)的基础上，利用各种机动模型和滤波算法实现对机动目标的跟踪，其中，Singer模型算法是一种常用的有效机动目标跟踪算法，主要由以下3个步骤实现(I)将雷达接收机输出的目标回波信号进行A/D变换，送雷达数据处理计算机执行以下步骤；(2)将机动加速度作为零均值时间相关色噪声建模。设机动加速度的概率密度函数在[_Amax，Amax]上近似服从零均值均匀分布，并假设其时间相关函数为指数衰减形式；(3)采用卡尔曼滤波算法对机动目标进行自适应跟踪，估计出目标的位置、速度和加速度状态。这种方法具有以下缺陷由于现有雷达在目标跟踪时只能利用目标距离和方位角测量信息，只有多普勒雷达可以利用径向速度和距离、方位角测量信息，因此当目标发生机动时，加速度发生剧变，雷达由于缺少加速度测量信息，会出现跟踪精度低甚至跟踪发散的问题。
技术实现思路
本专利技术的目的是提出一种基于径向加速度的RA-Singer-EKF机动目标跟踪算法，解决现有Singer模型算法对目标突发机动的自适应跟踪能力不强，以及速度和加速度估计误差较大的问题。本专利技术提出的基于径向加速度的RA-Singer-EKF方法的技术方案包括以下步骤步骤1:将雷达接收机接收到的线性调频信号s (t)通过采样器以采样间隔Ts进行采样，变为离散信号s(nTs)，其中η表示采样点序号；将S(nTs)送入雷达信号处理计算机；在雷达信号处理计算机中执行以下步骤步骤2 :初始化(设置分解参数)T设为雷达的脉冲宽度；λ为雷达波长；fs为采样频率；fu设为LFM信号的初始频率；kv，设为LFM信号的调频率；G((UXV)XN))设为过完备原子库，UXV为原子库中原子的个数，N = T/fs;λ 2设为判断OMP分解是否完成的相干比阈值；步骤3:形成过完备原子库(I)根据 LFM 信号回波 s(t)的特点，建立原子 gr = exp[j2 π (fun+kvn2/2N) ], r =1,2,...，N;(2)设定搜索精度和范围，假设搜索范围fu的取值为fu e [O, U] Δ fu, U = 1，2，…，U，U为起始频率的搜索个数，Afu为多普勒单元，搜索范围Akv的取值为kve [O,V] Akv, Akv为调频率单元，V = 1,2,…，V，V为调频率的搜索个数，构造的过完备原子库G为(UXV) XN的矩阵本文档来自技高网...

【技术保护点】
一种基于径向加速度的RA?Singer?EKF机动目标跟踪方法，是指在信号处理阶段采用正交匹配追踪(OMP)思想对目标径向加速度和径向速度进行提取，在数据处理阶段将径向加速度和径向速度估计值通过坐标转换引入量测方程和和状态方程中，并采用Singer模型和扩展卡尔曼滤波(EKF)方法实现机动目标的精确跟踪，其特征在于包括以下步骤：?步骤1：将雷达接收机接收到的线性调频信号s(t)通过采样器以采样间隔Ts进行采样，变为离散信号s(nTs)，其中n表示采样点序号；将s(nTs)送入雷达信号处理计算机；?在雷达信号处理计算机中执行以下步骤：?步骤2：初始化(设置分解参数)?T设为雷达的脉冲宽度；?λ为雷达波长；?fs为采样频率；?Ts＝1/fs为采样间隔；?fu设为LFM信号的初始频率；?kv设为LFM信号的调频率；?G((U×V)×N))设为过完备原子库，U×V为原子库中原子的个数，N＝T/fs；?λ2设为判断OMP分解是否完成的相干比阈值；?步骤3：形成过完备原子库?(1)根据LFM信号回波s(t)的特点，建立原子gr＝exp[j2π(fun+kvn2/2N)]，r＝1，2，…，N；?(2)设定搜索精度和范围，假设搜索范围fu的取值为fu∈[0，U]△fu，u＝1，2，…，U，U为起始频率的搜索个数，△fu为多普勒单元，搜索范围Δkv的取值为kv∈[0，V]Δkv，Δkv为调频率单元，v＝1，2，…，V，V为调频率的搜索个数，构造的过完备原子库G为(U×V)×N的矩阵：?G＝[g1，g2，…，gN]T；?步骤4：正交匹配追踪(OMP)?本专利技术采用正交匹配追踪(OMP)方法对信号进行参数提取，和标准匹配追踪(MP)方法相比，该方法在信号参数估计时精度更高，并且稀疏分解收敛速度更快；首先将过完备原子库中的原子与信号进行匹配程度比较，选择与信号最匹配的一组基gr，将所选的原子利用?Grant?Schmidt正交化方法进行正交化处理，然后将信号在这些正交原子构成的空间下进行分解，使其满足?||＝sup||?因此，信号可以分解为在最佳原子上的分量和残余两部分，即?s＝gr+R1s?其中，R1s是用最佳原子对原信号进行最佳匹配后的残余；然后从原子库中将最匹配的这组基删掉，接下来对最佳匹配后的残余可以不断进行上面同样的分解过程，即：?Ris＝gri+Ri+1s?其中gri满足：?＝sup?经过i次迭代后能量衰减程度△Ri+1为：ΔRi+1＝‖Ris‖2?‖Ri+1s‖2＝ω2‖Ris‖2，其中，‖Ris‖2和‖Ri+1s‖2分别是信号s的i阶和i+1阶残余能量；ω2＝λ2(Ris)表示信号衰减的速率；当λ2(Ris)≥λ2时，停止分解；假设经过L步分解后停止分解，信号被分解为：?用少量的原子L(相对于信号长度N而言，L《N)表示信号的主要成分，即：?步骤5：径向加速度和速度的确定?(1)经过上述对信号s(t)的OMP分解，计算出原子库中匹配的原子，得到信号的稀疏解能量图i∈(0，U×V)；?(2)在i∈(0，U×V)范围搜索最大的峰值；?(3)找到最大峰值的坐标Am(i)，并在原子库G中查找此坐标的位置的gr(fu，kv)，就可以得到目标的初始频率fu和调频率kv，i＝1，2，…，N，最后根据以下公式得到机动目标的径向加速度和速度：?α＝kλ/2，v＝fdλ/2?步骤6：将得到的径向加速度和径向速度估计值送到雷达数据处理计算机，进行极(球)坐标系到直角坐标系的转换?根据速度与状态向量几何关系图可以得到径向速度和直角坐标系下X轴方向速度vx和Y轴方向速度vy的关系：?同理可得径向加速度和直角坐标系下X轴方向速度αx和Y轴方向速度αy的关系：?步骤7：建立Singer模型?将机动加速度作为零均值时间相关色噪声建模。该模型中设机动加速度的概率密度函数在[?Amax，Amax]上近似服从零均值均匀分布，其时间相关函数为指数衰减形式：?式中，α是在区间(t，t+τ)内决定目标机动特性参数的待定参数；α为机动频率，即机动时间常数的倒数，其取值与具体机动状态有关；为目标的加速度方差，可由α(t)概率密度模型求取。?对时间相关函数R(τ)应用Wiener?Kolmogorov白化程序后，机动加速度可用输入为白噪声的一阶时间相关模型来表征，即?式中u(t)表示均值为零、方差为1的高斯白噪声；ω1(t)表示均值为零、方差为的高斯白噪声。?步骤8：确定量测方程：?设状态向量X表示为?X＝[x?vx?αx?y?vy?αy]T????(3)?式中x、vx、αx和y、vy、α...

【技术特征摘要】
1.一种基于径向加速度的RA-Singer-EKF机动目标跟踪方法，是指在信号处理阶段采用正交匹配追踪(OMP)思想对目标径向加速度和径向速度进行提取，在数据处理阶段将径向加速度和径向速度估计值通过坐标转换引入量测方程和和状态方程中，并采用Singer模型和扩展卡尔曼滤波(EKF)方法实现机动目标的精确跟踪，其特征在于包括以下步骤步骤1:将雷达接收机接收到的线性调频信号s (t)通过采样器以采样间隔Ts进行采样，变为离散信号s (nTs)，其中η表示采样点序号；将s(nTs)送入雷达信号处理计算机；在雷达信号处理计算机中执行以下步骤 步骤2 :初始化(设置分解参数) T设为雷达的脉冲宽度； λ为雷达波长； fs为采样频率； Ts = l/fs为采样间隔； f...

【专利技术属性】
技术研发人员：王国宏，贾舒宜，张磊，谭顺成，于洪波，
申请(专利权)人：中国人民解放军海军航空工程学院，
类型：发明
国别省市：