信号检测与估计仿真作业

1信号检测与估计计算机仿真作业一．实验目的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.2MUSIC原理MUSIC算法是一种基于矩阵特征空间分解的方法。从几何角度讲，信号处理的观测空间可以分解为信号子空间和噪声子空间，显然这两个空间是正交的。信号子空间由阵列接收到的数据协方差矩阵中与信号对应的特征向量组成，噪声子空间则由协方差矩阵中所有最小特征值（噪声方差）对应的特征向量组成。MUSIC算法就是利用这两个互补空间之间的正交特性来估计空间信号的方位。噪声子空间的所有向量被用来构造谱，所有空间方位谱中的峰值位置对应信号的来波方位。MUSIC算法大大提高了测向分辨率，同时适应于任意形状的天线阵列，但是原型MUSIC算法要求来波信号是不相干的。2.3ESPRIT算法原理ESPRIT算法估计信号参数时要求阵列的几何结构存在所谓的不变性，这个不变性可以通过两种手段得到：一是阵通过某些变换获得两个或两个以上的相同子阵。由于这种算法在有效性和方面都有非常突出的表现，已经被公认为空间谱估计的一种经典算法，随着ESPRIT算法的深入研究，ESPRIT算法进一步被广大学者接受并推广。2.4GEESE算法基本原理信号子空间特征向量的广义特征值法(GEESE)，可以在简化计算的情况下解决ESPRIT算法中实际噪声测量有误差的问题。它利用信号子空间的一个显著特征，那就是真实方向向量所张成的子空间与除了阵列输出互相关矩阵的最小多重特征值之外的所有相应特征向量所张成的子空间是一样的。22.5非相关源的数学模型图2.1.阵元接收信号与位置的关系阵元接收信号与位置的关系如图2-1所示。假设空间有M个阵元组成阵列，将阵元从1到M编号，并以阵元1作为参考点。由于各阵元无方向性，相对于基准点的位置向量分别为1(1,2,...,,0)iiNrr令信源信号为s(t),信号的载波为iwte，则基准点处的接收信号为()jwtste,各阵元上的接收信号的表达式为1()()exp[()]Tiiicststrjwtrk式中k为波数向量，为入射信号传播的传播方向，单位向量，1Ticr为信号相对于基准点的延迟时间，irk为信号传播到基准点r处的阵元相对于信号传播到基准点的滞后相位（弧度）。图中为入射信号传播方向角，k=k[cos,sin]。天线阵列中，信号的带宽B一般比载波频率ω小得多，所以就有1()()Ticstrst，即信号在各阵元上的差异可以忽略不计，称为窄带信号。因此，阵列信号用向量形式可以表示为：1212()[(),(),...,()]()[,,...,]TTTNTNjrkjrkjrkTstststststeee选定第一个阵元为基准点，则方向向量为2()[1,,...,]TTNjrkjrkTee式中，1iirrr当有P个信源时，波束的方向向量可分别表达为1()，2()，…，()p。这P个方向向量组成的矩阵12[(),(),...,()]pA称为阵列的方向矩阵或响应矩阵，它表3示所有信源的方向信息。当有N个窄带信号入射到空间M个阵元上时候，接收的信号可以写成如下的矢量形式：()()()XtAStNt式中，X(t)为阵列的M×1维快拍数据矢量，N(t)为阵列的M×1维噪声数据矢量，S(t)为空间信号的N×1维矢，A空间阵列的M×N维响应矩阵（导向矢量阵）。三．实验过程3.1实验1首先，当M=1，实验1就变成了一般二元随机振幅与随机相位信号的波形检测问题。这样，就有如下模型：𝐻0:𝑥(𝑡)=𝐵1cos(𝜔0𝑡+𝜃0)+𝑛(𝑡),0≤t≤T𝐻1:𝑥(𝑡)=𝐴1cos(𝜔1𝑡+𝜃0)+𝑛(𝑡),0≤t≤T其中，噪声n(t)是均值为零，功率谱密度为𝑃𝑛(𝜔)=𝑁02⁄的高斯白噪声；信号的振幅服从瑞丽分布，随机相位服从均匀分布。由于功率信噪比早0~60dB范围内，取𝐴0=100，𝑁0=2，T=1us。根据教材的推导，得到错误判决概率的表达式如下：P(𝐻0|𝐻1)=𝑁02𝑁0+𝜎𝑎2𝑇P(𝐻1|𝐻0)=𝑁02𝑁0+𝜎𝑎2𝑇从上面可以看到，当M=1时，模型很简单。当M1时，我们选取一个M值，用MATLAB软件进行仿真，得到仿真结果如下图所示（MATLAB仿真程序见附录1）：43.2实验23.2.1.MUSIC算法过程及仿真1)．收集阵列接收的数据样本,0,1,,1ixip得到数据协方差矩阵xxR1,,01pHxxiiiRxxP2)．对协方差矩阵xxR进行特征分解，分解为特征值和特征向量的形式xxRVV式中，0,1,1,[]NVqqq是xxR对应的特征向量所组成的矩阵，011011{,},NNVdiag为xxR的特征值。3)．利用最小特征值min的重数M来估计信号源数KKNM4)．计算MUSIC算法的功率谱()()()()HMUSICHHNNaaPaVVa5)．MUSICP中极大值对应的角度就是信号入射方向。在Matlab的命令窗口输入仿真程序，见附录2，得到如下结果：A=-102345DOA=44.998623.0019-10.05195图3.1MUSIC算法非相关源的模拟测向仿真3.2.2ESPRIT算法过程及仿真1)．由快拍数据X可以得到数据相关矩阵R的估计R2)．对相关矩阵的估计R做特征分解HRUU（33）特征值矩阵011{,}Ndiag，0,1,1,[]NUqqq。3)．利用小特征值的重数估计M得到信号源的估计数KNM。4)．将特征向量矩阵分解为子阵列矩阵得到：12ssUUU5).得到112()ssUU6).对Ψ进行特征分解，得到K个特征值，就可以得到对应K个信号的到达角。在Matlab的命令窗口输入仿真程序，见附录3，得到如下结果：A=-102345DOA=44.998623.0019-10.05196图3.2ESPRIT算法非相关源的模拟测向仿真3.2.3GEESE算法：1)．收集阵列接收的数据样本,0,1,...,1ixiP得到数据协方差矩阵xxR101PHxxiiiRxxP2).对xxR进行特征分解，得到sV和s：xxsssRVV3).选定J的值，将sV代入式中得到1E,2E:1,[]JJNJsEIOV2,1,1[]JJJNJsEOIOV4).将1SV，2SV代入式得到k；5).将k代入式中即得到信号源的波达方向:arg()kk1,2,...,kK在Matlab的命令窗口输入仿真程序，见附录4，得到如下结果：DOA=-30.024470.000030.0244图3.3GEESE算法非相关源的模拟测向仿真3.2.4三种算法性能比较MUSIC算法就是多重信号分类算法，它是一种信号参数估计算法，利用输入信号协方差矩阵的特征结构，给出的信息包括入射信号的数目、各个信号的波达方向、强度以及入射信号和噪声间的互相关。ESPRIT算法就是旋转不变子空间算法，也是一种基于子空间的波达方向估计技术，与MUSIC算法不同的是，ESPRIT算法不需要精确知道阵列的方向向量，仅需各子需各子阵列之间保持一致，因此降低了对阵列校准的严格性。GEESE算法是指信号子空间特征向量的广义特征值法，可以在简化计算的情况下解决ESPRIT算法中实际噪声测量有误差的问题。它利用信号子空间的一个显著特征，那就是真实方向向量所张成的子空间与除了阵列输出互相关矩阵的最小多重特征值之外的所有相应特征向量所张成的子空间是一样的。这三种算法是空间谱估计中最经典的算法。MUSIC算法估计值接近克拉美—罗界算法(CRB)，对参数的少量偏差不敏感，更接近实际应用，具有较好的应用前景,但需要对参数空间进行搜索，计算量大。随着信噪比的增加，MUSIC功率谱的峰值越高，估计精度越精确。在阵元数目不同，其他条件相同的情况下，阵元数目越大，旁瓣干扰越小，DOA估计越精确。在条件相同的情况下，相邻信号（以50为例）的MUSIC功率谱随着角度的增加而降低，信号源相关，MUSIC算法失效。色噪声下，MUSIC算法方位估计不准确。与MUSIC算法相比，ESPRIT算法还降低了计算量和存储量，且避免了参数空间的搜索，计算量小于MUSIC算法，但是算法数据协方差矩阵中提取噪声方差的8估计，有时会使估计结果变坏，当信号高度相关时估计性能同样会变坏，且对所设的参数有较高的要求，少量的误差也会导致算法的失败。在ESPRIT算法中随着信噪比的增加，均方误差越小，DOA估计效果越好。在阵元数目不同，其他条件相同的情况下，阵元数目越大，均方误差越小，ESPRIT算法的估计精度越高。在条件相同的情况下，相邻信号（以100为例）的均方差与信噪比关系随着角度的增加而性能降低。ESPRIT算法对相干信号的DOA估计失效。而GEESE算法，不仅计算量比较小，而且保证了算法的精度，但当信号高度相关时性能仍然变坏。在GEESE算法中随着信噪比的增加，均方误差越小，DOA估计效果越好。在阵元数目不同，其他条件相同的情况下，阵元数目越大，均方误差越小，GEESE算法的估计精度越高。在条件相同的情况下，相邻信号（以50为例）的均方差与信噪比关系随着角度的增加而性能降低，GEESE算法对相干信号的DOA估计失效。3.3实验33.3.1非平稳噪声下，三种算法对DOA估计影响的matlab仿真结果图3.4非平稳噪声下MUSIC算法对DOA估计的影响9图3.5非平稳噪声下ESPRIT算法对DOA估计的影响图3.6非平稳噪声下GEESE算法对DOA估计的影响103.3.2色噪声噪声下，三种算法对DOA估计影响的matlab仿真结果图3.7色噪声下MUSIC算法对DOA估计的影响图3.8色噪声下ESPRIT算法对DOA估计的影响11图3.9色噪声下GEESE算法对DOA估计的影响附录1第1题仿真程序clc;clear;closeallM=[1246810];%M为接收机个数N=[1510];%N为射频脉冲数fori=1:1Estimate=zeros(3,11);figureforj=1:3forsnr=0:60all_err=0;forcon=1:1100Error=mreceiver(M(5),N(j),snr);%计算j个接收机和i个射频脉冲数的最小错误概率all_err=all_err+sum(Error);endEstimate(j,snr+1)=all_err/2200.0;endendholdonplot(Estimate(1,:),'r');plot(Estimate(2,:),'b-');plot(Estimate(3,:),'m');xlabel('SNR=0~60dB');12ylabel('估计误差');title('16个接收单元最小误判概率仿真');legend('1个射频脉冲数','5个射频脉冲数','10个射频脉冲数');gridonholdoff;endfunctionda=anc_mac(x,y)if(size(x)~=size(y))e