无线信道参数估计算法
对通信系统分析之前,首先需要建立一个能够完整反映该系统的传输信道的模型,其模型中的参数包括接收端接收信号的时延,发射端和接收端的联合方向分布等,获取这些参数的过程就是无线信道参数估计的过程。
方法可以大致分为三类:
- 谱估计算法:多重信号分类法
- 参数子空间估计算法:旋转不变技术估计信号参数算法
- 确定性参数估计算法:期望最大算法以及演进的空间交替广义期望最大算法
0. 前言
阵列信号处理基础本质上属于参数估计问题,和信道估计知识基本上别无二致。末学在这里整理了阵列信号处理的基础知识。
1. 多重信号分类法
MUSIC 算法原理
MUSIC 算法,叫做多信号分类算法 (Multiple Signal classification)。MUSIC 算法的基本思想是将任意阵列输出数据的协方差矩阵进行特征分解,从而得到与信号分量相对应的信号子空间和信号分量相正交的噪声子空间,然后利用这两个子空间的正交性来估计信号的参数 (入射方向、极化信息和信号强度)。MUSIC 算法是一种高分辨率的、高精度的无线信道参数估计算法,而且该算法对于天线阵列的形状没有特殊要求,具有普遍的适用性,只要已知天线阵的布阵形式,无论是直线阵列还是圆阵列,不管阵元是否是等间距分布,都可以得到高分辨率的结果。但是该算法对于入射信号的要求非常高,应用该算法的前提是入射信号必须是互不相干的。
考虑由 N N N 个阵元组成的均匀线性阵列 (Uniform Linear Array),阵列之间的间距为半波长 d = λ / 2 d=\lambda/2 d=λ/2,假设有 M M M 个远场信号源,它们的功率相同且互不相干,分别从不同的方向入射到接收阵列,其中入射信号和噪声互不相关。以最左侧的阵元为参考点,各个阵元的位置相比前一个增加距离 d d d。信号入射方向角为 θ \theta θ,表示入射信号与线阵法线的夹角。一般假设信号源数 M M M 小于阵元数目 N N N。
在第 k k k 次快拍,得到的接收数据向量为:
X ( k ) = A ( θ ) S ( k ) + N ( k ) , k = 1 , 2 , ⋯ , K \boxed{ X(k)=A(\theta)S(k)+N(k),k=1,2,\cdots,K } X(k)=A(θ)S(k)+N(k),k=1,2,⋯,K
其中, A ( θ ) = [ a ( θ 1 ) , a ( θ 2 ) , ⋯ , a ( θ M ) ] A(\theta)=\left[ a(\theta_1),a(\theta_2),\cdots,a(\theta_M)\right] A(θ)=[a(θ1),a(θ2),⋯,a(θM)] 是 N × M N\times M N×M 的阵列矩阵,又称为阵列流形向量或者方向响应向量。 S ( k ) = [ s 1 ( k ) , s 2 ( k ) , ⋯ , s M ( k ) ] T S(k)=\left[ s_1(k),s_2(k),\cdots,s_M(k)\right]^T S(k)=[s1(k),s2(k),⋯,sM(k)]T 是 M M M 个入射信号矢量。
X ( k ) = [ x 1 ( k ) , x 2 ( k ) , ⋯ , x N ( k ) ] T X(k)=\left[ x_1(k),x_2(k),\cdots,x_N(k)\right]^T X(k)=[x1(k),x2(k),⋯,xN(k)]T 是 N N N 个输出信号矢量。
N ( k ) = [ n 1 ( k ) , n 2 ( k ) , ⋯ , n N ( k ) ] T N(k)=\left[ n_1(k),n_2(k),\cdots,n_N(k)\right]^T N(k)=[n1(k),n2(k),⋯,nN(k)]T 是阵列噪声矢量。
一般假设噪声为加性高斯即: E ( n i ( k ) ) = 0 , E ( n i ( k ) n i ( k ) H ) = σ 2 I \mathbb{E}(n_i(k))=0,\quad \mathbb{E}(n_i(k)n_i(k)^H)=\sigma^2I E(ni(k))=0,E(ni(k)ni(k)H)=σ2I。
阵列矩阵的元素满足:
a ( θ i ) = [ 1 , e − j 2 π d sin θ i / λ , ⋯ , e − j ( N − 1 ) 2 π d sin θ i / λ ] a(\theta_i)=\left[1,e^{-j2\pi d \sin{\theta_i}/ \lambda},\cdots,e^{-j(N-1)2\pi d \sin{\theta_i}/ \lambda} \right] a(θi)=[1,e−j2πdsinθi/λ,⋯,e−j(N−1)2πdsinθi/λ]
其中, λ \lambda λ 是中心波长, θ i \theta_i θi 为信号源的入射角。因此,问题转化为对 X X X 进行多点采样,通过对采样的输出信号 X ( k ) X(k) X(k) 估计出信号源的波达方向角 θ i \theta_i θi。
由上图中两束方向确定的信号发射到两个不同阵元,在同一时刻采样的信号值为随机值,待估的方向一定。于是,输出信号的协方差矩阵可以表示为:
R X X = E [ X ( k ) X ( k ) H ] = E [ A ( θ ) S ( k ) S ( k ) H A ( θ ) H ] + E [ N ( k ) N ( k ) H ] = A ( θ ) E [ S ( k ) S ( k ) H ] A ( θ ) H + R N = A R S A H + σ 2 I \begin{aligned} R_{XX}&=\mathbb{E}\left[X(k)X(k)^H \right] \\ &=\mathbb{E}\left[A(\theta)S(k)S(k)^HA(\theta)^H \right]+\mathbb{E}\left[N(k)N(k)^H \right]\\ &=A(\theta)\color{#F00}{\mathbb{E}\left[S(k)S(k)^H \right]}A(\theta)^H+R_N\\ &=A\color{#F00}{R_S}A^H+\sigma^2I \end{aligned} RXX=E[X(k)X(k)H]=E[A(θ)S(k)S(k)HA(θ)H]+E[N(k)N(k)H]=A(θ)E[S(k)S(k)H]A(θ)H+RN=ARSAH+σ2I
其中 R S R_S RS 是输入信号的协方差矩阵。由于信号和噪声相互独立,数据协方差可分解为信号子空间和噪声子空间。由于 R X X R_{XX} RXX 是对称矩阵,然后我们对其进行特征值分解,得到
R X X = U Σ U H = U S Σ S U S H + U N Σ N U N H R_{XX}=U \Sigma U^H = U_{S} \Sigma_{S}\ U_{S}^{H}+U_{N} \Sigma_{N} U_{N}^{H} RXX=UΣUH=USΣS USH+UNΣNUNH
其中 U = [ U S , U N ] U=\left[ U_S,U_N\right] U=[US,UN], Σ = diag ( σ 1 2 , ⋯ , σ N 2 ) \Sigma = \text{diag} \left( \sigma_{1}^2,\cdots,\sigma_{N}^2 \right) Σ=diag(σ12,⋯,σN2), U S U_S