单元激励不相同时平面阵方向图函数公式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函数绘制方向图。