申请/专利权人：哈尔滨工业大学

申请日：2018-05-17

公开（公告）日：2022-06-21

公开（公告）号：CN108919253B

主分类号：G01S13/72

分类号：G01S13/72

优先权：

专利状态码：有效-授权

法律状态：2022.06.21#授权;2018.12.25#实质审查的生效;2018.11.30#公开

摘要：一种基于帝国竞争算法的天波雷达多路径数据关联方法，本发明涉及一种基于帝国竞争算法的天波雷达多路径数据关联方法。本发明的目的是为了解决现有MPDA算法存在计算量大、算法执行时间过长的缺点。过程为：一，建立天波雷达跟踪系统模型；二，初始化，计算滤波预测值；三，随机生成初始国家并计算权力；四，初始帝国；五，殖民地同化作用；六，殖民地革命；七，帝国增强；八，殖民竞争；九，若迭代结束或只剩一个帝国，输出结果，否则，返回五；十，计算目标状态估计值，十一，判断目标跟踪是否结束，如果没有，使k＝k+1，执行二～十。本发明用于雷达数据处理领域。

主权项：1.一种基于帝国竞争算法的天波雷达多路径数据关联方法，其特征在于：步骤一，建立天波雷达运动目标的运动模型和多路径量测模型，将天波雷达运动目标地面坐标系下的运动状态标记为L维的运动状态向量x，天波雷达运动目标雷达坐标系下的量测值标记为M维的量测向量y，根据运动模型获得天波雷达状态误差协方差矩阵P和天波雷达过程噪声协方差矩阵Q，根据量测模型获得量测协方差矩阵R；步骤二，根据步骤一获得的第k时刻天波雷达运动目标在地面坐标系下的运动状态向量xk、天波雷达状态误差协方差矩阵Pk，计算天波雷达运动目标第k+1时刻的运动状态预测值和状态预测协方差矩阵Pk+1|k，并由此得到天波雷达运动目标第k+1时刻传播模式i＝1,...,nk+1在雷达坐标系下的量测预测值和量测预测协方差Sik+1；nk是时刻k传播模式的个数；步骤三，传感器获取当前时刻天波雷达运动目标在雷达坐标系下的量测值，随机生成N个初始量测-模式关联对，即初始国家，并计算每个国家的权力ct，t＝1，2，3，...，N，初始化当前迭代次数g＝1，预设算法迭代次数G；N取值为正整数；步骤四，依据国家权力大小从高到低排序，选择前Nimp个作为殖民国家，其余Ncol个作为殖民地，按殖民国家的权力大小将殖民地分配给初始殖民国家，形成初始帝国；步骤五，殖民地同化作用：殖民国家以自身量测-模式关联对为基础，依据一定同化率随机地替换和改变殖民地路径中的α个量测-模式关联对；步骤六，殖民地革命：殖民地按革命概率Pm随机地选取殖民地两个位置的关联量测进行交换，然后计算殖民地革命前后的权力，保留权力大的殖民国家，并重新评估殖民国家和殖民地的权力，若后者大于前者，则取代前者成为新的帝国；步骤七，帝国增强，根据每个殖民国家权力，对权力最小的国家进行增强，增强方式为：随机移除该殖民国家中β个关联量测，重新建立量测-模式关联对，重新计算殖民国家增强前后的权力，保留权力大的殖民国家；步骤八，殖民竞争，以帝国为单位，计算帝国内殖民国家和所有殖民地的标准化权力TCn；通过评估所有帝国权力，将最弱帝国的最弱殖民地分配给最强帝国；当某个帝国内没有任何殖民地时，帝国消失，该帝国的殖民国家也被分配给最强帝国作为殖民地；步骤九，若达到预设算法的迭代次数G或只剩一个帝国，算法结束，输出最终的量测-模式关联对作为数据关联的结果，否则，令g＝g+1，返回步骤五；步骤十，根据步骤九输出的关联结果计算当前时刻地面坐标系下的目标状态估计值；步骤十一，判断目标跟踪是否结束，如果没有结束，使k＝k+1，执行步骤二～步骤十，如果结束则完成数据关联方法。

