请用matlab实现直线阵列天线加上高斯白噪声的远场方向图

以下是一个简单的 Matlab 代码示例，用于生成具有高斯白噪声的直线阵列天线的远场方向图。假设阵列中有 $N$ 个天线，每个天线的位置为 $(x_i, y_i, z_i)$，其中 $i=1,2,\dots,N$。阵列的方向为 $\theta$ 和 $\phi$。噪声的方差为 $\sigma^2$。

% 阵列参数
N = 8; % 天线数目
d = 0.5; % 天线间距（米）
lambda = 0.1; % 波长（米）
k = 2*pi/lambda; % 波数

% 生成阵列中每个天线的位置坐标
x = 0:d:(N-1)*d;
y = zeros(1,N);
z = zeros(1,N);

% 生成远场方向图的角度范围
theta_range = linspace(0, pi, 181); % 0 到 180 度
phi_range = linspace(0, 2*pi, 361); % 0 到 360 度

% 生成远场方向图
F = zeros(length(theta_range), length(phi_range));
for t = 1:length(theta_range)
    for p = 1:length(phi_range)
        theta = theta_range(t);
        phi = phi_range(p);
        d_cos = cos(theta)*cos(phi)*x + cos(theta)*sin(phi)*y + sin(theta)*z;
        F(t,p) = abs(sum(exp(1j*k*d_cos)))^2/N;
    end
end

% 添加高斯白噪声
sigma = 0.1;
F_noisy = F + sigma*randn(size(F));

% 绘制远场方向图
figure;
imagesc(phi_range/pi*180, theta_range/pi*180, 10*log10(F_noisy));
xlabel('方位角（度）');
ylabel('俯仰角（度）');
title(['高斯白噪声远场方向图，\sigma = ', num2str(sigma)]);
colorbar;

这段代码将生成一个带有高斯白噪声的远场方向图。你可以根据需要调整阵列参数和噪声方差。