matlab 计算偶极子天线半功率波瓣宽度的循环函数代码

以下是 MATLAB 计算偶极子天线半功率波瓣宽度的循环函数代码：

function [HPBW,theta] = dipoleHPBW(l,lambda)

% 计算偶极子天线的半功率波瓣宽度(HPBW)
% l：天线长度，lambda：工作波长
% HPBW：半功率波瓣宽度，theta：方向角度数

N = 181; % 将180度分成181份
theta = linspace(0,pi,N); % 生成方向角度数向量
k = 2*pi/lambda; % 波数

% 计算方向角度数对应的辐射强度
U = zeros(1,N);
for i = 1:N
    U(i) = abs((cos(theta(i)*pi))*((sin(l*k*cos(theta(i)))/((l*k*cos(theta(i)))^2))));
end

% 计算最大辐射强度
Umax = max(U);

% 找到半功率波瓣宽度所对应的方向角度数
for i = 1:N
    if U(i) >= Umax/2
        index1 = i;
        break
    end
end
for i = N:-1:index1
    if U(i) <= Umax/2
        index2 = i;
        break
    end
end

% 计算半功率波瓣宽度
HPBW = (index2-index1)/(N-1)*180;
end

该函数输入天线长度 `l` 和工作波长 `lambda`，输出半功率波瓣宽度 `HPBW` 和方向角度数 `theta`。函数使用循环计算偶极子天线在不同方向上的辐射强度，并找到半功率波瓣宽度所对应的方向角度数。最后，将方向角度数转换为度数，并计算出半功率波瓣宽度。