matlab经典music,基于经典music的doa估计matlab仿真.doc

基于经典music的doa估计matlab仿真.doc

信息与通信工程学院阵列信号处理实验报告(基于经典MUSIC的DOA估计MATLAB仿真)学号XXXXXX专业XXXXXX学生姓名XXX任课教师XXX2015年X月题目基于经典MUSIC的DOA估计MATLAB仿真1算法简述基于天线阵列协方差矩阵的特征分解类DOA估计算法中，多重信号分类(MUSIC)算法具有普遍适用性，只要已知天线阵的分布形式，无论直线阵还是圆阵，不管阵元是否等间隔分布，都可以得到高分辨率的估计结果。阵列协方差矩阵R可以划分为两个空间，即。因导向矩阵的HSNRUΣ各矢量与噪声子空间正交，可得到阵列空间谱函数使变化，按照来搜索峰值来MUSICH1NPAMUSICP达到波达方向的估计值。2实验内容与结果实验使用8阵元均匀线阵，阵元间距为信号波长的一半，输入信号为1个BPSK信号，1个非相干的单频干扰，设置载波频率10MHZ、采样频率30MHZ、快拍数30、信干比0DB、信号方位角、干扰方位角，分析信噪比从051到20DB估计均方误差。实验结果见下图。100806040200204060801004020020406080估估估/DEGREE估估估/DB估估MUSIC估估MUSICSPECTRUM图1信噪比15DB情况下的波束图24681012141618200051X103SNRDB0°估估估估估估估24681012141618200051X103SNRDB5°估估估估估估估图2不同信噪比测得的均方误差3仿真分析由仿真结果可知，在实验前提条件下，经典MUSIC算法的DOA估计可有效估计出信号的来向。当信噪比较小时，估计得均方误差较大；随着信噪比的提高，估计均方误差逐渐减小；当信噪比大于10DB以后，角度估计成功率100。MUSIC算法在DOA估计中属于高精度算法，而算法的估计精度仍受阵元数量、信号角度、信号数量、信号间距、信噪比等影响，本实验仿真分析了一定情况下估计精度与信噪比的关系，得到了预期的实验结果。4程序CLEARALLCLOSEALLCLC线阵MUSICDOA参数设置SIGNAL_NO1INTERFERENCE_NO1S_NOSIGNAL_NOINTERFERENCE_NOSENSOR_NO8AZIMUTH0PI/1805PI/180FS30E6FC10E6F18E6干扰功率THETA1045529343PI/180RB2E6M2二进制WAVELENGTH3E8/FCDWAVELENGTH/2K300快拍数DBPS30SNR10SIR00M_NO40码速率等于信息速率DATA_NOM_NOFS/RBT1/FS1/FSDATA_NO/FSFORXUNHUAN11MEAN10MEAN20FORSNR120功率转换PS_L10PS/10PI_L10PSSIR/10PN_L10PSSNR/10信号生成BITRANDINT1,M_NO产生信息序列BITSTREAMFORI1M_NOIFBITI1BITSTREAMBITSTREAM,ONES1,FS/RBELSEBITSTREAMBITSTREAM,ONES1,FS/RBENDENDCARRIER_RCOS2PIFCTS_RCARRIER_RBITSTREAMCARRIER_ISIN2PIFCTS_ICARRIER_IBITSTREAMSIGNAL_RSQRTPS_LS_R/SQRTS_RS_R /LENGTHS_RSIGNAL_ISQRTPS_LS_I/SQRTS_IS_I /LENGTHS_I生成BPSK复信号SIGNAL1,COMPLEXSIGNAL_R,SIGNAL_I产生干扰I_RSQRT2PI_LCOS2PIF1TI_ISQRT2PI_LSIN2PIF1TFORI2INTERFERENCE_NO1SIGNALI,SQRT2PI_LCOMPLEXSIN2PIF1I1TTHETAI,SIN2PIF1I1TTHETAIEND模拟天线接收ASZEROSSENSOR_NO,S_NOFORI1SENSOR_NOFORII1S_NOASI,IIEXP1I2PISINAZIMUTHIID/WAVELENGTHI1ENDENDXASSIGNAL加噪声NOISE_RRANDNSENSOR_NO,DATA_NOSQRTPN_LNOISE_IRANDNSENSOR_NO,DATA_NOSQRTPN_LNOISECOMPLEXNOISE_R,NOISE_ISXNOISEX1S,1KRX1X1 /KV,DEIGRD1DIAGDUNV,1SENSOR_NOS_NOGNUNUN SEARCHING_DOA900190线阵的搜索范围为9090度FORI1LENGTHSEARCHING_DOAA_THETAEXPJ0SENSOR_NO1 2PIDSINPISEARCHING_DOAI/180/WAVELENGTHPMUSICI1/ABSA_THETA GNA_THETAENDPLOTSEARCHING_DOA,10LOGPMUSIC, R XLABEL 入射角/DEGREE YLABEL 空间谱/DB LEGEND MUSICSPECTRUM TITLE 经典MUSIC估计 GRIDON查最大值PMAX_DB10FORI1909IFPMUSICIPMAX_DB1PMAX_DB1PMUSICIPMAX_PITCHING1IENDENDPMAX_DB20FORI9101801IFPMUSICIPMAX_DB2PMAX_DB2PMUSICIPMAX_PITCHING2IENDENDERR1SNRPMAX_PITCHING11/10902ERR2SNRPMAX_PITCHING21/109052ENDMEAN1MEAN1ERR1MEAN2MEAN2ERR2ENDMSE1MEAN1/1000MSE2MEAN2/1000SUBPLOT211PLOTMSE1XLABEL SNRDB YLABEL 均方误差 SUBPLOT212PLOTMSE2XLABEL SNRDB YLABEL 均方误差 FIGURE2POLARPI/2PI/1800PI/2,10LOGPMUSIC33FIGUREGCF