%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%MUSIC算法%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear;clc;closeall;tic;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%M=8;L=1024;N=3;%%P信号个数M天线阵元数L采样点数thetaa=[-30-1020];snr=20;f0=10e6;C=3e8;fs=200;Ts=1/fs;lamda0=C/f0;d=0.5*lamda0;%%%%%%%%%%%%产生信号%%%%%%%%%%%%%%%%u0=-15;u1=0;u2=15;t=(0:L-1)*Ts;s0=exp(j*2*pi*(f0*t+0.5*u0*t.^2));s1=exp(j*2*pi*(f0*t+0.5*u1*t.^2));s2=exp(j*2*pi*(f0*t+0.5*u2*t.^2));S=[s0;s1;s2;];%%%%%%%%%%%生成导向矢量矩阵%%%%%%%%a=[0:M-1]';a_theta0=exp(j*2*pi*a*d*sin(thetaa(1)/180*pi)/lamda0);a_theta1=exp(j*2*pi*a*d*sin(thetaa(2)/180*pi)/lamda0);a_theta2=exp(j*2*pi*a*d*sin(thetaa(3)/180*pi)/lamda0);A=[a_theta0a_theta1a_theta2];X0=A*S;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%real_noise=randn(size(X0));imag_noise=randn(size(X0));noise=(real_noise+j*imag_noise)/2^0.5;noise=10^(-snr/20)*noise;%%%%%阵元接受数据%%%%%%%%%%%%X1=X0+noise;%%%%%%%%无噪声经典MUSIC算法%%%%%%%%%%Rx0=X0*X0'/L;[V0,D0]=eig(Rx0);%%%%%%%%求矩阵的特征值%%%%%%%E0=diag(D0);%%%%%%%%提取D矩阵的对角线构成一个一列矩阵%%%%%%%%Un0=V0(:,1:M-N);%%%%%%%%经典MUSIC算法%%%%%%%%%%Rx1=X1*X1'/L;[V1,D1]=eig(Rx1);%%%%%%%%求矩阵的特征值%%%%%%%E1=diag(D1);%%%%%%%%提取D矩阵的对角线构成一个一列矩阵%%%%%%%%Un1=V1(:,1:M-N);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%theta=-90:0.1:90;%%%%%从-90--90，每格取0.01%%%%forff=1:length(theta)a_theta=exp(j*2*pi*d/lamda0*a*sin(theta(ff)/180*pi));Pmusic0(ff)=100*log10(1/abs((a_theta'*Un0*Un0'*a_theta)));Pmusic1(ff)=100*log10(1/abs((a_theta'*Un1*Un1'*a_theta)));end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%figure(1)plot(theta,Pmusic0);xlabel('DOA(度)'),ylabel('Pmusic');title('无噪声MUSIC算法')gridon%%%%%%%%%figure(2);plot(theta,Pmusic1);xlabel('DOA(度)'),ylabel('Pmmusic');title('有噪声MUSIC算法')gridontoc%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
本文档为【music算法程序】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483