一种冲击噪声环境下的波束空间测向方法技术

本发明专利技术公开一种冲击噪声环境下的波束空间测向方法，包括步骤；获取信号采样数据，对数据进行波束空间处理，得到波束空间极大似然估计的目标方程；初始化海豚群群体和初始化信仰空间；计算适应度，记录海豚群群体全局最优位置；海豚群分成不同的等规模团队，团队里每个海豚依据信息共享获取临时位置；计算临时适应度，更每个海豚的局部最优位置；依据文化机制更新海豚的位置；计算适应度，更新每只海豚的局部最优位置，更新海豚群群体全局最优位置，根据部分局部最优位置更新信仰空间；判断是否达到最大迭代次数，并得出测向结果输出；本发明专利技术可在冲击噪声环境下对相干信源进行有效测向，具有收敛速度快和收敛精度高的优点。

A Beam Space Direction Finding Method in Shock Noise Environment

The invention discloses a beam space direction finding method in impulse noise environment, which includes steps: acquiring signal sampling data, processing data in beam space, obtaining the objective equation of maximum likelihood estimation of beam space; initializing dolphin population and belief space; calculating fitness, recording global optimal position of dolphin population; dividing dolphin population into different equal scales; Team, each dolphin in the team obtains temporary location based on information sharing; calculates temporary fitness, more local optimal location of each dolphin; updates dolphin location according to cultural mechanism; calculates fitness, updates local optimal location of each dolphin, updates global optimal location of dolphin population, and updates belief space according to partial optimal location; The method has the advantages of fast convergence speed and high convergence precision, and can effectively locate coherent sources in impulse noise environment.

