信号检测与估计仿真作业绪论.docxVIP

PAGE  PAGE 18 电 子 科 技 大 学 信号检测与估计计算机仿真作业 2013年12月6日  一． 实验目的 1.学习Matlab软件在信号检测与估计中的应用 2.学习MUSIC、ESPRIT、GEESE等的空间谱估计算法的原理，并通过仿真分析比较这三种算法的不同及性能特点 3.通过仿真分析了解非平稳噪声和色噪声对MUSIC、ESPRIT、GEESE方法性能的影响 二． 实验原理 2.1最小错误概率准则 出发点是如何使译码后的错误概率PE为最小。其基本思路为：收到yj后，对于所有的后验概率P(x1|yj),P(x2|yj), …,P(xi|yj),…，若其中P(x*|yj)具有最大值，则将x*判决为yj的估值。由于这种方法是通过寻找最大后验概率来进行译码的，故又常称之为最大后验概率准则。 最大后验概率译码方法是理论上最优的译码方法，但在实际译码时，既要知道先验概率又要知道后验概率，而后验概率的定量计算有时比较困难，需要寻找更为实际可行的译码准则。 2.2 MUSIC原理 MUSIC算法是一种基于矩阵特征空间分解的方法。从几何角度讲，信号处理的观测空间可以分解为信号子空间和噪声子空间，显然这两个空间是 HYPERLINK /view/935574.htm \t _blank 正交的。信号 HYPERLINK /view/1279998.htm \t _blank 子空间由阵列接收到的数据 HYPERLINK /view/1304423.htm \t _blank 协方差矩阵中与信号对应的 HYPERLINK /view/475996.htm \t _blank 特征向量组成，噪声子空间则由协方差矩阵中所有最小 HYPERLINK /view/689250.htm \t _blank 特征值（噪声方差）对应的特征向量组成。MUSIC算法就是利用这两个互补空间之间的 HYPERLINK /view/935574.htm \t _blank 正交特性来估计空间信号的方位。噪声子空间的所有 HYPERLINK /view/77260.htm \t _blank 向量被用来构造谱，所有空间方位谱中的峰值位置对应信号的来波方位。MUSIC算法大大提高了测向分辨率，同时适应于任意形状的天线阵列，但是原型MUSIC算法要求来波信号是不相干的。 2.3 ESPRIT算法原理 ESPRIT算法估计信号参数时要求阵列的几何结构存在所谓的不变性，这个不变性可以通过两种手段得到：一是阵通过某些变换获得两个或两个以上的相同子阵。由于这种算法在有效性和方面都有非常突出的表现，已经被公认为空间谱估计的一种经典算法，随着ESPRIT算法的深入研究，ESPRIT算法进一步被广大学者接受并推广。 2.4 GEESE算法基本原理 信号子空间特征向量的广义特征值法(GEESE)，可以在简化计算的情况下解决ESPRIT算法中实际噪声测量有误差的问题。它利用信号子空间的一个显著特征，那就是真实方向向量所张成的子空间与除了阵列输出互相关矩阵的最小多重特征值之外的所有相应特征向量所张成的子空间是一样的。 2.5非相关源的数学模型 图2.1 .阵元接收信号与位置的关系 阵元接收信号与位置的关系如图 2-1 所示。假设空间有M个阵元组成阵列，将阵元从1到M编号，并以阵元 1 作为参考点。由于各阵元无方向性，相对于基准点的位置向量分别为令信源信号为 s(t),信号的载波为，则基准点处的接收信号为,各阵元上的接收信号的表达式为式中k为波数向量，为入射信号传播的传播方向，单位向量，为信号相对于基准点的延迟时间，为信号传播到基准点r处的阵元相对于信号传播到基准点的滞后相位（弧度）。图中为入射信号传播方向角，k = k[,]。 天线阵列中，信号的带宽B一般比载波频率ω小得多，所以就有，即信号在各阵元上的差异可以忽略不计，称为窄带信号。 因此，阵列信号用向量形式可以表示为： 选定第一个阵元为基准点，则方向向量为 式中， 当有P个信源时，波束的方向向量可分别表达为，，…，。这P个方向向量组成的矩阵 称为阵列的方向矩阵或响应矩阵，它表示所有信源的方向信息。 当有N个窄带信号入射到空间 M 个阵元上时候，接收的信号可以写成如下的矢量形式： 式中，X(t)为阵列的M×1维快拍数据矢量，N(t)为阵列的 M ×1维噪声数据矢量， S(t)为空间信号的N×1维矢，A空间阵列的M×N维响应矩阵（导向矢量阵）。 三．实验过程 3.1 实验1 首先，当M=1，