classical_music.m

含有前向空间平滑、经典MUSIC算法、TOEPLITZ算法进行DOA估计的MATLAB源代码。

%MUSIC ALOGRITHM
%DOA ESTIMATION BY CLASSICAL_MUSIC
clear all;
close all;
clc;

source_number=2;%信元数
sensor_number=8;%阵元数
snapshot_number=256;%快拍数
w=[pi/4 pi/6].';%信号频率
l=[(2*pi*3e8)/w(1) (2*pi*3e8)/w(2)];%信号波长  
d=[0.5*l(1) 0.5*l(2)];%阵元间距
snr=0;%信噪比

theta1=-45;theta2=50;%两个信号的入射角度
A=[1 exp(-j*2*pi*d(1)*sin(theta1*pi/180)/l(1)) exp(-j*2*2*pi*d(1)*sin(theta1*pi/180)/l(1)) exp(-j*2*3*pi*d(1)*sin(theta1*pi/180)/l(1)) exp(-j*2*4*pi*d(1)*sin(theta1*pi/180)/l(1)) exp(-j*2*5*pi*d(1)*sin(theta1*pi/180)/l(1)) exp(-j*2*6*pi*d(1)*sin(theta1*pi/180)/l(1)) exp(-j*2*7*pi*d(1)*sin(theta1*pi/180)/l(1));...
   1 exp(-j*2*pi*d(2)*sin(theta2*pi/180)/l(2)) exp(-j*2*2*pi*d(2)*sin(theta2*pi/180)/l(2)) exp(-j*2*3*pi*d(2)*sin(theta2*pi/180)/l(2)) exp(-j*2*4*pi*d(2)*sin(theta2*pi/180)/l(2)) exp(-j*2*5*pi*d(2)*sin(theta2*pi/180)/l(2)) exp(-j*2*6*pi*d(2)*sin(theta2*pi/180)/l(2)) exp(-j*2*7*pi*d(2)*sin(theta2*pi/180)/l(2))].';%阵列流型
N_x=1024; %信号长度
s=exp(j*w*[0:N_x-1]);%仿真信号
%x=awgn(s,snr);
x=A*s+(1/sqrt(2))*(randn(sensor_number,N_x)+j*randn(sensor_number,N_x));%加了高斯白噪声后的阵列接收信号

R=(x*x')/snapshot_number;
disp(R);

%[V,D]=eig(R);
%Un=V(:,1:sensor_number-source_number);
%Gn=Un*Un';
[U,S,V]=svd(R);
Un=U(:,source_number+1:sensor_number);
Gn=Un*Un';

searching_doa=-90:0.1:90;%线阵的搜索范围为-90~90度
 for i=1:length(searching_doa)
   a_theta=exp(-j*(0:sensor_number-1)'*2*pi*d*sin(pi*searching_doa(i)/180)/l);
   Pmusic(i)=1./abs((a_theta)'*Gn*a_theta);
 end
plot(searching_doa,10*log(Pmusic));
%axis([-90 90 -90 90]);
xlabel('入射角/度');
ylabel('谱峰/dB');
legend('Music Spectrum');
title('经典MUSIC估计');
grid on;