MATLAB设计均匀线阵下MMSE波束形成器：N=16，θ0=-30°，d=λ/2， 干扰为 θ1=60°，θ2=19°，θ3=45°，干扰功率分别为40dB、 30dB和35dB

以下是MATLAB代码实现：

% 参数设置
N = 16; % 阵元数
d = 0.5; % 阵元间距，单位为波长
theta0 = -30; % 信源角度，单位为度
theta = (-90:0.1:90)'; % 探测角度范围，单位为度
theta1 = 60; % 干扰1角度，单位为度
theta2 = 19; % 干扰2角度，单位为度
theta3 = 45; % 干扰3角度，单位为度
P1 = 10^(-40/10); % 干扰1功率
P2 = 10^(-30/10); % 干扰2功率
P3 = 10^(-35/10); % 干扰3功率

% 构造阵列流形矩阵
A = zeros(N,length(theta));
for n = 1:N
A(n,:) = exp(-1j*2*pi*d*(n-1)*sind(theta))/sqrt(N);
end

% 构造信源方向矩阵
S = exp(-1j*2*pi*d*(0:N-1)'*sind(theta0))/sqrt(N);

% 构造干扰方向矩阵
I1 = exp(-1j*2*pi*d*(0:N-1)'*sind(theta1))/sqrt(N);
I2 = exp(-1j*2*pi*d*(0:N-1)'*sind(theta2))/sqrt(N);
I3 = exp(-1j*2*pi*d*(0:N-1)'*sind(theta3))/sqrt(N);

% 构造协方差矩阵
R = P1*I1*I1'+P2*I2*I2'+P3*I3*I3'+S*S';

% 计算权向量
w = inv(R)*S;

% 计算波束响应
B = abs(w'*A).^2;

% 绘制波束图
plot(theta,B);
xlabel('探测角度（度）');
ylabel('波束响应');
title('均匀线阵下MMSE波束形成器波束图');