给定两个方向角利用波束形成法进行波达方向角估计matlab代码

以下是给定两个方向角利用波束形成法进行波达方向角估计的 Matlab 代码，供您参考：

%% 参数设置
fc = 2.4e9;  % 信号频率
c = physconst('LightSpeed');  % 光速
lambda = c/fc;  % 波长
d = lambda/2;  % 天线间距
theta = [-30, 30];  % 给定两个方向角

%% 波束形成法
N = 4;  % 天线数
M = 181;  % 方向角分辨率
w = ones(N,1);  % 波束权重
A = zeros(N,M);  % 天线阵列流形矩阵
for k = 1:M
    A(:,k) = exp(-1j*2*pi*d*sin((k-1)*pi/180)*(0:N-1).'/lambda);
end
P = zeros(M,M);  % 信号的协方差矩阵
for k = 1:length(theta)
    a = exp(-1j*2*pi*d*sin(theta(k)*pi/180)*(0:N-1).'/lambda);
    P = P + a*a';
end
P = P/length(theta);  % 平均信号的协方差矩阵
P_noise = eye(N);  % 噪声的协方差矩阵
R = P + P_noise;  % 接收信号的协方差矩阵
w_opt = inv(R)*A*inv(A'*inv(R)*A)*w;  % 最优波束权重
theta_est = asin(-1j*w_opt(2)/(2*pi*d))*180/pi;  % 估计的波达方向角
disp(['估计的波达方向角为：', num2str(theta_est)]);

需要说明的是，以上代码仅是一个简单的示例，实际应用中需要根据具体场景进行参数的设置和算法的优化。另外，为了减小噪声的影响，实际应用中通常需要对信号进行平均处理。