写一下二维面阵天线方向图的matlab程序，要求能画出天线方向图，且能够求出3db波束宽度所对应的方位向和俯仰向的角度

二维面阵天线的方向图通常用于描述天线辐射功率在空间各个方向上的分布情况。在MATLAB中，我们可以利用`phased-array System Toolbox`来创建和绘制这个图形。下面是一个简单的示例，展示如何计算并绘制三维方向图，并找到3dB波束宽度对应的方位角和俯仰角：

% 导入必要的库
import phased.*
clear all; close all;

% 设定天线参数
numAntennas = 4;       % 数量为4的偶数阵列天线
lambda = 0.15;         % 波长 (假设厘米)
elementSpacing = lambda/2;   % 元件之间的间距

% 创建平面阵列天线系统
array = phased.ULA('NumElements', numAntennas, 'ElementSpacing', elementSpacing);

% 设置观测角度范围 (水平和垂直)
azimuths = linspace(-180, 180, 360);  % 方位角
elevations = linspace(-90, 90, 180);    % 俯仰角

% 计算方向图
[response, azEl] = pattern(array, azimuths, elevations, 'Type', 'power');

% 绘制方向图
surf(azimuths, elevations, 20*log10(response))
xlabel('Azimuth (degrees)');
ylabel('Elevation (degrees)');
zlabel('Power (dB)');
title('2D Array Antenna Directional Pattern');

% 找到3dB波束宽度
[~, idx3dB] = max(response);
azimuth3dB = azEl(1,idx3dB);
elevation3dB = azEl(2,idx3dB);

% 输出结果
fprintf('3dB beam width at:\n');
fprintf('Azimuth: %.2f degrees\n', azimuth3dB);
fprintf('Elevation: %.2f degrees\n', elevation3dB);