利用matlab写一段长短基线干涉仪测向的仿真程序。 要求： 1.共3个阵元，其中1号和2号阵元构成短基线，距离为0.125，1号和3号阵元构成长基线距离为0.625； 2.需要添加随机噪声，其中信噪比为10dB； 3.信源频率为1.2*10^9;

以下是一份简单的长短基线干涉仪测向的MATLAB仿真程序：

% 清空工作区和命令窗口
clear all;
clc;

% 定义参数
lambda = 3e8 / 1.2e9; % 波长
short_baseline = 0.125; % 短基线长度
long_baseline = 0.625; % 长基线长度
noise_snr = 10; % 噪声信噪比

% 生成信号
source_signal = exp(1j * 2 * pi * rand); % 信源信号
noise_signal = sqrt(0.5 / (10^(noise_snr / 10))) * (randn(3, 1) + 1j * randn(3, 1)); % 噪声信号
signal = source_signal * ones(3, 1) + noise_signal; % 总信号

% 计算阵元位置
array_pos = [0, 0, 0; short_baseline, 0, 0; 0, long_baseline, 0];

% 计算短基线干涉结果
short_baseline_phase = 2 * pi / lambda * (array_pos(2, :) - array_pos(1, :)) * [1; 0; 0]; % 短基线干涉相位
short_baseline_interference = exp(1j * short_baseline_phase); % 短基线干涉结果

% 计算长基线干涉结果
long_baseline_phase = 2 * pi / lambda * (array_pos(3, :) - array_pos(1, :)) * [1; 0; 0]; % 长基线干涉相位
long_baseline_interference = exp(1j * long_baseline_phase); % 长基线干涉结果

% 计算测向结果
bearing = atan2(imag(long_baseline_interference), imag(short_baseline_interference)) * 180 / pi; % 简单测向公式

% 输出结果
fprintf('信源方向：%f度\n', angle(source_signal) * 180 / pi);
fprintf('阵元1干涉相位：%f度\n', angle(short_baseline_interference(1)) * 180 / pi);
fprintf('阵元2干涉相位：%f度\n', angle(short_baseline_interference(2)) * 180 / pi);
fprintf('阵元3干涉相位：%f度\n', angle(long_baseline_interference(3)) * 180 / pi);
fprintf('测向结果：%f度\n', bearing);