• fanyi
    了解作者
  • matlab
    开发工具
  • 909B
    文件大小
  • rar
    文件格式
  • 0
    收藏次数
  • 5 积分
    下载积分
  • 0
    下载次数
  • 2022-08-18 21:50
    上传日期
MUSIC算法实现,可以跑的代码,优质代码。。。。
music.rar
  • music.m
    1.2KB
内容介绍
clc; clear; derad = pi/180; %角度->弧度 N = 15; % 阵元个数 M = 7; % 信源数目 theta = [-50 -30 -10 30 40 50 80]; % 待估计角度 snr = 10; % 信噪比 K = 1024; % 快拍数 dd = 0.5; % 阵元间距 d=lamda/2 d=0:dd:(N-1)*dd; A=exp(-1i*2*pi*d.'*sin(theta*derad));%方向相应矢量(N*M) S=randn(M,K); %信源信号,入射信号 X=A*S; % 构造接受信号 X1=awgn(X,snr,'measured'); %将白色高斯噪声添加到信号中 %% Rxx=X1*X1'/K; %协方差矩阵 [EV,D]=eig(Rxx) ; % 特征值分解 EVA=diag(D)'; %将特征值矩阵对角线提取并转为一行 [EVA,I]=sort(EVA); % 从小到大排列 返回索引 EV=fliplr(EV(:,I)); %对应特征矢量排序(左右翻转,即如今排序为从大到小特征值对应的特征向量) En=EV(:,M+1:N); % 取矩阵的第M+1到N列组成噪声子空间 for iang = 1:361 % 遍历 angle(iang)=(iang-181)/2; %从-90到90 phim=derad*angle(iang); a=exp(-1i*2*pi*d*sin(phim)).'; P(iang)=1/(a'*En*En'*a); end P=abs(P); Pmax=max(P); P=10*log10(P/Pmax); % 归一化处理 h=plot(angle,P); set(h,'Linewidth',2); xlabel('入射角/(degree)'); ylabel('空间谱/(dB)'); set(gca, 'XTick',[-90:20:90]); grid on;
评论
    相关推荐