CN109407045A - 一种非均匀传感器阵列宽带信号波达方向估计方法 - Google Patents
一种非均匀传感器阵列宽带信号波达方向估计方法 Download PDFInfo
- Publication number
- CN109407045A CN109407045A CN201811178910.3A CN201811178910A CN109407045A CN 109407045 A CN109407045 A CN 109407045A CN 201811178910 A CN201811178910 A CN 201811178910A CN 109407045 A CN109407045 A CN 109407045A
- Authority
- CN
- China
- Prior art keywords
- vector
- iteration
- band
- matrix
- broadband signal
- 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
-
- 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/80—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 ultrasonic, sonic or infrasonic waves
- G01S3/802—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)
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种非均匀传感器阵列宽带信号波达方向估计方法,包括以下步骤:S1、对非均匀传感器阵列中的样本序列作滑窗快速傅里叶变换,并构造实值加权样本协方差矢量yl;S2、构造对应的过完备阵列流型矩阵Φl;S3、利用实值加权样本协方差矢量yl和过完备阵列流型矩阵Φl设置迭代程序中空间联合稀疏表示向量γ的初始值γinit和正则化参数λ的初始值λinit;S4、利用实值加权样本协方差矢量yl和过完备阵列流型矩阵Φl更新空间联合稀疏表示向量γ和正则化参数λ;S5、当第i次迭代得到的γi与第i+1次迭代得到的γi+1之间满足||γi‑γi+1||2/||γi||2<α或者达到了最大迭代次数Pmax,终止迭代;否则返回S4;S6、空间联合稀疏表示向量γ的所有峰值,并输出宽带信号的波达方向估计值。
Description
技术领域
本发明涉及阵列信号处理领域,更具体的说涉及一种非均匀传感器阵列宽带信号波达方向估计方法。
背景技术
宽带信号具有抗干扰能力强、调制样式复杂、携带信息量大等特点,在雷达、通信、电子对抗侦察等领域具有十分广泛的应用,同时自然界中可接收的自然信号绝大多数都是宽带信号,因此对宽带信号的波达方向(Direction-of-Arrival,DOA)进行估计的问题正日益受到关注。
近年来,日益复杂多变的电磁环境导致来波信号数量增多,来波信号个数大于阵元数的欠定情况时有发生,而传统的传感器位置均匀排列的均匀阵列无法解决欠定DOA估计问题。对于欠定DOA估计问题,利用嵌套阵列或互质阵列等非均匀阵列进行DOA估计的研究正受到越来越多的关注。
利用非均匀阵列进行宽带信号的DOA估计,目前已有的研究都是先做多个窄频带划分,然后基于差分阵列的概念来构造虚拟扩展孔径阵列,而最终的DOA估计主要从两个不同的角度出发:1、合成各频带上虚拟阵列下的增广协方差矩阵,然后利用宽带信号DOA估计中经典的非相干子空间方法(Incoherent Signal Subspace Method,ISSM)实现DOA估计;2、基于稀疏表示思想,将该问题转化为多个频带上共享相同稀疏支撑集的非负稀疏信号重构问题并进行求解。
文献N.Hu,et al.Underdetermined DOA estimation method for widebandsignals using jointnonnegative sparse Bayesian learning,IEEE SPL,2017中公开了一种联合非负稀疏贝叶斯学习方法,该文献中利用各个频带上数据的联合非负稀疏性,建立了层级贝叶斯模型。利用期望最大化算法实现这种联合非负稀疏贝叶斯学习可以在迭代中给出波达方向的估值;该方法算法复杂度低、收敛速度快。然而在其层级贝叶斯模型中并未对基于凸优化的方法中常用于多频带联合稀疏约束的范数进行建模,因此其在低信噪比或者采样点数较少的情况下性能会受到影响。
发明内容
本发明要解决的技术问题是提供一种非均匀传感器阵列宽带信号波达方向估计方法,其能够提高波达方向的估计精度。
为了解决上述技术问题,本发明提供了一种非均匀传感器阵列宽带信号波达方向估计方法,其特征在于,包括以下步骤:
S1、对非均匀传感器阵列中各传感器接收通道采集的样本序列作滑窗快速傅里叶变换,并构造L个窄频带上的实值加权样本协方差矢量yl;
S2、构造与各窄频带上的实值加权样本协方差矢量yl相对应的过完备阵列流型矩阵Φl;
S3、利用L个窄频带上的实值加权样本协方差矢量yl和过完备阵列流型矩阵Φl设置迭代程序中空间联合稀疏表示向量γ的初始值γinit,并根据实值加权样本协方差矢量yl、过完备阵列流型矩阵Φl和空间联合稀疏表示向量的初始值γinit计算正则化参数λ的初始值λinit;
S4、利用L个窄频带上的实值加权样本协方差矢量yl和过完备阵列流型矩阵Φl更新空间联合稀疏表示向量γ和正则化参数λ;
S5、当第i次迭代得到的γi与第i+1次迭代得到的γi+1之间满足||γi-γi+1||2/||γi||2<α或者达到了最大迭代次数Pmax,终止迭代;否则继续返回S4迭代;其中α为一个任意趋于0的小数;
S6、当满足迭代终止条件时,搜索当前更新所得的空间联合稀疏表示向量γ的所有峰值,并根据空间联合稀疏表示向量γ的峰值位置输出宽带信号的波达方向估计值。
进一步的,步骤S1包括:
S11、对非均匀传感器阵列中各传感器接收通道采集的样本序列作滑窗快速傅里叶变换,得到输出快拍矢量序列xl(t),l=1,2,…,L;t=1,2,…,T;
S12、利用输出快拍矢量序列构造实值加权样本协方差矢量yl,l=1,2,…,L。
进一步的,步骤S4包括:
S41、利用实值加权样本协方差矢量yl、过完备阵列流型矩阵Φl以及上一轮迭代更新的空间联合稀疏表示向量γ,更新隐变量向量wl所服从的非负高斯分布的均值向量μl与协方差矩阵Σl,其中,若当前为第一轮迭代则空间联合稀疏表示向量γ=γinit;
S42、利用隐变量向量wl所服从的非负高斯分布的均值μl、协方差矩阵Σl和上一轮迭代更新的正则化参数λ,更新空间联合稀疏表示向量γ的值;其中,若当前为第一轮迭代则正则化参数λ=λinit;利用迭代更新的空间联合稀疏表示向量γ,更新正则化参数λ的值。
进一步的,步骤S11包括:
S111、利用窗函数对非均匀传感器阵列中各传感器接收通道采集的样本序列进行滑动数据选取,其中阵元个数为M;
S112、对各窗口中的加权非均匀阵列采样数据做快速傅里叶变换,其中,所加窗函数的窗长以及快速傅里叶变换点数为NFFT,相邻窗交迭Noverlap点;
S113、根据目标宽带信号所覆盖的频域范围,在第t个快速傅里叶变换窗口所得到的NFFT个频域值中选择L个作为待使用频带上的值,由于有M个传感器接收通道,可构造L个子频带上的M×1维向量:xl(t),l=1,2,…,L;t=1,2,…,T,其中每个向量是该频带上的“频域快拍”;
S114、对非均匀阵列采样数据共做T个滑窗快速傅里叶变换,可构造L个子频带上每个频带上T个“频域快拍”的数据向量:xl(t),l=1,2,…,L;t=1,2,…,T。
进一步的,步骤S12包括:
S121、计算L个窄频带上样本协方差矩阵其中(·)H表示共轭转置;
S122、对矩阵作向量化得到样本协方差矢量
S123、用加权矩阵J对进行加权处理以去除噪声项,得到加权样本协方差矢量其中权重矩阵J=[J1,J2,…JM-1]T,Jm=[em(M+1)-M+1,em(M+1)-M+2,…,em(M+1)],M2×1维向量em除第m个元素为1之外其他所有元素都为0,其中(·)T表示矩阵转置;
S124、利用J与构造加权样本协方差矢量的协方差矩阵
S125、构造实值加权样本协方差矢量其中Re(·)与Im(·)分别表示取实部与取虚部。
进一步的,步骤S2包括:
S21、以一定的角度网格间距构造可以覆盖阵列探测角度范围的密集角度抽样网格
S22、对于角度构造L个窄频带上与实值加权样本协方差矢量相对应的阵列导向矢量其中 为Kronecker积,原始阵列导向矢量fl为第l个频带的频率,dm为第m个传感器相对于参考传感器的间距,v为波速,(·)*表示求共轭;
S23、对于网格中的每一个角度计算并由此构造L个窄频带上与各实值加权样本协方差矢量相对应的过完备阵列流型矩阵Φl:
进一步的,步骤S3包括:
S31、计算空间联合稀疏表示向量γ的初值其中⊙为Hadamard积,||·||F为Frobenius范数;
S32、计算其中进而计算正则化参数λ的初值其中||·||2为范数;
S33、设置正则化参数λ的无信息超先验分布参数固定参数值c,d,其中c,d→0。
进一步的,步骤S41包括:
S421、计算L个窄频带上的隐变量向量wl中各元素的二阶矩:
其中特别当时,直接令
S422、利用计算所得的以及上一轮迭代所更新的正则化参数λ,更新空间联合稀疏表示向量γ的第n个元素
S423、利用上一轮迭代所更新的空间联合稀疏表示向量γ,更新正则化参数λ:
进一步的,步骤S6包括,当满足迭代终止条件时,搜索当前更新所得的空间联合稀疏表示向量γ的所有峰值,其峰值位置与网格中对应的角度即为宽带信号的波达方向估计值。
本发明的有益效果:
1、本发明将非均匀传感器阵列的宽带信号波达方向估计问题转化为多个窄频带上的具有相同非负稀疏支撑集的多向量联合非负稀疏重构问题,基于稀疏贝叶斯学习以及贝叶斯生成网络理论,建立层级贝叶斯模型并以此为基础进行求解,可实现欠定情况下的宽带信号波达方向估计;
2、本发明基于目标求解向量的非负性,将非均匀传感器阵列的宽带信号波达方向估计所对应的复数问题转化为实数问题,降低了算法的复杂度;
3、可确保欠定条件下信号源个数较多时的DOA估计精度;
4、可确保信噪比低、采样点数少的情况下DOA估计精度;
5、本发明所提出的非均匀传感器阵列的宽带信号波达方向估计方法,不需要预设任何经验参数,且其迭代算法的初值完全由传感器阵列输出数据得来,从而提高了DOA估计精度。
附图说明
图1是非均匀传感器阵列宽带信号波达方向估计方法的流程图;
图2是在SNR=-16dB,频域快拍数T=100条件下的空间谱图;
图3是在频域快拍数T=100条件下的RMSE随着SNR变化的趋势图;
图4是在SNR=-16dB条件下的RMSE随着频域快拍数变化的趋势图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明,以使本领域的技术人员可以更好地理解本发明并能予以实施,但所举实施例不作为对本发明的限定。
实施例
一种非均匀传感器阵列宽带信号波达方向估计方法,参照图1,仿真产生中心频率为100Hz,带宽40Hz宽带声信号源(声信号的传播速度为340m/s)。采用二级嵌套声传感器阵列作为非均匀阵列的例子,包括6个传感器,分别位于[0,d,2d,3d,7d,11d],其中d=1.7m。共生成9个覆盖相同频域范围且不相关的宽带信号源,从远场-69°、-45°、-30°、-15°、0°、12°、24°、42°、65°这9个不同的方向入射到该非均匀阵列上,各传感器接收信号时叠加与信号不相关的高斯白噪声,最后阵列输出采样率为240Hz的数据。
对各传感器通道输出数据样本序列作滑窗FFT以得到L=8个窄频带上的阵列输出快拍矢量序列:对一小段非均匀阵列采样数据做滑窗FFT,所加汉明窗的窗长以及FFT点数为NFFT=64,相邻窗交迭Noverlap=32点;根据目标宽带信号所覆盖的频域范围80~120Hz,在每个FFT窗口所得到的NFFT=64个频域值中选择86.25Hz、90Hz、93.75Hz、97.5Hz、101.25Hz、105Hz、108.75Hz和112.5Hz共L=8个频带作为目标子频带数据以供后续使用,对该段数据共做T个滑窗FFT,从而构造L=8个子频带上每个子频带T个“频域快拍”的数据向量:xl(t),l=1,2,…,L;t=1,2,…,T。
利用L个窄频带上的阵列输出快拍矢量序列构造L个窄频带上的实值加权样本协方差矢量:首先,计算每个频点对应的样本协方差矩阵其中(·)H表示共轭转置;然后,对作向量化得到样本协方差矢量其次,用加权矩阵J对进行加权处理以去除噪声项,得到加权样本协方差矢量其中权重矩阵J=[J1,J2,…JM-1]T,Jm=[em(M+1)-M+1,em(M+1)-M+2,…,em(M+1)],M2×1维向量em除第m个元素为1之外其他所有元素都为0,其中(·)T表示矩阵转置;利用J与构造加权样本协方差矢量的协方差矩阵最后,构造实值加权样本协方差矢量其中Re(·)与Im(·)分别表示取实部与取虚部。
构造L个窄频带上与各实值加权样本协方差矢量相对应的过完备阵列流型矩阵:首先,以1°的角度网格间距构造可以覆盖阵列探测角度范围的密集角度抽样网格然后,对于角度构造L个窄频带上与实值加权样本协方差矢量相对应的阵列导向矢量其中 为Kronecker积,原始阵列导向矢量fl为第l个频带的频率,dm为第m个传感器相对于参考传感器的间距,v为波速,(·)*表示求共轭;最后,对于网格中的每一个角度计算并由此构造L个窄频带上与各实值加权样本协方差矢量相对应的过完备阵列流型矩阵Φl:
设置迭代程序的初值以及固定参数:首先,计算空间联合稀疏表示向量γ的初值其中⊙为Hadamard积,||·||F为Frobenius范数;然后,计算其中进而计算空间联合稀疏表示向量λ的初值其中||·||2为范数;设置固定参数值c=d=10-6。
更新L个窄频带上的隐变量向量wl所服从的非负高斯分布的均值与协方差矩阵:首先,利用上一轮迭代所得的γ,更新L个窄频带上的隐变量向量wl所服从的非负高斯分布的协方差矩阵其中Γ=diag(γ);然后,更新L个窄频带上的隐变量向量wl所服从的非负高斯分布的均值
迭代更新空间联合稀疏表示向量γ与正则化参数λ的值:首先,计算L个窄频带上的隐变量向量wl中各元素的二阶矩:
其中特别当时,直接令然后,利用计算所得的以及上一轮迭代所更新的正则化参数λ,更新空间联合稀疏表示向量γ的第n个元素最后,利用上一轮迭代所更新的空间联合稀疏表示向量γ,更新正则化参数λ:
计算是否满足本实施例设置的迭代终止条件(最大迭代次数Pmax达到50),若不满足继续迭代;若满足,终止迭代,搜索当前更新所得的空间联合稀疏表示向量γ的所有峰值,其峰值位置与网格中对应的角度即为宽带信号的波达方向估计值。
参照图2,在信噪比(signal-to-noise ratio,SNR)为-16dB,频域快拍数T=100时本发明所提方法的空间谱图,其中各垂直虚线标出了真实DOA的位置,可见在欠定条件(信号源个数大于传感器个数)以及信噪比非常低的情况下该方法仍能通过搜索空间谱的峰值较为精确地对各信号进行定位,同时所提方法中所隐含的层级稀疏贝叶斯分布很好地抑制了空间谱伪峰的出现。
参照图3和图4,以信噪比和快拍数为变量,以均方根差(root-mean-suqare-error,RMSE)评估本发明所提宽带DOA估计方法的精度,图3为T=100条件下RMSE随着SNR变化的趋势图,图4为SNR=-5dB条件下RMSE随着频域快拍数变化的趋势图。在欠定条件下信号源个数较多时,以及信噪比低、采样点数少的情况下,本发明仍能达到接近于克拉美罗界(所有估计方法的误差都不会小于该下界)的估计效果。
此外,本发明所提方法所涉及到的每次迭代的计算复杂度约为O(L·min{M3(M-1)3,N3}),同时该方法能快速收敛,从而相比于其他稀疏重构宽带DOA估计方法在计算复杂度上要小得多。
本发明所涉及迭代计算的理论推导过程:
以估计一维波达方向的线型阵列为例,假设某非均匀分布线型传感器阵列共有M个全向传感器,分别位于0,d1,…,dM-1位置,其中dm-1为第m个传感器距参考传感器的距离。假设有K个不相关的远场宽带信号源,分别从K个方向θk,k=1,2,…,K入射。利用FFT或者窄带滤波器将阵列输出采样数据划分为L个窄带,每个窄带上的“频域快拍”数为T
xl(t)=Alsl(t)+nl(t),l=1,2,…,L;t=1,2,…,T, (1)
其中xl(t),sl(t),nl(t),分别表示第l个子频带上的阵列输出向量、源信号向量和加性噪声向量。Al为第l个子频带(fl频点)上的导向矢量矩阵:Al=[a(fl,θ1),a(fl,θ2),…,a(fl,θK)],其中阵列导向矢量v为波速。假设信号源与噪声不相关,则xl(t)的协方差矩阵可以写作
其中E{·}表示期望,是K个信号源在fl频点上的方差,Ql是fl频点上的对角的噪声协方差矩阵。
对Rl进行矢量化以构建协方差矢量,得到
式中vec(·)为矢量化运算,符号表示Khatri-Rao积,νl=vec(Ql)。的列数反映了阵列的自由度。由于νl的非零项仅存在于特定位置,加权矩阵J可用来消除噪声,并得到
其中J=[J1,J2,…JM-1]T,Jm=[em(M+1)-M+1,em(M+1)-M+2,…,em(M+1)]。
在实际应用中由于快拍数有限,我们只能得到样本协方差矩阵当源信号服从圆对称高斯分布时,样本协方差矩阵分布满足渐进复高斯分布,即
这里表示Kronecker积,利用加权矩阵J对进行加权,得到
其中又因为pl为非负实向量,(6)式可转换为如下的实值高斯分布形式
且
若以稀疏贝叶斯学习(SBL)思想进行DOA估计,设置密集角度抽样网格覆盖整个入射方向的范围,(7)式中的条件分布借助抽样网格可转化为稀疏表示形式
其中Φl是的过完备形式,列为其中wl是非负稀疏方差向量,它的非零项的位置与抽样网格中实际的DOA相对应。当l≠m,yl|wl与ym|wm相互独立,故可进一步得
其中wl,l=1,2,…,L具有相同的稀疏支撑集,因此假设这些向量独立同分布且满足n维非负高斯分布,其联合分布为
其中表示非负高斯概率密度函数。
空间联合稀疏表示向量γ决定了与向量wl,l=1,2,…,L相同的稀疏支撑集,这里我们认为空间联合稀疏表示向量γ是一个随机向量,并假定所有的γ[n]独立同分布且满足Gamma先验分布
将(11)式与(10)式结合,并利用广义逆高斯(GIG)概率密度函数的可积性和第二类修正贝塞尔函数可以得出:
由(12)式可知,该模型进行稀疏约束的方式与广泛应用的BPDN方法中的范数的作用相似。
最后假设λ具有伽马超先验分布
p(λ;c,d)=g(λ|c,d), (13)
其中超参数c,d被设置为常量,且应有c,d→0。
贝叶斯准则下,可推导wl的后验分布为非负高斯分布,且m≠l时wm|(ym,γ,λ)与wl|(yl,γ,λ)无关,则有
其中
基于期望最大化(EM)思想,将wl,l=1,2,,…,L视为隐藏变量,可推导空间联合稀疏表示向量γ和正则化参数λ的更新公式:
(1)推导γ[n]的更新公式,忽略p(y1,y2,…,yL,w1,w2,…,wL,γ,λ)中的γ[n],从而最大化
对(15)式中的γ[n]求导,令导数为0时可求得:
其中服从后验分布p(wl|yl,γ,λ)的wl的第n个元素的二阶矩为
(2)推导正则化参数λ的更新公式,忽略其他无关项,只需要将log p(γ|λ)p(λ)最大化,可得
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例,本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换,均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。
Claims (10)
1.一种非均匀传感器阵列宽带信号波达方向估计方法,其特征在于,包括以下步骤:
S1、对非均匀传感器阵列中各传感器接收通道采集的样本序列作滑窗快速傅里叶变换,并构造L个窄频带上的实值加权样本协方差矢量yl;
S2、构造与各窄频带上的实值加权样本协方差矢量yl相对应的过完备阵列流型矩阵Φl;
S3、利用L个窄频带上的实值加权样本协方差矢量yl和过完备阵列流型矩阵Φl设置迭代程序中空间联合稀疏表示向量γ的初始值γinit,并根据实值加权样本协方差矢量yl、过完备阵列流型矩阵Φl和空间联合稀疏表示向量的初始值γinit计算正则化参数λ的初始值λinit;
S4、利用L个窄频带上的实值加权样本协方差矢量yl和过完备阵列流型矩阵Φl更新空间联合稀疏表示向量γ和正则化参数λ;
S5、当第i次迭代得到的空间联合稀疏表示向量γi与第i+1次迭代得到的空间联合稀疏表示向量γi+1之间满足||γi-γi+1||2/||γi||2<α或者达到了最大迭代次数Pmax,终止迭代;否则继续返回S4迭代;其中α为一个任意趋于0的小数;
S6、当满足迭代终止条件时,搜索当前更新所得的空间联合稀疏表示向量γ的所有峰值,并根据空间联合稀疏表示向量γ的峰值位置输出宽带信号的波达方向估计值。
2.如权利要求1所述的非均匀传感器阵列宽带信号波达方向估计方法,其特征在于,步骤S1包括:
S11、对非均匀传感器阵列中各传感器接收通道采集的样本序列作滑窗快速傅里叶变换,得到输出快拍矢量序列xl(t),l=1,2,…,L;t=1,2,…,T;
S12、利用输出快拍矢量序列构造实值加权样本协方差矢量yl,l=1,2,…,L。
3.如权利要求1所述的非均匀传感器阵列宽带信号波达方向估计方法,其特征在于,步骤S4包括:
S41、利用实值加权样本协方差矢量yl、过完备阵列流型矩阵Φl以及上一轮迭代更新的空间联合稀疏表示向量γ,更新隐变量向量wl所服从的非负高斯分布的均值向量μl与协方差矩阵Σl,其中,若当前为第一轮迭代则空间联合稀疏表示向量γ=γinit;
S42、利用隐变量向量wl所服从的非负高斯分布的均值μl、协方差矩阵Σl和上一轮迭代更新的正则化参数λ,更新空间联合稀疏表示向量γ的值;其中,若当前为第一轮迭代则正则化参数λ=λinit;利用迭代更新的空间联合稀疏表示向量γ,更新正则化参数λ的值。
4.如权利要求2所述的非均匀传感器阵列宽带信号波达方向估计方法,其特征在于,步骤S11包括:
S111、利用窗函数对非均匀传感器阵列中各传感器接收通道采集的样本序列进行滑动数据选取,其中阵元个数为M;
S112、对各窗口中的加权非均匀阵列采样数据做快速傅里叶变换,其中,所加窗函数的窗长以及快速傅里叶变换点数为NFFT,相邻窗交迭Noverlap点;
S113、根据目标宽带信号所覆盖的频域范围,在第t个快速傅里叶变换窗口所得到的NFFT个频域值中选择L个作为待使用频带上的值,由于有M个传感器接收通道,可构造L个子频带上M×1维向量:xl(t),l=1,2,…,L;t=1,2,…,T,其中每个向量是该频带上的“频域快拍”;
S114、对非均匀阵列采样数据共做T个滑窗快速傅里叶变换,可构造L个子频带上每个频带上T个“频域快拍”的数据向量:xl(t),l=1,2,…,L;t=1,2,…,T。
5.如权利要求2所述的非均匀传感器阵列宽带信号波达方向估计方法,其特征在于,步骤S12包括:
S121、计算L个窄频带上样本协方差矩阵其中(·)H表示共轭转置;
S122、对矩阵作向量化得到样本协方差矢量
S123、用加权矩阵J对进行加权处理以去除噪声项,得到加权样本协方差矢量其中权重矩阵J=[J1,J2,…JM-1]T,Jm=[em(M+1)-M+1,em(M+1)-M+2,…,em(M+1)],M2×1维向量em除第m个元素为1之外其他所有元素都为0,其中(·)T表示矩阵转置;
S124、利用J与构造加权样本协方差矢量的协方差矩阵
S125、构造实值加权样本协方差矢量其中Re(·)与Im(·)分别表示取实部与取虚部。
6.如权利要求1所述的非均匀传感器阵列宽带信号波达方向估计方法,其特征在于,步骤S2包括:
S21、以一定的角度网格间距构造可以覆盖阵列探测角度范围的密集角度抽样网格
S22、对于角度构造L个窄频带上与实值加权样本协方差矢量相对应的阵列导向矢量其中 为Kronecker积,原始阵列导向矢量fl为第l个频带的频率,dm为第m个传感器相对于参考传感器的间距,v为波速,(·)*表示求共轭;
S23、对于网格中的每一个角度计算并由此构造L个窄频带上与各实值加权样本协方差矢量相对应的过完备阵列流型矩阵Φl:
7.如权利要求1所述的非均匀传感器阵列宽带信号波达方向估计方法,其特征在于,步骤S3包括:
S31、计算空间联合稀疏表示向量γ的初值其中⊙为Hadamard积,||·||F为Frobenius范数;
S32、计算其中进而计算正则化参数λ的初值其中||·||2为范数;
S33、设置正则化参数λ的无信息超先验分布参数固定参数值c,d,其中c,d→0。
8.如权利要求3所述的非均匀传感器阵列宽带信号波达方向估计方法,其特征在于,步骤S41包括:
S411、利用上一轮迭代所得的空间联合稀疏表示向量γ,更新L个窄频带上的隐变量向量wl所服从的非负高斯分布的协方差矩阵其中Γ=diag(γ);
S412、更新L个窄频带上的隐变量向量wl所服从的非负高斯分布的均值
9.如权利要求3所述的非均匀传感器阵列宽带信号波达方向估计方法,其特征在于,步骤S42包括:
S421、计算L个窄频带上的隐变量向量wl中各元素的二阶矩:
其中特别当时,直接令
S422、利用计算所得的以及上一轮迭代所更新的正则化参数λ,更新空间联合稀疏表示向量γ的第n个元素
S423、利用上一轮迭代所更新的空间联合稀疏表示向量γ,更新正则化参数λ:
10.如权利要求6所述的非均匀传感器阵列宽带信号波达方向估计方法,其特征在于,步骤S6包括,当满足迭代终止条件时,搜索当前更新所得的空间联合稀疏表示向量γ的所有峰值,其峰值位置与网格中对应的角度即为宽带信号的波达方向估计值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811178910.3A CN109407045B (zh) | 2018-10-10 | 2018-10-10 | 一种非均匀传感器阵列宽带信号波达方向估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811178910.3A CN109407045B (zh) | 2018-10-10 | 2018-10-10 | 一种非均匀传感器阵列宽带信号波达方向估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109407045A true CN109407045A (zh) | 2019-03-01 |
CN109407045B CN109407045B (zh) | 2020-05-22 |
Family
ID=65466930
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811178910.3A Active CN109407045B (zh) | 2018-10-10 | 2018-10-10 | 一种非均匀传感器阵列宽带信号波达方向估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109407045B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110045323A (zh) * | 2019-03-14 | 2019-07-23 | 电子科技大学 | 一种基于矩阵填充的互质阵稳健自适应波束形成算法 |
CN110109050A (zh) * | 2019-04-08 | 2019-08-09 | 电子科技大学 | 嵌套阵列下基于稀疏贝叶斯的未知互耦的doa估计方法 |
CN110457748A (zh) * | 2019-07-04 | 2019-11-15 | 中国人民解放军63892部队 | 一种等水平的二部覆盖阵的试验设计方法 |
CN110579737A (zh) * | 2019-07-17 | 2019-12-17 | 电子科技大学 | 一种杂波环境中基于稀疏阵列的mimo雷达宽带doa计算方法 |
CN111505567A (zh) * | 2020-03-25 | 2020-08-07 | 山东省科学院海洋仪器仪表研究所 | 基于单矢量水听器波达方向估计的多目标跟踪方法 |
CN111929638A (zh) * | 2020-09-24 | 2020-11-13 | 成都启英泰伦科技有限公司 | 一种语音波达方向估计方法及装置 |
CN111983553A (zh) * | 2020-08-20 | 2020-11-24 | 上海无线电设备研究所 | 一种基于互质多载频稀疏阵列的无网格doa估计方法 |
CN112198508A (zh) * | 2020-10-29 | 2021-01-08 | 中国人民武装警察部队工程大学 | 一种基于支撑集约束的雷达目标成像与识别方法 |
CN112487703A (zh) * | 2020-11-09 | 2021-03-12 | 南京信息工程大学滨江学院 | 基于稀疏贝叶斯在未知噪声场的欠定宽带信号doa估计方法 |
CN113406560A (zh) * | 2021-05-20 | 2021-09-17 | 华南理工大学 | 一种非相干分布宽带源的角度和频率参数估计方法 |
CN113411142A (zh) * | 2021-06-29 | 2021-09-17 | 齐鲁工业大学 | 宽带确定性最大似然测向方法、***、终端及存储介质 |
CN113589255A (zh) * | 2021-08-23 | 2021-11-02 | 武汉大学 | 一种基于多频联合稀疏贝叶斯学习的到达角估计方法 |
CN115656998A (zh) * | 2022-11-22 | 2023-01-31 | 中国人民解放军空军预警学院 | 一种低样本数下阵列信号自适应检测方法与*** |
CN115792792A (zh) * | 2022-11-22 | 2023-03-14 | 中国人民解放军空军预警学院 | 极低样本下基于协方差矩阵迭代估计的测向方法与*** |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103954950A (zh) * | 2014-04-25 | 2014-07-30 | 西安电子科技大学 | 一种基于样本协方差矩阵稀疏性的波达方向估计方法 |
-
2018
- 2018-10-10 CN CN201811178910.3A patent/CN109407045B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103954950A (zh) * | 2014-04-25 | 2014-07-30 | 西安电子科技大学 | 一种基于样本协方差矩阵稀疏性的波达方向估计方法 |
Non-Patent Citations (5)
Title |
---|
LIN ZHANG ET AL.: "Wideband DOA estination in spherical harmonic domain using sparse bayesian learning", 《PROCEEDINGS OF THE 36TH CHINESE CONTROL CONFERENCE》 * |
NAN HU ET AL.: "Wideband DOA estimation from the sparse recovery perspective for the spatial-only modeling of array data", 《SIGNALPROCESSING》 * |
何振清 等: "宽带协方差矩阵的多字典联合稀疏表示DOA估计", 《信号处理》 * |
殷允杰: "基于非均匀阵列的波达方向估计技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
赵永红 等: "一种新的基于稀疏表示的宽带信号DOA估计方法", 《电子与信息学报》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110045323A (zh) * | 2019-03-14 | 2019-07-23 | 电子科技大学 | 一种基于矩阵填充的互质阵稳健自适应波束形成算法 |
CN110109050B (zh) * | 2019-04-08 | 2022-05-03 | 电子科技大学 | 嵌套阵列下基于稀疏贝叶斯的未知互耦的doa估计方法 |
CN110109050A (zh) * | 2019-04-08 | 2019-08-09 | 电子科技大学 | 嵌套阵列下基于稀疏贝叶斯的未知互耦的doa估计方法 |
CN110457748A (zh) * | 2019-07-04 | 2019-11-15 | 中国人民解放军63892部队 | 一种等水平的二部覆盖阵的试验设计方法 |
CN110579737A (zh) * | 2019-07-17 | 2019-12-17 | 电子科技大学 | 一种杂波环境中基于稀疏阵列的mimo雷达宽带doa计算方法 |
CN110579737B (zh) * | 2019-07-17 | 2022-10-11 | 电子科技大学 | 一种杂波环境中基于稀疏阵列的mimo雷达宽带doa计算方法 |
CN111505567A (zh) * | 2020-03-25 | 2020-08-07 | 山东省科学院海洋仪器仪表研究所 | 基于单矢量水听器波达方向估计的多目标跟踪方法 |
CN111505567B (zh) * | 2020-03-25 | 2022-09-06 | 山东省科学院海洋仪器仪表研究所 | 基于单矢量水听器波达方向估计的多目标跟踪方法 |
CN111983553A (zh) * | 2020-08-20 | 2020-11-24 | 上海无线电设备研究所 | 一种基于互质多载频稀疏阵列的无网格doa估计方法 |
CN111983553B (zh) * | 2020-08-20 | 2024-02-20 | 上海无线电设备研究所 | 一种基于互质多载频稀疏阵列的无网格doa估计方法 |
CN111929638A (zh) * | 2020-09-24 | 2020-11-13 | 成都启英泰伦科技有限公司 | 一种语音波达方向估计方法及装置 |
CN112198508A (zh) * | 2020-10-29 | 2021-01-08 | 中国人民武装警察部队工程大学 | 一种基于支撑集约束的雷达目标成像与识别方法 |
CN112487703A (zh) * | 2020-11-09 | 2021-03-12 | 南京信息工程大学滨江学院 | 基于稀疏贝叶斯在未知噪声场的欠定宽带信号doa估计方法 |
CN112487703B (zh) * | 2020-11-09 | 2024-05-28 | 南京信息工程大学滨江学院 | 基于稀疏贝叶斯在未知噪声场的欠定宽带信号doa估计方法 |
CN113406560A (zh) * | 2021-05-20 | 2021-09-17 | 华南理工大学 | 一种非相干分布宽带源的角度和频率参数估计方法 |
CN113411142A (zh) * | 2021-06-29 | 2021-09-17 | 齐鲁工业大学 | 宽带确定性最大似然测向方法、***、终端及存储介质 |
CN113411142B (zh) * | 2021-06-29 | 2022-09-06 | 齐鲁工业大学 | 宽带确定性最大似然测向方法、***、终端及存储介质 |
CN113589255A (zh) * | 2021-08-23 | 2021-11-02 | 武汉大学 | 一种基于多频联合稀疏贝叶斯学习的到达角估计方法 |
CN113589255B (zh) * | 2021-08-23 | 2023-08-01 | 武汉大学 | 一种基于多频联合稀疏贝叶斯学习的到达角估计方法 |
CN115656998A (zh) * | 2022-11-22 | 2023-01-31 | 中国人民解放军空军预警学院 | 一种低样本数下阵列信号自适应检测方法与*** |
CN115792792A (zh) * | 2022-11-22 | 2023-03-14 | 中国人民解放军空军预警学院 | 极低样本下基于协方差矩阵迭代估计的测向方法与*** |
CN115656998B (zh) * | 2022-11-22 | 2023-06-13 | 中国人民解放军空军预警学院 | 一种低样本数下阵列信号自适应检测方法与*** |
CN115792792B (zh) * | 2022-11-22 | 2023-10-27 | 中国人民解放军空军预警学院 | 极低样本下基于协方差矩阵迭代估计的测向方法与*** |
Also Published As
Publication number | Publication date |
---|---|
CN109407045B (zh) | 2020-05-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109407045A (zh) | 一种非均匀传感器阵列宽带信号波达方向估计方法 | |
Gu et al. | Joint 2-D DOA estimation via sparse L-shaped array | |
CN106021637B (zh) | 互质阵列中基于迭代稀疏重构的doa估计方法 | |
CN108375763B (zh) | 一种应用于多声源环境的分频定位方法 | |
CN110045323B (zh) | 一种基于矩阵填充的互质阵稳健自适应波束形成算法 | |
CN110113085B (zh) | 一种基于协方差矩阵重构的波束形成方法及*** | |
CN109298383B (zh) | 一种基于变分贝叶斯推断的互质阵波达方向角估计方法 | |
CN109407046A (zh) | 一种基于变分贝叶斯推断的嵌套阵列波达方向角估计方法 | |
CN109143153A (zh) | 一种基于稀疏重构的超嵌套阵列波达方向估计方法 | |
CN109490850A (zh) | 主瓣干扰下宽带阵列自适应波束形成方法 | |
CN111273229B (zh) | 基于低秩矩阵重建的水声宽频散射源的定位方法 | |
CN110109058A (zh) | 一种平面阵列反卷积声源识别方法 | |
CN112487703B (zh) | 基于稀疏贝叶斯在未知噪声场的欠定宽带信号doa估计方法 | |
CN106093920B (zh) | 一种基于对角加载的自适应波束形成算法 | |
CN104950297A (zh) | 基于矩阵1范数拟合的阵元误差估计方法 | |
CN109600152A (zh) | 一种基于子空间基变换的自适应波束形成方法 | |
CN110221249A (zh) | 基于压缩感知的宽带声源定位方法 | |
CN107290732A (zh) | 一种量子大***的单基地mimo雷达测向方法 | |
CN106788655A (zh) | 互耦条件下未知互耦信息的干扰相干稳健波束形成方法 | |
Poursheikhali et al. | Source localization in inhomogeneous underwater medium using sensor arrays: Received signal strength approach | |
CN103837858A (zh) | 一种用于平面阵列的远场波达角估计方法及*** | |
CN108416105B (zh) | 脉冲和高斯噪声下稳健的自适应波束形成算法 | |
CN104407319A (zh) | 阵列信号的目标源测向方法和*** | |
Tan et al. | Efficient FFT based multi source DOA estimation for ULA | |
Suo et al. | Single acoustic vector sensor DOA enhanced by unsaturated bistable stochastic resonance with linear amplitude response constrained |
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 |