发明名称 一种宽带雷达目标离散特征状态估计方法
摘要 本发明公开了一种宽带雷达目标离散特征状态估计方法,属于雷达跟踪领域。本方法包括如下步骤:1、M个目标具有n<sub>t</sub>个观测变量,令第r个观测<img file="DDA00002828198900011.GIF" wi="111" he="66" />来源于第m个目标;第m个目标的离散特征状态变量为<img file="DDA00002828198900012.GIF" wi="292" he="66" />具有N个取值;初始化目标m在t=1时刻的离散特征状态为第j种的代价函数;2、在2≤t≤T区间,m在t时刻离散特征状态为j的最大代价函数;将取得最大代价函数δ′<sub>t,m</sub>(j)值时的u记录为最优状态路径y<sub>t,m</sub>(j);3、提取T时刻代价函数值最大的离散特征状态变量;4、由T经由最优状态路径回溯,提取1≤t≤T-1的离散特征状态序列;则1≤t≤T的离散特征状态序列作为第m个目标的离散特征状态估计值。本发明适用于在宽带雷达关联、滤波一体化框架下对目标离散特征状态估计。
申请公布号 CN103176177A 申请公布日期 2013.06.26
申请号 CN201310048877.3 申请日期 2013.02.06
申请人 北京理工大学 发明人 刘泉华;陈新亮;郑乐;曾涛;李阳;龙腾;毛二可
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 付雷杰;高燕燕
主权项 1.一种宽带雷达目标离散特征状态估计方法,其特征在于,包括如下步骤:步骤1、M个目标在观测时间段[1,T]内具有n<sub>t</sub>个观测变量,已知时刻t中第r个观测<img file="FDA00002828198600011.GIF" wi="111" he="65" />来源于第m个目标的关联权重为w<sub>m,r</sub>(t),1≤r≤nt;第m个目标的离散特征状态变量为<img file="FDA00002828198600012.GIF" wi="139" he="66" />假设第m个目标具有N种离散特征状态,将N种离散特征状态由1~N进行编号,1≤m≤M;使用编号j代表第j种离散特征状态,1≤j≤N;初始化目标m在t=1时刻的离散特征状态为第j种离散特征状态时的代价函数δ(<sub>t=1</sub>)<sub>,m</sub>(j),其中j取遍区间[1,N]的所有整数:<img file="FDA00002828198600013.GIF" wi="1031" he="125" /><img file="FDA00002828198600014.GIF" wi="123" he="73" />为第m个目标的离散特征观测概率:<img file="FDA00002828198600015.GIF" wi="591" he="83" />步骤2、对于2≤t≤T,计算目标m在t时刻离散特征状态j=1时的最大代价函数δ<sub>t,m</sub>(j=1)值:<img file="FDA00002828198600016.GIF" wi="1497" he="145" />其中,将t时刻的各离散特征状态用j表示,将t-1时刻的各离散特征状态用u表示,u的取值范围与j相同,1≤u≤N;上述公式表示令u取遍1~N,计算t-1时刻的每个状态u转移到t时刻第1个离散特征状态的代价函数,并取最大值;<img file="FDA00002828198600017.GIF" wi="71" he="64" />为第m个目标由t-1时刻的状态u转移至t时刻的状态j的特征状态转移概率;获取最大代价函数δ<sub>t,m</sub>(j=1)对应的u值,记为u<sub>max</sub>,将u<sub>max</sub>赋值给y<sub>t,m</sub>(j=1),y<sub>t,m</sub>(j=1)=u<sub>max</sub>表示对于第m个目标来说t-1时刻的第u<sub>max</sub>个离散特征状态是转移到t时刻第1个离散特征状态的最优状态路径;依次计算目标m在t时刻的离散特征状态j=2~N时的最大代价函数值δ<sub>t,m</sub>(j=2)~δ<sub>t,m</sub>(j=N),并记录对应的最优状态路径y<sub>t,m</sub>(j=2)~y<sub>t,m</sub>(j=N);针对2≤t≤T范围内的每个时刻,计算δ<sub>t,m</sub>(j=1)~δ<sub>t,m</sub>(j=N)以及对应的y<sub>t,m</sub>(j=1)~y<sub>t,m</sub>(j=N);步骤3、在最终时刻T的N个最大代价函数值δ<sub>T,m</sub>(j=1)~δ<sub>T,m</sub>(j=N)中找到一个最大值,获取所述最大值对应的j,记为j<sub>max</sub>,将j<sub>max</sub>赋值给x'<sub>m</sub>(T);x'<sub>m</sub>(T)为T时刻最有可能的离散特征状态;步骤4、将最优状态路径y<sub>T,m</sub>(j=x′<sub>m</sub>(T))作为T-1时刻最有可能转移到T时刻第x′<sub>m</sub>(T)个离散特征状态的状态,并赋值给x′<sub>m</sub>(T-1);再将最优状态路径y<sub>T-1,m</sub>(j=x′<sub>m</sub>(T-1))作为T-2时刻最有可能转移到T-1时刻第x′<sub>m</sub>(T-1)个离散特征状态的状态,并赋值给x′<sub>m</sub>(T-2),依次类推,获得所有时刻的x′<sub>m</sub>(t),1≤t≤T,将x′<sub>m</sub>(1)~x′<sub>m</sub>(T)作为第m个目标的离散特征状态估计值。
地址 100081 北京市海淀区中关村南大街5号