capon算法Doa估计测角精度仿真matlab

capon算法DOA估计合集，包含随SNR，阵元个数的测角精度变化，以及与其他算法的对比，程序全部可用

clear all;
close all;
clc;
%%设置信号为窄带信号，幅度为1
S=[111]
imag1=sqrt（-1）;
element_num=8;%阵元数为8
d_lamda=1/2;%阵元间距d与波长lamda的关系
theta=linspace（-pi/2pi/22000）;
theta0=0/180*pi;%来波方向 （我觉得应该是天线阵的指向）
%theta0=0;%来波方向
theta1=12/180*pi;
theta2=30/180*pi;
L=20;
kk_sum=0;
SNR=8;
w_l0=exp（-imag1*2*pi*d_lamda*sin（theta0）*[0:element_num-1]‘）;
w_l1=exp（-imag1*2*pi*d_lamda*sin（theta1）*[0:element_num-1]‘）;
w_l2=exp（-imag1*2*pi*d_lamda*sin（theta2）*[0:element_num-1]‘）;
A=[w_l0w_l1w_l2];
S=randn（length（S）L）
x=A*S;
X=awgn（xSNR‘measured‘）;
R=X*X‘/L;


[VD]=eig（R）;
Ci=inv（R）; 

for  j=1:length（theta）  %（我认为是入射角度，即来波方向计算阵列流形矩阵A）
     a=exp（-imag1*2*pi*d_lamda*sin（theta（j））*[0:element_num-1].‘）;
     
     p（j）=1/（a‘*Ci*a）;  %（matlab中的‘默认为共轭转置，如果要计算转置为w.‘*a）
end

figure（1）
plot（theta*180/pidb（abs（p）/max（abs（p））））;
xlabel（‘角度/°‘）;
ylabel（‘归一化增益/dB‘）;
title（‘0°，12°，30°时capon算法天线方向图‘）
figure（2）                             %绘出应用MUSIC法估计的波达方向图            
polar（thetaabs（p））;
%axis（[-90 90 -60 0]）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1152  2017-12-15 20:43  capon.m
     文件        1703  2017-12-15 21:19  capon_L.m
     文件        1420  2017-12-15 22:21  capon_MUSIC.m
     文件        1634  2017-12-15 21:05  capon_SNR.m