写出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。