基于被动时反的 ac-MVDR 水下声源定位matlab代码

以下是基于被动时反的 ac-MVDR 水下声源定位的 MATLAB 代码：

% 声明变量
N = 8; % 阵元个数
fs = 8000; % 采样率
c = 1500; % 水中声速
d = 0.05; % 阵元间距
theta = -90:1:90; % 目标角度范围
theta = theta .* pi / 180; % 转换为弧度制
freq = 2000; % 信号频率
lambda = c / freq; % 信号波长

% 生成阵列流型矩阵A
A = zeros(N, length(theta));
for n = 1:N
    A(n,:) = exp(-1i * 2 * pi * (n-1) * d / lambda * cos(theta));
end

% 加载数据
load('data.mat'); % data为包含8个信道数据的矩阵

% 计算协方差矩阵
X = data.';
Rxx = X * X' / size(X, 2);

% ac-MVDR算法
theta_est = zeros(1, length(theta));
for i = 1:length(theta)
    w = (Rxx + 0.01 * eye(size(Rxx))) \ A(:,i) / (A(:,i)' * (Rxx + 0.01 * eye(size(Rxx))) \ A(:,i));
    theta_est(i) = atan2(imag(w'*X), real(w'*X));
end

% 结果可视化
figure;
plot(theta*180/pi, abs(theta_est)*180/pi, '*');
xlabel('True angle (degree)');
ylabel('Estimated angle (degree)');
title('ac-MVDR underwater source localization');

在代码中，首先定义了阵元个数、采样率、水中声速、阵元间距、目标角度范围、信号频率、信号波长等变量。然后生成阵列流型矩阵A，加载数据并计算协方差矩阵Rxx。接着使用ac-MVDR算法估计目标角度，并将结果可视化。