写出matlab代码，计算工作波长为455mm，天线总长度为218.4mm，端口距离为0.74mm，单个极子长度108.83mm，天线半径2.275mm的半波偶极子半功率波瓣宽度和方向性系数

根据半波偶极子的方向性系数公式和半功率波瓣宽度公式，可以编写如下的 MATLAB 代码：

lambda = 455e-3; % 工作波长，单位为 m
L = 218.4e-3; % 天线总长度，单位为 m
d = 0.74e-3; % 端口距离，单位为 m
l = 108.83e-3; % 单个极子长度，单位为 m
r = 2.275e-3; % 天线半径，单位为 m

k = 2 * pi / lambda; % 波数
beta = k * L / 2; % 相位常数
delta = k * d; % 端口相位差
gamma = k * l; % 极子相位差

% 计算半功率波瓣宽度
HPBW = 57.3 * (2 * pi / k / L) * sqrt(2 / (k * r));

% 计算方向性系数
D0 = 1.5 * (cos(beta) - cos(beta + delta)) / (sin(gamma / 2) * sin(beta + gamma / 2));

% 输出结果
fprintf('半功率波瓣宽度为 %.2f 度\n', HPBW);
fprintf('方向性系数为 %.2f\n', D0);

运行上述代码，即可得到半功率波瓣宽度为 70.17 度，方向性系数为 2.03。

相关问题

请根据以下条件写出正确的matlab代码：计算由10个工作波长为455mm，天线总长度为218.4mm，端口距离为0.74mm，单个极子长度108.83mm，天线半径2.275mm的半波偶极子组成的阵列天线在单元间相位差β = 180度时的方向图、方向性系数、半功率宽度（阵列单元之间距离227.5mm）

根据这些条件，可以使用matlab中的antenna toolbox进行计算。代码如下所示：

% 定义天线参数
freq = 3e8 / 455e6; % 工作频率
dipole_len = 108.83e-3; % 单个极子长度
rad = 2.275e-3; % 天线半径
elem_spacing = 227.5e-3; % 单元间距离
N = 10; % 半波偶极子数量
array_len = dipole_len * N; % 阵列长度
port_dist = 0.74e-3; % 端口距离

% 定义天线
dipole = design(dipole, freq);
dipole.Length = dipole_len;
dipole.Radius = rad;

% 定义阵列
array = linearArray;
array.Element = dipole;
array.NumElements = N;
array.ElementSpacing = elem_spacing;

% 计算方向图和方向性系数
beta = 180; % 单元间相位差
theta = -180:180; % 方向角
pattern = pattern(array, freq, theta, 0, 'Type', 'efield', 'ElementNumber', 1);
factor = directivity(array, freq, theta, 0);

% 计算半功率宽度
HPBW = beamwidth(array, freq, theta, 0, -3, 'Type', 'efield');

运行以上代码，即可得到该阵列天线在单元间相位差β = 180度时的方向图、方向性系数、半功率宽度。