CN102621527A - 基于数据重构的宽带相干源的方位估计方法 - Google Patents
基于数据重构的宽带相干源的方位估计方法 Download PDFInfo
- Publication number
- CN102621527A CN102621527A CN201210073516XA CN201210073516A CN102621527A CN 102621527 A CN102621527 A CN 102621527A CN 201210073516X A CN201210073516X A CN 201210073516XA CN 201210073516 A CN201210073516 A CN 201210073516A CN 102621527 A CN102621527 A CN 102621527A
- Authority
- CN
- China
- Prior art keywords
- data
- subband
- covariance matrix
- frequency
- coherent source
- 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
Images
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明提供的是一种基于数据重构的宽带相干源的方位估计方法。(1)对宽带阵列接收数据进行子带分解;(2)选取宽带信号的最低频率作为聚焦频率,利用空间数据重构公式将各子带数据聚焦到同一频率;(3)计算聚焦后的各子带数据的协方差矩阵,利用矩阵共轭重排原理改进协方差矩阵,将改进后的各子带数据的协方差矩阵取均值,作为最终的协方差矩阵;(4)将最终的协方差矩阵利用MVDR算法进行空间谱估计,得到宽带相干源的方位信息。本发明不仅在低信噪比和少快拍数下有更高的角度分辨率,并且无需对角加载等处理即可克服高信噪比条件下MVDR算法的不稳定性,具有更优越的方位估计性能。
Description
技术领域
本发明涉及的是一种阵列信号处理方法。具体地说是一种目标方位估计方法。
背景技术
基于高分辨算法的阵列信号处理技术一直是科研的热点,而由于海洋环境的变化性以及水声信道的复杂性,许多经典的子空间类高分辨算法失去了其原有的优势,可靠性下降,实用性变差。最小方差无畸变波束形成(MVDR)算法能够提供最佳的信号保护、干扰消除和噪声降低能力,更适用于海洋环境。在实际的海洋环境中,声纳接收宽带信号的情况越来越多,宽带信号因其具有目标回波携带信息量大、混响背景相关性弱、有利于目标检测、参量估计和目标特征提取等特点,也在阵列信号处理领域中扮演日益重要的角色。随着科技的发展,宽带信号在移动通信、无线网路等和我们生活密切相关的多个方面中发挥着重要的作用。
传统的宽带相干源解相干技术需要进行方向预估,其性能受预估方位的影响,因此空间重采样等方法被陆续提出。文献[1][2]介绍了无需方向预估的宽带信号处理方法,文献[3]给出了一种适用于MVDR算法的更优越的子带分解方法。本发明结合文献[1][2][3]的优点并予以改进,可以获得更好的宽带相干源方位估计性能。
与本发明相关的参考文献包括:
[1]朱维杰,孙进才.基于阵列接收数据重采样的频率不变波束形成器[J].自然科学进展,2002,12(6):669-672;
[2]朱维杰,孙进才.基于快速傅里叶变换插值的宽带信号方位估计[J].声学学报,2002,27(6):513-517;
[3]Biao Jiang,Ye Zhu,Changyu Sun.DOA estimation of wideband sources usinga constant-Q subband coherent MVDR method[C]//7th International Conference on SignalProcessing Proceedings.Beijing,China,2004:2437-2440。
发明内容
本发明的目的在于提供一种能够在高信噪比和小快拍数下具有稳定性能的基于数据重构的宽带相干源的方位估计方法。
本发明的目的是这样实现的:
(1)利用恒定相对带宽方法对宽带阵列接收数据进行子带分解,并与FFT分段叠加子带分解进行比较;
(2)选取宽带信号的最低频率作为聚焦频率,利用空间数据重构公式将各子带数据聚焦到同一频率;
(3)计算聚焦后的各子带数据的协方差矩阵,利用矩阵共轭重排原理改进协方差矩阵,将改进后的各子带数据的协方差矩阵取均值,作为最终的协方差矩阵;
(4)将最终的协方差矩阵利用MVDR算法进行空间谱估计,得到宽带相干源的方位信息。
本发明运用恒定相对带宽方法将接收数据进行子带分解,然后将时域的数据重构引入空间阵元域数据,结合矩阵共轭重排方法,提出一种基于数据重构的宽带相干源MVDR算法,不仅在低信噪比和少快拍数下有更高的角度分辨率,并且无需对角加载等处理即可克服高信噪比条件下MVDR算法的不稳定性,具有更优越的方位估计性能。
附图说明
图1为16元均匀线列阵,两个相干入射信号的中心频率f0=800Hz,相对带宽为50%,入射角度分别为0°、8°,采样频率Fs=10f0,快拍数为1000,信噪比为15dB时,采用FFT分段叠加子带分解方法,选取20个子带,基于FFT插值法获得的MVDR方位谱图。
图2为16元均匀线列阵,两个相干入射信号的中心频率f0=800Hz,相对带宽为50%,入射角度分别为0°、8°,采样频率Fs=10f0,快拍数为1000,信噪比为15dB时,采用恒定相对带宽子带分解方法,选取12个子带,基于FFT插值法获得的MVDR方位谱图。
图3为16元均匀线列阵,两个相干入射信号的中心频率f0=800Hz,相对带宽为50%,入射角度分别为0°、8°,采样频率Fs=10f0,快拍数为1000,信噪比为10dB时,基于FFT插值法获得的MVDR方位谱图。
图4为相同条件下,信噪比为50dB时,基于FFT插值法获得的MVDR方位谱图。
图5为入射角度分别为0°、6°,信噪比为10dB,快拍数为1000时,基于FFT插值法获得的MVDR方位谱图。
图6为入射角度分别为0°、6°,信噪比为10dB,快拍数为1000时,基于数据重构法获得的MVDR方位谱图。
图7为入射角度分别为0°、6°,信噪比为50dB,快拍数为1000时,基于FFT插值法获得的MVDR方位谱图。
图8为入射角度分别为0°、6°,信噪比为50dB,快拍数为1000时,基于数据重构法获得的MVDR方位谱图。
图9为入射角度分别为0°、6°,信噪比为20dB,快拍数为500时,基于FFT插值法获得的MVDR方位谱图。
图10为入射角度分别为0°、6°,信噪比为20dB,快拍数为500时,基于数据重构法获得的MVDR方位谱图。
图11为该发明方法的整体流程框图。
具体实施方式
下面结合附图举例对本发明做更详细地描述:
1.设有M个阵元组成的阵列,接收P个位于远场的宽带平稳随机信号,当存在加性噪声时,第m个阵元接收到的信号可以表示为:
其中m=1,2·,M,sp(t)表示第p个信号源,nm(t)表示第m个阵元上的加性噪声,θp表示第p个信号的方位角,τm(θp)表示第m个阵元接收的第p个信号(相对于参考点)的传播延迟.在均匀线列阵中τm(θp)=(m-1)d·sin(θp)/c,d为阵元间距,c为声速。
设接收数据x(t)=[x1(t),x2(t),·,xM(t)]T经离散化后变成M×N的矩阵,N为数据长度,采样频率为fs,将其进行N点FFT得到M×N频域数据X(f)。可表示为:
X(f)=A(f,Θ)S(f)+N(f) (2)
式中:
X(f)=[X1(f),X2(f),·,XM(f)]T,
A(f,Θ)=[a(f,θ1),a(f,θ2),·,a(f,θp)],
S(f)=[S1(f),S2(f),·,SP(f)]T,
N(f)=[N1(f),N2(f),·,NM(f)]T,
其中X(f),S(f),N(f)分别为观测矢量、信号矢量和噪声矢量,Xm(f),Sm(f),Nm(f)分别为xm(t),sp(t),nm(t)的频域形式,A(f,Θ)为阵列流形,a(f,θp)为方向矢量,T表示转置。
2宽带信号子带分解
2.1FFT分段叠加子带分解法
实际应用中,由于一次FFT的结果存在很大误差,因此不能将接收数据直接进行FFT,而是采用分段叠加的方法,以减少FFT运算引起的误差,具体方法如下:
①设接收数据x(t)=[x1(t),x2(t),·,xM(t)]T经离散化后变成M×N的矩阵,M×N表示M行N列,N为数据长度,采样频率为fs,将其平均分为K段,每段数据长度为N1点,相邻两段数据重叠50%,则K即为需要独立处理的频段数目,并满足:
其中round[.]表示取整运算。也就是说,接收数据最终需要进行K次的循环叠加过程。
②对每段数据进行N1点FFT,则所得到的频域数据为M×N1的矩阵,M×N1表示M行N1列,设信号最低频率为fl,最高频率为fh,因此,fl在频域数据中对应的位置为第fL个数据点,fh在频域数据中对应的位置为第fH个数据点。fL和fH分别满足:
其中round[.]表示取整运算,fs为采样频率。
③在宽带信号处理中,只有(fl,fh)区间的频域数据包含有用信息,因此,只需要在每段获得的频域数据中选取(fL,fH)区间的频点进行处理即可,且每一个频点对应的M×1数据代表了一个子带,M×1表示M行1列,则分解的子带数为:
Bn=fH-fL+1 (5)
④每一次循环过程中,求解Bn个子带数据的协方差矩阵并相加取均值,K次循环后再对K次的处理结果进行叠加取均值,作为最终的协方差矩阵.利用FFT分段叠加子带分解法得到的MVDR方位谱如图1所示。
2.2恒定相对带宽子带分解方法
利用FFT进行子带分解时,为了减少FFT运算引起的误差,需要进行频域平滑,所以增加了计算量,而利用恒定相对带宽的子带分解方法则无需频域平滑,有效地减少了计算量。
具体方法如下:
①设接收数据x(t)=[x1(t),x2(t),·,xM(t)]T经离散化后变成M×N的矩阵,N为数据长度,采样频率为fs,将其进行N点FFT得到M×N频域数据X(f)。
②设K为分解的子带数目,选取相对带宽值Q,则有:
Q=2(f(k+1)-f(k))/(f(k+1)+f(k)) (6)
从而,
f(k+1)=(2+Q)·f(k)/(2-Q) (7)
其中k=1,2,·,K,f(k)为第k个子带的下限频率,且令f(1)=fl,f(K+1)=fh。
③f(k)在频域数据中对应数据点为:
因此第k个子带对应于频域数据X(f)中M×[F(k),F(k+1)]的数据,round[.]表示取整运算。
④求解K个子带数据的协方差矩阵并相加取均值,作为最终的协方差矩阵。利用恒定相对带宽子带分解方法得到的MVDR方位谱如图2所示。
显然,恒定相对带宽子带分解方法的计算量远小于FFT分段叠加子带分解法,并且利用恒定相对带宽子带分解方法可以获得和FFT分段叠加子带分解法相当的方位估计性能。
3.选取信号的最低频率fl作为聚焦频率,则为使每个子带数据聚焦到同一频率fl上,Xk经重构后得到数据矩阵Xk′的采样间隔为:
T′s=Ts·fl/f(k)=d·fl/f(k) (9)
其中,f(k)为第k个子带的中心频率,k=1,2,·K,K为子带总数.所以Xk′的数据长度为:
M′=(M-1)d/Ts′=(M-1)f(k)/fl (10)
则空间数据重构后得到:
其中E1M×1=[1,1,·1]T,T表示转置,E21×N=[1,1,·1],T0=[0,d,·,(M-1)d],T=[0,Ts′,·,M′Ts′];ωc满足ωm<ωc<ωs-ωm,并且有ωs=2π/d,ωm为奈奎斯特频率。
4.选取Xk′中间的M行数据矩阵作为插值后的矩阵,然后计算其协方差矩阵R。
利用共轭重排的原理将协方差矩阵R进行如下处理,得到新的协方差矩阵Rx:
其中ii是除副对角线元素为1,其余元素均为0的方阵,R表示协方差矩阵R的共轭矩阵。
5.按照以上步骤求得K个子带的协方差矩阵,求其平均值,然后利用MVDR算法进行DOA估计,获得目标的方位谱。
其中fj为聚焦频率, 为聚焦频率下的方向矢量,T表示转置,H表示共轭转置。
按照以上5个步骤,即可求出不同条件下宽带相干源的MVDR方位谱,实现目标的方位估计,对比图3、图4,可以看出高信噪比下传统FFT插值法具有不稳定性;对比图5、图6,可以看出较低信噪比下本发明方法角度分辨率优于FFT插值法;对比图7、图8,可以看出高信噪比下本发明方法具有稳定性;对比图9、图10,可以看出少快拍数下本发明方法角度分辨率优于FFT插值法。
Claims (2)
1.一种基于数据重构的宽带相干源的方位估计方法,其特征是包括如下步骤:
(1)对宽带阵列接收数据进行子带分解;
(2)选取宽带信号的最低频率作为聚焦频率,利用空间数据重构公式将各子带数据聚焦到同一频率;
(3)计算聚焦后的各子带数据的协方差矩阵,利用矩阵共轭重排原理改进协方差矩阵,将改进后的各子带数据的协方差矩阵取均值,作为最终的协方差矩阵;
(4)将最终的协方差矩阵利用MVDR算法进行空间谱估计,得到宽带相干源的方位信息。
2.根据权利要求1所述的基于数据重构的宽带相干源的方位估计方法,其特征是所述对宽带阵列接收数据进行子带分解是利用恒定相对带宽子带分解方法,具体步骤如下:
①接收数据x(t)=[x1(t),x2(t),·,xM(t)]T经离散化后变成M×N的矩阵,N为数据长度,采样频率为fs,将其进行N点FFT得到M×N频域数据X(f);
②K为分解的子带数目,选取相对带宽值Q,则:
Q=2(f(k+1)-f(k))/(f(k+1)+f(k))
从而,
f(k+1)=(2+Q)·f(k)/(2-Q)
其中k=1,2,·,K,f(k)为第k个子带的下限频率,且令f(1)=fl,f(K+1)=fh;
③f(k)在频域数据中对应数据点为:
因此第k个子带对应于频域数据X(f)中M×[F(k),F(k+1)]的数据,round[.]表示取整运算;④求解K个子带数据的协方差矩阵并相加取均值,作为最终的协方差矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210073516.XA CN102621527B (zh) | 2012-03-20 | 2012-03-20 | 基于数据重构的宽带相干源的方位估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210073516.XA CN102621527B (zh) | 2012-03-20 | 2012-03-20 | 基于数据重构的宽带相干源的方位估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102621527A true CN102621527A (zh) | 2012-08-01 |
CN102621527B CN102621527B (zh) | 2014-06-11 |
Family
ID=46561556
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210073516.XA Expired - Fee Related CN102621527B (zh) | 2012-03-20 | 2012-03-20 | 基于数据重构的宽带相干源的方位估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102621527B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102882617A (zh) * | 2012-10-10 | 2013-01-16 | 上海师范大学 | 一种基于谱相关特征的频谱检测方法 |
CN102932034A (zh) * | 2012-10-31 | 2013-02-13 | 哈尔滨工程大学 | 一种快速的宽带相干源方位估计方法 |
CN103576137A (zh) * | 2013-09-27 | 2014-02-12 | 电子科技大学 | 一种基于成像策略的多传感器多目标定位方法 |
CN104537202A (zh) * | 2014-10-31 | 2015-04-22 | 哈尔滨工业大学深圳研究生院 | 基于卫星编队协作的空间天线阵列合成方法 |
CN105445723A (zh) * | 2015-12-22 | 2016-03-30 | 中国船舶重工集团公司第七一五研究所 | 一种基于背景噪声频谱起伏特性的宽带阵列信号处理方法 |
WO2016119388A1 (zh) * | 2015-01-30 | 2016-08-04 | 华为技术有限公司 | 一种基于语音信号构造聚焦协方差矩阵的方法及装置 |
CN107577872A (zh) * | 2017-09-05 | 2018-01-12 | 湖南时变通讯科技有限公司 | 一种频率不变波束形成器设计方法及装置 |
CN109725198A (zh) * | 2017-10-31 | 2019-05-07 | 北京普源精电科技有限公司 | 一种频谱仪的余辉生成方法及频谱仪 |
CN109799495A (zh) * | 2019-01-02 | 2019-05-24 | 东南大学 | 一种用于高保真阵列处理的宽带时延估计方法 |
CN113030983A (zh) * | 2021-03-17 | 2021-06-25 | 中国科学院声学研究所 | 一种基于测深侧扫声纳的近场逐点聚焦doa方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101470187A (zh) * | 2007-12-26 | 2009-07-01 | 中国科学院声学研究所 | 一种用于线列阵的高精度测向方法 |
-
2012
- 2012-03-20 CN CN201210073516.XA patent/CN102621527B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101470187A (zh) * | 2007-12-26 | 2009-07-01 | 中国科学院声学研究所 | 一种用于线列阵的高精度测向方法 |
Non-Patent Citations (6)
Title |
---|
BIAO JIANG 等: "DOA Estimation of Wideband Sources Using a Constant-Q Subband Coherent MVDR Method", 《IEEE》, 31 December 2004 (2004-12-31), pages 1 - 4 * |
丁婷 等: "降低多径信号相关性的修正DOA估计方法", 《通信技术》, 29 February 2008 (2008-02-29), pages 53 - 55 * |
何心怡 等: "基于子带分解的宽带波束域最小方差无畸变响应高分辨方位估计方法研究", 《声学学报》, vol. 29, no. 6, 30 November 2004 (2004-11-30), pages 533 - 538 * |
卞红雨 等: "宽带相干源MVDR算法性能研究", 《声学技术》, vol. 30, no. 3, 30 June 2011 (2011-06-30), pages 271 - 273 * |
蒋飚 等: "一种宽带高分辨MVDR有效算法研究", 《***工程与电子技术》, vol. 27, no. 7, 31 July 2005 (2005-07-31), pages 1186 - 1188 * |
郭娜 等: "基于矩阵重构算法的宽带相干信号方位估计", 《火力与指挥控制》, vol. 35, no. 7, 31 July 2010 (2010-07-31), pages 50 - 52 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102882617A (zh) * | 2012-10-10 | 2013-01-16 | 上海师范大学 | 一种基于谱相关特征的频谱检测方法 |
CN102932034A (zh) * | 2012-10-31 | 2013-02-13 | 哈尔滨工程大学 | 一种快速的宽带相干源方位估计方法 |
CN102932034B (zh) * | 2012-10-31 | 2014-09-17 | 哈尔滨工程大学 | 一种快速的宽带相干源方位估计方法 |
CN103576137A (zh) * | 2013-09-27 | 2014-02-12 | 电子科技大学 | 一种基于成像策略的多传感器多目标定位方法 |
CN104537202A (zh) * | 2014-10-31 | 2015-04-22 | 哈尔滨工业大学深圳研究生院 | 基于卫星编队协作的空间天线阵列合成方法 |
CN104537202B (zh) * | 2014-10-31 | 2017-12-22 | 哈尔滨工业大学深圳研究生院 | 基于卫星编队协作的空间天线阵列合成方法 |
WO2016119388A1 (zh) * | 2015-01-30 | 2016-08-04 | 华为技术有限公司 | 一种基于语音信号构造聚焦协方差矩阵的方法及装置 |
CN105445723A (zh) * | 2015-12-22 | 2016-03-30 | 中国船舶重工集团公司第七一五研究所 | 一种基于背景噪声频谱起伏特性的宽带阵列信号处理方法 |
CN107577872A (zh) * | 2017-09-05 | 2018-01-12 | 湖南时变通讯科技有限公司 | 一种频率不变波束形成器设计方法及装置 |
CN107577872B (zh) * | 2017-09-05 | 2021-02-19 | 湖南时变通讯科技有限公司 | 一种时域频率不变波束形成器设计方法及装置 |
CN109725198A (zh) * | 2017-10-31 | 2019-05-07 | 北京普源精电科技有限公司 | 一种频谱仪的余辉生成方法及频谱仪 |
CN109725198B (zh) * | 2017-10-31 | 2021-03-12 | 北京普源精电科技有限公司 | 一种频谱仪的余辉生成方法及频谱仪 |
CN109799495A (zh) * | 2019-01-02 | 2019-05-24 | 东南大学 | 一种用于高保真阵列处理的宽带时延估计方法 |
CN109799495B (zh) * | 2019-01-02 | 2022-06-17 | 东南大学 | 一种用于高保真阵列处理的宽带时延估计方法 |
CN113030983A (zh) * | 2021-03-17 | 2021-06-25 | 中国科学院声学研究所 | 一种基于测深侧扫声纳的近场逐点聚焦doa方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102621527B (zh) | 2014-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102621527B (zh) | 基于数据重构的宽带相干源的方位估计方法 | |
CN102608565B (zh) | 一种基于均匀圆阵列的波达方向估计方法 | |
CN103592642B (zh) | Mimo雷达波形的设计方法 | |
CN100466061C (zh) | 一种宽带波束形成方法和装置 | |
CN106788653A (zh) | 一种基于协方差矩阵重构的自适应波束形成方法 | |
CN102445691B (zh) | 一种多通道星载合成孔径雷达方位频谱稀疏重建方法 | |
CN103353596A (zh) | 基于压缩感知的波束空间域米波雷达测高方法 | |
CN102175989B (zh) | 一种非相干分布式信号二维波达角的测定方法 | |
CN109407055A (zh) | 基于多径利用的波束形成方法 | |
CN107966677B (zh) | 一种基于空间稀疏约束的圆阵模态域方位估计方法 | |
CN110515038B (zh) | 一种基于无人机-阵列的自适应无源定位装置及实现方法 | |
CN103546221B (zh) | 一种宽带相干信号波达角估计方法 | |
CN103713276B (zh) | 基于最小互熵谱分析的波达方向估计方法 | |
CN109490850A (zh) | 主瓣干扰下宽带阵列自适应波束形成方法 | |
CN104730513A (zh) | 一种分级子阵聚焦mvdr波束形成方法 | |
CN107064901A (zh) | 一种舰载高频地波超视距雷达的目标方位估计方法 | |
CN104811867A (zh) | 基于阵列虚拟扩展的麦克风阵列空域滤波方法 | |
CN103902830A (zh) | 一种圆环形阵列稳健旁瓣控制超指向性波束形成方法 | |
CN109521393A (zh) | 一种基于信号子空间旋转特性的波达方向估计算法 | |
Chi | Underwater Real-Time 3D Acoustical Imaging: Theory, Algorithm and System Design | |
CN102932034B (zh) | 一种快速的宽带相干源方位估计方法 | |
CN109491009A (zh) | 一种光纤组合阵及基于光纤组合阵的栅瓣抑制方法 | |
CN109061597B (zh) | 基于盲源分离与时频脊波域滤波的电离层杂波抑制方法 | |
CN115236589B (zh) | 一种基于协方差矩阵修正的极地冰下doa估计方法 | |
CN114563760B (zh) | 一种基于sca阵型的二阶超波束形成方法、设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140611 Termination date: 20200320 |
|
CF01 | Termination of patent right due to non-payment of annual fee |