【技术实现步骤摘要】
一种冲击噪声环境下的波束空间测向方法
本专利技术涉及阵列信号处理领域，具体涉及一种冲击噪声环境下的波束空间测向方法。
技术介绍
测向是阵列信号处理领域一个重要的研究方向，得到很多研究人员的关注。波束空间处理技术可以将阵元合成一定数量的波束通道作为数据接收通道，而波束空间的空间谱估计方法具有降低计算量、提高稳健性和降低系统的复杂性的优势，具有重要的研究意义。但现有的特征分解类波束空间测向方法，由于某些相干源的导向矢量与噪声子空间不完全正交，若不做特殊处理，无法直接对相干信源进行有效估计，空间平滑技术虽然可以对相干源进行波达方向估计，但其要求阵列具有平移不变性，且需要具有特殊的阵列结构。为了测相干信源方向对阵列有效孔径有损失，使其分辨力和测向精度都会下降。而冲击噪声更是一种复杂的噪声背景，不同特征指数对应不同噪声环境，特征指数越小，拖尾越严重，传统测向方法失效的可能性就越大，若在冲击噪声环境下可以设计新的低阶矩阵，必然会促进波束空间测向方法的近一步应用。在测向中传统的极大似然估计虽然是一种最佳估计方法，但是在阵列较大的情况下其运算量太大，且在冲击噪声环境下算法性能恶化。波束空间处理技术与阵元空间处理中每个阵元对应一个接收通道不同，波束空间方法能在大阵列和小信号数的场合有效地降低计算量。同时，就大规模天线测向系统而言，数据接收通道的减少可大幅减少测向系统的复杂性。但是，传统的极大似然测向方程的轮换投影变换方法不仅具有量化误差，而且具有计算复杂并易陷入局部收敛的缺点。鉴于上述缺陷，本专利技术创作者经过长时间的研究和实践终于获得了本专利技术。
技术实现思路
为解决上述技术缺陷，本专利技术采用的技术方案在于，提供一种冲击噪声环境下的波束空间测向方法，包括步骤：S1，获取信号采样数据，对所述信号采样数据进行波束空间处理，得到波束空间极大似然估计的目标方程；S2，初始化海豚群群体和初始化信仰空间；S3，计算适应度，记录海豚群群体全局最优位置；S4，将所述海豚群群体分成不同的等规模团队，所述团队里每个海豚依据信息共享获取个体海豚临时位置；根据所述个体海豚临时位置计算适应度，更新每个海豚的局部最优位置；S5，依据所述信仰空间更新海豚位置；S6，计算所述步骤S5更新后海豚位置的适应度，更新每只海豚的局部最优位置，更新海豚群群体全局最优位置，根据部分局部最优位置更新信仰空间；S7，判断是否达到最大迭代次数，并得出测向结果输出。较佳的，N个远场信号入射到设置有M个阵元的均匀线阵上，将第k次快拍数据写成矢量形式；所述第k次快拍数据的表达式为：其中，E(k)为快拍数据矢量，为冲击噪声数据矢量，S(k)为空间信号矢量，A(θ)为M×N维阵列流型矩阵；θ为到达角度矢量；所述到达角度矢量θ＝[θ1，θ2，...，θN]，其中，θi为第i个信号的入射方向与线阵法线的夹角，i＝1，2，...，N；所述M×N维阵列流型矩阵A(θ)的表达式为：A(θ)＝[a1(ω0)，a2(ω0)，...，aN(ω0)]，其中，第i个信号的导向矢量：其中，τli表示第i个信号到达第l个阵元时相对于参考阵元的时延，l＝1，2，...，M，ω0是信号的频率，j为复数单位；exp表示以e为底的指数函数；所述信号频率ω0的表达式为：其中，c为光速，λ为波长；以第一个阵元位置为参考点，第l个阵元的坐标为yl，时延其中，l＝1，2，...，M，i＝1，2，...，N。较佳的，所述步骤S1对空间阵元通过变换合成B个波束，再利用合成的波束数据进行波达方向估计；为形成的主瓣波束指向，T代表转置运算；通过阵列的接收数据形成B个波束所需要的归一化加权矩阵经过波束空间变换的波束输出为：其中，H代表共轭转置运算；所述归一化加权矩阵的表达式为：其中，m为波束覆盖空间中起始波束的指向角度影响参数；第次快拍采样获得的波束输出数据表达式为：V(k)＝[V1(k)，V2(k)，...，VB(k))]T，通过表达式变换得到分数低阶化向量为：其中，b＝1，2，...，B，为权系数，*代表共轭，T代表转置运算。较佳的，构造所述信号采样数据的低阶协方差矩阵，所述低阶协方差矩阵的表达式为：其中，K为最大快拍数，E为求期望运算；则波束空间导向矩阵的表达式：正交投影矩阵的表达式：角度估计的波束空间极大似然方程：其中，tr代表矩阵的迹运算。较佳的，海豚群体大小设置为为海豚团队个数，为每个海豚团队所包含的海豚个数，第h只海豚的位置为：其中，z代表迭代次数；在定义区间形成包含个维度为N的海豚位置其中，n＝1，2，...，N，z代表迭代次数；将所述第h只海豚至第z代所搜索到的最优位置记作所述迭代次数z的初始值为1。较佳的，所述信仰空间包含规范知识和形式知识，所述形式知识为所有海豚群中具有最优适应度的海豚位置即全局最优位置；所述信仰空间中所述规范知识为用<In，Ln，Un>定义其中In表示第n个变量的变化范围，In＝[ln，un]，ln和un分别为第n维变量的下限和上限；初始化信仰空间中的所述规范知识，海豚根据问题所给定的变量取值范围来初始化规范知识的下限和上限表示第n维变量的下限所对应的评价值，表示第n维变量的上限所对应的评价值，均初始化为-∞，初始时令所述迭代次数z为1。较佳的，所述步骤S3中，将所述第h只海豚位置对应于波达方向的角度，计算适应度，第h只海豚的位置的适应度函数为：其中，C为低阶协方差矩阵，为正交投影矩阵，tr代表矩阵的迹运算，记录所述整体适应度最大的位置定义为所述海豚群群体全局最优位置较佳的，所述步骤S4中，把整个海豚群等分成个团队，对于第h只海豚，其所属的团队标号为对应团队局部最优位置为：生成第h只海豚的临时位置为：其中，r1和r2是[0，1]之间均匀随机数，z代表迭代次数，计算第h只海豚临时位置的适应度值若则否则从而更新海豚局部最优位置。较佳的，所述步骤S5中，使用文化机制规范知识和海豚的局部最优位置调整位置变量变化步长及前进方向的影响函数；所述影响函数的表达式为：其中，η为缩放比例因子，为标准正态分布的随机数，z代表迭代次数，n＝1，2，...，N。较佳的，在所述步骤S6中，计算所述步骤S5更新后海豚位置的适应度若则z代表迭代次数，从而更新每个海豚的局部最优位置；更新所述海豚群群体全局最优位置根据接受函数挑选标号顺序前20％的海豚局部最优位置更新规范知识，所述接受函数用于选择能够直接影响当前信仰空间的知识经验的个体；设置第个位置影响规范知识的下限，第个位置影响规范知识的上限，更新所述规范知识；所述规范知识的更新方程为：判断是否达到最大迭代次数，若未达到，则返回所述步骤S5，并令所述迭代次数z加1；若达到最大迭代次数，输出所述海豚群群体全局最优位置，并将所述海豚群群体全局最优位置对应的波达方向角度作为测向结果输出。与现有技术比较本专利技术的有益效果在于：本专利技术有效解决了波束空间测向和极大似然测向中所遇到的难题，可在冲击噪声环境下对相干信源进行有效测向，具有收敛速度快和收敛精度高的优点，促进了波束空间测向方法的应用与发展。附图说明图1为本专利技术所述冲击噪声环境下的波束空间测向方法流程图；图2为特征指数为1.7时独立信源波达方向估计对比图；图3为特征指数为1.7时相干信源波达方向估计对比图；图4为特本文档来自技高网...

