matlab实现平面阵列波束成形，请给出具体代码

以下是MATLAB实现平面阵列波束成形的代码示例：

clear all;
clc;
close all;

% 设置阵列参数
N = 16; % 阵元数
d = 0.5; % 阵元间距，单位：波长
theta0 = 30; % 信源入射角，单位：度
theta0 = theta0*pi/180; % 转化为弧度制
lambda = 1; % 波长

% 生成阵列导向矢量
theta = 0:0.1:90; % 采样角度范围
theta = theta*pi/180; % 转化为弧度制
a = exp(-1j*pi*d/lambda*(0:N-1)'*sin(theta)); % 阵列导向矢量

% 生成接收信号
SNR = 10; % 信噪比
sigma2 = 10^(-SNR/10); % 噪声方差
s = exp(1j*2*pi/lambda*d*(0:N-1)'*sin(theta0)); % 信源信号
n = sqrt(sigma2/2)*(randn(N,length(theta))+1j*randn(N,length(theta))); % 噪声信号
x = a*s+n; % 接收信号

% 理论波束形成
w0 = exp(-1j*pi*d/lambda*(0:N-1)'*sin(theta0)); % 理论加权系数
P_theory = abs(w0'*a).^2; % 理论波束形成

% 最大信噪比波束形成
Rxx = x*x'/length(theta); % 采样自相关矩阵
w_MVDR = Rxx\a/(a'/Rxx*a); % MVDR加权系数
P_MVDR = abs(w_MVDR'*a).^2; % MVDR波束形成

% 平均信噪比波束形成
w_MWF = Rxx\a; % MWF加权系数
P_MWF = abs(w_MWF'*a).^2; % MWF波束形成

% 绘制波束形成图像
figure;
plot(theta*180/pi,P_theory,'r',theta*180/pi,P_MVDR,'g',theta*180/pi,P_MWF,'b');
xlabel('入射角度（度）');
ylabel('输出功率（dB）');
title('波束形成');
legend('理论波束形成','MVDR波束形成','MWF波束形成');

代码中，首先设置了阵列参数：阵元数、阵元间距、信源入射角、波长。然后，通过阵列导向矢量的计算得到不同入射角下的导向矢量。接着，生成接收信号，其中包括信源信号和噪声信号。在理论波束形成中，通过计算加权系数得到理论波束形成；在最大信噪比波束形成和平均信噪比波束形成中，分别采用MVDR和MWF算法计算加权系数，并计算波束形成。最后，绘制三种波束形成的图像。