求进行10次蒙特卡洛实验后的二维的music算法的方位角和俯仰角的角度均值，用matlab实现完整代码，检查无误后再发给我

这里给出使用MATLAB实现10次蒙特卡洛实验的代码，并求出二维music算法的方位角和俯仰角的角度均值。

首先，生成随机信号源和阵列位置：

% 生成随机信号源和阵列位置
M = 6;  % 阵元数
N = 3;  % 信号源数
d = 0.5;  % 阵元间距
lambda = 1;  % 波长
fc = lambda/2;  % 信号频率
theta_s = sort(rand(1, N)*180-90);  % 信号源方位角
phi_s = sort(rand(1, N)*180-90);  % 信号源俯仰角
theta_a = rand(1, M)*180-90;  % 阵列方位角
phi_a = rand(1, M)*180-90;  % 阵列俯仰角

然后，生成接收信号矩阵：

% 生成接收信号矩阵
A = zeros(M, N);
for i = 1:N
    A(:, i) = exp(-1j*2*pi*d/lambda*[0:M-1]'*sind(theta_a)*cosd(phi_a(i)));
end
S = exp(-1j*2*pi*fc*[0:N-1]'*cosd(phi_s));
X = A*S;

接下来，进行蒙特卡洛实验并求解方位角和俯仰角均值：

% 进行蒙特卡洛实验并求解均值
num_exp = 10;  % 实验次数
theta_hat = zeros(1, num_exp);
phi_hat = zeros(1, num_exp);
for i = 1:num_exp
    % 加入高斯噪声
    Xn = X + randn(M, N)*0.1;
    % 二维MUSIC算法
    Rxx = Xn*Xn'/N;
    [V, D] = eig(Rxx);
    [D, idx] = sort(diag(D), 'descend');
    V = V(:, idx);
    En = V(:, N+1:M);
    theta_range = -90:0.1:90;
    phi_range = -90:0.1:90;
    Pmusic = zeros(length(theta_range), length(phi_range));
    for j = 1:length(theta_range)
        for k = 1:length(phi_range)
            a = exp(-1j*2*pi*d/lambda*[0:M-1]'*sind(theta_a)*cosd(phi_a(k)));
            Pmusic(j, k) = 1/norm(En'*a)^2;
        end
    end
    [max_P, max_idx] = max(Pmusic(:));
    [theta_idx, phi_idx] = ind2sub(size(Pmusic), max_idx);
    theta_hat(i) = theta_range(theta_idx);
    phi_hat(i) = phi_range(phi_idx);
end
mean_theta_hat = mean(theta_hat)
mean_phi_hat = mean(phi_hat)

其中，`num_exp`为实验次数，可以根据需要进行修改。在每次实验中，我们加入高斯噪声，然后使用二维MUSIC算法求解方位角和俯仰角的估计值。最后，取所有实验的估计值的均值作为最终的估计结果。

完整代码如下：