CN102866383B - 一种基于空域自适应滤波的波达方向估计方法 - Google Patents
一种基于空域自适应滤波的波达方向估计方法 Download PDFInfo
- Publication number
- CN102866383B CN102866383B CN201210404468.8A CN201210404468A CN102866383B CN 102866383 B CN102866383 B CN 102866383B CN 201210404468 A CN201210404468 A CN 201210404468A CN 102866383 B CN102866383 B CN 102866383B
- Authority
- CN
- China
- Prior art keywords
- array element
- signal
- array
- sup
- vector
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Radio Transmission System (AREA)
Abstract
Description
技术领域
本发明涉及阵列信号处理中的波达方向估计方法,具体涉及一种基于空域自适应滤波的窄带信号的波达方向估计方法。
技术背景
信号的波达方向估计是阵列信号处理中的重要研究内容,广泛应用于雷达、声纳、通信等***中。目前,针对大信噪比情况下的窄带信号的各种波达方向(Direction Of Arrival,DOA)估计方法中,多重信号分类法(Multiple SignalsClass,MUSIC)具有超分辨特性,受到广泛关注。但由于MUSIC算法在实现过程中,需要估计信源个数,同时也需要对接收信号协方差矩阵进行特征分解,一旦接收阵元的数量较大,对信号进行处理的计算量大大增加,这些都不利于在FPGA或者DSP器件中实现。由于MUSIC算法的这些缺点,导致该方法工程实现困难。
发明内容
本发明专利要解决的技术问题是针对窄带信号设计新的波达方向估计方法,以克服使用经典MUSIC算法中需要已知信源个数,计算过程复杂,计算量大的不足。
解决上述问题的技术方案是一种基于空域自适应滤波的波达方向估计方法,其实现步骤为:
第一步,由N+1个阵元构成的均匀线阵中,阵元间距l为半个波长,第一个 阵元为参考阵元,其余N个阵元为辅助阵元,对阵元接收的窄带信号进行K次采样,采样频率为fs,则参考阵元接收信号为d=[d(1)...d(K)],辅助阵元接收信号为 ,......, 其中N和K都为自然数,K>2N;第二步,对所有接收数据矢量进行正交变换,得到相应的解析信号 ......,
第三步,对解析信号采用基于功率倒置的空域自适应滤波算法,进行M次迭代计算求权值,在权值w(k)收敛后,最后一次迭代计算出最优权值w(M),该权矢量为N行列矢量,其中,迭代次数满足k<M<K;
第四步,根据得到的最优权值矢量w(M),把所有N+1个阵元视为一个新的阵列,构造一个新的权值矢量为 及方向矢量 v(θ)为入射信号进入N个辅助阵元的方向矢量;
上述一种基于空域自适应滤波的波达方向估计方法,其第三步采用的基于功率倒置的空域自适应滤波算法,具体的实现步骤为:
第一步,取初值w(0)=0,k=1,收敛判决门限ε,迭代步长μ;
第三步,迭代计算权值
第四步,计算|e(k)|2,进行判断,若|e(k)|2≥ε,则k=k+1,转第二步;若|e(k)|2<ε,权值w(k)收敛,此时取M=k,停止计算。
本发明的有益效果是:新方法实现了窄带信号的DOA估计,不仅无需进行 信源数估计和矩阵特征分解,计算量大大减小,而且该方法在保持高分辨率情况下,各个计算步骤都适合FPGA或者DSP设计,便于工程实现。
附图说明
图1均匀线阵下阵列接收信号的示意图
图2本发明方法的流程示意图
图3空域自适应滤波结构框图
图4本发明方法的窄带信号DOA估计空间谱图
具体实施方法
本发明方法是一种基于空域自适应滤波的窄带信号的波达方向估计方法,该方法采用一种基于功率倒置的自适应空域滤波算法来构造阵列权值矢量及方向矢量,然后定义空间谱函数,从而实现信号DOA估计。本发明无需进行信源数估计和矩阵特征分解,减少了运算量,便于工程实现。下面对本发明方法进行详细地描述。
本发明方法针对远场不相干的窄带信号,图1为均匀线阵下阵列接收信号的示意图。图中,阵列由N+1个相同的全向阵元构成,阵元间距l为半个波长。空间具有J个入射信号,这些信号DOA分别为θj,j=1,2,...,J,每个DOA位于区间J<N+1。在***阵列为N+1阵元构成的均匀线阵中,设置第一个阵元为参考阵元,位于坐标原点,其余N个为辅助阵元,则第j个信号到达该点时为sj(t)。阵列接收信号可以表示为y(t)=[d(t),x1(t),...,xN(t)]T,符号[]T表示对矢量和矩阵求转置,其中i0(t)为参考阵元接收的噪声信号,第n个 辅助阵元接收信号为τnj为第n个辅助阵元接收到的第j个信号相对于原点信号的延时,in(t)为第n个辅助阵元接收的噪声信号,n=1,2,...,N。
图2为本发明方法的流程示意图,具体方法分为五个步骤:
第一步,对阵元接收的窄带信号进行K次采样,采样频率为fs,则参考阵元接收信号为d=[(d(1)...d(K)],辅助阵元接收信号为 ,....., 其中N和K都为自然数,K>2N;
第三步,对解析信号采用基于功率倒置的空域自适应滤波算法,进行M次迭代计算求权值,在权值w(k)收敛后,最后一次迭代计算出最优权值w(M),该权矢量为Ⅳ行列矢量,其中,迭代次数满足k<M<K;
第四步,根据得到的最优权值矢量w(M),把所有N+1个阵元视为一个新的阵列,构造一个新的权值矢量为 及方向矢量 v(θ)为入射信号进入N个辅助阵元的方向矢量;
上述一种基于空域自适应滤波的波达方向估计方法,其第三步采用的基于功率倒置的空域自适应滤波算法,具体的实现步骤为:
第一步,取初值w(o)=0,k=1,收敛判决门限ε,迭代步长μ;
第三步,迭代计算权值
第四步,计算|e(e)|2,进行判断,若|e(k)|2≥ε,则k=k+1,转第二步;若|e(k)|2<ε,权值w(k)收敛,此时取M=k,停止计算。
图4是按照本发明方法对非相干窄带信号进行DOA估计仿真,得到的空间谱图,其中阵列个数N+1=8,l为半个波长,空间信号J=2,DOA分别为和采样频率fs=10MHz,快拍数K=2000,信噪比分别为10dB和20dB。图中,两个谱峰对应的角度等于信号入射的角度,说明本发明方法的正确性。
本发明利用功率倒置自适应空间滤波产生权矢量,计算空间谱,从而估计信号DOA。相对于传统MUSIC算法,保持高分辨率情况下,运算量大大减小,满足***实时性要求,实现更为简单。
Claims (2)
1.一种基于空域自适应滤波的波达方向估计方法,其实现步骤为:
第一步,由N+1个阵元构成的均匀线阵中,阵元间距l为半个波长,第一个阵元为参考阵元,其余N个阵元为辅助阵元,对阵元接收的窄带信号进行K次采样,采样频率为fs,则参考阵元接收信号为d=[d(1) ... d(K)],辅助阵元接收信号为 ......, 其中N和K都为自然数,K>2N;
第二步,对所有接收数据矢量进行正交变换,得到相应的解析信号 ......,
第三步,对解析信号采用基于功率倒置的空域自适应滤波算法,进行M次迭代计算求权值,第k次迭代权值为w(k),收敛后的最后一次迭代计算出最优权值w(M),该权矢量为N行列矢量,其中,迭代次数满足k<M<K,k=1,2,…,M-1;
第四步,根据得到的最优权值矢量w(M),把所有N+1个阵元视为一个新的阵列,构造一个新的权值矢量为 及方向矢量 v(θ)为入射信号进入N个辅助阵元的方向矢量;
第五步,定义空间谱在可见区间内,对空间谱进行谱峰搜索,谱峰位置即为估计出的信号波达方向,谱峰个数即为信号个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210404468.8A CN102866383B (zh) | 2012-10-22 | 2012-10-22 | 一种基于空域自适应滤波的波达方向估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210404468.8A CN102866383B (zh) | 2012-10-22 | 2012-10-22 | 一种基于空域自适应滤波的波达方向估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102866383A CN102866383A (zh) | 2013-01-09 |
CN102866383B true CN102866383B (zh) | 2014-01-22 |
Family
ID=47445366
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210404468.8A Expired - Fee Related CN102866383B (zh) | 2012-10-22 | 2012-10-22 | 一种基于空域自适应滤波的波达方向估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102866383B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104023397B (zh) * | 2014-06-20 | 2017-06-13 | 深圳大学 | 分布式网络中基于gossip算法的多目标DOA估计***及估计方法 |
CN104777449B (zh) * | 2014-10-31 | 2017-06-06 | 哈尔滨工业大学深圳研究生院 | 基于贝叶斯信息准则的信源数估计方法 |
CN104569920B (zh) * | 2014-11-14 | 2017-11-28 | 哈尔滨工业大学深圳研究生院 | 无需信源数信息的联合估计到达角‑频率的方法 |
CN104502885A (zh) * | 2015-01-04 | 2015-04-08 | 哈尔滨工程大学 | 基于变换矩阵的特征值差分信源数估计方法 |
CN105262530B (zh) * | 2015-09-21 | 2019-03-19 | 梁海浪 | 一种来波方向快速估计的方法 |
CN108802669B (zh) * | 2018-07-13 | 2020-08-25 | 中国人民解放军陆军工程大学 | 二维波达方向估计方法、二维波达方向估计装置及终端 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2713761B2 (ja) * | 1989-04-20 | 1998-02-16 | 沖電気工業株式会社 | 時間空間信号のスペクトル推定方法 |
CN100380135C (zh) * | 2005-07-26 | 2008-04-09 | 武汉大学 | 基于天线子阵的高频地波雷达抗干扰方法 |
CN101887129B (zh) * | 2010-04-23 | 2013-11-06 | 西安航空电子科技有限公司 | 北斗卫星用户设备接收抗干扰方法 |
CN101908918B (zh) * | 2010-07-26 | 2015-04-01 | 重庆大学 | 无线通信接收机中的一种波束合成方法 |
CN102621539B (zh) * | 2012-04-24 | 2013-12-11 | 重庆大学 | 单次搜索宽带信号波达方向估计方法 |
-
2012
- 2012-10-22 CN CN201210404468.8A patent/CN102866383B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102866383A (zh) | 2013-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102866383B (zh) | 一种基于空域自适应滤波的波达方向估计方法 | |
Liu et al. | Sparsity-inducing direction finding for narrowband and wideband signals based on array covariance vectors | |
CN103091661B (zh) | 基于迭代谱重构的宽带信号波达方向估计方法 | |
CN103353596B (zh) | 基于压缩感知的波束空间域米波雷达测高方法 | |
CN103344940B (zh) | 低复杂度的doa估计方法及*** | |
CN103901395B (zh) | 一种冲击噪声环境下相干信号波达方向动态跟踪方法 | |
CN103605108B (zh) | 声矢量阵高精度远程方位估计方法 | |
CN105403856A (zh) | 基于嵌套式最小冗余阵列的波达方向估计方法 | |
CN104698433A (zh) | 基于单快拍数据的相干信号doa估计方法 | |
CN104977558A (zh) | 一种基于贝叶斯压缩感知的分布源中心波达方向估计方法 | |
CN105093185A (zh) | 基于稀疏表示的单基地多输入多输出雷达目标波达方向估计方法 | |
CN103744061A (zh) | 基于迭代最小二乘方法的mimo雷达doa估计方法 | |
CN103713276B (zh) | 基于最小互熵谱分析的波达方向估计方法 | |
CN105259550A (zh) | 基于压缩感知的多输入多输出雷达二维角度估计方法 | |
CN103364772A (zh) | 基于实数域广义多重信号分类算法的目标低仰角估计方法 | |
CN106019214A (zh) | 宽带相干信号源doa估计方法 | |
CN105403874A (zh) | 非均匀阵列欠定波达方向估计方法 | |
CN104678372A (zh) | 正交频分复用雷达超分辨距离与角度值联合估计方法 | |
CN104865556A (zh) | 基于实域加权最小化l1范数方法的MIMO雷达***DOA估计方法 | |
CN103616661A (zh) | 一种稳健的远场窄带信号源个数估计方法 | |
CN101483280B (zh) | 一种稳健波束合成器权值求解方法 | |
Wu et al. | Subspace-based method for joint range and DOA estimation of multiple near-field sources | |
CN102621539B (zh) | 单次搜索宽带信号波达方向估计方法 | |
CN104330766A (zh) | 一种稳健的波达方向估计方法 | |
CN106483193B (zh) | 一种基于高阶累计量的波达快速估计方法 |
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: 20140122 Termination date: 20141022 |
|
EXPY | Termination of patent right or utility model |