一种多鸽群信息交互的单快拍测向方法技术

本发明专利技术提供一种多鸽群信息交互的单快拍测向方法，建立均匀线阵单快拍采样信号模型；获得单快拍极大似然方程；初始化鸽群，并将其划分为三个子鸽群；计算鸽群中鸽子位置的适应度值，确定每个子种群的局部最优位置和整个鸽群的全局最优位置；更新基本鸽群的速度和位置，产生混沌权重；反向鸽群中的鸽子根据跳转操作更新位置；更新鸽子位置；确定子鸽群中鸽子的局部最优位置和整个鸽群的全局最优位置；更新信仰空间；最终输出的鸽群全局最优位置即为来波方向估计值。本发明专利技术实现了仅对单个快拍的数据进行处理从而得到对阵列接收信号的波达方向估计，降低了DOA估计的运算量，同时提高了系统的实时性，实现了对目标来波的高精度测向。

【技术实现步骤摘要】
一种多鸽群信息交互的单快拍测向方法
本专利技术涉及一种多鸽群信息交互的单快拍测向方法，属于阵列信号处理领域。
技术介绍
波达方向角(DirectionofArrival,DOA)估计一直是阵列信号处理的热点内容，在通信、雷达和声纳等系统中有着广泛应用。传统的多重信号分类法(MUSIC)和信号参数估计的旋转不变子空间技术(ESPRIT)已经具有较高的估计性能，但是这些算法都是基于特征值分解运算的基础上进行的，而且往往为了获得良好的估计性能需要成百上千的快拍数，不仅实时性低，计算量也很大。为了降低DOA估计的运算量，提高系统的实时性，单快拍DOA估计受到了学者的广泛关注。单快拍DOA估计即只对单个快拍的数据进行处理，实现输入信号的波达方向估计。根据已有的技术文献发现，王凌等在《系统工程与电子技术》(2012,Vol.34,No.7,pp.1323–1328)发表的“单次快拍波达方向矩阵法”中，利用接收信号重排构造伪协方差矩阵，但该方法构造的伪协方差矩阵由于利用的信息有限，这种方法存在阵列自由度的损失。伍逸枫等在《电光与控制》(2010.Vol.17,No.3,pp.60–63)发表的“基于Toeplitz矩阵重构的相干信源波达方向估计研究”提出的重构协方差矩阵，在一定程度上提高了分辨率，但都是建立在较大的快拍数情况下，并不利于实时处理。已有的文献表明，单快拍DOA估计可以提高系统的实时性，减少运算量，但快拍数量的减少将导致估计性能不准确，因此需要设计一种高性能的单快拍测向方法。本专利技术设计了一种多鸽群信息交互搜索机制的单快拍测向方法，结合极大似然估计对信源进行精准测向。仿真结果表明，该方法可以实现单快拍测向，并且无论是独立源还是相干源都具有较高的测向精度。
技术实现思路
本专利技术的目的是为了提供一种多鸽群信息交互的单快拍测向方法，是一种在高斯噪声环境下基于文化、反向和混沌的多鸽群信息交互搜索机制的单快拍测向方法。本专利技术的目的是这样实现的：步骤如下：步骤一：建立均匀线阵单快拍采样信号模型；给出一个阵元数为N的均匀线阵，阵元间距有M个远场窄带信号分别从θi方向入射到该阵列，且入射信号和噪声信号不相关，i＝1,2,…,M；选取第一个阵元为参考阵元，则在时刻第k个阵元接收到的信号为：其中：为入射信号，为第k个阵元的噪声信号，k＝1,2,…,N；则阵列接收的一次快拍信号模型为：y(1)＝A(θ)s(1)+n(1)式中：y(1)＝[y1(1),y2(1),…,yN(1)]T，A(θ)＝[a(θ1),a(θ2),…,a(θM)]N×M为导向矩阵，其中第i个导向矢量为θ＝[θ1,θ2,…,θM]为来波方位矢量；s(1)＝[s1(1),s2(1),…,sM(1)]T为信号矢量，n(1)＝[n1(1),n2(1),…,nN(1)]T为阵列噪声矢量；步骤二：利用阵列接收到的单快拍数据构造伪协方差矩阵，再利用均匀线阵的导向矩阵构造正交投影矩阵，获得极大似然方程；利用阵列接收的单快拍数据构造矩阵：并对所构造的矩阵Ry(1)进行协方差处理，即正交投影矩阵为PA(θ)＝A(θ)(AH(θ)A(θ))-1AH(θ)，极大似然方程的角度估计值为：其中H代表共轭转置，tr()为矩阵求迹函数；步骤三：初始化鸽群，计算适应度值；产生3个子鸽群，将鸽群中的只鸽子划分为3个子鸽群，基本鸽群、反向鸽群和文化鸽群的鸽子数量分别为N1、N2和N3，分别按照鸽群基本算子、反向机制和文化机制进行演化，其中第h个鸽群随机产生Nh只鸽子，h＝1,2,3；每只鸽子搜索空间的维数定义为M维，第h个鸽群第i只鸽子的位置为：其中：t为迭代次数，初始设t＝1，定义第h个鸽群第i只鸽子位置的适应度函数为：其中：第1个鸽群中第i只鸽子的速度为其中1≤i≤N1；步骤四：计算鸽群中每只鸽子位置的适应度值，确定第h个鸽群第i只鸽子的局部最优位置和整个鸽群的全局最优位置其中，h＝1,2,3；鸽群根据知识空间的产生规则生成信仰空间；步骤五：基本鸽群中，每只鸽子按照鸽群基本算子更新其速度和位置；第i只鸽子第m维的速度更新方程为第i只鸽子第m维的位置更新方程为其中i＝1,2,…,N1，m＝1,2,…,M，为鸽群全局最优位置的第m维，为地图和指南针因数，为混沌方程所产生的混沌权重；具体方法为：初始的为[0,1]间的均匀随机数，初始值不能等于1、0.75、0.5、0.25和0，按照如下规则更新：步骤六：反向鸽群的鸽子通过跳转操作更新位置；反向鸽群中设置跳转概率J，依概率J决定是否对更新个体进行种群跳转；步骤七：文化鸽群根据影响函数生成鸽子新位置；根据文化机制的规范知识和局部最优位置来调整位置变量变化步长及前进方向的影响函数为：其中：η为缩放比例因子，为标准正态分布的随机数，1≤i≤N3，1≤m≤M；步骤八：计算每只鸽子位置的适应度值，进行鸽群间信息交互，确定第h个子鸽群中第i只鸽子的局部最优位置和整个鸽群的全局最优位置其中，h＝1,2,3；利用来计算第h个子鸽群中第i个鸽子位置的适应度值，将鸽群中鸽子位置的适应度值与其相应的局部最优位置的适应度值比较，若鸽子位置的适应度值更大，则将其更新为鸽群的局部最优位置；将鸽群中鸽子位置的适应度值与全局最优位置的适应度值比较，若鸽子位置的适应度值更大，则将其更新为鸽群的全局最优位置；步骤九：鸽群根据接受函数，按照文化演进规则更新信仰空间；鸽群设定接受函数，根据接受函数挑选鸽群中局部极值最优前20％的鸽子位置更新规范知识，更新信仰空间，和分别为影响规范函数下界和上界更新的鸽子的局部最优位置，其更新方程为：步骤十：判断是否达到最大迭代次数，若没有达到，令t＝t+1，返回步骤五继续进行；否则，输出鸽群全局最优位置即来波方向估计值。本专利技术还包括这样一些结构特征：1.步骤四还包括：生成初始信仰空间，初始化规范知识，规范知识的下限和上限根据待解决问题的变量取值范围来初始化；表示第m维变量的下限所对应的适应度值，表示第m维变量的下限所对应的评价值，均初始化为-∞。2.步骤六中在计算反向个体时用来代替初始变量范围([am,bm])，即其中是均匀分布在[0,1]间的随机数，1≤i≤N2，1≤m≤M。与现有技术相比，本专利技术的有益效果是：(1)本专利技术仅对单个快拍的数据进行处理，实现对阵列接收信号的波达方向估计，降低了DOA估计的运算量，同时提高了系统的实时性。(2)本专利技术对接收到的单快拍数据重新构造了矩阵，并对其进行协方差处理，再结合极大似然测向方法，实现了对目标的高精度测向。(3)所设计的多鸽群信息交互的搜索机制的单快拍测向方法同时适用于相干信源、独立信源和混合信源，并且具有较优秀的抗噪声能力，和较高的测向成功概率。附图说明图1：基于多鸽群信息交互的单快拍测向方法流程图。图2：三个相干信源时所提方法与运用MUSIC算法测向对比图。图3：三个独立信源时所提方法与运用MUSIC算法测向对比图。图4：四个信源混合有独立信源和相干信源，信源1、2和3独立，信源4与信源1相干时所提方法与运用MUSIC算法测向对比图。具体实施方式下面结合附图与具体实施方式对本专利技术作进一步详细描述。为便于叙述，将基于多鸽群信息交互的单快拍测向方法简记为MPIO。如图1所示，本专利技术技术方案包括如下步骤：步本文档来自技高网...

