CAPON算法估计DOA

%%CAPON算法
clc
clear
close all

%% 基本参数设置
M = 10;
N = 100;
theta = [10 20 30];
K = length(theta);
d = 0.5;
dd = 0:M-1;
SNR = 20;

%% 构造信号接收模型
s = randn(K,N);
a = exp(1i*2*pi*d*dd'*sind(theta));
x = a*s;
x = awgn(x,SNR,'measured');%  measured  表示使用实际测量的信噪比
R = x*x'/N;

%% 空间谱搜索
A = zeros(M,1);
for i = 1:361
    angle(i) = (i-181)/2;%俯仰角范围[-90,90]
    A = exp(1i*pi*dd'*sind(angle(i)));
    P(i) = 1/(A'*inv(R)*A);
end
P = abs(P);
P_max = max(P);
P = 10*log10(P/P_max);

%% 绘图
figure
h_1 = plot(angle,P,'b','LineWidth',2);
axis([-90 90 -60 0])
set(gca, 'XTick',-90:30:90)
grid on
xlabel('angle(degree)')
ylabel('magnitude spectr(dB)')
legend('CAPON')