专利名称：基于龙伯格观测器与子空间更新的波达方向跟踪方法及装置的制作方法
技术领域：
本发明属于无线电波信号分析处理技术领域，涉及一种波达方向的自适应估计与跟踪。具体涉及一种基于龙伯格观测器与子空间更新的波达方向跟踪方法及装置，用于估计和跟踪自适应均勻线性阵列接受到的远场窄带无线电波的到达方向，该远场窄带无线电波包括相干和非相干电波。
背景技术：
阵列信号处理将多个天线阵子设置在空间的不同位置组成天线阵列，并利用这一天线阵列对空间无线电波进行接受和处理，提取阵列所接收的无线电波信号(从信号处理的角度，无线电波成为信号)及其特征信息。阵列信号处理具有灵活的波束控制、高的信号增益、很强的抗干扰能力和高的空间超分辨能力等特点，受到人们的高度关注。阵列信号处理中一个典型的问题是基于自适应阵列天线估计噪声环境中入射信号的到达方向，由于其计算简单和高分辨能力，基于子空间的方向估计方法成为公认的用来估计到达方向的方法得到很多的关注，此方法中，利用了信号子空间和噪声子空间之间的正交性。基于子空间的方法中多重信号分类(MUSIC)技术是一种典型的方向估计方法。在用于估计非相关信号的到达方向的基于子空间的方法中，首先根据噪声接收阵列数据来估计阵列协方差矩阵，然后通过该阵列协方差矩阵的特征值分解(EVD)来获得信号子空间和噪声子空间。此后，采用该信号子空间和噪声子空间之间的正交性，来估计入射信号的到达方向。另一方面，对于相关信号(包括完全相关信号，即相干信号)的到达方向的估计，为了抑制入射信号的相关性，将具有沿直线以相同的间隔设置在不同的空间位置的M个天线阵子的天线(下文中，将这种天线称为均勻线性阵列(ULA))分为多个重叠的子阵列，然后对这些子阵列的协方差矩阵求平均，以将空间平均的协方差矩阵的信号子空间的维数恢复为入射信号的数量。由此，可以采用信号子空间和噪声子空间之间的正交关系， 按照与用于估计非相关信号的到达方向的基于子空间的方法相同的方式，来估计相干信号的到达方向。在许多实际应用中，需要自适应估计确定方向信号或者即使跟踪随时间变化的到达方向。但是，由于信号子空间或者噪声子空间计算中批处理特征分解(例如特征值分解或者奇异值分解)较大的计算负荷和时间耗费，大多数的基于子空间的方法很难在实时操作中实施。因此，为了克服传统的需要特征分解的子空间方法的计算缺陷，子空间更新或者跟踪算法引起了极大的关注，并且在很多文献中出现了大量的更新算法。除了引进各种有效技术和一些替换矩阵分解法，子空间估计也可以被表达成有限制或无限制的优化问题。 同时，人们也提出许多自适应算法来解决子空间更新问题，例如高斯-牛顿迭代法，共轭梯度法，递推最小二乘法以及它们的改进算法。虽然这些子空间更新算法已经被成功应用在特殊的非稳定状态信号环境的到达方向估计中，但是大部分方法并不是适合各种反射导致的多径传播中的相干信号。在现有技术中提出的无需特征分解的基于子空间的方法(SUMWE)的到达方向估计方法，其中利用阵列数据计算得一线性算子组成零空间的基，该方法不需要特征分解，并且计算效率高(例如，参见J. Xin and A. Sano, "Computationally efficient subspace-basedmethod for direction-of-arrival estimation without eigendecomposition,，，IEEE Trans. Signal Process. , vol. 52, no. 4, pp. 876-893, Apr. 2004.)。然而该方法没有考虑在线到达方向估计以及对随时间而变化的到达方向进行跟踪问题。
发明内容
本发明的目的在于克服上述现有技术的缺点，提供一种基于龙伯格观测器与子空间更新的波达方向跟踪方法，该方法通过运用RLS算法来对噪声子空间进行实时更新，进而使用龙伯格观测器来实现对无线电波的方向估计与跟踪。本发明的目的是通过以下技术方案来解决的这种基于龙伯格观测器与子空间更新的波达方向跟踪方法，包括以下步骤1)初始化龙伯格观测器的增益以及可检测的最大来波数；2)根据瞬时相关性矩阵估计出来波数量；3)根据来波数量和瞬时相关性矩阵估计出初始波达方向；4)用龙伯格观测器进行波达方向的一步预测；5)更新瞬时噪声子空间；6)运用瞬时的测量值和龙伯格观测器来估计该瞬时的波达方向。以上如果以M表示天线阵子的数量，则通过第M天线阵子所接收的数据与第一、第
二.....第(M-I)天线阵子中的每一个所接收的数据之间的相关性，来获得所述瞬时相关
矩阵，并且在时间和空间非相关白噪声环境中对所述无线电波的波数和方向进行自适应估
计和跟踪；则通过第1天线阵子所接收的数据与第二、第三.....第M天线阵子中的每一个
所接收的数据之间的相关性，来获得所述瞬时相关矩阵，并且在时间和空间非相关白噪声环境中对所述无线电波的波数和方向进行自适应估计和跟踪。进一步，以上根据所述瞬时相关性来生成1至4个所述瞬时相关矩阵，并使用该1 至4个所述瞬时相关性矩阵来估计所述无线电波的来波个数；并使用该1至4个所述瞬时相关性矩阵来估计所述瞬时噪声子空间。以上步骤3)中，根据所述瞬时噪声子空间估计无线电波的波达方向。以上步骤4)中，使用递归最小二乘来实时更新加权矩阵。以上步骤6)中，根据所述的波达方向的一步预测值和所述的瞬时噪声子空间来估计当前来波方向的测量值。以上所述方法不仅对时间和空间非相关白噪声环境中的相干信号进行自适应估计和跟踪；而且对时间和空间非相关白噪声环境中的部分相干信号进行自适应估计和跟踪；另外所述方法对时间和空间非相关白噪声环境中的非相干信号进行自适应估计和跟
S示ο本发明还提出一种实现上述方法的方向估计与跟踪装置，其包括
瞬时相关性计算模块用于计算一天线阵子所接收的数据与另一天线阵子所接收的数据之间的瞬时相关性；瞬时相关性矩阵计算模块用于根据所述瞬时相关性来计算瞬时相关性矩阵；噪声子空间估计模块用于利用所述瞬时相关性矩阵，通过线性运算来估计噪声子空间；到达方向估计模块用于根据所述噪声子空间来计算所述到达方向的估计值。本发明具备以下有益效果本发明针对基于均勻线阵(ULA，uniform 1 inear array)的多个相干窄带信号的波达方向(DOA，directions-of-arrival)问题，提出新的计算简单而有效的基于子空间的自适应波达方向估计方法(AMEND，adaptive method for estimating directions-of-arrival)。其中，利用改进的基于QR分解的MENSE方法确定入射信号的个数，而且提出一个新的递归最小二乘(RLS，recursive least-squares)算法来实现零空间更新，同时采用动态模型和Luenberger状态观测器来解决方向估计值的自动关联。理论分析和数值仿真结果表明该方法具有良好的在线估计和实时跟踪效果。
图1表示根据本发明实施例的方向估计与跟踪算法的概要流程图；图2表示方向更新间隔与采样间隔之间的关系；图3表示线性阵列天线和发送源的结构；图4表示非相关白噪声环境下的阵列协方差矩阵；图5表示均勻线性阵列中的子阵列；图6表示使用龙伯格观测器的方向估计与跟踪示意图；图7表示方向估计与跟踪算法的方框图。
具体实施例方式本发明的该种基于龙伯格观测器与子空间更新的波达方向跟踪方法，具体包括以下步骤1)初始化龙伯格观测器的增益以及可检测的最大来波数；2)根据瞬时相关性矩阵估计出来波数量；如果以M表示天线阵子的数量，则通过第M天线阵子所接收的数据与第一、第
二.....第(M-I)天线阵子中的每一个所接收的数据之间的相关性，来获得所述瞬时相关
矩阵，并且在时间和空间非相关白噪声环境中对所述无线电波的波数和方向进行自适应估
计和跟踪；通过第1天线阵子所接收的数据与第二、第三.....第M天线阵子中的每一个所
接收的数据之间的相关性，来获得所述瞬时相关矩阵，并且在时间和空间非相关白噪声环境中对所述无线电波的波数和方向进行自适应估计和跟踪。并且根据所述瞬时相关性来生成1至4个所述瞬时相关矩阵，使用该1至4个所述瞬时相关性矩阵来估计所述无线电波的来波个数；并使用该1至4个所述瞬时相关性矩阵来估计所述瞬时噪声子空间。3)根据来波数量和瞬时相关性矩阵估计出初始波达方向；具体为根据所述瞬时噪声子空间估计无线电波的波达方向。
4)用龙伯格观测器进行波达方向的一步预测；在龙伯格观测器中，使用递归最小二乘来实时更新加权矩阵。5)更新瞬时噪声子空间；6)运用瞬时的测量值和龙伯格观测器来估计该瞬时的波达方向。具体为根据波达方向的一步预测值和所述的瞬时噪声子空间来估计当前来波方向的测量值。本发明所述的方法不仅能够对时间和空间非相关白噪声环境中的相干信号进行自适应估计和跟踪，而且能对时间和空间非相关白噪声环境中的部分相干信号进行自适应估计和跟踪，还可以对时间和空间非相关白噪声环境中的非相干信号进行自适应估计和跟
S示ο为了实现本发明的上述方法，本发明提出一种方向估计与跟踪装置，该装置用于利用阵列天线来估计和跟踪无线电波的到达方向，该阵列天线具有沿直线以相同的间隔设置在不同空间位置上的多个天线阵子，该装置包括瞬时相关性计算模块、瞬时相关性矩阵计算模块、噪声子空间估计模块和到达方向估计模块。其中瞬时相关性计算模块用于计算一天线阵子所接收的数据与另一天线阵子所接收的数据之间的瞬时相关性；瞬时相关性矩阵计算模块用于根据所述瞬时相关性来计算瞬时相关性矩阵；噪声子空间估计模块用于利用所述瞬时相关性矩阵，通过线性运算来估计噪声子空间；到达方向估计模块用于根据所述噪声子空间来计算所述到达方向的估计值。下面将参照附图来详细说明本发明的实施例。下面概要描述根据本发明实施例的无线电波到达方向估计和跟踪方法。根据本发明实施例的无线电波到达方向估计和跟踪方法尤其适用下述的情况在线到达方向估计以及对随时间而变化的到达方向进行跟踪。当前的描述假定构成阵列天线的M个天线阵子接收ρ个到达的无线电波(其中M > 2p)。本发明适用于前向子阵列、后向子阵列，以及前向和后向子阵列。现假定ρ个窄带信号Isi (t)}沿角度{ θ i(t)}入射到M个阵元组成的均勻线性阵列(ULA)上。各个阵元所接收到的信号可以表示为公式(1)。y(k) = [y1(k),y2(k),...,yM(k)]T = A(6 (η)) s (k)+w (k)A 全(η)), α(θ2 (η)),…,α{θρ ( ))]αψχη)) = ReM 明("《，…，…《。(财-耿物)^s(k) = [S1 (k), S2(k), . . . , sp(k)]T…(1)w (k) = [W1 (k), w2(k), ... ,wM(k)]Tω0 = Inf0,糊⑷)kd/c sin O1 (k)其中k = nNs+l，nNs+2，· · ·，(n+l)Ns。η = 0,1,2,…。f。、c 和 d 分别表示载波频率和载波的传播速度，以及阵元间隔(半波长)。(0工表示转置(transposition)，而 B(O^n))和A分别与阵列响应向量和矩阵相对应。Wi(k)表示具有零均值和方差Q2的空间和时间上的不相关的复白噪声。现假定无线电波的到达角度θ i (t)相对于采样频率1/TS来说是缓慢变化的。即在te (nT, (η+1)Τ]时间间隔内θ Jt) θ , (ηΤ) 0 和殘⑷分别表示Mn)变化的速度和加速度。在不存在过程噪声和量测噪声的情况下使用确定状态模型来进行状态估计
其中状态向量X 全[没 Ql(Ji)殘⑷]“。图4示出了非相关白噪声环境下的阵列协方差矩阵。考虑信号到达方向对于时间而言是恒定的情况，即ejn) = 0iO在这种情况下， 阵列协方差矩阵表示为公式(3).R ^ E{y(k)yH (k)} = ARsAh + σ Μ■■■ (3)其中Ε{ · }和(·)H分别表示期望运算和复共轭转置，Rs = E{s(k)sH(k)}是入射信号的协方差矩阵，而Lm是MXM单位矩阵。此外如果将所接收到的数据yi(k)之间的相关性Am定义为^ =E{yi(k)ym(k)}，则关系^ =《(丫成立，其中( 广表示复共轭。公
式(3)中的阵列协方差矩阵R可以表示为公式G)。
图5示出了均勻线性阵列中的子阵列。将整个阵列分为L个叠加的前向和后向子阵列，每个子阵列都包含ρ个阵元。则第1个前向和后向子阵接收到的信号为yfl(k) ^ [yrikly^kl-^-ml 和^ “#)，>^⑷，"·，Λ—ω⑷广， 则可以得到k时刻的四个瞬时相关性Hankel矩阵(具体参见J. Xin and A. Sano,"Computationally efficient subspace-based method for direction-of-arrivalestimation without eigendecomposition,” IEEE Trans. Signal Process.，vol. 52，no. 4，pp. 876-893，Apr. 2004.)如公式(5)禾口(6)，Of(Zc) = Yf{k)yM{k\ = ff{k)y;{k)…(5)Φ, (k) = Yb (k)yi (k\ Φ, (k) = Yb (k)yM (k)…(6)其中((k)^ [^1 (k\yf2(k\···,y^(k)f，Ff (k) ^ [yf2(k\yf,(k\(k)f，
Yb(k) = [yM(k),yb2(kl···,ybL^⑷]r，ζ⑷=U2⑷,yb3⑷,…，ybL(k)f，并且L = M-P+i。把公式 (5)禾Π (6)中的每个(Μ-ρ) XpHankel矩阵分割为一个ρΧρ前向子矩阵和一个(M_2p) Xp
后向子矩阵如公式(8)， 「 ⑷
权利要求
1.一种基于龙伯格观测器与子空间更新的波达方向跟踪方法，其特征在于，包括以下步骤1)初始化龙伯格观测器的增益以及可检测的最大来波数；2)根据瞬时相关性矩阵估计出来波数量；3)根据来波数量和瞬时相关性矩阵估计出初始波达方向；4)用龙伯格观测器进行波达方向的一步预测；5)更新瞬时噪声子空间；6)运用瞬时的测量值和龙伯格观测器来估计该瞬时的波达方向。
2.根据权利要求1所述的基于龙伯格观测器与子空间更新的波达方向跟踪方法，其特征在于，如果以M表示天线阵子的数量，则通过第M天线阵子所接收的数据与第一、第二.....第(M-I)天线阵子中的每一个所接收的数据之间的相关性，来获得所述瞬时相关矩阵，并且在时间和空间非相关白噪声环境中对所述无线电波的波数和方向进行自适应估计和跟踪；则通过第1天线阵子所接收的数据与第二、第三.....第M天线阵子中的每一个所接收的数据之间的相关性，来获得所述瞬时相关矩阵，并且在时间和空间非相关白噪声环境中对所述无线电波的波数和方向进行自适应估计和跟踪。
3.根据权利要求1所述的基于龙伯格观测器与子空间更新的波达方向跟踪方法，其特征在于，根据所述瞬时相关性来生成1至4个所述瞬时相关矩阵，并使用该1至4个所述瞬时相关性矩阵来估计所述无线电波的来波个数；并使用该1至4个所述瞬时相关性矩阵来估计所述瞬时噪声子空间。
4.根据权利要求1所述的基于龙伯格观测器与子空间更新的波达方向跟踪方法，其特征在于，步骤3)中，根据所述瞬时噪声子空间估计无线电波的波达方向。
5.根据权利要求1所述的基于龙伯格观测器与子空间更新的波达方向跟踪方法，其特征在于，步骤4)中，使用递归最小二乘来实时更新加权矩阵。
6.根据权利要求1所述的基于龙伯格观测器与子空间更新的波达方向跟踪方法，其特征在于，步骤6)中，根据所述的波达方向的一步预测值和所述的瞬时噪声子空间来估计当前来波方向的测量值。
7.根据权利要求1所述的基于龙伯格观测器与子空间更新的波达方向跟踪方法，其特征在于，所述方法对时间和空间非相关白噪声环境中的相干信号进行自适应估计和跟踪。
8.根据权利要求1所述的基于龙伯格观测器与子空间更新的波达方向跟踪方法，其特征在于，所述方法对时间和空间非相关白噪声环境中的部分相干信号进行自适应估计和跟S示ο
9.根据权利要求1所述的基于龙伯格观测器与子空间更新的波达方向跟踪方法，其特征在于，所述方法对时间和空间非相关白噪声环境中的非相干信号进行自适应估计和跟S示ο
10.一种实现权利要求1所述方法的方向估计与跟踪装置，其特征在于，包括 瞬时相关性计算模块用于计算一天线阵子所接收的数据与另一天线阵子所接收的数据之间的瞬时相关性；瞬时相关性矩阵计算模块用于根据所述瞬时相关性来计算瞬时相关性矩阵； 噪声子空间估计模块用于利用所述瞬时相关性矩阵，通过线性运算来估计噪声子空间；到达方向估计模块用于根据所述噪声子空间来计算所述到达方向的估计值。
全文摘要
本发明公开了一种基于龙伯格观测器与子空间更新的波达方向跟踪方法及装置，本发明针对基于均匀线阵(ULA，uniform linear array)的多个相干窄带信号的波达方向(DOA，directions-of-arrival)问题，提出一种新的计算简单而有效的基于子空间的自适应波达方向估计方法(AMEND，adaptive method for estimating directions-of-arrival)。其中，利用改进的基于QR分解的MENSE方法确定入射信号的个数，而且提出一个新的递归最小二乘(RLS，recursive least-squares)算法来实现零空间更新，同时采用动态模型和Luenberger状态观测器来解决方向估计值的自动关联。理论分析和数值仿真结果表明该方法具有良好的在线估计和实时跟踪效果。
文档编号G01S3/12GK102175990SQ20111002996
公开日2011年9月7日 申请日期2011年1月27日 优先权日2011年1月27日
发明者曹祥, 辛景民, 郑南宁 申请人:西安交通大学