基于相关干涉仪测向方法的MATLAB仿真实现

一、核心代码实现

1. 参数配置

%% 系统参数
c = 3e8;        % 光速(m/s)
f0 = 2e9;       % 信号频率(Hz)
lambda = c/f0;  % 波长(m)
fs = 10*f0;     % 采样频率
T = 1e-6;       % 信号持续时间(s)
N = fs*T;       % 采样点数
d1 = lambda/2;  % 基线1长度
d2 = 3*lambda/2;% 基线2长度
theta_true = 35;% 真实来波方向(度)

2. 信号生成

%% 生成阵列信号
t = 0:1/fs:T-1/fs;
s0 = exp(1j*2*pi*f0*t); % 参考信号

% 阵元1信号
tau1 = d1*sin(deg2rad(theta_true))/c;
s1 = exp(1j*2*pi*f0*(t - tau1));

% 阵元2信号
tau2 = d2*sin(deg2rad(theta_true))/c;
s2 = exp(1j*2*pi*f0*