二维DOA估计matlab代码

二维DOA估计是指通过接收到的信号来估计信号源在二维平面上的方向。下面是一个简单的基于协方差矩阵的二维DOA估计的Matlab代码，供参考：

% 二维DOA估计的Matlab代码
clc;
clear all;
close all;

% 参数设置
M = 8;                  % 阵元数
d = 0.5;                % 阵元间距
fs = 1000;              % 采样频率
f = [100 200];          % 信号频率
theta = [-30 20];       % 信号入射角度
SNR = 10;               % 信噪比

% 构造阵列流型矩阵
A = zeros(M, length(theta));
for i = 1:length(theta)
    for j = 1:M
        A(j,i) = exp(-1j*2*pi*d*(j-1)*sind(theta(i)));
    end
end

% 生成噪声和信号
N = 1024;
t = (0:N-1)/fs;
noise = randn(M,N);
signal = 0;
for i = 1:length(theta)
    signal = signal + exp(1j*2*pi*f(i)*t).*A(:,i);
end
signal = signal./norm(signal);
x = signal + 10^(-SNR/20)*norm(signal).*noise;

% 协方差矩阵
Rxx = x*x'/N;

% MUSIC算法
theta_range = -90:0.1:90;
Pmusic = zeros(length(theta_range),1);
for i = 1:length(theta_range)
    A_range = zeros(M, 1);
    for j = 1:M
        A_range(j) = exp(-1j*2*pi*d*(j-1)*sind(theta_range(i)));
    end
    Pmusic(i) = abs(1/((A_range')*inv(Rxx)*A_range));
end

% 画图
figure;
subplot(2,1,1);
plot(t,real(x));
xlabel('时间 t (秒)');
ylabel('幅值');
title('接收信号');

subplot(2,1,2);
plot(theta_range,10*log10(Pmusic/max(Pmusic)));
xlabel('入射角度 (度)');
ylabel('功率谱 (dB)');
title('MUSIC算法估计的入射角度');

该代码主要包括以下几个步骤：

1. 构造阵列流型矩阵；
2. 生成噪声和信号；
3. 计算协方差矩阵；
4. 使用MUSIC算法估计信号入射角度；
5. 画图展示结果。

需要注意的是，该代码仅适用于二维平面上的单频信号，如果需要估计多频信号或三维信号，则需要进行相应的修改。