二维ROOT-MUSIC算法matlab代码

二维ROOT-MUSIC DOA估计算法的MATLAB实现代码，可以参考学习。

%ROOT_MUSIC ALOGRITHM
%DOA ESTIMATION BY ROOT_MUSIC
clear all;
close all;
clc;
tic
N=100;
source_number=2;%信元数
sensor_number=8;%阵元数
N_x=1024; %信号长度
snapshot_number=N_x;%快拍数
w=[pi/4 pi/6].‘;%信号频率
l=（（2*pi*3e8）/w（1）+（2*pi*3e8）/w（2））/2;%信号波长  
d=0.5*l;%阵元间距
snr=0;%信噪比

theta1=50;theta2=60;%两个信号的入射角度
fai1=60;fai2=50;%方位角

Ax=[exp（-j*（0:sensor_number-1）*d*2*pi*cos（fai1*pi/180）*sin（theta1*pi/180）/l）;exp（-j*（0:sensor_number-1）*d*2*pi*cos（fai2*pi/180）*sin（theta2*pi/180）/l）].‘;%阵列流型

s=10.^（（snr/2）/10）*exp（j*w*[0:N_x-1]）;%仿真信号
%x=awgn（ssnr）;
for n=1:N
x=Ax*s+（1/sqrt（2））*（randn（sensor_numberN_x）+j*randn（sensor_numberN_x））;%加了高斯白噪声后的阵列接收信号

R=（x*x‘）/N_x;
%对协方差矩阵进行特征分解
[VD]=eig（R）;
Un=V（:1:sensor_number-source_number）;
Gn=Un*Un‘;

%找出多项式的系数，并按阶数从高至低排列
a = zeros（2*

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        3470  2019-04-27 22:55  二维Root_MUSIC.m