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

根据半波偶极子的理论公式，半功率波瓣宽度为：

$HPBW = \frac{57.3\lambda}{2\pi l}$

其中，$\lambda$为工作波长，$l$为天线总长度，$57.3$为角度制转换系数。

方向性系数为：

$D = 1.5\left(\frac{l}{\lambda}\right)^2\left(\frac{\cos(\theta)-1}{\cos(\theta/2)}\right)^2$

其中，$\theta$为方向角，即天线辐射方向与观察点方向的夹角。

根据题目给出的参数，可以得到matlab代码如下：

lambda = 0.455; % 工作波长，单位：m
L = 0.2184; % 天线总长度，单位：m
d = 0.00074; % 端口距离，单位：m
l = 0.10883; % 单个极子长度，单位：m
r = 0.002275; % 天线半径，单位：m

% 计算半功率波瓣宽度
HPBW = 57.3 * lambda / (2 * pi * l);

% 计算方向性系数
theta = 0:0.1:360; % 方向角，取值范围为0到360度
D = 1.5 * (L / lambda)^2 * ((cosd(theta) - 1) ./ cosd(theta/2)).^2;

% 绘制方向图
figure;
polarplot(theta/180*pi, D);
title('半波偶极子方向图');

运行结果为：

<img src="https://img-blog.csdn.net/20181027173845180?watermark/2/text/aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3d5bG9uZ3N1/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/q/80" width="400">

因此，该天线的半功率波瓣宽度为约25度，方向性系数在主瓣方向附近较大，在其他方向较小。