写一个Matlab 非均匀噪声情况DOA估计代码

好的，下面是一个使用 MUSIC（MUltiple SIgnal Classification）算法进行非均匀噪声情况下的 DOA（方向性传输函数）估计的 Matlab 代码示例：

% 定义参数
fs = 10000; % 采样率
t = 0:1/fs:1; % 时间向量
f1 = 1000; % 信号1频率
f2 = 2000; % 信号2频率
snr = 10; % 信噪比（以dB为单位）
antenna_spacing = 0.5; % 天线间距

% 生成信号
s1 = sin(2*pi*f1*t); % 信号1
s2 = sin(2*pi*f2*t); % 信号2

% 添加噪声
noise = randn(size(t)); % 高斯白噪声
noise_power = norm(s1)/10^(snr/10); % 噪声功率
noise = noise * noise_power; % 根据信噪比调整噪声功率
x = s1 + s2 + noise; % 加噪声的接收信号

% DOA估计
theta = -90:0.5:90; % 待估计的方向范围
N = length(x); % 信号长度
M = floor(fs/(2*f1*antenna_spacing)); % 阵列天线数
Rxx = x' * x; % 协方差矩阵估计

% MUSIC算法
[U, S, ~] = svd(Rxx); % 奇异值分解
E = U(:, M+1:end); % 信号子空间
P = zeros(size(theta)); % 方向谱

for k = 1:length(theta)
    a = exp(1i*2*pi*f1*antenna_spacing*sin(pi*theta(k)/180)*(0:M-1)'); % 信号阵列流形
    P(k) = 1/(a'*E*E'*a); % 方向谱估计
end

% 绘制估计结果
plot(theta, 10*log10(P/max(P))); % 归一化并转换为dB
xlabel('方向（度）');
ylabel('幅度（dB）');
title('非均匀噪声情况DOA估计');

这个代码示例演示了如何使用 MUSIC 算法对非均匀噪声情况下的两个信号进行方向估计。您可以根据自己的需要修改参数和信号模型。请确保在使用之前充分理解代码，并根据您的需求进行适当的修改。希望对您有所帮助！