CN108663654A - 一种基于连续量子鸽群的360度全方位动态测向方法 - Google Patents
一种基于连续量子鸽群的360度全方位动态测向方法 Download PDFInfo
- Publication number
- CN108663654A CN108663654A CN201810668563.6A CN201810668563A CN108663654A CN 108663654 A CN108663654 A CN 108663654A CN 201810668563 A CN201810668563 A CN 201810668563A CN 108663654 A CN108663654 A CN 108663654A
- Authority
- CN
- China
- Prior art keywords
- quantum
- pigeon
- array
- dove group
- dove
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/14—Systems for determining direction or deviation from predetermined direction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种基于连续量子鸽群的360度全方位动态测向方法,属于阵列信号处理技术领域。实现步骤如下:获取阵列接收的快拍数据,计算阵列信号经模式空间变换后的协方差矩阵;计算适应度值,更新全局最优量子位置;判断是否达到要求迭代次数:若未达到继续迭代;进入地标算子操作,选择适应度值较优的前一半作为当前鸽群并计算当前鸽群的量子位置中心;更新每只鸽子的量子位置及其相应映射态,计算适应度值,更新全局最优量子位置;判断是否达到要求的迭代次数,若未达到继续迭代;否则输出鸽群全局最优位置;更新搜索区间;判断是否达到最大快拍数,输出动态测向结果。本方法搜索速度快,跟踪精度高,可360度全方位动态测向,应用前景广泛。
Description
技术领域
本发明属于阵列信号处理技术领域,具体涉及一种基于连续量子鸽群的360度全方位动态测向方法。
背景技术
测向又称为空间谱估计或者波达方向(Direction of Arrival,DOA)估计,是阵列信号处理的一个重要研究领域,在卫星和移动通信系统、雷达、地震学等方面有着广泛的应用。传统的DOA估计往往针对的是固定信源,然而在实际情况中,信源的角度是随时变化的,因此,对入射角度随时间变化的信号源进行动态DOA跟踪是DOA估计理论应用中的一个重要课题,如何在非相干信源情况下以及相干信源情况下实现360度全方位动态DOA估计更是工程应用中遇到的难点之一。
经对已有文献的检索发现,刁鸣等在《系统与电子技术》(2009,Vol.29,No.12,pp.2046-2049)
上发表的“一种新的基于粒子群算法的DOA跟踪方法”中利用均匀线阵进行多移动目标的动态DOA估计,然而均匀线阵只能提供180度的方位角测向。为了实现360度全方位的动态测向,本发明选择均匀圆阵的阵列结构,均匀圆阵以其特殊的阵列结构和良好的测向性能,在测向应用中备受青睐,但是由于圆阵的阵列流形矩阵是非Vandermonde矩阵,在数学处理上比较困难,所以要对均匀圆形阵列的输出信号进行相位模式激励,使其成为模式空间内的虚拟线阵,其阵列流形矩阵才具有Vandermonde结构,一些只适用于等距均匀线阵的算法才可以运用于等距均匀圆阵。现有的圆阵模式空间处理方法没有应用到动态目标的DOA跟踪研究方向上。
此外,已有文献还表明测向方法大多使用基于子空间的测向方法,这些方法不仅在低信噪比时性能较差,而且在跟踪相干信源的情况下还需要解相干处理,步骤繁琐。本发明设计的是圆阵模式空间下基于极大似然的动态测向方法,不仅在信噪比低,快拍数小以及相干信源的情况下具有优越的DOA估计性能,还能直接处理阵列数据协方差矩阵,避免了子空间类测向方法不断进行协方差矩阵分解的问题,实现了非相干源和相干源360度的全方位动态测向。但是,极大似然方法的一个主要缺点是它涉及一个多维非线性优化问题,计算量巨大,耗时长。
因此,本发明设计了一种圆阵模式空间下基于连续量子鸽群的360度全方位动态测向方法。该方法利用连续量子鸽群机制在搜索区间内搜索极大似然方程估计的最优角度,通过逐步缩小搜索区间减少运算量,同时依据量子编码和模拟量子演化方程设计的鸽群计算方法还可以加快算法的收敛速度,增强全局搜索能力,提高跟踪精度。
发明内容
本发明的目的在于提供可减少搜索时间,提高搜索速度和跟踪速度,具有良好的实时性的一种基于连续量子鸽群的360度全方位动态测向方法。
本发明的目的通过如下技术方案来实现:
一种基于连续量子鸽群的360度全方位动态测向方法,包括如下步骤:
步骤一:获取阵列接收的快拍数据,定义变换矩阵T,计算阵列信号经模式变换后的协方差矩阵。
假设N个窄带远场信源以方位角θn,俯仰角方向入射到一个由M个阵元构成的半径为r的均匀圆阵上,n=1,2,…,N,入射波长为λ,那么阵列接收的第k次快拍数据的数学模型为式中,K为最大快拍数,x(k)=[x1(k),x2(k),…,xM(k)]T为M×1维的阵列接收的快拍数据矢量,为阵列M×N维的流形矩阵,θ=[θ1,θ2,…,θN]和为信源的方向矢量,是流形矩阵的第n个导向矢量,其中s(k)为N×1维的信号矢量,n(k)为M×1维的噪声矢量,阵列噪声假定为空间和时间均独立的高斯白噪声。这里只讨论所有信源都与阵列共面的情形,即所有的俯仰角均为90°。
均匀圆阵可激发的最大模式为向下取整函数。变换矩阵T定义为式中J=diag{j-LJ-L(-β),…,jLJL(-β)},其中Jl(·)为l阶第一类贝塞尔函数,l=-L,…,0,…,L。式中FH=[w-L,w-L+1,…,wL],其中 用T左乘x(k)可得y(k)=Tx(k)=Bs(k)+Tn(k),其中阵列流形其具有Vandermonde矩阵结构,此时均匀圆阵已经转换为虚拟线阵。第1次快拍采样数据经过模式空间变换后的协方差矩阵为
步骤二:初始化搜索区间。
在第k次快拍,N个角度的搜索区间定义为其中un(k)和gn(k)分别为第k次快拍第n维角度搜索区间的上限和下限,其初始值分别取搜索区间定义域的上限和下限,n=1,2,…,N。连续量子鸽群的最大迭代次数取第k次快拍所有维搜索区间上限与下限之差整倍数的最大值,即其中取整数,为向下取整函数。
步骤三:初始化鸽群每只鸽子的量子位置和速度,计算适应度值,更新全局最优量子位置。
鸽群中共有P只鸽子,每只鸽子的量子位置在量子域[0,1]内随机初始化,速度在[-vmax,vmax]内随机初始化。第i只鸽子量子位置定义为其中t为当前迭代次数,为第i只鸽子的第n维量子位置,且第i只鸽子的位置为量子位置的相应映射态,映射关系为第i只鸽子的速度为且
评估第i只鸽子位置状态的适应度函数为其中为的映射矩阵,trace()为求矩阵的迹。利用此适应度函数计算每只鸽子位置的适应度值,选择适应度值最大的解的量子位置作为全局最优量子位置,记作其中为到第t次迭代为止鸽群最优量子位置的第n维,n=1,2,…,N,其相应映射态为全局最优鸽子位置
步骤四:根据地图和指南针算子更新每只鸽子的速度和量子位置及其相应映射态,计算适应度值,更新全局最优量子位置。
第i只鸽子速度的更新方程为其中i=1,2,…,P,n=1,2,…,N,R是地图和指南针因数,取值范围为(0,1),是(0,1)之间的随机数。如果速度超出边界值,则将其限制在边界,即若若
第i只鸽子量子位置更新方程为 同时按照映射方程获得量子位置的相应映射态。
利用适应度函数计算第i只鸽子位置的适应度值,更新全局最优量子位置,即对于第i只鸽子,若否则,
步骤五:判断地图和指南针算子操作的循环次数是否达到所要求的迭代次数:若未达到,令t=t+1,返回步骤四继续循环;否则,令t=t+1,停止地图和指南针算子操作,进入地标算子操作。
步骤六:进入地标算子操作,选择适应度值较优的前一半作为当前鸽群并计算当前鸽群的量子位置中心。
按照适应度值从大到小的规则对鸽子进行排序,当前代鸽子数量为Pt+1=Pt/2,选择适应度值较优的前一半作为当前鸽群,计算当前鸽群的量子位置中心的第n维,即并将其作为参考方向,i=1,2,…,P,n=1,2,…,N。
步骤七:更新每只鸽子的量子位置及其相应映射态,计算适应度值,更新全局最优量子位置。
按照方程更新每只鸽子的量子位置,其中i=1,2,…,P,n=1,2,…,N,是(0,1)之间的随机数。按照映射方程获得量子位置的相应映射态。利用适应度函数计算第i只鸽子位置的适应度值,更新全局最优量子位置,即对于第i只鸽子,若否则,
步骤八:判断地标算子操作的循环次数是否达到要求的迭代次数:若未达到,令t=t+1,返回步骤六;否则,输出鸽群全局最优位置bt+1。
步骤九:更新N个角度的搜索区间。
在第(k+1)次快拍,N个角度搜索区间更新为其中为第(k+1)个快拍第n维的中心值,即δ为遗传因子,ω为搜索区间中影响收敛速度的收敛因子,为第k个快拍第n维的估计值,d为搜索空间的搜索半径,n=1,2,…,N。
步骤十:判断是否达到最大快拍数:若未达到,继续获取下一个快拍采样经模式空间变换后的数据y(k+1)=[y1(k+1),…,yM(k+1)]T,进而协方差矩阵的更新方程为其中令k=k+1,返回步骤三;否则,根据得到的所有快拍采样下动态目标的估计值,输出动态测向的结果。
本发明的效果在于:
基于圆阵模式空间实现了适用于非相干信源和相干信源的360度全方位动态测向,具有更广泛的适用性。
设计了连续量子鸽群计算方法进而能对阵列的极大似然方程进行高精度求解,具有收敛速度快,收敛精度高的优点。
通过更新搜索中心、搜索区间和迭代次数,可减少搜索时间,提高搜索速度和跟踪速度,具有良好的实时性。
附图说明
图1为基于连续量子鸽群的360度全方位动态测向方法的示意图;
图2为连续量子鸽群计算方法的流程图;
图3为均匀圆阵示意图;
图4为信噪比SNR=10dB时两个独立信源的动态测向结果;
图5为信噪比SNR=10dB时两个相干信源的动态测向结果;
图6为信噪比SNR=10dB时三个独立信源的动态测向结果;
图7为信噪比SNR=10dB时两个相干信源一个独立信源的动态测向结果;
图8为信噪比SNR=10dB时三个相干信源的动态测向结果。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明:
图1、图2分别为圆阵模式空间下基于连续量子鸽群的360度全方位动态测向方法的示意图和连续量子鸽群计算方法的流程图。本发明能够完成高斯噪声环境下的全方位360度测向,设计了连续量子鸽群计算方法来计算阵列的极大似然方程,通过不断更新搜索中心、搜索区间和迭代次数,减少搜索时间,提高搜索速度和跟踪速度,以满足性能要求。本发明解决问题所采用的方案步骤如下:
步骤一:获取阵列接收的快拍数据,定义变换矩阵T,计算阵列信号经模式变换后的协方差矩阵。
假设N个窄带远场信源以方位角θn,俯仰角方向入射到一个由M个阵元构成的半径为r的均匀圆阵上,n=1,2,…,N,入射波长为λ,那么阵列接收的第k次快拍数据的数学模型为式中,K为最大快拍数,x(k)=[x1(k),x2(k),…,xM(k)]T为M×1维的阵列接收的快拍数据矢量,为阵列M×N维的流形矩阵,θ=[θ1,θ2,…,θN]和为信源的方向矢量,是流形矩阵的第n个导向矢量,其中s(k)为N×1维的信号矢量,n(k)为M×1维的噪声矢量,阵列噪声假定为空间和时间均独立的高斯白噪声。这里只讨论所有信源都与阵列共面的情形,即所有的俯仰角均为90°。
均匀圆阵可激发的最大模式为向下取整函数。变换矩阵T定义为式中J=diag{j-LJ-L(-β),…,jLJL(-β)},其中Jl(·)为l阶第一类贝塞尔函数,l=-L,…,0,…,L。式中FH=[w-L,w-L+1,…,wL],其中m=1,2,…,M-1。用T左乘x(k)可得y(k)=Tx(k)=Bs(k)+Tn(k),其中阵列流形其具有Vandermonde矩阵结构,此时均匀圆阵已经转换为虚拟线阵。第1次快拍采样数据经过模式空间变换后的协方差矩阵为
步骤二:初始化搜索区间。
在第k次快拍,N个角度的搜索区间定义为其中un(k)和gn(k)分别为第k次快拍第n维角度搜索区间的上限和下限,其初始值分别取搜索区间定义域的上限和下限,n=1,2,…,N。连续量子鸽群的最大迭代次数取第k次快拍所有维搜索区间上限与下限之差整倍数的最大值,即其中取整数,为向下取整函数。
步骤三:初始化鸽群每只鸽子的量子位置和速度,计算适应度值,更新全局最优量子位置。
鸽群中共有P只鸽子,每只鸽子的量子位置在量子域[0,1]内随机初始化,速度在[-vmax,vmax]内随机初始化。第i只鸽子量子位置定义为其中t为当前迭代次数,为第i只鸽子的第n维量子位置,且第i只鸽子的位置为量子位置的相应映射态,映射关系为第i只鸽子的速度为且
评估第i只鸽子位置状态的适应度函数为其中为的映射矩阵,trace()为求矩阵的迹。利用此适应度函数计算每只鸽子位置的适应度值,选择适应度值最大的解的量子位置作为全局最优量子位置,记作其中为到第t次迭代为止鸽群最优量子位置的第n维,n=1,2,…,N,其相应映射态为全局最优鸽子位置
步骤四:根据地图和指南针算子更新每只鸽子的速度和量子位置及其相应映射态,计算适应度值,更新全局最优量子位置。
第i只鸽子速度的更新方程为其中i=1,2,…,P,n=1,2,…,N,R是地图和指南针因数,取值范围为(0,1),是(0,1)之间的随机数。如果速度超出边界值,则将其限制在边界,即若若
第i只鸽子量子位置更新方程为 同时按照映射方程获得量子位置的相应映射态。
利用适应度函数计算第i只鸽子位置的适应度值,更新全局最优量子位置,即对于第i只鸽子,若否则,
步骤五:判断地图和指南针算子操作的循环次数是否达到所要求的迭代次数:若未达到,令t=t+1,返回步骤四继续循环;否则,令t=t+1,停止地图和指南针算子操作,进入地标算子操作。
步骤六:进入地标算子操作,选择适应度值较优的前一半作为当前鸽群并计算当前鸽群的量子位置中心。
按照适应度值从大到小的规则对鸽子进行排序,当前代鸽子数量为Pt+1=Pt/2,选择适应度值较优的前一半作为当前鸽群,计算当前鸽群的量子位置中心的第n维,即并将其作为参考方向,i=1,2,…,P,n=1,2,…,N。
步骤七:更新每只鸽子的量子位置及其相应映射态,计算适应度值,更新全局最优量子位置。
按照方程更新每只鸽子的量子位置,其中i=1,2,…,P,n=1,2,…,N,是(0,1)之间的随机数。按照映射方程获得量子位置的相应映射态。利用适应度函数计算第i只鸽子位置的适应度值,更新全局最优量子位置,即对于第i只鸽子,若否则,
步骤八:判断地标算子操作的循环次数是否达到要求的迭代次数:若未达到,令t=t+1,返回步骤六;否则,输出鸽群全局最优位置bt+1。
步骤九:更新N个角度的搜索区间。
在第(k+1)次快拍,N个角度搜索区间更新为其中为第(k+1)个快拍第n维的中心值,即δ为遗传因子,ω为搜索区间中影响收敛速度的收敛因子,为第k个快拍第n维的估计值,d为搜索空间的搜索半径,n=1,2,…,N。
步骤十:判断是否达到最大快拍数:若未达到,继续获取下一个快拍采样经模式空间变换后的数据y(k+1)=[y1(k+1),…,yM(k+1)]T,进而协方差矩阵的更新方程为其中令k=k+1,返回步骤三;否则,根据得到的所有快拍采样下动态目标的估计值,输出动态测向的结果。
图3中,均匀圆阵的M个阵元均匀分布在半径为r的圆周上,采用球坐标系来表示来波方向,坐标原点位于均匀圆阵的圆心位置上,来波的俯仰角是来波方向与z轴的夹角且方位角θ是来波方向在xoy平面内的投影与x轴的夹角且θ∈[0,360°]。
图4至图8所示的仿真中,均匀圆阵的阵元数M=16,阵元间距d=0.3λ,半径为噪声为空间和时间均独立的高斯白噪声。波达方向估计时的参数如下:信噪比SNR=10dB,最大快拍数K=300,每个方向的初始搜索区间均为[0°,360°],收敛因子ω=0.995,搜索半径d=3,遗传因子δ=0.8,协方差矩阵更新方程中的μ=0.95。连续量子鸽群计算方法中鸽群的鸽子数量P=40,地图和指南针算子R=0.3,鸽子速度最大值vmax=0.15。最大迭代次数取第k次快拍所有维搜索区间上限与下限之差整数倍的最大值,倍数地图和指南针算子操作迭代次数取最大迭代次数的地标算子操作的迭代次数取最大迭代次数的假设有2个信号源分别从θ1(k)=[210+5sin(2πk/200)]°,θ2(k)=[100+5sin(2πk/200)]°方向入射到均匀圆阵上,这两个独立信号源和相干信号源的动态测向结果分别如图4和图5所示。假设有3个信号源分别从θ1(k)=[210+5sin(2πk/200)]°,θ2(k)=[150+5sin(2πk/200)]°,θ3(k)=[100+5sin(2πk/200)]°方向入射到均匀圆阵上,这三个独立信号源,两个相干源一个独立源和三个相干信号源的动态测向结果分别如图6、图7、图8所示。通过仿真图可以看出,仿真环境下本发明所设计的方法可有效地对动态来波方向进行估计和跟踪,解决了动态目标360度全方位测向难题,同时在高斯噪声环境下无论是非相干信源还是相干信源都能够被准确的跟踪也说明了本发明所设计的方法具有广泛的适用性。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于连续量子鸽群的360度全方位动态测向方法,其特征在于,包括以下步骤:
步骤一:获取阵列接收的快拍数据,定义变换矩阵T,计算阵列信号经模式变换后的协方差矩阵;
步骤二:初始化搜索区间;
步骤三:初始化鸽群每只鸽子的量子位置和速度,计算适应度值,更新全局最优量子位置;
步骤四:根据地图和指南针算子更新每只鸽子的速度和量子位置及其相应映射态,计算适应度值,更新全局最优量子位置;
步骤五:判断地图和指南针算子操作的循环次数是否达到所要求的迭代次数:若未达到,令t=t+1,返回步骤四继续循环;否则,令t=t+1,停止地图和指南针算子操作,进入地标算子操作;
步骤六:进入地标算子操作,选择适应度值较优的前一半作为当前鸽群并计算当前鸽群的量子位置中心;
步骤七:更新每只鸽子的量子位置及其相应映射态,计算适应度值,更新全局最优量子位置;
步骤八:判断地标算子操作的循环次数是否达到要求的迭代次数:若未达到,令t=t+1,返回步骤六;否则,输出鸽群全局最优位置;
步骤九:更新N个角度的搜索区间;
步骤十:判断是否达到最大快拍数:若未达到,继续获取下一个快拍采样经模式空间变换后的数据,更新协方差矩阵,令k=k+1,返回步骤三;否则,根据得到的所有快拍采样下动态目标的估计值,输出动态测向的结果。
2.根据权利要求1所述的一种基于连续量子鸽群的360度全方位动态测向方法,其特征是步骤一具体包括:
假设N个窄带远场信源以方位角θn,俯仰角方向入射到一个由M个阵元构成的半径为r的均匀圆阵上,n=1,2,…,N,入射波长为λ,那么阵列接收的第k次快拍数据的数学模型为
式中,K为最大快拍数,x(k)=[x1(k),x2(k),…,xM(k)]T为M×1维的阵列接收的快拍数据矢量,为阵列M×N维的流形矩阵,θ=[θ1,θ2,…,θN]和为信源的方向矢量,是流形矩阵的第n个导向矢量,其中n=1,2,…,N,m=1,2,…,M-1,s(k)为N×1维的信号矢量,n(k)为M×1维的噪声矢量;
均匀圆阵可激发的最大模式 为向下取整函数,变换矩阵T定义为
式中J=diag{j-LJ-L(-β),…,jLJL(-β)},其中Jl(·)为l阶第一类贝塞尔函数,l=-L,…,0,…,L,式中FH=[w-L,w-L+1,…,wL],其中m=1,2,…,M-1,用T左乘x(k)可得y(k)=Tx(k)=Bs(k)+Tn(k),其中阵列流形第1次快拍采样数据经过模式空间变换后的协方差矩阵为
3.根据权利1所述的一种基于连续量子鸽群的360度全方位动态测向方法,其特征是步骤二具体包括:
在第k次快拍,N个角度的搜索区间定义为
其中un(k)和gn(k)分别为第k次快拍第n维角度搜索区间的上限和下限,其初始值分别取搜索区间定义域的上限和下限,n=1,2,…,N;
连续量子鸽群的最大迭代次数取第k次快拍所有维搜索区间上限与下限之差整倍数的最大值,即其中取整数,为向下取整函数。
4.根据权利要求1所述的一种基于连续量子鸽群的360度全方位动态测向方法,其特征是步骤四具体包括:
第i只鸽子速度的更新方程为
其中i=1,2,…,P,n=1,2,…,N,R是地图和指南针因数,取值范围为(0,1),是(0,1)之间的随机数,如果速度超出边界值,则将其限制在边界,即若若
第i只鸽子量子位置更新方程为
i=1,2,…,P,n=1,2,…,N,同时按照映射方程获得量子位置的相应映射态;
利用适应度函数计算第i只鸽子位置的适应度值,更新全局最优量子位置,即对于第i只鸽子,若否则,
5.根据权利要求1所述的一种基于连续量子鸽群的360度全方位动态测向方法,其特征是步骤六具体包括:
按照适应度值从大到小的规则对鸽子进行排序,当前代鸽子数量为Pt+1=Pt/2,选择适应度值较优的前一半作为当前鸽群,计算当前鸽群的量子位置中心的第n维,即并将其作为参考方向,i=1,2,…,P,n=1,2,…,N。
6.根据权利要求1所述的一种基于连续量子鸽群的360度全方位动态测向方法,其特征是步骤十具体包括:
判断是否达到最大快拍数:若未达到,继续获取下一个快拍采样经模式空间变换后的数据y(k+1)=[y1(k+1),…,yM(k+1)]T,进而协方差矩阵的更新方程为
其中
令k=k+1,返回步骤三;否则,根据得到的所有快拍采样下动态目标的估计值,输出动态测向的结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810668563.6A CN108663654B (zh) | 2018-06-26 | 2018-06-26 | 一种基于连续量子鸽群的360度全方位动态测向方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810668563.6A CN108663654B (zh) | 2018-06-26 | 2018-06-26 | 一种基于连续量子鸽群的360度全方位动态测向方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108663654A true CN108663654A (zh) | 2018-10-16 |
CN108663654B CN108663654B (zh) | 2022-03-18 |
Family
ID=63772335
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810668563.6A Active CN108663654B (zh) | 2018-06-26 | 2018-06-26 | 一种基于连续量子鸽群的360度全方位动态测向方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108663654B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109669155A (zh) * | 2018-11-16 | 2019-04-23 | 中国电子科技集团公司第三十八研究所 | 一种冲击噪声环境下的波束空间测向方法 |
CN109683125A (zh) * | 2018-11-16 | 2019-04-26 | 中国电子科技集团公司第三十八研究所 | 一种免疫飞蛾扑火机制的特殊阵列测向方法 |
CN110940949A (zh) * | 2019-12-11 | 2020-03-31 | 哈尔滨工程大学 | 强冲击噪声环境下基于量子企鹅搜索机制的互质阵列doa估计方法 |
CN111487594A (zh) * | 2020-04-23 | 2020-08-04 | 中国民航大学 | 一种基于粒子群优化的圆阵列波束形成方法 |
CN116520300A (zh) * | 2023-07-03 | 2023-08-01 | 天津知海科技有限公司 | 解算设备配置方法、装置、电子设备及可读存储介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101119147A (zh) * | 2006-08-01 | 2008-02-06 | 中兴通讯股份有限公司 | 一种空间到达方向的估计方法及装置 |
CN102608565A (zh) * | 2012-03-23 | 2012-07-25 | 哈尔滨工程大学 | 一种基于均匀圆阵列的波达方向估计方法 |
KR20120097296A (ko) * | 2011-02-24 | 2012-09-03 | 곽근창 | 다중화자의 다채널 음성신호로부터 음원분리를 통한 로봇청각시스템 및 방법 |
CN103901394A (zh) * | 2014-03-21 | 2014-07-02 | 哈尔滨工程大学 | 一种冲击噪声环境下的量子万有引力搜索动态doa估计方法 |
CN104049234A (zh) * | 2014-03-18 | 2014-09-17 | 电子科技大学 | 采用均匀圆阵快速测定空间谱的方法 |
CN105066998A (zh) * | 2015-08-03 | 2015-11-18 | 北京航空航天大学 | 一种基于量子鸽群优化的无人机自主空中加油目标检测方法 |
CN107238812A (zh) * | 2017-05-16 | 2017-10-10 | 哈尔滨工程大学 | 一种基于最小间隙阵列的鲁棒动态测向方法 |
-
2018
- 2018-06-26 CN CN201810668563.6A patent/CN108663654B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101119147A (zh) * | 2006-08-01 | 2008-02-06 | 中兴通讯股份有限公司 | 一种空间到达方向的估计方法及装置 |
KR20120097296A (ko) * | 2011-02-24 | 2012-09-03 | 곽근창 | 다중화자의 다채널 음성신호로부터 음원분리를 통한 로봇청각시스템 및 방법 |
CN102608565A (zh) * | 2012-03-23 | 2012-07-25 | 哈尔滨工程大学 | 一种基于均匀圆阵列的波达方向估计方法 |
CN104049234A (zh) * | 2014-03-18 | 2014-09-17 | 电子科技大学 | 采用均匀圆阵快速测定空间谱的方法 |
CN103901394A (zh) * | 2014-03-21 | 2014-07-02 | 哈尔滨工程大学 | 一种冲击噪声环境下的量子万有引力搜索动态doa估计方法 |
CN105066998A (zh) * | 2015-08-03 | 2015-11-18 | 北京航空航天大学 | 一种基于量子鸽群优化的无人机自主空中加油目标检测方法 |
CN107238812A (zh) * | 2017-05-16 | 2017-10-10 | 哈尔滨工程大学 | 一种基于最小间隙阵列的鲁棒动态测向方法 |
Non-Patent Citations (6)
Title |
---|
BIN LIAO, ET AL: "DOA Estimation and Tracking of ULAs with Mutual Coupling", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 * |
DAYONG ZHAO, ET AL: "Direction Finding of Maximum Likelihood Algorithm Using Artificial Bee Colony in the Impulsive Noise", 《2010 INTERNATIONAL CONFERENCE ON ARTIFICIAL INTELLIGENCE AND COMPUTATIONAL INTELLIGENCE》 * |
MING DIAO, ET AL: "DOA TRACKING BASED ON MIMO RADAR IN IMPULSIVE NOISE BACKGROUND", 《ICSP2014 PROCEEDINGS》 * |
刁鸣等: "一种新的基于粒子群算法的 DOA 跟踪方法", 《系统工程与电子技术》 * |
张亚平: "鸽群智能算法的改进及其在高超声速飞行控制中的应用", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
邱新建: "基于进化策略的DOA估计与跟踪算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109669155A (zh) * | 2018-11-16 | 2019-04-23 | 中国电子科技集团公司第三十八研究所 | 一种冲击噪声环境下的波束空间测向方法 |
CN109683125A (zh) * | 2018-11-16 | 2019-04-26 | 中国电子科技集团公司第三十八研究所 | 一种免疫飞蛾扑火机制的特殊阵列测向方法 |
CN110940949A (zh) * | 2019-12-11 | 2020-03-31 | 哈尔滨工程大学 | 强冲击噪声环境下基于量子企鹅搜索机制的互质阵列doa估计方法 |
CN110940949B (zh) * | 2019-12-11 | 2023-04-28 | 哈尔滨工程大学 | 强冲击噪声环境下基于量子企鹅搜索机制的互质阵列doa估计方法 |
CN111487594A (zh) * | 2020-04-23 | 2020-08-04 | 中国民航大学 | 一种基于粒子群优化的圆阵列波束形成方法 |
CN116520300A (zh) * | 2023-07-03 | 2023-08-01 | 天津知海科技有限公司 | 解算设备配置方法、装置、电子设备及可读存储介质 |
CN116520300B (zh) * | 2023-07-03 | 2023-09-01 | 天津知海科技有限公司 | 解算设备配置方法、装置、电子设备及可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108663654B (zh) | 2022-03-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108663654A (zh) | 一种基于连续量子鸽群的360度全方位动态测向方法 | |
CN107015191B (zh) | 一种在多径干扰环境下单偶极子极化敏感阵列降维doa估计方法 | |
CN103383450B (zh) | 共形阵列雷达幅相误差校正快速实现方法 | |
CN106597364B (zh) | 一种单天线单站无源定位的目标辐射源初始位置估计方法 | |
CN109239646A (zh) | 一种冲击噪声环境下连续量子水蒸发的二维动态测向方法 | |
CN103902826B (zh) | 一种冲击噪声环境下的多移动目标跟踪方法 | |
CN109633525A (zh) | 一种快速的均匀圆阵空间谱测向实现方法 | |
CN110687500B (zh) | 智能天线获取无线信号到达角的识别和定位方法及系统 | |
CN112130111B (zh) | 一种大规模均匀十字阵列中单快拍二维doa估计方法 | |
CN106501801B (zh) | 一种基于混沌多种群共生进化的双基地mimo雷达跟踪方法 | |
CN105445718B (zh) | 一种基于阵列重构的分布式多载舰超视距雷达的doa估计方法 | |
CN110007266A (zh) | 一种冲击噪声下的任意阵列相干源测向方法 | |
CN109669156B (zh) | 冲击噪声下基于量子帝王蝶的圆阵模式空间动态测向方法 | |
CN109782218A (zh) | 一种基于双平行天线阵的相干分布式非圆信号doa估计方法 | |
CN109633522A (zh) | 基于改进的music算法的波达方向估计方法 | |
CN107238812B (zh) | 一种基于最小间隙阵列的鲁棒动态测向方法 | |
CN104793177B (zh) | 基于最小二乘法的麦克风阵列测向方法 | |
CN109116295A (zh) | 基于相控阵选取基线的无源测向算法 | |
CN104931920A (zh) | Iesprit,一种基于任意阵列的空间信号doa的快速估计算法 | |
CN103096465B (zh) | 一种环境自适应的多目标直接定位方法 | |
CN109212466A (zh) | 一种基于量子蜻蜓演化机制的宽带测向方法 | |
CN108387868A (zh) | 一种基于信号到达角度差测量的扫描式辐射源定位方法 | |
CN109884582B (zh) | 利用一维测向快速确定目标三维坐标的方法 | |
CN111382830A (zh) | 一种基于混合粒子群算法的信号到达角度估计方法和装置 | |
CN110412501A (zh) | 基于改进极坐标表示方法模型下的雷达信号测向方法、装置及计算机存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |