• 大小: 972B
    文件类型: .rar
    金币: 1
    下载: 0 次
    发布日期: 2021-05-28
  • 语言: Matlab
  • 标签: music  算法  

资源简介

经典MUSIC算法MATLAB仿真,带详细注释,仿真通过。可以直接使用。MATLAB阵列信号处理,DOA估计,MATLAB信号处理,信息与信号处理 专业,通信专业等学科学习的例程,经典仿真案例,学习入门必备,程序正确,可运行。

资源截图

代码片段和文件信息

% DOA estimation by MUSIC 
% Developed by xiaofei zhang (南京航空航天大学 电子工程系 张小飞)
% EMAIL:zhangxiaofei@nuaa.edu.cn
clear all
close all
derad = pi/180;        % deg -> rad
radeg = 180/pi;
twpi = 2*pi;
kelm = 8;               % 阵列数量
dd = 0.5;               % space 
d=0:dd:(kelm-1)*dd;     % 
iwave = 3;              % number of DOA
theta = [10 30 60];     % 角度
snr = 10;               % input SNR (dB)
n = 500;                 % 
A=exp(-j*twpi*d.‘*sin(theta*derad));%%%% direction matrix
S=randn(iwaven);
X=A*S;
X1=awgn(Xsnr‘measured‘);
Rxx=X1*X1‘/n;
InvS=inv(Rxx); %%%%
[EVD]=eig(Rxx);%%%% 
EVA=diag(D)‘;
[EVAI]=sort(EVA);
EVA=fliplr(EVA);
EV=fliplr(EV(:I));

% MUSIC
for iang = 1:361
        angle(iang)=(iang-181)/2;
        phim=derad*angle(iang);
        a=exp(-j*twpi*d*sin(phim)).‘;
        L=iwave;    
        En=EV(:L+1:kelm);
        SP(iang)=(a‘*a)/(a‘*En*En‘*a);
end
   
% 
SP=abs(SP);
SPmax=max(SP);
SP=10*log10(SP/SPmax);
h=plot(angleSP);
set(h‘Linewidth‘2)
xlabel(‘angle (degree)‘)
ylabel(‘magnitude (dB)‘)
axis([-90 90 -60 0])
set(gca ‘XTick‘[-90:30:90])
grid on  





 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1122  2010-08-16 09:41  张小飞-MUSIC算法MATLAB程序\doa_music.m

     目录          0  2017-11-02 16:19  张小飞-MUSIC算法MATLAB程序

----------- ---------  ---------- -----  ----

                 1122                    2


评论

共有 条评论