发明名称 一种基于G0杂波背景及恒定目标幅度的检测前跟踪方法
摘要 本发明提供一种基于G0杂波背景及恒定目标幅度的检测前跟踪方法,针对恒定幅度的目标,将杂波背景建模为G0分布,由于量测值幅度似然比并不存在解析解,本发明采用加权gamma函数和对其进行近似,因此该方法适用于G0杂波背景及恒定目标幅度情况。本发明采用加权不完全gamma函数和作为值函数,由于采用的加权不完全gamma函数和量测值幅度似然比利用了杂波参数和目标参数,所以该方法较之现有动态规划检测前跟踪算法,充分利用了背景与杂波信息,能更好地体现目标与杂波的差异。
申请公布号 CN104062651B 申请公布日期 2016.06.29
申请号 CN201410305868.2 申请日期 2014.06.30
申请人 电子科技大学 发明人 孔令讲;姜海超;刘睿;李小龙;易伟;崔国龙;陈建;杨晓波
分类号 G01S13/66(2006.01)I;G06F19/00(2011.01)I 主分类号 G01S13/66(2006.01)I
代理机构 电子科技大学专利中心 51203 代理人 邹裕蓉
主权项 一种基于G0杂波背景及恒定目标幅度的检测前跟踪方法,其特征在于,包括以下步骤:步骤1、初始化系统参数:G0杂波的形状参数α、G0杂波的尺度参数β和目标回波幅度A,动态规划检测前跟踪算法处理帧数K,初始化变量k=1;步骤2、从雷达接收机中读取第k帧回波数据的量测值Z<sub>k</sub>,<img file="FDA0000873471790000011.GIF" wi="668" he="95" /><img file="FDA0000873471790000012.GIF" wi="69" he="70" />为第k帧回波数据在量测单元(i,j)中的幅度值,N<sub>x</sub>为X轴量化的单元总数,N<sub>y</sub>为Y轴量化单元的总数;步骤3、利用G0杂波的形状参数α和尺度参数β以及目标回波幅度A,对每个量测单元(i,j)计算加权不完全gamma函数和s<sub>ij</sub>:s<sub>ij</sub>=u<sub>1</sub>v<sub>1</sub>+u<sub>2</sub>v<sub>2</sub>+u<sub>3</sub>v<sub>3</sub>其中,u<sub>1</sub>、u<sub>2</sub>、u<sub>3</sub>均为不完全gamma函数,u<sub>1</sub>=γ(p<sub>1</sub>,q<sub>1</sub>u),u<sub>2</sub>=γ(p<sub>2</sub>,q<sub>2</sub>u);u<sub>3</sub>=γ(p<sub>3</sub>,q<sub>3</sub>u);v<sub>1</sub>、v<sub>2</sub>、v<sub>3</sub>均为权重系数,<img file="FDA0000873471790000013.GIF" wi="1197" he="181" />函数因子p<sub>1</sub>、p<sub>2</sub>、p<sub>3</sub>、q<sub>1</sub>、q<sub>2</sub>、q<sub>3</sub>、u分别为p<sub>1</sub>=α+1、p<sub>2</sub>=α+3、p<sub>3</sub>=α+0.5、<img file="FDA0000873471790000014.GIF" wi="527" he="103" /><img file="FDA0000873471790000015.GIF" wi="725" he="142" />不完全gamma函数<img file="FDA0000873471790000016.GIF" wi="652" he="110" />Γ(a)=Γ(a,b)+γ(a,b),<img file="FDA0000873471790000017.GIF" wi="429" he="102" />步骤4、利用加权不完全gamma函数和s<sub>ij</sub>对每个量测单元(i,j)计算量测值幅度似然比<img file="FDA0000873471790000018.GIF" wi="205" he="71" /><img file="FDA0000873471790000019.GIF" wi="605" he="212" />其中,|·|表示取模;步骤5、值函数累积:当k=1,则第1帧中每个状态<img file="FDA00008734717900000110.GIF" wi="61" he="77" />的值函数<img file="FDA00008734717900000111.GIF" wi="134" he="95" />为该状态对应的量测单元(i,j)的量测值幅度似然比<img file="FDA00008734717900000112.GIF" wi="622" he="95" />一帧中一个量测单元(i,j)一一对应一个状态, f=1,…,N<sub>x</sub>N<sub>y</sub>;如果k≠1,则第k帧中的每个状态<img file="FDA0000873471790000021.GIF" wi="69" he="78" />的值函数<img file="FDA0000873471790000022.GIF" wi="138" he="94" />为第k帧第f个状态在第k‑1帧的有效状态转移区域<img file="FDA0000873471790000023.GIF" wi="166" he="86" />内值函数的最大值与第f个状态对应的量测值幅度似然比之和,即<img file="FDA0000873471790000024.GIF" wi="710" he="125" />同时记录第k帧第f个状态<img file="FDA0000873471790000025.GIF" wi="70" he="75" />与其在第k‑1帧的有效状态转移区域<img file="FDA0000873471790000026.GIF" wi="167" he="87" />内值函数最大值对应的状态之间对应关系;步骤6、如果k&lt;K,更新k=k+1返回步骤2;如果k=K,执行步骤7;步骤7、用目标候选状态集合D存储第K帧中超过门限V<sub>T</sub>的值函数所对应的状态;步骤8、航迹恢复:如果目标候选状态集合D为空,则宣布没有目标,算法结束;如果目标候选状态集合D不为空,对于集合D中每个元素作为候选目标状态,进行航迹回溯,根据记录的第k帧第f个状态<img file="FDA0000873471790000027.GIF" wi="61" he="77" />与其在第k‑1帧的有效状态转移区域<img file="FDA0000873471790000028.GIF" wi="166" he="87" />内值函数最大值对应的状态之间对应关系,得到集合D中每个候选目标状态对应恢复的航迹;步骤9、经虚假航迹删除后输出检测结果与目标航迹。
地址 611731 四川省成都市高新区(西区)西源大道2006号