MUSIC算法.doc

专业综合
课程设计报告

空间谱估计算法设计任务
实现空间谱估计算法，并考察算法性能。
方案设计
由均匀线阵形式，确定阵列的导向矢量；
由阵列导向矢量，对接收信号进行建模仿真；
根据多重信号分类算法实现空间谱估计；
考察算法性能与信噪比，采样率，观测时间等参数的关系。
设计原理
3.1空间谱估计数学模型
空间谱估计就是利用空间阵列实现空间信号的参数估计的一项专门技术。整个空间谱估计系统应该由三部分组成：空间信号入射、空间阵列接收及参数估计。相应地可分为三个空间，即目标空间、观察空间及估计空间，也就是说空间谱估计系统由这三个空间组成，其框图见图1。

图1空间谱估计的系统结构
对于上述的系统结构，作以下几点说明。
(1)目标空间是一个由信号源的参数与复杂环境参数张成的空间。对于空间谱估计系统，就是利用特定的一些方法从这个复杂的目标空间中估计出信号的未知参数。
(2)观察空间是利用空间按一定方式排列的阵元，来接收目标空间的辐射信号。由于环境的复杂性，所以接收数据中包括信号特征(方位、距离、极化等)和空间环境特征(噪声、杂波、干扰等)。另外由于空间阵元的影响，接收数据中同样也含有空间阵列的某些特征(互耦、通道不一致、频带不一致等)。这里的观察空间是一个多维空间，即系统的接收数据是由多个通道组成，而传统的时域处理方法通常只有一个通道。特别需要指出的是：通道与阵元并不是一一对应，通道是由空间的一个、几个或所有阵元合成的(可用加权或不加权)，当然空间某个特定的阵元可包含在不同的通道内。
(3)估计空间是利用空间谱估计技术(包括阵列信号处理中的一些技术，如阵列校正、空域滤波等技术)从复杂的观察数据中提取信号的特征参数。
从系统框图中可以清晰的看出，估计空间相当于是对目标空间的一个重构过程，这个重构的精度由众多因素决定，如环境的复杂性、空间阵元间的互耦、通道不一致、频带不一致等。
3.2阵列信号处理
首先，考虑N个远场的窄带信号入射到空间某阵列上，阵列天线由M个阵元组成，这里假设阵元数等于通道数，即各阵元接收到信号后经过各自的传输信道送到处理器，也就是说处理器接收来自M个通道的数据。

(3.2-1)
式中，是接受信号的幅度，是接收信号的相位，是接收信号的频率。在窄带远场信号源的假设下，有
(3.2-2)
根据式(3.2-1)和式(3.2-2)，显然有下式成立：
(3.2-3)
则可以得到第L个阵元接收信号为
(3.2-4)
式中，为第L个阵元对第i个信号的增益，表示第L个阵元在t时刻的噪声，表示第i个信号到达第L个阵元时相对参考阵元的时延。
将M个阵元在特定时刻接收的信号排列成一个列矢量，可得
(3.2-5)
在理想情况下，假设阵列中各阵元是各向同性的且不存在通道不一致、互耦等因素的影响，则式(3.2-4)中的增益可以省略(即归一化1)，在此假设下式(3.2-5)可以简化为
(3.2-6)
将式(3.2-6)写成矢量形式如下：
(3.2-7)
式中，为阵列的维快拍数据矢量，为阵列的维噪声数据矢量，为空间信号的维矢量，A为空间阵列的维流型矩阵(导向矢量阵)，且
(3.2-8)
其中导向矢量
(3.2-9)
式中，c为光速，为波长。
由上述的知识可知，一旦知道阵元间的延迟表达式τ，就很容易得出待定空间阵列的导向矢量或阵列流型。下面推导一下空间阵元间的延迟表达式。假设空间任意两个阵元，其中一个为参考阵元(位于原点)，另一个阵元的坐标为(x,y,z)，两阵元的几何关系见图，图中“×”表示阵元。

图2空间任意两阵元的几何关系
由几何关系可以推导出两阵元的波程差为
(3.2-10)
这里的波程差其实就是位于x轴上两阵元间的延迟、位于y轴上两阵元间的延迟和位于z轴上两阵元间的延迟之和。
根据式(3.2-10)的结论，下面给出实际环境中常用的几种阵列及阵元间的相互延迟表达式。
(1)平面阵设阵元的位置为，以原点为参考点，另假设信号入射参数为，分别表示方位角与俯仰角，其中方位角表示与x轴的夹角。
(2)线阵设阵元的位置为，以原点为参考点，另假设信号入射参数为，表示方位角，其中方位角表示与y轴的夹角（即与线阵法线的夹角），则有
(3.2-11)
(3)均匀圆阵设以均匀圆阵的圆心为参考点，则有
(3.2-12)
其中方位角表示与x轴的夹角，r为圆半径。
3.3多重信号分类算法
多重信号分类(MUSIC)算法是SchmidtRO等人在1979年提出的。这一算法的提出促进了特征结构分类算法的兴起和发展，该算法已成为空间谱估计理论体系中的标志性算法。此算法提出之前的有关算法都是针对阵列接收数据协方差矩阵进行直接处理，而MUSIC算法的基本思想则是将任意阵列输出数据的协方差矩阵进行特征分解，从而得到与信号分量相对应的信号子空间和与信号分量相正交的噪声子空间，