单元激励不相同时平面阵方向图函数公式matlab程序实现 不调用array_pattern

以下是使用MATLAB实现平面阵方向图函数公式的程序：

% 定义常量
c = 3e8; % 光速

% 定义阵列参数
N = 8; % 阵元数量
d = 0.5 * c / fc; % 阵元间距

% 定义角度范围
theta = -90:0.1:90; % 角度范围

% 计算方向图函数
AF = zeros(size(theta)); % 初始化方向图函数

for n = 1:N
    % 计算阵元位置
    xn = (n-1) * d;
    
    % 计算阵元相位
    phase = exp(-1j * 2 * pi * xn * sind(theta) / lambda);
    
    % 累加阵元激励
    AF = AF + phase;
end

% 绘制方向图
plot(theta, abs(AF));
xlabel('角度 (度)');
ylabel('方向图函数');
title('平面阵方向图');

请注意，此程序假设已经定义了变量fc（载波频率）和lambda（波长）。您需要在程序中根据您的实际情况进行相应的更改。