相干信号空间谱估计测向Matlab仿真研究

相干信号空间谱估计测向Matlab 仿真研究1. 引言由于多径传播、电磁干扰等因素的影响，相干信源存在的电磁环境是经常碰到的。

当空间存在相干源时，经典的超分辨DOA 估计方法：MUSIC 算法和ESPRIT 算法，已经失去了其高分辨性能优势，有时甚至不能正确地估计出信源的真实方位。

新MUSIC 算法在空间不存在相干源时，其估计性能基本上是和MUSIC 算法是接近的，但若有相干源存在时，其估计性能也是大大降低。

因此，若将其用于相干源，必须和经典的MUSIC 算法一样，首先对阵列输出的协方差矩阵进行各种去相干处理，然后再采用新MUSIC 算法实现对相干信源的DOA 估计。

基于加权空间平滑的MUSIC 算法，该算法充分利用了子阵输出的自相关信息和互相关信息，将阵列协方差矩阵的所有子阵阵元数阶子矩阵进行加权平均，而权矩阵的选取以平滑后等价的信源协方差矩阵与对角阵的逼进为约束条件，以期对相干信源最大限度地去相干，改进常规空间平滑算法对相干源的分辨力。

基于此本文提出的基于加权空间平滑的新MUSIC 算法，以实现对相干源最大限度的去相干，实现相干源的高分辨DOA 估计。

2. 窄带阵列相干源的数学模型和空间平滑算法 2.1 窄带阵列相干源的数学模型对于M 元均匀线阵，阵源间距为d ，且假设均为各向同性阵元。

阵列远场中在以线阵轴线法线为参考的(1,2,...)k P =k θ处有 P 个窄带点源以平面波入射，以阵列第一阵元为参考点，某一特定信号到达线阵时，各阵元接收信号间仅仅存在因波程差引起的相位差。

阵列接收的快拍数据可由下式表示为:t t t θX()=A()S()+N() (2-1)(2-1)式中()t X 为M 1⨯快拍数据矢量，是阵元输出信号，1()[(),......()]T M t x t x t =X 。

()t N 为M 1⨯阵列阵元噪声矢量，且是与信号源不相关的白噪声，均值为0，方差为2n σ，1()[(),......()]T M t n t n t =N 。

()t S 为输入信号矢量1()[(),......()]T M t s t s t =S 。

()θA 为阵列的流形矩阵,如式(2-2)所示。

1()[(),......()]P θθθ=A a a =11()(1)()1...j j M e e ϕθϕθ---⎡⎢⎢⎢⎢⎢⎣ 22()(1)()1...j j M e e ϕθϕθ---.........()(1)()1...p p j j M e e ϕθϕθ---⎤⎥⎥⎥⎥⎥⎦(2-2) （2()sin()i i d πϕθθλ=, 1,2,...i M =）向量()i θa 称为第i 个信号源的方向向量。

矩阵A 列向量即是阵列方向向量的集合，它表示所有信源的方向，称为阵列的方向矩阵。

阵列的协方差矩阵R 定义为2[()()]H H E t t σ==+S R X X AR A I ,式中[()()]H E t t =S R S S 为信源的协方差矩阵。

当空间信源互不相干时，对协方差矩阵R 进行特征分解，构造信号子空间和噪声子空间，利用其正交性直接采用MUSIC 或ESPRIT 算法进行DOA估计，令其特征值分解为H =∑s R U U ,式中，221(,......,)Mdiag σσ=∑，如果信号不相干，则A 为列满秩，()rank P =H S AR A ，于是221(,...,,0,...0)P diag αα=H H S U AR A U ，将R 同时左乘HU 和右乘U ，可以得到：22221(,...,,0,...0)H H Pdiag σαασ++H H S U RU =U AR A U U U =I ，所以自相关矩阵R 的特征值为2222221,...,,...Pασασσσ++，前P 个特征值为信号特征值, 后M-P 个特征值为噪声特征值, 由信号特征值对应的所有特征向量形成的矢量空间称为信号子空间, 信号子空间同时可以用信号的方向矢量()i θa 表示。

由噪声特征值对应的特征向量形成的矢量空间称为噪声子空间, 二者组成信号自相关矩阵的信息空间,利用信号矢量与噪声空间的所有矢量都是正交的这一特性来估计来波信号的到达角（如式2-3所示）, 即MUSIC 法。

1=()H HPθθθ（）（）a UU a （2-3） 代码如下：仿真结果如下：图1 无相干信号源时的 MUSIC 测向MUSIC 测向波达方向M U S I C 谱但若存在相干信源时，阵列输出信号协方差的秩()rank P <R ，对信号协方差矩阵进行特征值分解后，得到的较大的特征值个数小于P ，而特征值为2σ 的个数将大于M − P 。

与此相对应的信号子空间的向量也少于P ，即特征向量展开的信号子空间的维数少于1()[(),......()]P θθθ=A a a 的列数。

