构建一个八个阵元的均匀线阵，三个远场信号源，入射角度分别为-30 30 60度，构建信号模型，并加入方差为0.5，均值为0的随机分布的幅度误差，方差50，均值为0随机分布的相位误差，并使用单辅助信号源的有源校正法对幅相误差校正，并输出校正前后的对比图和数值比较，使用matlab实现的代码

以下是使用MATLAB实现的代码：

% 定义信号源的入射角度和波长
theta = [-30 30 60]; % 入射角度
lambda = 1; % 波长

% 定义均匀线阵参数
N = 8; % 阵元个数
d = lambda/2; % 阵元间距

% 构建均匀线阵的位置向量
pos = zeros(2,N);
for n = 1:N
    pos(1,n) = (n-1)*d;
end

% 构建信号模型
A = ones(N,length(theta)); % 幅度矩阵
for k = 1:length(theta)
    phase = (2*pi/lambda)*d*sin(theta(k)*(pi/180))*(0:N-1)';
    A(:,k) = exp(1i*phase); % 构建相位矩阵
end

% 加入幅度误差和相位误差
A_error = sqrt(0.5)*randn(N,length(theta)) + 1i*sqrt(0.5)*randn(N,length(theta)); % 幅度误差
p_error = sqrt(50)*randn(N,length(theta)); % 相位误差
A = A.*A_error;
A = A.*exp(1i*p_error);

% 构建接收信号矩阵
R = zeros(N,length(theta));
for k = 1:length(theta)
    R(:,k) = A(:,k).*exp(-1i*(2*pi/lambda)*d*sin(theta(k)*(pi/180))*(0:N-1)');
end

% 有源校正法进行幅相误差校正
s_ref = exp(1i*(2*pi/lambda)*d*sin(theta(1)*(pi/180))*(0:N-1)'); % 辅助信号源
s_ref = s_ref/norm(s_ref); % 归一化
W = zeros(N,length(theta));
for k = 1:length(theta)
    W(:,k) = s_ref.*conj(A(:,k));
end

% 构建校正后的接收信号矩阵
R_corrected = R - s_ref*sum(W,2)';

% 绘制校正前后的对比