muscialgorythm.txt

相关信号源的MUSIC估计

%MUSIC ALGORITHM 
%DOA Estimation by ULA 
clear all; 
close all; 
clc 
N_x=1024; % Length of Signal 
N=8; % Size of Rx Matrix 
A=[4 3 3]; 
l=1.8;%波长 
d=0.5*l;%阵元间距 
w=[pi/6 pi/10 pi/6]';%信号频率 
xx=pi/4;yy=0*pi/4;zz=-pi/4; %三个信号的入射角 
B=[1 exp(-j*2*pi*d*sin(xx)/l) exp(-j*2*2*pi*d*sin(xx)/l) exp(-j*2*3*pi*d*sin(xx)/l) exp(-j*2*4*pi*d*sin(xx)/l) exp(-j*2*5*pi*d*sin(xx)/l) exp(-j*2*6*pi*d*sin(xx)/l) exp(-j*2*7*pi*d*sin(xx)/l);... 
1 exp(-j*2*pi*d*sin(yy)/l) exp(-j*2*2*pi*d*sin(yy)/l) exp(-j*2*3*pi*d*sin(yy)/l) exp(-j*2*4*pi*d*sin(yy)/l) exp(-j*2*5*pi*d*sin(yy)/l) exp(-j*2*6*pi*d*sin(yy)/l) exp(-j*2*7*pi*d*sin(yy)/l);... 
1 exp(-j*2*pi*d*sin(zz)/l) exp(-j*2*2*pi*d*sin(zz)/l) exp(-j*2*3*pi*d*sin(zz)/l) exp(-j*2*4*pi*d*sin(zz)/l) exp(-j*2*5*pi*d*sin(zz)/l) exp(-j*2*6*pi*d*sin(zz)/l) exp(-j*2*7*pi*d*sin(zz)/l)]'; %8*3阵列流型 
M=3; % Number of Signals 
xxx=2*exp(j*w*[0:N_x-1]);%仿真信号 3*1024
x=B*xxx+randn(8,N_x)+j*randn(8,N_x); %加噪声 
R=x*x'; 
[V D]=eig(R); 
[lambda,index] = sort((diag(D))); 
UU=V(:,index(1:5)); 
theta=-pi/2:.001:pi/2;%ULA估计角度变化的范围和频率选择 
for i = 1:length(theta) 
AA=[1 exp(-j*2*pi*d*sin(theta(i))/l) exp(-j*2*2*pi*d*sin(theta(i))/l) exp(-j*2*3*pi*d*sin(theta(i))/l) exp(-j*2*4*pi*d*sin(theta(i))/l) exp(-j*2*5*pi*d*sin(theta(i))/l) exp(-j*2*6*pi*d*sin(theta(i))/l) exp(-j*2*7*pi*d*sin(theta(i))/l)];%方向矢量 
WW=AA*UU*UU'*AA'; 
Pmusic(i)=abs(1/WW);%角谱 
% PCap(i) = abs(1/(AA*inv(R)*AA')); 
end 
Pmusic = 10*log10(Pmusic/max(Pmusic)); 
plot(theta*180/pi,Pmusic); 
grid 
xlabel('rad')