CN113111304B - 强冲击噪声下基于量子射线机理的相干分布源测向方法 - Google Patents
强冲击噪声下基于量子射线机理的相干分布源测向方法 Download PDFInfo
- Publication number
- CN113111304B CN113111304B CN202110357999.5A CN202110357999A CN113111304B CN 113111304 B CN113111304 B CN 113111304B CN 202110357999 A CN202110357999 A CN 202110357999A CN 113111304 B CN113111304 B CN 113111304B
- Authority
- CN
- China
- Prior art keywords
- quantum
- ray
- rotation angle
- phi
- rays
- 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.)
- Active
Links
- 230000001427 coherent effect Effects 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 47
- 239000011159 matrix material Substances 0.000 claims abstract description 37
- 238000013507 mapping Methods 0.000 claims abstract description 19
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 15
- 238000009792 diffusion process Methods 0.000 claims abstract description 13
- 239000013598 vector Substances 0.000 claims description 86
- 238000005315 distribution function Methods 0.000 claims description 9
- 101100187346 Aspergillus sp. (strain MF297-2) notP gene Proteins 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 description 10
- 238000005457 optimization Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- LVNGJLRDBYCPGB-LDLOPFEMSA-N (R)-1,2-distearoylphosphatidylethanolamine Chemical compound CCCCCCCCCCCCCCCCCC(=O)OC[C@H](COP([O-])(=O)OCC[NH3+])OC(=O)CCCCCCCCCCCCCCCCC LVNGJLRDBYCPGB-LDLOPFEMSA-N 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Image Generation (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明提供一种强冲击噪声下基于量子射线机理的相干分布源测向方法,包括:建立相干分布源的广义阵列流型,构造基于加权无穷范数低阶协方差矩阵的极大似然测向方程;计算更新后所有射线的适应度函数值,更新全局最优量子位置和局部最优量子位置;每条射线依概率从斯涅尔折射定律演化和随机演化两种演化规则中选择一种更新其量子位置;计算更新后所有射线的适应度函数值,更新全局最优量子位置和局部最优量子位置;判断是否达到最大迭代次数,若未达到,返回步骤三;若达到则终止循环迭代,输出全局最优量子位置,经过映射变换为全局最优位置对应中心方位角和角度扩散的极大似然估计值。本发明在强冲击噪声环境下具有鲁棒性,突破现有应用局限。
Description
技术领域
本发明涉及一种强冲击噪声环境下的基于量子射线机理的相干分布源测向方法,属于阵列信号处理领域。
背景技术
近年来,高分辨率的阵列信号处理技术备受关注,尤其是阵列测向技术,在地震探测、雷达、无源声纳、无线通信等诸多领域得到广泛应用,并形成了许多高分辨率估计算法。这些估计算法大多假设信号源是远场的点目标,点目标模型是实际环境的一种近似,且这些算法往往假设背景噪声为高斯噪声,利用二阶或高阶累积量进行分析可以获得理想的结果。然而,由于目标的尺寸取决于与阵列接受天线的距离,加之多径效应、散射等现象的影响,无线移动通信中信号波达角常在一定范围内被展开,呈现出分布信号源,而且有非高斯噪声的存在,如海杂波噪声、大气噪声、无线信道噪声等,这些噪声的模型可表示成SαS随机过程,它与高斯噪声模型失配,使得传统的基于二阶或高阶累积量的算法失效,根据分布源的散射特性,可以将其分为相干分布源和非相干分布源两类,因此研究在冲击噪声下基于量子射线机理的相干分布源测向方法具有重要的意义和价值。
使用极大似然算法进行相干分布源测向,可以获得高精度高分辨的测向性能,且可以分辨相干信源,但是需要对多维非线性优化问题进行全局最优值搜索,如何快速高精度得到搜索结果是极大似然测向方法应用的经典难题,使用智能优化算法对其求解是一种具有潜力的解决方案,但现有的智能优化算法在求解复杂测向问题时又存在一些不足,如收敛速度慢、易陷入局部极值等,因此需要针对具体问题设计新的智能优化算法进行求解。
通过对现有技术文献的检索发现,马培锋等在《声学技术》(2007,26(05):817-821)上发表的“相干分布源的方位估计”中提出Toeplitz化的相干分布源波达方向估计方法,此方法把点目标的Toeplitz化去相干算法引入到分布式目标的方位估计中,并对噪声子空间进行重构,在高斯噪声环境下表现出一定的鲁棒性,但在冲击噪声下效果不理想。
已有文献的检索结果表明,现有的相干分布源测向方法多采用广义Capon波束形成方法、DSPE方法和广义ESPRIT方法等,此类算法实时性好且计算量小,但大多不能直接对相干信源进行求解,并且在低信噪比和强冲击噪声的环境下性能较差,因此提出一种在强冲击噪声下基于量子射线机理的相干分布源测向方法,其具体是在强冲击噪声下设计了一种加权无穷范数低阶协方差矩阵,设计基于加权无穷范数低阶协方差矩阵的极大似然测向方程进行分布源测向,并通过量子射线机理快速得到测向结果,解决现有的相干分布源测向方法在强冲击噪声和低信噪比背景下性能恶化的技术难题。
发明内容
针对现有相干分布源测向方法的缺点和不足,本发明设计了一种强冲击噪声下基于量子射线机理的相干分布源测向方法,在强冲击噪声环境下具有鲁棒性,并设计了量子射线机理进行高效求解,突破了现有相干分布源测向方法的一些应用局限。
本发明的目的是这样实现的:步骤如下:
步骤一:建立相干分布源的广义阵列流型,获取阵列接收信号的快拍采样数据,构造基于加权无穷范数低阶协方差矩阵的极大似然测向方程;
步骤二:初始化每条射线的量子位置并设定相关参数,构造适应度函数,计算所有射线的适应度函数值,保存全局最优量子位置和局部最优量子位置,初始化每条射线的量子旋转角矢量并更新其量子位置,计算更新后所有射线的适应度函数值,更新全局最优量子位置和局部最优量子位置;
步骤三:每条射线依概率从斯涅尔折射定律演化和随机演化两种演化规则中选择一种更新其量子位置;
步骤四:将每条射线更新后的量子位置映射为射线位置,计算更新后所有射线的适应度函数值,更新全局最优量子位置和局部最优量子位置;
步骤五:判断是否达到最大迭代次数G,若未达到,令g=g+1,返回步骤三;若达到则终止循环迭代,输出全局最优量子位置,经过映射变换为全局最优位置对应中心方位角和角度扩散的极大似然估计值。
本发明还包括这样一些结构特征:
1.步骤一具体包括:分布源的阵列接收信号模型可以描述为角信号密度函数在分布空间积分的形式,即:
其中,x(t)为分布源的阵列接收信号,a(θ)为阵列导向矢量,θ表示分布源信号所处的方向,si(θ-θi,t)为第i个分布源的角信号密度函数,n(t)为噪声矢量,P为分布源个数;对于相干分布源模型,角信号密度函数可以表示为随机信号幅度与分布源的空间分布函数乘积的形式,即:
si(θ-θi,t)=si(t)gi(θ-θi)
其中,ψi=(θi,Δi)代表第i个相干分布源的角度参数,θi和Δi分别表示第i个相干分布源的中心方位角度和角度扩散,建立相干分布源的广义阵列流型及阵列接收信号模型,即:
x(t)=B(ψ)s(t)+n(t)
定义最大快拍数为Kp,假设阵列有P个波长为λ的相干分布源入射到由M个阵元组成的均匀线阵上,则均匀线阵接收的第k次快拍数据可以表示为x(k)=B(ψ)s(k)+n(k),其中,k=1,2,…,Kp,B(ψ)=[b(ψ1),b(ψ2),…,b(ψP)]为M×P维广义导向矩阵,第p个广义导向矢量为p=1,2…,P,x(k)=[x1(k),x2(k),…,xM(k)]T为M×1维阵列快拍数据,其中k为快拍次数,s(k)=[s1(k),s2(k),…,sP(k)]T为P×1维信号,n(k)为M×1维服从特征指数为α的标准SαS分布的复冲击噪声,j为复数单位,T为矩阵转置;第k次采样数据的加权无穷范数归一化信号可以表示为其中,β∈[0.8,1]为加权常数,则阵元接收数据的加权无穷范数低阶协方差矩阵为其中,Kp为最大快拍数,上标H表示共轭转置;构造基于加权无穷范数低阶协方差矩阵的极大似然测向方程为:
其中,PB(ψ)=B(ψ)[BH(ψ)B(ψ)]-1BH(ψ)为广义阵列流形矩阵B(ψ)的投影矩阵,上标H为矩阵共轭转置,R为阵元接收数据的加权无穷范数低阶协方差矩阵,argmax()表示寻找具有最大函数值的变量,trace表示对矩阵求迹。
2.步骤二具体包括:首先设定射线数目为N,最大迭代次数为G;在第g次迭代中,第条射线在2P维搜索空间中的量子位置为其中,表示第条射线第φ组量子位置矢量,φ=1,2,…,P,当g=1时,初代所有射线的量子位置的每一维初始化为[0,1]之间的均匀随机数,根据量子位置映射成位置的映射规则将所有射线的量子位置的每一维映射到搜索空间范围内,其中,X1,max,X1,min分别为中心方位角搜索空间上下限,X2,max,X2,min分别为角度扩散搜索空间上下限,φ=1,2,…,P,得到第条射线的位置令与分布源中心方位角矢量θ=[θ1,θ2,…,θP]相对应,与角度扩散矢量Δ=[Δ1,Δ2,…,ΔP]相对应,第条射线的适应度函数为根据适应度函数计算每条射线的适应度函数值,将拥有最优适应度函数值射线的量子位置保存为全局最优量子位置将当前所有射线的量子位置保存为局部最优量子位置集合其中,表示第条射线到目前代为止所搜索到的局部最优量子位置;
初始化第条射线的量子旋转角矢量其中,q为[0,1]之间均匀分布的随机数,使用模拟量子旋转门更新量子位置:其中代表第g+1代中第条射线第p维量子旋转角,其中,p=1,2…,2P,之后将更新后的量子位置映射为位置然后对所有射线进行适应度函数的计算和评价;更新到当前代为止全局最优量子位置和每条射线的局部最优量子位置,令g=g+1。
3.步骤三具体步骤为:对于第条射线,产生一个[0,1]之间均匀分布的随机数如果PS为选择概率,则对这条射线进行随机演化,该射线的量子旋转角矢量为其中,q为[0,1]之间均匀分布的随机数,将该射线的量子旋转角矢量进行分组,即:其中表示第g+1代中第条射线第φ组量子旋转角矢量,根据更新量子旋转角矢量,其中,X1,max,X1,min分别为中心方位角搜索空间上下限,X2,max,X2,min分别为角度扩散搜索空间上下限,为分裂因子,为[0,1]之间均匀分布的随机数,norm表示求矢量的长度;使用模拟量子旋转门更新量子位置:其中代表第g+1代中第条射线第p维量子旋转角,其中,p=1,2…,2P;
如果则对这条射线根据斯涅尔折射定律进行演化,首先根据生成量子原点,其中表示第g代中第条射线第φ组量子原点矢量,φ=1,2,…,P;将该射线的量子旋转角矢量进行分组,即:其中表示第g+1代中第条射线第φ组量子旋转角矢量,分两种情况:
(1)如果量子原点与当前量子位置相同,则根据更新量子旋转角矢量,其中,为[0,1]之间均匀分布的随机数,ζ为常数,φ=1,2,…,P;使用模拟量子旋转门更新量子位置:其中代表第g+1代中第条射线第p维量子旋转角,其中,p=1,2…,2P;
(2)如果量子原点与当前量子位置不相同,根据斯涅尔折射定律,生成量子法线矢量其中,表示第g代中第条射线第φ组量子法线矢量,其中,φ=1,2…,P;然后根据和分别对量子旋转角矢量和量子法线矢量进行归一化,φ=1,2,…,P,之后根据更新量子旋转角矢量,其中,nr为折射率,dot表示两矢量点积,φ=1,2…,P;使用模拟量子旋转门更新量子位置:其中代表第g+1代中第条射线第p维量子旋转角,其中,p=1,2…,2P。
与现有技术相比,本发明的有益效果是:
(1)针对现有相干分布源测向方法在强冲击噪声环境下性能恶化的问题,设计了更具鲁棒性的基于量子射线机理的相干分布源测向方法,在冲击噪声下设计了加权无穷范数低阶协方差矩阵,并利用极大似然测向方法实现了相干分布源测向。
(2)本发明设计的相干分布源测向方法设计了加权无穷范数低阶协方差矩阵,能够分辨相干信源,在强冲击噪声下实现了对目标的有效测向,设计的量子射线机理可以对加权无穷范数低阶协方差矩阵极大似然方程进行高精度求解,快速准确的得到测向结果。
(3)仿真实验证明了强冲击噪声下基于量子射线机理的相干分布源测向方法的有效性,突破了传统方法在强冲击噪声下性能恶化甚至失效的应用局限,且相对于传统的求解方法速度更快、精度更高。
附图说明
图1是本发明所设计的强冲击噪声下基于量子射线机理的相干分布源测向方法示意图。
图2是α=0.8,GSNR=20dB时对两个相干分布源的测向结果。
图3是α=0.8时测向角度均方根误差与广义信噪比关系曲线。
图4是α=0.8时估计成功概率与广义信噪比关系曲线。
图5是α=1.8时测向角度均方根误差与广义信噪比关系曲线。
图6是α=1.8估计成功概率与广义信噪比关系曲线。
图7是GSNR=20dB时估计成功概率与特征指数关系曲线。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述。
结合图1至图7,本发明的步骤如下:
步骤一,建立相干分布源的广义阵列流型,获取阵列接收信号的快拍采样数据,构造基于加权无穷范数低阶协方差矩阵的极大似然测向方程。
分布源的阵列接收信号模型可以描述为角信号密度函数在分布空间积分的形式,即:其中,x(t)为分布源的阵列接收信号,a(θ)为阵列导向矢量,θ表示分布源信号所处的方向,si(θ-θi,t)为第i个分布源的角信号密度函数,n(t)为噪声矢量,P为分布源个数。对于相干分布源模型,角信号密度函数可以表示为随机信号幅度与分布源的空间分布函数乘积的形式,即:si(θ-θi,t)=si(t)gi(θ-θi),其中,si(t)为随机信号幅度,gi(θ-θi)为分布源的空间分布函数。假设相干分布源的空间分布函数为高斯分布,即:其中,Δi为未知角分布参数,则有,其中,ψi=(θi,Δi)代表第i个相干分布源的角度参数,θi和Δi分别表示第i个相干分布源的中心方位角度和角度扩散,则可以建立相干分布源的广义阵列流型及阵列接收信号模型,即:x(t)=B(ψ)s(t)+n(t),其中,B(ψ)=[b(ψ1),b(ψ2),…,b(ψP)]为相干分布源的广义阵列流型矩阵,d为均匀线阵阵元间距,λ为分布源传播波长,p=1,2…,P,M为接收阵元数。
定义最大快拍数为Kp,假设阵列有P个波长为λ的相干分布源入射到由M个阵元组成的均匀线阵上,则均匀线阵接收的第k次快拍数据可以表示为x(k)=B(ψ)s(k)+n(k),其中,k=1,2,…,Kp,B(ψ)=[b(ψ1),b(ψ2),…,b(ψP)]为M×P维广义导向矩阵,第p个广义导向矢量为p=1,2…,P,x(k)=[x1(k),x2(k),…,xM(k)]T为M×1维阵列快拍数据,其中k为快拍次数,s(k)=[s1(k),s2(k),…,sP(k)]T为P×1维信号,n(k)为M×1维服从特征指数为α的标准SαS分布的复冲击噪声,j为复数单位,T为矩阵转置。
第k次采样数据的加权无穷范数归一化信号可以表示为其中,β∈[0.8,1]为加权常数,则阵元接收数据的加权无穷范数低阶协方差矩阵可以表示为其中,Kp为最大快拍数,上标H表示共轭转置。构造基于加权无穷范数低阶协方差矩阵的极大似然测向方程其中,PB(ψ)=B(ψ)[BH(ψ)B(ψ)]-1BH(ψ)为广义阵列流形矩阵B(ψ)的投影矩阵,上标H为矩阵共轭转置,R为阵元接收数据的加权无穷范数低阶协方差矩阵,argmax()表示寻找具有最大函数值的变量,trace表示对矩阵求迹。
步骤二,初始化每条射线的量子位置并设定相关参数,构造适应度函数,计算所有射线的适应度函数值,保存全局最优量子位置和局部最优量子位置,初始化每条射线的量子旋转角矢量并更新其量子位置,计算更新后所有射线的适应度函数值,更新全局最优量子位置和局部最优量子位置。
首先设定射线数目为N,最大迭代次数为G。在第g次迭代中,第条射线在2P维搜索空间中的量子位置为其中,表示第条射线第φ组量子位置矢量,φ=1,2,…,P,当g=1时,初代所有射线的量子位置的每一维初始化为[0,1]之间的均匀随机数,根据量子位置映射成位置的映射规则将所有射线的量子位置的每一维映射到搜索空间范围内,其中,X1,max,X1,min分别为中心方位角搜索空间上下限,X2,max,X2,min分别为角度扩散搜索空间上下限,φ=1,2,…,P,得到第条射线的位置令与分布源中心方位角矢量θ=[θ1,θ2,…,θP]相对应,与角度扩散矢量Δ=[Δ1,Δ2,…,ΔP]相对应,第条射线的适应度函数为根据适应度函数计算每条射线的适应度函数值,将拥有最优适应度函数值射线的量子位置保存为全局最优量子位置将当前所有射线的量子位置保存为局部最优量子位置集合其中,表示第条射线到目前代为止所搜索到的局部最优量子位置。
初始化第条射线的量子旋转角矢量其中,q为[0,1]之间均匀分布的随机数,使用模拟量子旋转门更新量子位置:其中代表第g+1代中第条射线第p维量子旋转角,其中,p=1,2…,2P,之后将更新后的量子位置映射为位置然后对所有射线进行适应度函数的计算和评价。更新到当前代为止全局最优量子位置和每条射线的局部最优量子位置,令g=g+1。
步骤三,每条射线依概率从斯涅尔折射定律演化和随机演化两种演化规则中选择一种更新其量子位置,具体步骤为:
对于第条射线,产生一个[0,1]之间均匀分布的随机数如果其中,PS为选择概率,则对这条射线进行随机演化,该射线的量子旋转角矢量为其中,q为[0,1]之间均匀分布的随机数,将该射线的量子旋转角矢量进行分组,即:其中表示第g+1代中第条射线第φ组量子旋转角矢量,根据更新量子旋转角矢量,其中,X1,max,X1,min分别为中心方位角搜索空间上下限,X2,max,X2,min分别为角度扩散搜索空间上下限,为分裂因子,为[0,1]之间均匀分布的随机数,norm表示求矢量的长度。使用模拟量子旋转门更新量子位置:其中代表第g+1代中第条射线第p维量子旋转角,其中,p=1,2…,2P。
如果则对这条射线根据斯涅尔折射定律进行演化,首先根据生成量子原点,其中表示第g代中第条射线第φ组量子原点矢量,φ=1,2,…,P。将该射线的量子旋转角矢量进行分组,即:其中表示第g+1代中第条射线第φ组量子旋转角矢量,分两种情况:
(1)如果量子原点与当前量子位置相同,则根据更新量子旋转角矢量,其中,为[0,1]之间均匀分布的随机数,ζ为常数,φ=1,2,…,P。使用模拟量子旋转门更新量子位置:其中代表第g+1代中第条射线第p维量子旋转角,其中,p=1,2…,2P。
(2)如果量子原点与当前量子位置不相同,根据斯涅尔折射定律,生成量子法线矢量其中,表示第g代中第条射线第φ组量子法线矢量,其中,φ=1,2…,P。然后根据和分别对量子旋转角矢量和量子法线矢量进行归一化,φ=1,2,…,P,之后根据更新量子旋转角矢量,其中,nr为折射率,dot表示两矢量点积,φ=1,2…,P。使用模拟量子旋转门更新量子位置:其中代表第g+1代中第条射线第p维量子旋转角,其中,p=1,2…,2P。
步骤四,将每条射线更新后的量子位置映射为射线位置,计算更新后所有射线的适应度函数值,更新全局最优量子位置和局部最优量子位置。
步骤五,判断是否达到最大迭代次数G,若未达到,令g=g+1,返回步骤三;若达到则终止循环迭代,输出全局最优量子位置,经过映射变换为全局最优位置对应中心方位角和角度扩散的极大似然估计值。
本发明所设计的基于量子射线机理的相干分布源测向方法记作QRO-INF-ML;基于粒子群搜索机制的分数低阶协方差极大似然测向方法记作PSO-FLOC-ML,所涉及的粒子群搜索机制参考《山东大学学报(工学版)》(2010,40(01):133-138)上发表的“冲击噪声背景下的动态DOA跟踪”。
在仿真实验中两个中心方位角矢量为θ=[30°,-30°],角度扩散矢量为Δ=[5°,1°]的相干分布源入射到间距为载波波长一半的均匀线阵上。
仿真实验参数设置如下:Kp=1000,G=1000,M=8,N=30,PS=0.35,nr=0.5,[X1,min,X1,max]=[-90,90],[X2,min,X2,max]=[0,10],s=7.5,ζ=0.001,β=1。Monte Carlo实验次数为100,并规定对单个相干分布源估计中心方位角和实际中心方位角差值小于1°为估计成功一次。分别在特征指数α=0.8和α=1.8两种不同强度的冲击噪声下,对比测向角度均方根误差随广义信噪比变化曲线和估计成功概率随广义信噪比变化曲线,传统方法在强冲击噪声下完全失效,体现了本发明所设计的方法相对于传统方法的在冲击噪声环境下的优势。
从仿真图中可以看出本发明所设计的基于量子射线机理的相干分布源测向方法在不同强度的冲击噪声下都表现出优于传统方法的测向性能。
Claims (5)
1.强冲击噪声下基于量子射线机理的相干分布源测向方法,其特征在于:步骤如下:
步骤一:建立相干分布源的广义阵列流型,获取阵列接收信号的快拍采样数据,构造基于加权无穷范数低阶协方差矩阵的极大似然测向方程;
步骤二:初始化每条射线的量子位置并设定相关参数,构造适应度函数,计算所有射线的适应度函数值,保存全局最优量子位置和局部最优量子位置,初始化每条射线的量子旋转角矢量并更新其量子位置,计算更新后所有射线的适应度函数值,更新全局最优量子位置和局部最优量子位置;
步骤三:每条射线依概率从斯涅尔折射定律演化和随机演化两种演化规则中选择一种更新其量子位置;
步骤四:将每条射线更新后的量子位置映射为射线位置,计算更新后所有射线的适应度函数值,更新全局最优量子位置和局部最优量子位置;
步骤五:判断是否达到最大迭代次数G,若未达到,令g=g+1,返回步骤三;若达到则终止循环迭代,输出全局最优量子位置,经过映射变换为全局最优位置对应中心方位角和角度扩散的极大似然估计值。
2.根据权利要求1所述的强冲击噪声下基于量子射线机理的相干分布源测向方法,其特征在于:步骤一具体包括:分布源的阵列接收信号模型可以描述为角信号密度函数在分布空间积分的形式,即:
其中,x(t)为分布源的阵列接收信号,a(θ)为阵列导向矢量,θ表示分布源信号所处的方向,si(θ-θi,t)为第i个分布源的角信号密度函数,n(t)为噪声矢量,P为分布源个数;对于相干分布源模型,角信号密度函数可以表示为随机信号幅度与分布源的空间分布函数乘积的形式,即:
si(θ-θi,t)=si(t)gi(θ-θi)
其中,ψi=(θi,Δi)代表第i个相干分布源的角度参数,θi和Δi分别表示第i个相干分布源的中心方位角度和角度扩散,建立相干分布源的广义阵列流型及阵列接收信号模型,即:
x(t)=B(ψ)s(t)+n(t)
定义最大快拍数为Kp,假设阵列有P个波长为λ的相干分布源入射到由M个阵元组成的均匀线阵上,则均匀线阵接收的第k次快拍数据可以表示为x(k)=B(ψ)s(k)+n(k),其中,k=1,2,…,Kp,B(ψ)=[b(ψ1),b(ψ2),…,b(ψP)]为M×P维广义导向矩阵,第p个广义导向矢量为x(k)=[x1(k),x2(k),…,xM(k)]T为M×1维阵列快拍数据,其中k为快拍次数,s(k)=[s1(k),s2(k),…,sP(k)]T为P×1维信号,n(k)为M×1维服从特征指数为α的标准SαS分布的复冲击噪声,j为复数单位,T为矩阵转置;第k次采样数据的加权无穷范数归一化信号可以表示为其中,β∈[0.8,1]为加权常数,则阵元接收数据的加权无穷范数低阶协方差矩阵为其中,Kp为最大快拍数,上标H表示共轭转置;构造基于加权无穷范数低阶协方差矩阵的极大似然测向方程为:
其中,PB(ψ)=B(ψ)[BH(ψ)B(ψ)]-1BH(ψ)为广义阵列流形矩阵B(ψ)的投影矩阵,上标H为矩阵共轭转置,R为阵元接收数据的加权无穷范数低阶协方差矩阵,argmax()表示寻找具有最大函数值的变量,trace表示对矩阵求迹。
3.根据权利要求1或2所述的强冲击噪声下基于量子射线机理的相干分布源测向方法,其特征在于:步骤二具体包括:首先设定射线数目为N,最大迭代次数为G;在第g次迭代中,第条射线在2P维搜索空间中的量子位置为其中,表示第条射线第φ组量子位置矢量,φ=1,2,…,P,当g=1时,初代所有射线的量子位置的每一维初始化为[0,1]之间的均匀随机数,根据量子位置映射成位置的映射规则将所有射线的量子位置的每一维映射到搜索空间范围内,其中,X1,max,X1,min分别为中心方位角搜索空间上下限,X2,max,X2,min分别为角度扩散搜索空间上下限,φ=1,2,…,P,得到第条射线的位置令与分布源中心方位角矢量θ=[θ1,θ2,…,θP]相对应,与角度扩散矢量Δ=[Δ1,Δ2,…,ΔP]相对应,第条射线的适应度函数为根据适应度函数计算每条射线的适应度函数值,将拥有最优适应度函数值射线的量子位置保存为全局最优量子位置将当前所有射线的量子位置保存为局部最优量子位置集合其中,表示第条射线到目前代为止所搜索到的局部最优量子位置;
4.根据权利要求1或2所述的强冲击噪声下基于量子射线机理的相干分布源测向方法,其特征在于:步骤三具体步骤为:对于第条射线,产生一个[0,1]之间均匀分布的随机数如果PS为选择概率,则对这条射线进行随机演化,该射线的量子旋转角矢量为其中,q为[0,1]之间均匀分布的随机数,将该射线的量子旋转角矢量进行分组,即:其中表示第g+1代中第条射线第φ组量子旋转角矢量,根据更新量子旋转角矢量,其中,X1,max,X1,min分别为中心方位角搜索空间上下限,X2,max,X2,min分别为角度扩散搜索空间上下限,为分裂因子,为[0,1]之间均匀分布的随机数,norm表示求矢量的长度;使用模拟量子旋转门更新量子位置:其中代表第g+1代中第条射线第p维量子旋转角,其中,p=1,2…,2P;
如果则对这条射线根据斯涅尔折射定律进行演化,首先根据生成量子原点,其中表示第g代中第条射线第φ组量子原点矢量,φ=1,2,…,P;将该射线的量子旋转角矢量进行分组,即:其中表示第g+1代中第条射线第φ组量子旋转角矢量,分两种情况:
(1)如果量子原点与当前量子位置相同,则根据更新量子旋转角矢量,其中,为[0,1]之间均匀分布的随机数,ζ为常数,φ=1,2,…,P;使用模拟量子旋转门更新量子位置:其中代表第g+1代中第条射线第p维量子旋转角,其中,p=1,2…,2P;
5.根据权利要求3所述的强冲击噪声下基于量子射线机理的相干分布源测向方法,其特征在于:步骤三具体步骤为:对于第条射线,产生一个[0,1]之间均匀分布的随机数如果PS为选择概率,则对这条射线进行随机演化,该射线的量子旋转角矢量为其中,q为[0,1]之间均匀分布的随机数,将该射线的量子旋转角矢量进行分组,即:其中表示第g+1代中第条射线第φ组量子旋转角矢量,根据更新量子旋转角矢量,其中,X1,max,X1,min分别为中心方位角搜索空间上下限,X2,max,X2,min分别为角度扩散搜索空间上下限,为分裂因子,为[0,1]之间均匀分布的随机数,norm表示求矢量的长度;使用模拟量子旋转门更新量子位置:其中代表第g+1代中第条射线第p维量子旋转角,其中,p=1,2…,2P;
如果则对这条射线根据斯涅尔折射定律进行演化,首先根据生成量子原点,其中表示第g代中第条射线第φ组量子原点矢量,φ=1,2,…,P;将该射线的量子旋转角矢量进行分组,即:其中表示第g+1代中第条射线第φ组量子旋转角矢量,分两种情况:
(1)如果量子原点与当前量子位置相同,则根据更新量子旋转角矢量,其中,为[0,1]之间均匀分布的随机数,ζ为常数,φ=1,2,…,P;使用模拟量子旋转门更新量子位置:其中代表第g+1代中第条射线第p维量子旋转角,其中,p=1,2…,2P;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110357999.5A CN113111304B (zh) | 2021-04-01 | 2021-04-01 | 强冲击噪声下基于量子射线机理的相干分布源测向方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110357999.5A CN113111304B (zh) | 2021-04-01 | 2021-04-01 | 强冲击噪声下基于量子射线机理的相干分布源测向方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113111304A CN113111304A (zh) | 2021-07-13 |
CN113111304B true CN113111304B (zh) | 2022-09-27 |
Family
ID=76713700
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110357999.5A Active CN113111304B (zh) | 2021-04-01 | 2021-04-01 | 强冲击噪声下基于量子射线机理的相干分布源测向方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113111304B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117970227B (zh) * | 2024-02-04 | 2024-07-19 | 哈尔滨工程大学 | 强冲击噪声下基于相干分布源的幅相误差和角度参数联合估计方法及系统 |
CN117993038A (zh) * | 2024-02-07 | 2024-05-07 | 哈尔滨工程大学 | 一种基于离散量子电磁场优化机制的分布式阵列布局结构优化方法、系统及存储介质 |
CN118191763A (zh) * | 2024-03-28 | 2024-06-14 | 哈尔滨工程大学 | 一种强冲击噪声下非相干分布源幅相误差校正和测向方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101982782A (zh) * | 2010-09-08 | 2011-03-02 | 西安电子科技大学 | 类微光子源导航体系中的混沌测频系统与方法 |
CN103901394A (zh) * | 2014-03-21 | 2014-07-02 | 哈尔滨工程大学 | 一种冲击噪声环境下的量子万有引力搜索动态doa估计方法 |
CN103901395A (zh) * | 2014-03-28 | 2014-07-02 | 哈尔滨工程大学 | 一种冲击噪声环境下相干信号波达方向动态跟踪方法 |
CN107290732A (zh) * | 2017-07-11 | 2017-10-24 | 哈尔滨工程大学 | 一种量子大爆炸的单基地mimo雷达测向方法 |
CN107658573A (zh) * | 2017-08-22 | 2018-02-02 | 哈尔滨工程大学 | 一种基于量子搜寻者搜索机制的圆环阵方向图综合方法 |
CN109376329A (zh) * | 2018-09-05 | 2019-02-22 | 哈尔滨工程大学 | 一种基于量子鸟群演化机制的阵列幅相误差校正方法 |
CN110046326A (zh) * | 2019-04-28 | 2019-07-23 | 哈尔滨工程大学 | 一种时频doa估计方法 |
CN112255629A (zh) * | 2020-09-30 | 2021-01-22 | 中国人民解放军空军工程大学 | 基于联合uca阵列的序贯esprit二维不相干分布源参数估计方法 |
-
2021
- 2021-04-01 CN CN202110357999.5A patent/CN113111304B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101982782A (zh) * | 2010-09-08 | 2011-03-02 | 西安电子科技大学 | 类微光子源导航体系中的混沌测频系统与方法 |
CN103901394A (zh) * | 2014-03-21 | 2014-07-02 | 哈尔滨工程大学 | 一种冲击噪声环境下的量子万有引力搜索动态doa估计方法 |
CN103901395A (zh) * | 2014-03-28 | 2014-07-02 | 哈尔滨工程大学 | 一种冲击噪声环境下相干信号波达方向动态跟踪方法 |
CN107290732A (zh) * | 2017-07-11 | 2017-10-24 | 哈尔滨工程大学 | 一种量子大爆炸的单基地mimo雷达测向方法 |
CN107658573A (zh) * | 2017-08-22 | 2018-02-02 | 哈尔滨工程大学 | 一种基于量子搜寻者搜索机制的圆环阵方向图综合方法 |
CN109376329A (zh) * | 2018-09-05 | 2019-02-22 | 哈尔滨工程大学 | 一种基于量子鸟群演化机制的阵列幅相误差校正方法 |
CN110046326A (zh) * | 2019-04-28 | 2019-07-23 | 哈尔滨工程大学 | 一种时频doa估计方法 |
CN112255629A (zh) * | 2020-09-30 | 2021-01-22 | 中国人民解放军空军工程大学 | 基于联合uca阵列的序贯esprit二维不相干分布源参数估计方法 |
Non-Patent Citations (4)
Title |
---|
An Efficient Approximation for Nakagami-m Quantile Function Based on Generalized Opposition-Based Quantum Salp Swarm Algorithm;Hongyuan Gao 等;《Mathmatical Problems in Engineering》;20190814;8-30 * |
Quantum fireworks algorithm for optimal cooperation mechanism of energy harvesting cognitive radio;Hongyuan Gao 等;《Journal of Systems Engineering and Electronics》;20180228;1-13 * |
在冲击噪声环境下基于子空间的测向算法研究;吕泽均 等;《航空学报》;20030331;174-177 * |
基于时间调制阵列的宽带信号测向技术;姜艳;《舰船电子工程》;20190720;60-62 * |
Also Published As
Publication number | Publication date |
---|---|
CN113111304A (zh) | 2021-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113111304B (zh) | 强冲击噪声下基于量子射线机理的相干分布源测向方法 | |
Xiang et al. | Improved de-multipath neural network models with self-paced feature-to-feature learning for DOA estimation in multipath environment | |
CN111123192B (zh) | 一种基于圆形阵列和虚拟扩展的二维doa定位方法 | |
CN112698264B (zh) | 增广互质阵列脉冲噪声环境下相干信源的doa估计方法 | |
CN109375154B (zh) | 一种冲击噪声环境下基于均匀圆阵的相干信号参数估计方法 | |
CN109669156B (zh) | 冲击噪声下基于量子帝王蝶的圆阵模式空间动态测向方法 | |
CN109239646B (zh) | 一种冲击噪声环境下连续量子水蒸发的二维动态测向方法 | |
CN109212466B (zh) | 一种基于量子蜻蜓演化机制的宽带测向方法 | |
CN113296049B (zh) | 互质阵列脉冲环境下非圆信号的共轭增广doa估计方法 | |
Ahmad et al. | Wideband DOA estimation based on incoherent signal subspace method | |
CN111273269B (zh) | 基于ipso-bp的频率分集阵列的雷达目标定位方法 | |
CN113759303B (zh) | 一种基于粒子群算法的无网格波达角估计方法 | |
CN108594165B (zh) | 一种基于期望最大化算法的窄带信号波达方向估计方法 | |
CN115236584A (zh) | 基于深度学习的米波雷达低仰角估计方法 | |
Ma et al. | A novel DOA estimation for low-elevation target method based on multiscattering center equivalent model | |
CN112800596B (zh) | 强冲击噪声下基于嵌套阵列的鲁棒动态测向方法 | |
Yang et al. | A correlation-aware sparse Bayesian perspective for DOA estimation with off-grid sources | |
Reaz et al. | A comprehensive analysis and performance evaluation of different direction of arrival estimation algorithms | |
CN116559768A (zh) | 一种强冲击噪声下非相干分布源测向方法 | |
CN115932714A (zh) | 一种基于双曲正切核相关熵的单快拍波达方向估计方法 | |
CN112363106B (zh) | 基于量子粒子群的信号子空间波达方向检测方法及系统 | |
CN114184999B (zh) | 一种互耦小孔径阵列的生成式模型处理方法 | |
Zhang et al. | DOA estimation method for underwater acoustic signals based on two-dimensional power distribution (TPD) for few element array | |
CN113378103B (zh) | 一种强冲击噪声下相干分布源动态跟踪方法 | |
CN117970227B (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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230224 Address after: 314000 rooms 103 and 105, building 1, No. 79, Jinsui Road, Jiaxing Economic and Technological Development Zone, Zhejiang Province Patentee after: JIAXING NUOAIDI COMMUNICATION TECHNOLOGY CO.,LTD. Address before: 150001 Intellectual Property Office, Harbin Engineering University science and technology office, 145 Nantong Avenue, Nangang District, Harbin, Heilongjiang Patentee before: HARBIN ENGINEERING University |
|
TR01 | Transfer of patent right |