【技术保护点】
1.一种多鸽群信息交互的单快拍测向方法，其特征在于：步骤如下：步骤一：建立均匀线阵单快拍采样信号模型；给出一个阵元数为N的均匀线阵，阵元间距

【技术特征摘要】
1.一种多鸽群信息交互的单快拍测向方法，其特征在于：步骤如下：步骤一：建立均匀线阵单快拍采样信号模型；给出一个阵元数为N的均匀线阵，阵元间距有M个远场窄带信号分别从θi方向入射到该阵列，且入射信号和噪声信号不相关，i＝1,2,…,M；选取第一个阵元为参考阵元，则在时刻第k个阵元接收到的信号为：其中：为入射信号，为第k个阵元的噪声信号，k＝1,2,…,N；则阵列接收的一次快拍信号模型为：y(1)＝A(θ)s(1)+n(1)式中：y(1)＝[y1(1),y2(1),…,yN(1)]T，A(θ)＝[a(θ1),a(θ2),…,a(θM)]N×M为导向矩阵，其中第i个导向矢量为θ＝[θ1,θ2,…,θM]为来波方位矢量；s(1)＝[s1(1),s2(1),…,sM(1)]T为信号矢量，n(1)＝[n1(1),n2(1),…,nN(1)]T为阵列噪声矢量；步骤二：利用阵列接收到的单快拍数据构造伪协方差矩阵，再利用均匀线阵的导向矩阵构造正交投影矩阵，获得极大似然方程；利用阵列接收的单快拍数据构造矩阵：并对所构造的矩阵Ry(1)进行协方差处理，即正交投影矩阵为PA(θ)＝A(θ)(AH(θ)A(θ))-1AH(θ)，极大似然方程的角度估计值为：其中H代表共轭转置，tr()为矩阵求迹函数；步骤三：初始化鸽群，计算适应度值；产生3个子鸽群，将鸽群中的只鸽子划分为3个子鸽群，基本鸽群、反向鸽群和文化鸽群的鸽子数量分别为N1、N2和N3，分别按照鸽群基本算子、反向机制和文化机制进行演化，其中第h个鸽群随机产生Nh只鸽子，h＝1,2,3；每只鸽子搜索空间的维数定义为M维，第h个鸽群第i只鸽子的位置为：其中：t为迭代次数，初始设t＝1，定义第h个鸽群第i只鸽子位置的适应度函数为：其中：第1个鸽群中第i只鸽子的速度为其中1≤i≤N1；步骤四：计算鸽群中每只鸽子位置的适应度值，确定第h个鸽群第i只鸽子的局部最优位置和整个鸽群的全局最优位置其中，h＝1,2,3；鸽群根据知识空间的产生规则生成信仰空间；步骤五：基本鸽群中，每只鸽子按照鸽群基本算子更新其速度和位置；第i只鸽子第m维的速度更新方程为第i只鸽子第m维...

【专利技术属性】
技术研发人员：高洪元，马雨微，谢婉婷，刁鸣，张世铂，池鹏飞，苏雨萌，刘子奇，吕阔，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：黑龙江,23