对某些相干源的方向矢量()i θa ，1,2,...i P = 将不正交于噪声子空间，不出现零点，所以，有些源在空间谱曲线中将不呈现峰值，造成谱估计的漏报。

在这里我们不妨将上面程序的信号源是S1和S2设置成相干信号，S2=2S1； 仿真结果如下：图2 有相干信号源时的MUSIC 测向因此，我们要对阵列输出的协方差矩阵首先进行预处理，使其阵列协方差矩阵的秩恢复为信号元数P ，然后再采用MUSIC 或ESPRIT 算法。

MUSIC 测向波达方向M U S I C 谱2.2 空间平滑算法（spatial smoothing ）2.2.1 前向空间平滑算法将M 个阵元的均匀线阵，分成相互交错的P 个子阵，每个子阵包含的阵元数为m 个，即满足M ＝p+m-1。

信号源数为N 。

图3 前向空间算法原理图如图 3所示，取第一个子阵(最左边的子阵)为参考子阵，那么各个子阵的输出矢量分别为:11222311[,,...,][,,...,]...[,,...,]f m f m f p p p M x x x x x x x x x ++===X X X （2-4）对于第k 个子阵有：(1)11()[,,...,]()()()f k k k k k m m k t x x x t t θ-++-==+X A D s n （2-5）其中：12sin()0...0dj e D πθλ⎡⎢⎢=⎢⎢⎢⎣ 22sin()0...0d j e πθλ ............ 2sin()00...N dj e πθλ⎤⎥⎥⎥⎥⎥⎦（2-6） 那么该子阵的数据协方差矩阵为：(1)(1)2()(())k k H k m s m R θθσ--=+R A DA D I （2-7） 其中，m A 是一个m ×p 的参考子阵(通常取第一个子阵)的导向矢量矩阵， 1()[(),......()]m m m N θθθ=A a a ，22sin()sin()()[1,,....,]k k ddj jT m k eeππθθλλθ--=a ，s R 为信号的协方差矩阵，}{Hs E ss =R 。

前向空间平滑技术是通过求各个子阵协方差矩阵的均值来实现的，即取前向平滑修正的协方差矩阵为:11pfk p ==∑k R R （2-8）可以证明，当满足m>N, p>N 时，前向空间平滑数据协方差矩阵f R 是满秩的。

即可以通过特征分解求得相应的信号子空间和噪声子空间。

2.2.2 前后向空间平滑算法如果按照图4划分阵列，即称为后向平滑的方法划分子阵，那么各个子阵的输出矢量为:图4 后向空间平滑算法原理图11121211[,,...,][,,...,]...[,,...,]b M M M m b M M M m b p m m x x x x x x x x x --+----===X X X (2-9)那么，第k 个子阵的数据矢量为：*12()[,,...,]bk M k M k M m k t x x x -+---+=X （2-10）比较前向平滑和后向平滑的数据矢量，可以得到前向平滑中第k 个子阵与后向平滑中第p-k+1个子阵之间存在如下关系：**(1)**1()(())()()bfk p k k m k t J X t J s t Jn t ---+==+X A D（2-11）其中J 为m 的交换矩阵。

00...1⎡⎢⎢=⎢⎢⎣J 0 0...1...010...0⎤⎥⎥⎥⎥⎦，所以后向平滑第p-k+1个子阵的数据协方差矩阵为:*(1)*(1)*21()()b k k H Hp k m m σ-----+=+s R JA DR D A J I (2-12) 那么后向空间平滑修正的数据矩阵为:111p bbp k k p -+==∑R R (2-13)取前向平滑和后向平滑数据协方差矩阵的平均，即前后向空间平滑的数据矩阵，即2b ffb+=R R R (2-14)同样可以证明，当满足m>N, p>N 时，后向空间平滑数据协方差矩阵b R 是满秩的。

3.Matlab 仿真程序代码4.Matlab 仿真结果分析图5 MUSIC 算法与改进MUSIC 算法的比较由此可见，空间平滑算法有效的达到了去相干的目的，较为准确的测出了方向。

下面将研究将子阵列划分的个数L 变化时的情形,仿真波形如下：-100-80-60-40-20020*********MUSIC 测向波达方向M U S I C 谱图6 测向随L 的变化由此可见，当空间中存在两个相干信号时，L=2时已经不能正确测向，因此L>信号源数。

5.参考文献[1]通信原理与MATLAB 仿真教程 [2]MATLAB 及其在理工课程中的应用指南[3]论文：一种基于加权空间平滑的新MUSIC 算法 [4] 谱估计测向的MUSIC 算法研究MUSIC 测向波达方向M U S I C 谱-1001020304050MUSIC 测向波达方向M U S I C 谱-100-50050100MUSIC 测向波达方向M U S I C 谱-100-5005010000.20.40.60.81MUSIC 测向波达方向M U S I C 谱。