【技术保护点】
1.一种冲击噪声环境下的波束空间测向方法，其特征在于，包括步骤：S1，获取信号采样数据，对所述信号采样数据进行波束空间处理，得到波束空间极大似然估计的目标方程；S2，初始化海豚群群体和初始化信仰空间；S3，计算适应度，记录海豚群群体全局最优位置；S4，将所述海豚群群体分成不同的等规模团队，所述团队里每个海豚依据信息共享获取个体海豚临时位置；根据所述个体海豚临时位置计算适应度，更新每个海豚的局部最优位置；S5，依据所述信仰空间更新海豚位置；S6，计算所述步骤S5更新后海豚位置的适应度，更新每只海豚的局部最优位置，更新海豚群群体全局最优位置，根据部分局部最优位置更新信仰空间；S7，判断是否达到最大迭代次数，并得出测向结果输出。

【技术特征摘要】
1.一种冲击噪声环境下的波束空间测向方法，其特征在于，包括步骤：S1，获取信号采样数据，对所述信号采样数据进行波束空间处理，得到波束空间极大似然估计的目标方程；S2，初始化海豚群群体和初始化信仰空间；S3，计算适应度，记录海豚群群体全局最优位置；S4，将所述海豚群群体分成不同的等规模团队，所述团队里每个海豚依据信息共享获取个体海豚临时位置；根据所述个体海豚临时位置计算适应度，更新每个海豚的局部最优位置；S5，依据所述信仰空间更新海豚位置；S6，计算所述步骤S5更新后海豚位置的适应度，更新每只海豚的局部最优位置，更新海豚群群体全局最优位置，根据部分局部最优位置更新信仰空间；S7，判断是否达到最大迭代次数，并得出测向结果输出。2.如权利要求1所述的冲击噪声环境下的波束空间测向方法，其特征在于，N个远场信号入射到设置有M个阵元的均匀线阵上，将第k次快拍数据写成矢量形式；所述第k次快拍数据的表达式为：其中，E(k)为快拍数据矢量，为冲击噪声数据矢量，S(k)为空间信号矢量，A(θ)为M×N维阵列流型矩阵；θ为到达角度矢量；所述到达角度矢量θ＝[θ1，θ2，...，θN]，其中，θi为第i个信号的入射方向与线阵法线的夹角，i＝1，2，...，N；所述M×N维阵列流型矩阵A(θ)的表达式为：A(θ)＝[a1(ω0)，a2(ω0)，...，aN(ω0)]，其中，第i个信号的导向矢量：其中，τli表示第i个信号到达第l个阵元时相对于参考阵元的时延，l＝1，2，...，M，ω0是信号的频率，j为复数单位；exp表示以e为底的指数函数；所述信号频率ω0的表达式为：其中，c为光速，λ为波长；以第一个阵元位置为参考点，第l个阵元的坐标为yl，时延其中，l＝1，2，...，M，i＝1，2，...，N。3.如权利要求2所述的冲击噪声环境下的波束空间测向方法，其特征在于，所述步骤S1对空间阵元通过变换合成B个波束，再利用合成的波束数据进行波达方向估计；为形成的主瓣波束指向，T代表转置运算；通过阵列的接收数据形成B个波束所需要的归一化加权矩阵经过波束空间变换的波束输出为：其中，H代表共轭转置运算；所述归一化加权矩阵的表达式为：其中，m为波束覆盖空间中起始波束的指向角度影响参数；第次快拍采样获得的波束输出数据表达式为：V(k)＝[V1(k)，V2(k)，...，VB(k))]T，通过表达式变换得到分数低阶化向量为：其中，b＝1，2，...，B，为权系数，*代表共轭，T代表转置运算。4.如权利要求3所述的冲击噪声环境下的波束空间测向方法，其特征在于，构造所述信号采样数据的低阶协方差矩阵，所述低阶协方差矩阵的表达式为：其中，K为最大快拍数，E为求期望运算；则波束空间导向矩阵的表达式：正交投影矩阵的表达式：角度估计的波束空间极大似然方程：其中，tr代表矩阵的迹运算。5.如权利要求4...

【专利技术属性】
技术研发人员：杜亚男，李武，查月波，
申请(专利权)人：中国电子科技集团公司第三十八研究所，
类型：发明
国别省市：安徽,34