写一段matlab程序，模拟仿真4个阵元的阵列天线的随机噪声

clc;
clear all;

% 设置仿真参数
n = 10000; % 仿真次数
fs = 1e6; % 采样频率
fc = 2.4e9; % 中心频率
bw = 20e6; % 带宽
t = 0:1/fs:(1-1/fs); % 时间序列
N = length(t); % 采样点数
dt = 1/fs; % 采样时间间隔

% 生成随机噪声
noise = randn(4,N);

% 生成信号
signal = sin(2*pi*fc*t);

% 配置阵列天线
d = 0.5; % 阵元间距
theta = [0 30 60 90]; % 信号入射角度
lambda = 3e8/fc; % 波长
k = 2*pi/lambda; % 波数
phase = k*d*sind(theta); % 阵元间相位差

% 计算阵列输出
array_out = zeros(1,N);
for i = 1:4
array_out = array_out + exp(1j*phase(i))*signal + noise(i,:);
end

% 绘制结果
figure;
subplot(2,1,1);
plot(t,signal);
xlabel('Time (s)');
ylabel('Signal');
title('Input Signal');

subplot(2,1,2);
plot(t,abs(array_out));
xlabel('Time (s)');
ylabel('Amplitude');
title('Array Output');