全文数据：一种基于帝国竞争算法的天波雷达多路径数据关联方法技术领域[0001]本发明涉及一种基于帝国竞争算法的天波雷达多路径数据关联方法，属于雷达数据处理领域。背景技术[0002]数据关联模块是雷达数据处理系统的关键组成部分，数据关联的作用是在杂波中将量测数据与真实目标航迹进行关联匹配，从而实现对目标的精准跟踪。数据关联算法性能的优劣直接影响到跟踪系统整体的性能。[0003]与常规雷达不同，天波雷达面临的环境更为复杂，且存在着低检测概率、低数据率、低测量精度、高虚警率和多径效应等问题，使其数据处理变得非常困难。其中，多径效应的存在降低了单一传播模式下的检测概率，提高了航迹失跟率，是天波雷达迫切需要解决的问题。由此可见，传统的数据关联算法已经不能满足天波雷达迅速建立航迹形成早期预警的要求，研究、改进和发展新的天波超视距雷达数据关联算法是十分迫切和必要的。[0004]随着电离层探测技术的进步，学者开始探讨利用电离层信息解决多径传播问题，提出了多路径概率数据关联算法MultipathProbabilityDataAssociation，MPDA等算法。其中，MPDA算法将各种传播模式的回波作为备选的竞争量测，在最大似然意义下选择最优的回波，进行传播模式与目标的关联以及状态更新，是应用最为广泛的一种算法，但由于该算法需要穷举所有的关联事件，计算量会随着杂波密度的增加而指数增加，使算法成为NP难题，存在计算量大、算法执行时间过长的缺点。发明内容[0005]本发明的目的是为了解决现有MPDA算法存在计算量大、算法执行时间过长的缺点，不适用于实际工程，从而提出一种基于帝国竞争算法的天波雷达多路径数据关联方法。[0006]—种基于帝国竞争算法的天波雷达多路径数据关联方法过程为：[0007]步骤一，建立天波雷达运动目标的运动模型和多路径量测模型，将天波雷达运动目标地面坐标系下的运动状态标记为L维的运动状态向量X，天波雷达运动目标雷达坐标系下的量测值标记为M维的量测向量y，根据运动模型获得天波雷达状态误差协方差矩阵P和天波雷达过程噪声协方差矩阵Q，根据量测模型获得量测协方差矩阵R;[0008]步骤二，根据步骤一获得的第k时刻天波雷达运动目标在地面坐标系下的运动状态向量X⑹、天波雷达状态误差协方差矩阵P⑹，计算天波雷达运动目标第k+Ι时刻的运动状态预测值+11幻和状态预测协方差矩阵Pk+1Ik，并由此得到天波雷达运动目标第k+1时刻传播模式i=l，...，1^+1在雷达坐标系下的量测预测值f,伏+1|幻和量测预测协方差Sik+1;nk是时亥Ijk传播模式的个数；[0009]步骤三，传感器获取当前时刻天波雷达运动目标在雷达坐标系下的量测值，随机生成N个初始量测-模式关联对，即初始国家，并计算每个国家的权力Ct，t=1，2，3，...，N，初始化当前迭代次数g=1，预设算法迭代次数G;N取值为正整数；[0010]步骤四，依据国家权力大小从高到低排序，选择前化_个作为殖民国家，其余Ν™ι个作为殖民地，按殖民国家的权力大小将殖民地分配给初始殖民国家，形成初始帝国；[0011]步骤五，殖民地同化作用:殖民国家以自身量测-模式关联对步骤三得到的）为基础，依据一定同化率随机地替换和改变殖民地路径中的α个量测-模式关联对；[0012]步骤六，殖民地革命:殖民地按革命概率?„随机地选取殖民地两个位置的关联量测进行交换，然后计算殖民地革命前后的权力，保留权力大的殖民国家，并重新评估殖民国家和殖民地的权力，若后者大于前者，则取代前者成为新的帝国；[0013]步骤七，帝国增强，根据每个殖民国家权力，对权力最小的国家进行增强，增强方式为:随机移除该殖民国家中β个关联量测，重新建立量测-模式关联对，重新计算殖民国家增强前后的权力，保留权力大的殖民国家；[0014]步骤八，殖民竞争，以帝国为单位，计算帝国内殖民国家和所有殖民地的标准化权力TCn;通过评估所有帝国权力，将最弱帝国的最弱殖民地分配给最强帝国；当某个帝国内没有任何殖民地时，帝国消失，该帝国的殖民国家也被分配给最强帝国作为殖民地；[0015]步骤九，帝国竞争算法不断迭代运算，若达到预设算法的迭代次数G或只剩一个帝国，算法结束，输出最终的量测-模式关联对作为数据关联的结果，否则，令g=g+l，返回步骤五；[0016]步骤十，根据步骤九输出的关联结果计算当前时刻地面坐标系下的目标状态估计值；[0017]步骤十一，判断目标跟踪是否结束，如果没有结束，使k=k+l，执行步骤二〜步骤十，如果结束则完成数据关联方法。[0018]本发明的有益效果为：[0019]本发明将天波雷达的多路径数据关联问题建模为组合优化问题，并将帝国竞争算法引入其中，为该领域提供了一个基于社会群智能算法的新思路，在保证算法跟踪精度的前提下大大缩短了算法执行时间和计算量;如图8到图11给出50次蒙特卡洛仿真后目标各维度目标各维度的均方根误差，从图中可以看出，随着跟踪的增加，目标径向距离、径向距离变化率、方位角、方位角变化率的均方根误差均逐渐减少并趋近于〇,跟踪精度提高10倍左右，同样证明本发明可以有效解决天波雷达多路径数据关联问题，达到了收敛的效果。附图说明[0020]图1为基于帝国竞争算法的天波雷达多路径数据关联方法流程图；[0021]图2为天波雷达运动目标在雷达坐标系下的径向距量测图；[0022]图3为天波雷达运动目标在雷达坐标系下的多普勒量测图；[0023]图4为天波雷达运动目标在雷达坐标系下的方位角量测图；[0024]图5为天波雷达运动目标在地面坐标系下的目标跟踪航迹图；[0025]图6为天波雷达运动目标在地面坐标系下的目标径向距离变化率跟踪图；[0026]图7为天波雷达运动目标在地面坐标系下的目标方位角变化率跟踪图；[0027]图8为天波雷达运动目标的径向距离均方根误差随时间变化曲线图；[0028]图9为天波雷达运动目标的径向距离变化率均方根误差随时间变化曲线图；[0029]图10为天波雷达运动目标的方位角均方根误差随时间变化曲线图；[0030]图11为天波雷达运动目标的方位角变化率均方根误差随时间变化曲线图。具体实施方式[0031]具体实施方式一:结合图1说明本实施方式，本实施方式的一种基于帝国竞争算法的天波雷达多路径数据关联方法具体过程为：[0032]步骤一，建立天波雷达运动目标的运动模型和多路径量测模型（多路径量测模型就是多径效应下天波雷达的系统量测模型，包括量测向量和量测方程y⑹），将天波雷达运动目标地面坐标系下的运动状态标记为L维的运动状态向量X，天波雷达运动目标雷达坐标系下的量测值标记为M维的量测向量y，根据运动模型包括运动状态向量和运动状态方程）获得天波雷达状态误差协方差矩阵P和天波雷达过程噪声协方差矩阵Q，根据量测模型获得量测协方差矩阵R;[0033]步骤二，根据步骤一获得的第k时刻天波雷达运动目标在地面坐标系下的运动状态向量X⑹、天波雷达状态误差协方差矩阵P⑹，计算天波雷达运动目标第k+Ι时刻的运动状态预测值+II幻和状态预测协方差矩阵Pk+1Ik，并由此得到天波雷达运动目标第k+1时刻传播模式i=l，...，nk+1在雷达坐标系下的量测预测值尤供+ΐμ和量测预测协方差Sik+1;nk是时亥Ijk传播模式的个数；[0034]步骤三，传感器获取当前时刻天波雷达运动目标在雷达坐标系下的量测值，随机生成N个初始量测-模式关联对，即初始国家，并计算每个国家的权力Ct，t=1，2，3，...，N，初始化当前迭代次数g=1，预设算法迭代次数G;N取值为正整数；[0035]步骤四，依据国家权力大小从高到低排序，选择前化_个作为殖民国家，其余Ν™ι个作为殖民地，按殖民国家的权力大小将殖民地分配给初始殖民国家，形成初始帝国；[0036]步骤五，殖民地同化作用:殖民国家以自身量测-模式关联对步骤三得到的）为基础，依据一定同化率随机地替换和改变殖民地路径中的α个量测-模式关联对；[0037]步骤六，殖民地革命:殖民地按革命概率?„随机地选取殖民地两个位置的关联量测进行交换，然后计算殖民地革命前后的权力，保留权力大的殖民国家，并重新评估殖民国家和殖民地的权力，若后者大于前者，则取代前者成为新的帝国；[0038]步骤七，帝国增强，根据每个殖民国家权力，对权力最小的国家进行增强，增强方式为:随机移除该殖民国家中β个关联量测，重新建立量测-模式关联对，重新计算殖民国家增强前后的权力，保留权力大的殖民国家；[0039]步骤八，殖民竞争，以帝国为单位，计算帝国内殖民国家和所有殖民地的标准化权力TCn;通过评估所有帝国权力，将最弱帝国的最弱殖民地分配给最强帝国；当某个帝国内没有任何殖民地时，帝国消失，该帝国的殖民国家也被分配给最强帝国作为殖民地；[0040]步骤九，帝国竞争算法不断迭代运算，若达到预设算法的迭代次数G或只剩一个帝国，算法结束，输出最终的量测-模式关联对作为数据关联的结果，否则，令g=g+l，返回步骤五；[0041]步骤十，根据步骤九输出的关联结果计算当前时刻地面坐标系下的目标状态估计值；[0042]步骤十一，判断目标跟踪是否结束，如果没有结束，使k=k+l，执行步骤二〜步骤十，如果结束则完成数据关联方法。[0043]具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤一中建立天波雷达运动目标的运动模型和多路径量测模型（多路径量测模型就是多径效应下天波雷达的系统量测模型，包括量测向量和量测方程y⑹），将天波雷达运动目标地面坐标系下的运动状态标记为L维的运动状态向量X，天波雷达运动目标雷达坐标系下的量测值标记为M维的量测向量y，根据运动模型（包括运动状态向量和运动状态方程获得天波雷达状态误差协方差矩阵P和天波雷达过程噪声协方差矩阵Q，根据量测模型获得量测协方差矩阵R;具体过程为：[0044]在采样时刻k，一个勾速天波雷达运动目标在地面坐标系下的运动状态向量表示为：[0046]其中，Pk表示天波雷达运动目标在地面坐标系下的径向距离，PA表示天波雷达运动目标在地面坐标系下的径向距离变化率径向速度），b⑹表示天波雷达运动目标在地面坐标系下的方位角，翁幻表示天波雷达运动目标在地面坐标系下的方位角变化率，L=4,分别为的转置；[0047]则天波雷达运动目标运动状态方程表示为：[0048]Xk+1=FkXk+vk[0049]其中，FGO是LXL阶的状态转移矩阵，L维过程噪声V⑹是一个零均值的高斯白噪声，V⑹协方差矩阵为Q⑹；[0050]在采样时刻k，天波雷达运动目标一次量测值在雷达坐标系下的量测向量表示为：[0051]y⑹=[Rg⑹，Rr⑹，Az⑹]'[0052]式中，心为天波雷达运动目标在在雷达坐标系下的径向距，Rr为天波雷达运动目标在在雷达坐标系下的多普勒，Az为天波雷达运动目标在在雷达坐标系下的面方位角，M=3，分别为Rg、Rr、Az;[Rg⑹，Rrk，AZk]为求[Rg⑹，Rrk，AZ⑹]的转置；[0053]则考虑多径效应多径效应指天波雷达的多路径传播效应，多径效应导致了有多个传播模式，两者在概念上等价下天波雷达的系统量测模型表示为：[0054]y⑹=HiX⑹）+Wi⑹，i=1，···，nk[0055]其中，H1X是非线性量测方程;Wlk是零均值的高斯白噪声，量测协方差矩阵为Ri⑹；nk是时亥Uk传播模式的个数；[0056]天波雷达过程噪声协方差矩阵Q为：[0057]Qk=E[VkVk’][0058]式中，E[v⑹V⑹’]为求V⑹V⑹’的期望；[0059]量测协方差矩阵R为：[0060]Rik=E[WikWik’][0061]其它步骤及参数与具体实施方式一相同。[0062]具体实施方式三:本实施方式与具体实施方式一或二不同的是:所述步骤二中根据步骤一获得的第k时刻天波雷达运动目标在地面坐标系下的运动状态向量xk、天波雷达状态误差协方差矩阵Pk，计算天波雷达运动目标第k+Ι时刻的运动状态预测值啦+1μ和状态预测协方差矩阵pk+iIk，并由此得到天波雷达运动目标第k+i时刻传播模式i=l，...，nk+1在雷达坐标系下的量测预测值和量测预测协方差S1k+Ι;具体为：[0063]当k=0时，采用传感器获得天波雷达运动目标在雷达坐标系下的测量值向量，初始化天波雷达运动目标在地面坐标系下的运动状态向量X〇和天波雷达状态误差协方差矩阵P⑼；[0064]当k0时，用k时刻的滤波估计值和滤波协方差Pk|k初始化天波雷达运动目标在地面坐标系下的运动状态向量Xk和天波雷达状态误差协方差矩阵Pk，汁算第k+Ι时刻的运动状态预测值和状态预测协方差矩阵Pk+1Ik，并由此得到天波雷达运动目标第k+Ι时刻传播模式i=I，...，nk+1的量测预测值：和量测预测协方差S1k+1;计算方法具体为：[0069]其中，J1k+1为量测函数H1X在当前估计值+lI:下的雅可比矩阵，J1k+1’为心k+1的转置，nk+1是时刻k+1传播模式的个数。[0070]其它步骤及参数与具体实施方式一或二相同。[0071]具体实施方式四：本实施方式与具体实施方式一至三之一不同的是:所述步骤三中传感器获取当前时刻天波雷达运动目标在雷达坐标系下的量测值，随机生成N个初始量测-模式关联对，即初始国家，并计算每个国家的权力ct，t=l，2,3,...，N，初始化当前迭代次数g=1，预设算法迭代次数G;N取值为正整数；[0072]具体过程为：[0073]随机形成N个初始国家C={Ci，C2,...，Cn}，并计算每个国家的权力C={ci，C2,...，cN}，将国家权力c与似然函数d建立相关性cn=lDn;[0074]其中Cn，Dn分别表示第η个国家的权力和对应的似然函数，n=l，2,3，...，N;[0075]似然函数的计算公式为、[0076]式中，dzl表示第z个量测与第i个传播模式关联的似然函数，yzCO表示第k时刻第z个量测，表示第k时刻第i个传播模式的量测预测值，S1k表示第k时刻第i个传播模式的新息协方差，为的转置。[0077]在本发明中，每个国家代表一组量测-传播模式关联对。假设有三条传播模式A、B、C，另有三个量测a、b、c，则初始化形成的量测-传播模式关联对表示为A，c-B，a-C，b，则该国家的权力c=ID，其中，D=dcA+daB+dbc;[0078]cU表示第c个量测与第A个传播模式关联的似然函数，daB表示第a个量测与第B个传播模式关联的似然函数，dbC表示第b个量测与第C个传播模式关联的似然函数。[0079]其它步骤及参数与具体实施方式一至三之一相同。[0080]具体实施方式五:本实施方式与具体实施方式一至四之一不同的是:所述步骤四中依据国家权力大小从高到低排序，选择前NimA作为殖民国家，其余PW个作为殖民地，按殖民国家的权力大小将殖民地分配给初始殖民国家，形成初始帝国；具体为：[0081]根据国家权力大小从高到低排序，选择前化_个作为殖民国家，殖民国家的权力定义为c_impf，f=1，2，…,Nimp;[0082]其余Ncoi个作为殖民地，殖民地的权力定义Sc_colj，j=l，2，···,Ncoi;[0083]然后对殖民国家权力进行标准化：[0085]式中，c_imph为第h个殖民国家的权力，h=1，2，…，Nimp;[0086]根据公SN_impf=IOUncHpjmpfXU将所有殖民地分配给殖民国家，第f个殖民国家分个殖民地；[0087]其中，round{p_impfXNc〇i}为对p_impfXNcoi取整;N_impf为殖民地个数。[0088]其它步骤及参数与具体实施方式一至四之一相同。[0089]具体实施方式六:本实施方式与具体实施方式一至五之一不同的是:所述步骤五中殖民地同化作用:殖民国家以自身量测-模式关联对步骤三得到的为基础，依据一定同化率随机地替换和改变殖民地路径中的α个量测-模式关联对;具体过程为：[0090]步骤五一、定义ρ为同化率，并满足00时[0168]计算第k+Ι时刻的运动状态预测值和状态预测协方差Pk+1Ik[0171]第k+1时刻传播模式i=l，2,3,4的量测预测值和量测预测协方差S1Qc+D[0174]式中，J1k+Ι量测函数出X在当前估计值iAr+11幻下的雅可比矩阵[0176]忽略时间k可以表示为：[0178]式中，[0182]然后，在步骤三，产生初始帝国，随机形成N个初始国家C={，，...，CN}，在本实施例中，在本发明中，每个国家代表一组量测-传播模式关联对。假设有四个量测a、b、c，d，则初始化形成的量测-传播模式关联对可表示为EE，c-EF，a-FE，b-FF，d，即一个国家。其中，（EE，c表示量测c与传播模式EE相关联，计算每个国家的权力C={ci，C2,...，cn}，cn=1Dn，其中Dn表示第η个国家对应的似然函数之和：[0184]式中，表不第j个量测与第q个传播模式关联的似然函数，初始化当前迭代次数g，预设算法迭代次数G。[0185]之后，在步骤四，依据国家权力大小从高到低排序，选择Nimpf作为殖民国家，它们的权力定义为c_impi，i=1，2，…,Nimp其余NcxjI个作为殖民地，它们的权力定义为C_C〇lj，j=I，2，…，Ν™ι;然后对殖民国家权力进行标准化：L0187J根据公SN_impi=IOuncUpjmpiXU将所有殖民地分配给殖民国家，第η个殖民国家分得1；[11^个殖民地；[0188]其中，round{p_impiXNc〇i}为对p_impiXNcoi取整;N_impi为殖民地个数；[0189]然后，在步骤五，发生殖民地同化作用，定义同化率P=O.5,并在殖民地每个传播模式的位置生成一个[0，1]之间的随机数，若随机数大于等于P，则发生同化，殖民国家中与该传播模式关联的量测作为同化后殖民地相同传播模式的关联量测，若随机数小于P，且在殖民地与这些传播模式的关联量测在同化时未出现过，则视作被保持，保持方式为殖民地在这些位置处的量测-传播模式关联对直接作为同化后殖民地的量测-传播模式关联对，若随机数小于P，但殖民地与这些传播模式的关联量测已在同化时出现过，则此位置也发生同化，将殖民国家在这些位置处的关联量测替换此位置的关联量测，最后，将同化后的殖民地替换原殖民地。[0190]然后，在步骤六，发生殖民地革命，定义pm为革命概率：[0192]其中，p」mpi为殖民地所属殖民国家的标准化权力，p」mpmax所有殖民国家中最大的标准化权力，g为当前运算迭代数，G为总共迭代次数，Po为初始革命概率，定义Po=O.3,每个殖民地生成一个[0，1]之间的随机数，若该随机数大于pm，则该殖民地发生革命，革命方式为随机将殖民地的两个传播模式的关联量测交换，每个殖民地完成革命操作后，重新评估帝国内的殖民国家和殖民地权力并进行排序，若殖民地的最高权力大于殖民国家，则该殖民地和殖民国家交换身份。[0193]接着，在步骤七，进行帝国增强，根据殖民国家权力表，对权力表中权力最小的殖民国家进行增强操作，增强方式为随机移除m个关联量测，从剩余量测集中选取相同数量的量测插入移除的位置，重新计算该殖民国家的权力，若替换后的权力大于原来的权力，则增强成功，否则，保留原殖民国家。[0194]然后，在步骤八，发生市国克争，克争过程为：[0195]1计算帝国权力，帝国的权力由组成该帝国的殖民国家的权力与殖民地平均权力加权和构成：[0197]其中，TCf为第f个帝国的权力，c_C〇lj为第f个帝国第j个殖民地的权力，N_impf为当前该帝国殖民地个数为权系数;c_impf为第f个殖民国家的权力，f=1，2，…，Nimp;[0198]⑵对帝国权力进行标准化，第h个帝国的标准化权力TPh表示为：[0201]且满足’建立帝国权力向量TP:[0202]式中，NTCh、NTCf为中间变量;TCh为第h个帝国的权力，h=1，2，…，Nimp;[0203]3建立一个维数与帝国数量相同的，由均匀分布随机数rh组成的向量：[0206]其中，，NIh表示第h个帝国内拥有国家的数量，U为均匀分布，并定义概率向量：[0208]假设dh为D中所有元素的最大值，hN_，则把最弱帝国的最弱殖民地分配给第h个帝国，更新最弱和最强帝国的殖民地数量；[0209]第h个帝国即是最强帝国；[0210]4判断是否存在无殖民地的帝国，若是，则该帝国的殖民国家作为殖民地分配给最强帝国，更新帝国和殖民国家的数量Nimp=Nimp-I。[0211]然后，在步骤九，判断当前帝国数量是否为1或当前算法迭代次数g多G，若是，算法迭代结束，输出最终的量测-模式关联对作为数据关联的结果，若否，返回步骤五，并令g=g+1〇[0212]然后，在步骤十，根据步骤九输出的关联结果计算当前时刻地面坐标系下的目标状态估计值，计算公式为：[0215]式中，y;幻表示第Z个量测与第i个传播模式相关联;i々+lπ+1为第k+1时刻天波雷达运动目标在地面坐标系下的滤波估计值，Pk+11k+1为第k+1时刻天波雷达运动目标在地面坐标系下的滤波协方差。[0216]最后，在步骤十一，返回执行所述步骤二，并使k=k+l，直到目标跟踪结束。[0217]选取均方根误差RMSE作为性能的量度，进行100次蒙特卡罗仿真实验。[0218]仿真结果如图2至晒11所示[0219]图2到图4给出单次仿真的径向距、多普勒和方向角典型量测图，由图可见，对于同一目标，其四种传播模式的量测非常接近，并存在相互干扰的情况，而且每条传播模式的检测概率都很低，这也尽一步增大了数据关联的难度。[0220]图5给出地面坐标系下目标跟踪航迹图，由图3可见，随着跟踪时间的增加，跟踪结果向着真实值逐渐趋于稳定，消除了量测噪声和杂波的影响，达到了收敛的效果，证明了本发明的有效性。[0221]图6图7分别给出地面坐标系下目标径向距离变化率和方位角变化率的跟踪结果图，由图中可见，同样地，随着跟踪时间的增加，跟踪结果均向着真实值逐渐趋于稳定，消除了量测噪声和杂波的影响，达到了收敛的效果，证明了本发明的有效性。[0222]图8到图11给出50次蒙特卡洛仿真后目标各维度目标各维度的均方根误差，从图中可以看出，随着跟踪的增加，目标径向距离、径向距离变化率、方位角、方位角变化率的均方根误差均逐渐减少并趋近于〇,跟踪精度提高10倍左右，同样证明本发明可以有效解决天波雷达多路径数据关联问题，达到了收敛的效果。[0223]本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

