用matlab画出惠更斯元极坐标归一化方向图，并求出其半功率宽度HP

在MATLAB中，你可以使用`radar`函数和`anglevec`函数来绘制惠更斯元极坐标归一化的方向图。首先，你需要设置信号的频率、辐射角度和相关的参数。以下是简单的步骤：

1. 定义信号参数：

[theta, r] = meshgrid(anglevec(0:180)); % 创建角度向量
freq = 1; % 频率（通常用赫兹表示）

2. 使用惠更斯原理计算每个点的幅度：

H = exp(-1i * 2*pi * freq * r.*sind(theta)); % 惠更斯元计算

3. 归一化并绘制方向图：

H_normalized = H / max(abs(H(:))); % 归一化到最大值为1
imagesc(theta, r, abs(H_normalized)); % 绘制幅值图像
colorbar; % 显示颜色条
title('惠更斯元极坐标归一化方向图');
xlabel('方位角 (度)');
ylabel('距离 (单位)');

4. 计算半功率宽度（Half-Power Beamwidth, HPBW）：

% 一般通过找到幅值下降到峰值一半的两个临界角度来计算HPBW
half_power_threshold = 0.707; % 选取70.7%作为阈值
[H_half_power, idx] = min(abs(H_normalized), [], 2); % 找到最小的幅值点
HPBW_degrees = [theta(idx(1)) theta(idx(2))] - 180; % 考虑对称性减去180度