1. 系统架构与核心模块
ADS-B接收机的MATLAB实现需包含以下核心模块：
- 卫星信号源初始化：模拟或接收真实卫星ADS-B信号。
- 接收天线参数配置：设置天线类型、频率响应及方向性。
- 射频前端配置：AD9361芯片参数设置（LO频率、增益、滤波器）。
- 信号处理链路：下变频、滤波、解调与解码。
2. 卫星信号源初始化
2.1 信号生成参数设置
%% 卫星信号参数
fc = 1090e6; % ADS-B中心频率 (Hz)
fs = 2.5e6; % 采样率 (Hz)
bw = 10e6; % 信号带宽 (Hz)
modType = 'PPM'; % 脉冲位置调制
prf = 1e6; % 脉冲重复频率 (Hz)
2.2 信号生成代码
%% 生成ADS-B信号（示例）
t = 0:1/fs:1e-3; % 1ms信号时长
carrier = exp(1j*2*pi*fc*t); % 载波信号
modSignal = pammod(randi([0 1],1,fs/2), 1); % PPM调制
txSignal = carrier .* modSignal; % 调制后信号
2.3 卫星轨道参数（可选）
%% 卫星轨道模型（TLE数据）
tle = read_tle('satellite.tle'); % 读取TLE文件
satellite = satelliteScenario(tle); % 构建卫星场景
[az, el, r] = satellitePosition(satellite, datetime('now')); % 实时位置计算
3. 接收天线初始化
3.1 天线类型选择
- 双半波偶极子天线：适用于1090MHz频段，需设置方向图和增益。
- 螺旋天线：宽频带特性，需调整螺距和圈数。
3.2 天线参数配置
%% 天线参数设置
antenna = phased.IsotropicAntennaElement('FrequencyRange',[70e6 6e9]); % 全向天线
antenna.Gain = 5; % 天线增益 (dBi)
antenna.Polarization = 'RHCP'; % 右旋圆极化
**3.3 天线阵列配置（可选）
%% 天线阵列（提升方向性）
array = phased.ULA('NumElements',4,'ElementSpacing',0.5*lambda); % 4元均匀线阵
array.Element = antenna; % 设置单元天线
4. 射频前端配置（AD9361）
4.1 IIO系统对象初始化
%% 创建IIO系统对象
s = iio_sys_obj_matlab;
s.ip_address = '192.168.10.2'; % SDR板IP地址
s.dev_name = 'ad9361'; % 设备名称
s.in_ch_no = 1; % 输入通道数
s.out_ch_no = 1; % 输出通道数
4.2 AD9361关键参数配置
%% 设置AD9361参数
input_content = struct();
input_content.RX_LO_FREQ = 1090e6; % 接收本振频率
input_content.RX_SAMPLING_FREQ = 2.5e6; % 采样率
input_content.RX_RF_BANDWIDTH = 10e6; % 射频带宽
input_content.RX1_GAIN_MODE = 'fast_attack'; % 快速攻击增益模式
input_content.RX1_GAIN = 20; % 增益值 (dB)% 写入配置
for i = 1:numel(fieldnames(input_content))s.setAttr(input_content.(fieldnames(input_content)){i});
end
4.3 滤波器配置
%% FIR滤波器设计（抑制带外噪声）
fdesign = fdesign.lowpass('Fp,Fst,Ap,Ast',1.6e6,2e6,-1,-40,fs);
d = design(fdesign,'equiripple');
coeff = d.Coefficients; % 获取滤波器系数
s.writeFirData('adsb_fir', coeff); % 写入AD9361 FIR滤波器
5. 信号接收与处理流程
5.1 信号捕获
%% 接收信号
[rxData, len] = s.readTxRxBuffer(); % 读取IQ数据
rxData = rxData(1:len); % 截取有效数据
5.2 下变频与滤波
%% 下变频至基带
mixer = comm.Mixer('LocalOscillatorFrequency', fc, 'PhaseOffset', 0);
baseband = mixer(rxData, complex(ones(size(rxData)), zeros(size(rxData))));%% 应用FIR滤波器
b = coeff.';
filtered = filter(b, 1, baseband);
5.3 解调与解码
%% PPM解调
threshold = max(abs(filtered))/2;
bits = (abs(filtered) > threshold) .* 2 - 1; % 二进制判决%% 帧同步与解析
[header, payload] = parseModeSFrame(bits); % 解析Mode-S帧
flightID = parseFlightID(payload); % 提取航班号
6. 可视化与结果输出
6.1 信号频谱分析
%% 绘制频谱
Pxx = pwelch(filtered,1024,512,1024,fs);
plot(10*log10(Pxx));
xlabel('Frequency (Hz)');
ylabel('Power/Frequency (dB/Hz)');
title('ADS-B Signal Spectrum');
6.2 解码信息显示
%% 显示解码结果
disp('--- Decoded ADS-B Information ---');
disp(['Flight ID: ', num2str(flightID)]);
disp(['Latitude: ', num2str(lat)]);
disp(['Longitude: ', num2str(lon)]);
disp(['Altitude: ', num2str(alt), ' meters']);
7. 卫星信号模拟扩展
7.1 多普勒频移补偿
%% 多普勒效应模拟
[~, v] = satelliteVelocity(satellite, datetime('now')); % 计算卫星速度
dopplerShift = (v * fc) / speedOfLight; % 多普勒频移
rxData = rxData .* exp(1j*2*pi*dopplerShift*t); # 补偿频移
7.2 信道建模
%% 添加大气衰减
attenuation = 10^(-0.003 * pathLoss); # 路径损耗模型
rxData = rxData * attenuation;
8. 关键问题与优化建议
-
同步问题：
-
使用早迟门同步算法补偿传播延迟。
-
代码示例：
[early, prompt, late] = earlyLateGate(baseband, 0.5e-6, 1e-6); syncIndex = findCorrelationPeak(early, prompt); # 寻找同步点
-
-
抗干扰措施：
-
采用自适应滤波（LMS算法）抑制窄带干扰。
-
代码示例：
[y, e] = adaptfilt.lms(32, 0.01); filtered = filter(y, 1, noisySignal);
-
-
实时性优化：
-
使用MATLAB Coder生成C代码加速信号处理。
-
示例命令：
codegen processADSBSignal -args {rxData};
-
9. 完整代码框架
%% 主程序
clear; clc;% 初始化卫星与天线
initSatellite();
initAntenna();% 配置AD9361
configureAD9361();% 信号捕获与处理
[rxData, len] = captureSignal();
[baseband, filtered] = processSignal(rxData);% 解调与解码
[header, payload] = parseModeSFrame(filtered);
displayFlightInfo(payload);
10. 参考文献与资源
- AD9361官方文档：AD9361 User Guide www.analog.com/media/en/technical-documentation/data-sheets/ad9361.pdf
- 代码： 自动化程序(ADS-B端) www.youwenfan.com/contentcni/65533.html
- ADS-B协议标准：EUROCONTROL ASTERIX www.eurocontrol.int/publications/standardisation/asterix