权利要求：1.一种基于帝国竞争算法的天波雷达多路径数据关联方法，其特征在于：步骤一，建立天波雷达运动目标的运动模型和多路径量测模型，将天波雷达运动目标地面坐标系下的运动状态标记为L维的运动状态向量X，天波雷达运动目标雷达坐标系下的量测值标记为M维的量测向量y，根据运动模型获得天波雷达状态误差协方差矩阵P和天波雷达过程噪声协方差矩阵Q，根据量测模型获得量测协方差矩阵R;步骤二，根据步骤一获得的第k时刻天波雷达运动目标在地面坐标系下的运动状态向量X⑹、天波雷达状态误差协方差矩阵Pk，计算天波雷达运动目标第k+Ι时刻的运动状态预测值+1i幻和状态预测协方差矩阵Pk+1Ik，并由此得到天波雷达运动目标第k+Ι时刻传播模式i=l，...，取+1在雷达坐标系下的量测预测值tA：+l|幻和量测预测协方差S1k+1;nk是时刻k传播模式的个数；步骤三，传感器获取当前时刻天波雷达运动目标在雷达坐标系下的量测值，随机生成N个初始量测-模式关联对，即初始国家，并计算每个国家的权力ct，t=1，2，3，...，N，初始化当前迭代次数g=1，预设算法迭代次数G;N取值为正整数；步骤四，依据国家权力大小从高到低排序，选择前NimA作为殖民国家，其余Ν™ι个作为殖民地，按殖民国家的权力大小将殖民地分配给初始殖民国家，形成初始帝国；步骤五，殖民地同化作用：殖民国家以自身量测-模式关联对为基础，依据一定同化率随机地替换和改变殖民地路径中的α个量测-模式关联对；步骤六，殖民地革命:殖民地按革命概率匕随机地选取殖民地两个位置的关联量测进行交换，然后计算殖民地革命前后的权力，保留权力大的殖民国家，并重新评估殖民国家和殖民地的权力，若后者大于前者，则取代前者成为新的帝国；步骤七，帝国增强，根据每个殖民国家权力，对权力最小的国家进行增强，增强方式为：随机移除该殖民国家中β个关联量测，重新建立量测-模式关联对，重新计算殖民国家增强前后的权力，保留权力大的殖民国家；步骤八，殖民竞争，以帝国为单位，计算帝国内殖民国家和所有殖民地的标准化权力TCn;通过评估所有帝国权力，将最弱帝国的最弱殖民地分配给最强帝国；当某个帝国内没有任何殖民地时，帝国消失，该帝国的殖民国家也被分配给最强帝国作为殖民地；步骤九，若达到预设算法的迭代次数G或只剩一个帝国，算法结束，输出最终的量测-模式关联对作为数据关联的结果，否则，令g=g+l，返回步骤五；步骤十，根据步骤九输出的关联结果计算当前时刻地面坐标系下的目标状态估计值；步骤十一，判断目标跟踪是否结束，如果没有结束，使k=k+l，执行步骤二〜步骤十，如果结束则完成数据关联方法。2.根据权利要求1所述一种基于帝国竞争算法的天波雷达多路径数据关联方法，其特征在于:所述步骤一中建立天波雷达运动目标的运动模型和多路径量测模型，将天波雷达运动目标地面坐标系下的运动状态标记为L维的运动状态向量X，天波雷达运动目标雷达坐标系下的量测值标记为M维的量测向量y，根据运动模型获得天波雷达状态误差协方差矩阵P和天波雷达过程噪声协方差矩阵Q，根据量测模型获得量测协方差矩阵R;具体过程为：在采样时刻k，天波雷达运动目标在地面坐标系下的运动状态向量表示为：其中，Pk表示天波雷达运动目标在地面坐标系下的径向距离，Hfc表示天波雷达运动目标在地面坐标系下的径向距离变化率，b⑹表示天波雷达运动目标在地面坐标系下的方位角，以幻表示天波雷达运动目标在地面坐标系下的方位角变化率，L=4,分别为Pk、的转置；则天波雷达运动目标运动状态方程表示为：其中，F⑹是LXL阶的状态转移矩阵，L维过程噪声V⑹是一个零均值的高斯白噪声，Vk协方差矩阵为Qk;在采样时刻k，天波雷达运动目标一次量测值在雷达坐标系下的量测向量表示为：式中，Rg为天波雷达运动目标在在雷达坐标系下的径向距，Rr为天波雷达运动目标在在雷达坐标系下的多普勒，Az为天波雷达运动目标在在雷达坐标系下的面方位角，M=3,分别为Rg、Rr、Az;[Rg⑹，Rrk，AZ⑹]为求[Rg⑹，Rrk，AZ⑹]的转置；则考虑多径效应下天波雷达的系统量测模型表示为：其中，H1X是非线性量测方程;Wlk是零均值的高斯白噪声，量测协方差矩阵为R1⑹；nk是时亥Ijk传播模式的个数；天波雷达过程噪声协方差矩阵Q为：式中，E[v⑹V⑹’]为求V⑹V⑹’的期望；量测协方差矩阵R为：3.根据权利要求2所述一种基于帝国竞争算法的天波雷达多路径数据关联方法，其特征在于:所述步骤二中根据步骤一获得的第k时刻天波雷达运动目标在地面坐标系下的运动状态向量X⑹、天波雷达状态误差协方差矩阵P⑹，计算天波雷达运动目标第k+Ι时刻的运动状态预测值+lI幻和状态预测协方差矩阵Pk+11k，并由此得到天波雷达运动目标第k+Ι时刻传播模式i=l，...，nk+1在雷达坐标系下的量测预测值tA+l|幻和量测预测协方差Sik+Ι;具体为：当k=0时，采用传感器获得天波雷达运动目标在雷达坐标系下的测量值向量，初始化天波雷达运动目标在地面坐标系下的运动状态向量X⑼和天波雷达状态误差协方差矩阵PO；当ko时，用k时刻的滤波估计值i到幻和滤波协方差PkIk初始化天波雷达运动目标在地面坐标系下的运动状态向量Xk和天波雷达状态误差协方差矩阵Pk，计算第k+i时刻的运动状态预测值ia+i㈨和状态预测协方差矩阵Pk+11k，并由此得到天波雷达运动目标第k+Ι时刻传播模式i=I，...，nk+1的量测预测值+1μ和量测预测协方差S1k+Ι;计算方法具体为：其中，J1k+1为量测函数H1X在当前估计值维-i-1μ下的雅可比矩阵，J1k+ι’为心k+l的转置，nk+i是时刻k+l传播模式的个数。4.根据权利要求3所述一种基于帝国竞争算法的天波雷达多路径数据关联方法，其特征在于:所述步骤三中传感器获取当前时刻天波雷达运动目标在雷达坐标系下的量测值，随机生成N个初始量测-模式关联对，即初始国家，并计算每个国家的权力Ct，t=1，2，3,...，N，初始化当前迭代次数g=l，预设算法迭代次数G;N取值为正整数；具体过程为：随机形成N个初始国家C={Cl，C2,...，Cn}，并计算每个国家的权力C={C1，C2,...，cn}，将国家权力C与似然函数d建立相关性cn=lDn;其中cn，Dn分别表示第η个国家的权力和对应的似然函数，n=l，2,3,...，N;似然函数的计算公式为式中，dzl表示第z个量测与第i个传播模式关联的似然函数，yzGO表示第k时刻第z个量测，表示第k时刻第i个传播模式的量测预测值，S1QO表示第k时刻第i个传播模式的新息协方差，ΛβI々-1'为λΘ-1的转置。5.根据权利要求4所述一种基于帝国竞争算法的天波雷达多路径数据关联方法，其特征在于:所述步骤四中依据国家权力大小从高到低排序，选择前Nimp个作为殖民国家，其余Ν™ι个作为殖民地，按殖民国家的权力大小将殖民地分配给初始殖民国家，形成初始帝国；具体为：根据国家权力大小从高到低排序，选择前化_个作为殖民国家，殖民国家的权力定义为c_impf，f=1，2，…，Nimp;其余NcxjI个作为殖民地，殖民地的权力定义为3_3〇1j，j=1，2，…，NcxjI;然后对殖民国家权力进行标准化：式中，c_imph为第h个殖民国家的权力，h=1，2，…，Nimp;根据公SN_impf=round{p_impfXNc〇i}将所有殖民地分配给殖民国家，第f个殖民国家分得N_impf个殖民地；其中，round{p_impfXNc〇i}为对p_impf\仏。1取整;1'1_;[1]^为殖民地个数。6.根据权利要求5所述一种基于帝国竞争算法的天波雷达多路径数据关联方法，其特征在于:所述步骤五中殖民地同化作用:殖民国家以自身量测-模式关联对为基础，依据一定同化率随机地替换和改变殖民地路径中的α个量测-模式关联对;具体过程为：步骤五一、定义P为同化率，并满足〇Ρ1，在殖民地每个传播模式的位置生成一个[〇，1]之间的随机数，分别执行步骤五二、步骤五三和步骤五四；步骤五二、随机数大于等于P的殖民地传播模式的位置处发生同化；同化方式为：用殖民国家中与随机数大于等于P的殖民地传播模式的位置处殖民地传播模式关联的量测作为同化后殖民地相同传播模式的关联量测；步骤五三、随机数小于P的传播模式的位置处，若殖民地与随机数小于P的传播模式的关联量测在同化时未出现过，则视作被保持；保持方式为:殖民地在随机数小于P的传播模式的位置处的量测-传播模式关联对直接作为同化后殖民地的量测-传播模式关联对；步骤五四、随机数小于P的传播模式的位置处，若殖民地与随机数小于P的传播模式的关联量测已在同化时出现过，则将殖民国家在随机数小于P的传播模式的位置处的关联量测替换此位置的关联量测；步骤五五、将同化后的殖民地替换原殖民地。7.根据权利要求6所述一种基于帝国竞争算法的天波雷达多路径数据关联方法，其特征在于:所述步骤六中殖民地革命:殖民地按革命概率Pm随机地选取殖民地两个位置的关联量测进行交换，然后计算殖民地革命前后的权力，保留权力大的殖民国家，并重新评估殖民国家和殖民地的权力，若后者大于前者，则取代前者成为新的帝国；具体为：步骤六一、定义Pm为革命概率：其中，p_impf为殖民地所属殖民国家的标准化权力，p_impmax为所有殖民国家中最大的标准化权力，g为当前运算迭代数，G为预设算法迭代次数，Po为初始革命概率，OPo1;步骤六二、每个殖民地生成一个[〇，1]之间的随机数，若该随机数大于pm，则该殖民地发生革命；革命方式为：随机将殖民地的两个传播模式的关联量测交换，直至每个殖民地的量测交换完，即完成革命操作；步骤六三、每个殖民地完成革命操作后，重新评估帝国内的殖民国家和殖民地权力并进行排序，若殖民地的权力大于殖民国家，则该殖民地和殖民国家交换身份。8.根据权利要求7所述一种基于帝国竞争算法的天波雷达多路径数据关联方法，其特征在于:所述步骤七中帝国增强，根据每个殖民国家权力，对权力最小的国家进行增强，增强方式为:随机移除该殖民国家中β个关联量测，重新建立量测-模式关联对，重新计算殖民国家增强前后的权力，保留权力大的殖民国家;具体过程为：根据每个殖民国家权力，对权力中权力最小的殖民国家进行增强操作，增强方式为：随机移除β个关联量测，从剩余量测集中选取β个量测插入移除的位置，重新计算该殖民国家的权力，若替换后的殖民国家的权力大于原来的殖民国家的权力，则增强成功，否则，保留原殖民国家。9.根据权利要求8所述一种基于帝国竞争算法的天波雷达多路径数据关联方法，其特征在于:所述步骤八中殖民竞争，以帝国为单位，计算帝国内殖民国家和所有殖民地的标准化权力1'〇1;通过评估所有帝国权力，将最弱帝国的最弱殖民地分配给最强帝国；当某个帝国内没有任何殖民地时，帝国消失，该帝国的殖民国家也被分配给最强帝国作为殖民地;具体为：步骤八一、计算帝国权力，帝国的权力由组成该帝国的殖民国家的权力与殖民地平均权力加敕的和构成，其中，TCf为第f个帝国的权力，c_C〇lj为第f个帝国第j个殖民地的权力，N_impf为当前该帝国殖民地个数为权系数;c_impf为第f个殖民国家的权力，f=1，2，…，Nimp;步骤八二、对帝国权力进行标准化，第h个帝国的标准化权力TPh表示为：且满月，建立帝国权力向量TP式中，NTCh、NTCf为中间变量;TCh为第h个帝国的权力，h=1，2，…，Nimp;步骤八三、建立一个维数与帝国数量相同的，由均匀分布随机数η组成的向量：其中：，NIh表示第h个帝国内拥有国家的数量，U为均匀分布，并定义概率向量：假设dh为D中所有元素的最大值，hN_，则把最弱帝国的最弱殖民地分配给第h个帝国，更新最弱和最强帝国的殖民地数量；第h个帝国即是最强帝国；步骤八四、判断是否存在无殖民地的帝国，若是，则该帝国的殖民国家作为殖民地分配给最强帝国，更新帝国即殖民国家的数量Nimp=Nimp-1。1〇.根据权利要求9所述一种基于帝国竞争算法的天波雷达多路径数据关联方法，其特征在于:所述步骤十中根据步骤九输出的关联结果计算当前时刻地面坐标系下的目标状态估计值;具体为：按照步骤九输出的量测-模式关联对更新目标的状态：式中，y;幻表示第Z个量测与第i个传播模式相关联；iA：+i|A+i为第k+1时刻天波雷达运动目标在地面坐标系下的滤波估计值，Pk+11k+Ι为第k+Ι时刻天波雷达运动目标在地面坐标系下的滤波协方差。

百度查询： 哈尔滨工业大学 一种基于帝国竞争算法的天波雷达多路径数据关联方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD