基于差分进化狼群算法的欺骗干扰识别方法
【技术领域】
[0001] 本发明属于通信技术领域，具体涉及一种基于差分进化狼群算法的欺骗干扰识别 方法。
【背景技术】
[0002] 卫星导航定位系统易受到多种形式的干扰。其中，最主要的干扰样式为欺骗干扰， 该干扰具有与真实信号相同的增益，并在时、频、空等多域与真实信号完全或者部分重叠， 识别效果往往不太理想。卫星导航系统接收机如何抗欺骗干扰已经成为当前迫切需要解决 的技术问题。
[0003] 抗欺骗干扰的前提是接收机能够检测识别出欺骗干扰。目前国内外针对欺骗干 扰识别方法的研究主要集中在基于特征提取的干扰识别算法，提取的特征参数包括时钟偏 移、多普勒频移和接收机位置等。近年来，对发射机建模并基于模型参数来识别无线发射机 成为了一个新的研究热点。研究表明，基于模型参数的识别方法即使在样本数不多，SNR较 低时也有效。然而，其估计模型参数的方法为最小二乘法和迭代法。这些方法存在收敛性 不能被严格证明的不足。
【发明内容】
[0004] 为了识别欺骗干扰和真实信号，本发明提出了一种基于发射系统建模以及模型系 数辨识的欺骗干扰识别方法。在对信号发射机和干扰机的建模中主要考虑功率放大器的非 线性，因此，将发射机和干扰机简化为无记忆多项式模型，在对信道进行建模中考虑其多径 效应，将其简化为FIR滤波器模型。整个系统可以等效为一个Hammerstein模型，利用差分 进化狼群算法根据输入输出信号来进行系统辨识，确定系统模型参数。以此模型参数为基 础，通过统计检测方法进行欺骗干扰的识别。
[0005] 本发明采取以下技术方案：
[0006] 基于差分进化狼群算法的欺骗干扰识别方法，按如下步骤进行：
[0007] 第一步：考虑到非线性，将发射机和干扰机等效为无记忆多项式模型。
[0009] 其中，M是多项式系数个数，d(n)是输入信号，bk是多项式系数。另外，在不改变 系统特性的情况下，为了确定系统的唯一性，令h = 1。
[0010] 第二步：考虑到多径响应，将无线传输信道等效为FIR滤波器模型。
[0012] 上式中，hk是信道响应系数，N是FIR滤波器的阶数，w(n)~N(0, 〇2)是加性高斯 白噪声，x(n)是输入信号，y (η)为接收机接收到的信道输出信号。
[0013] 第三步：对整个通信传输系统进行建模。
[0015] 第四步：接收机根据接收到的数据码元y (η)，估计得到发射端的C/A码，根据输入 输出利用差分进化狼群算法进行系统辨识，得到模型参数的估计值。
[0016] 本发明采用差分进化算法的差分策略，通过该算法中的变异、交叉和选择的过程， 使得算法能跳出局部最优，同时也充分利用差分进化算法在局部寻优方面的优点，提高算 法的寻优精度。
[0017] DE-WPA算法的基本思想是：种群先按照基本狼群算法进行，即依次执行游走、奔 袭、围攻行为，之后进行差分进化，即变异、交叉及贪婪选择操作，之后再对头狼位置进行更 新。总体上可以认为是依次进行WPA算法搜索和DE算法变异，是一种两层的串行结构。
[0018] DE-WPA算法的具体步骤如下：
[0019] 步骤1，数值初始化。初始化狼群中人工狼位置&及其数目Ν，最大迭代次数1^_， 探狼比例因子α，最大游走次数Tniax，距离判定因子w，步长因子S，更新比例因子R。
[0020] 步骤2,选取最优人工狼为头狼，除头狼外最佳的S_num匹人工狼为探狼并执行游 走行为，直到某只探狼侦察到的猎物气味浓度Y1大于头狼所感知的猎物气味浓度Y lMd或达 到最大游走次数Tniax，则转步骤3。
[0021] 步骤3,人工猛狼据式（4)向猎物奔袭，若途中猛狼感知的猎物气味浓度Y^Ywd, 则Ylead= Y i，替代头狼并发起召唤行为；若Y lMd，则人工猛狼继续奔袭直到dls< dn_， 转步骤4。
[0023] 式中，Stepb是召唤行为的步长，G k是第k次迭代中头狼的位置，Jif是第k次迭代 中第i头狼的位置。
[0024] 步骤4,按式（5)对参与围攻行为的人工狼的位置进行更新，执行围攻行为。
[0026] 式中，step。是围攻行为的步长，λ是一个随机数，且λ~U[-l, 1]。
[0027] 步骤5,对这狼群进行差分变异。对种群中每个个体X1，随机生成三个互不相同的 整数ri，r2, r3e {1，2,…，N}，且要求i，r i，r2, r3这四个数互不相同，然后按（6)式生成变异 个体
[0029] 式中，F是变异尺度系数，V1是变异个体，心，义是根据ri，r 2，r3所选的三个 个体，若所生成的变异个体超出了解空间范围，则随机产生一个变异个体来代替。
[0030] 步骤6,对这狼群进行交叉操作。交叉操作按每个个体向量的分量进行。将变 异产生的变异个体与目标个体按照下面的公式进行交叉操作，具体执行过程如下：首先生 成一个随机整数rancK，然后对变异个体V1和种群的目标个体X i按式（7)得到试验个体 οχ^···,^,/ = ^···, V，为了保证个体的进化，首先通过随机选择使U1至少有一位 由V1贡献，其它位由交叉概率CR来决定具体哪位由X i贡献哪位由V i贡献：
[0031]
[0032] 式中randj是位于[0, 1]间的均勾分布的随机实数，而randn ;是属于{1，2,…D} 内随机产生的维数索引号，其保证了至少有一位由变异向量￥贡献。此处CR为交叉概 率因子，也是位于[0，1]间的一个常数。差分进化算法引入交叉操作是为了增加种群的多 样性。
[0033] 步骤7,对狼群进行选择操作：采用的是"贪婪选择"策略，经变异以及交叉操作后 生成的候选个体U1与目标个体X 1进行竞争：
[0035] 上式中，f是适应度函数，在UJP X i中选择适应度函数值最优者作为新的个体，替 换原来的个体。需要指出的是式（8)适合最大化。
[0036] 步骤8,按"胜者为王"的头狼产生规则对头狼位置进行更新。
[0037] 步骤9,判断是否达到优化精度要求或最大迭代次数Tniax，若达到则输出头狼的位 置，即所求问题的最优解，否则转步骤2。
[0038] 第五步：在利用差分进化狼群算法获得系统参数的估计值之后，采用直观的欧氏 距离检测法进行识别，判决准则为：
[0040] 其中，假设H。表示接收机接收到的是欺骗干扰信号，H i表示接收机接收到的信号 是真实信号。式中，G1SGPS发射机系统参数向量，Θ 2为干扰机系统参数向量，§是根据 系统的输入输出估计得到的参数向量。
[0041] 本发明基于差分进化狼群算法的欺骗干扰识别方法综合考虑发射机的非线性特 性和信道的多径响应，进行数学统计建模，并利用系统辨识算法进行系统辨识，得到了参数 向量6，于该参数向量，统计检测中最直观的欧氏距离法实现对欺骗干扰的检测。本发明 方法运算简单，识别性能好。
【附图说明】
[0042] 图1是通信系统结构图。
[0043] 图2是干扰识别总体方案图。
[0044] 图3是差分进化狼群算法流程图。
[0045] 图4是分别用差分进化狼群算法和狼群算法进行参数估计并以此进行干扰识别 的效果的比较。
【具体实施方式】
[0046] 本发明针对GPS欺骗干扰检测问题展开研究，通过综合考虑发射机的非线性和无 线通信信道的多径效应，将整个通信系统等效为Hammerstein模型，然后利用差分进化狼 群算法得到模型参数的估计值，以此进行欺骗干扰检测。
[0047] 以下结合附图对本发明作进一步说明。
[0048] 本实施例基于差分进化狼群算法的欺骗干扰识别方法，按如下步骤进行：
[0049] 第一步：考虑到非线性，将发射机和干扰机等效为无记忆多项式模型。
[0051] 仿真中 M = 4。
[0052] 第二步：考虑到多径响应，将无线传输信道等效为FIR滤波器模型。
[0054] 仿真中 N = 3。
[0055] 第三步：对整个通信传输系统进行建模。
[0057] 第四步：接收机根据接收到的数据码元y (η)，估计得到发射端的C/A码，根据输入 输出利用差分进化狼群算法进行系统辨识，得到模型参数的估计值。
[0058] 本发明采用差分进化算法的差分策略，通过该算法中的变异、交叉和选择的过程， 使得算法能跳出局部最优，同时也充分利用差分进化算法在局部寻优方面的优点，提高算 法的寻优精度。
[0059] DE-WPA算法的基本思想是：种群先按照基本狼群算法进行，即依次执行游走、奔 袭、围攻行为，之后进行差分进化，即变异、交叉及贪婪选择操作，之后再对头狼位置进行更 新，并按照"强者生存"的狼群更新机制进行群体更新。总