MATLAB实现平面阵方向图函数公式 - 单元激励不相同时

在MATLAB中，可以使用以下代码实现平面阵方向图函数公式的计算，适用于单元激励不相等的情况：

% 定义单元激励不相同时的平面阵方向图函数公式
function pattern = array_pattern(theta, phi, weights)
    % theta: 仰角
    % phi: 方位角
    % weights: 单元激励权重
    
    % 计算方向图函数公式
    pattern = zeros(size(theta));
    for i = 1:length(weights)
        pattern = pattern + weights(i) * exp(1i * 2 * pi * (i-1) * cos(theta) * sin(phi));
    end
    
    % 归一化
    pattern = abs(pattern) / max(abs(pattern));
end

使用示例：

% 定义仰角和方位角的范围
theta = linspace(0, pi, 100);
phi = linspace(0, 2*pi, 100);

% 定义单元激励权重
weights = [1, 0.8, 0.6, 0.4, 0.2];

% 计算方向图函数
pattern = array_pattern(theta, phi, weights);

% 绘制方向图
figure;
polarplot3d(theta, phi, pattern);

这段代码定义了一个名为'array_pattern'的函数，用于计算单元激励不相同时的平面阵的方向图函数。函数的输入参数包括仰角'theta'、方位角'phi'和单元激励权重'weights'。函数通过对每个单元的激励权重进行加权求和来计算方向图函数。最后，使用'polarplot3d'函数绘制方向图。