二维music算法matlab

二维MUSIC算法（Multiple Signal Classification）是一种用于估计二维平面上信号源位置的算法，MATLAB中也提供了相应的函数实现。以下是一个简单的二维MUSIC算法的MATLAB实现：

% 二维MUSIC算法实现
% Rx: 接收信号的协方差矩阵
% K: 需要估计的信号源个数
% d: 阵元间距
% theta_range: 估计角度的范围
function [theta, Pmusic] = music_2d(Rx, K, d, theta_range)
    M = size(Rx, 1);  % 阵列元素个数
    Ntheta = length(theta_range);  % 角度个数
    [U, S, ~] = svd(Rx);  % 对协方差矩阵进行奇异值分解
    En = U(:, K+1:M);  % 估计噪声子空间
    Pmusic = zeros(Ntheta, 1);  % 初始化MUSIC谱
    for ii = 1:Ntheta
        a = exp(-1j*2*pi*d*sin(theta_range(ii))*(0:M-1)');  % 构造阵列流形
        Pmusic(ii) = 1/(a'*En*En'*a);  % 计算MUSIC谱
    end
    Pmusic = 10*log10(abs(Pmusic)/max(abs(Pmusic)));  % 对MUSIC谱进行归一化
    [~, idx] = sort(Pmusic, 'descend');  % 对MUSIC谱进行排序
    theta = theta_range(idx(1:K));  % 提取前K个最大的MUSIC峰值对应的角度
end

其中，输入参数Rx为接收信号的协方差矩阵，K为需要估计的信号源个数，d为阵元间距，theta_range为估计角度的范围。输出参数theta为估计得到的信号源位置，Pmusic为对应的MUSIC谱。