基于Simulink的GPS信号与导航电文生成仿真系统

资 源 简 介

本项目利用MATLAB/Simulink平台构建高精度的GPS卫星信号仿真模型，重点实现GPS L1频点信号的完整生成流程。系统功能涵盖核心的扩频码生成与信号调制，具体包括：1. C/A码生成模块，采用两个10级线性反馈移位寄存器（LFSR）组合生成Gold码序列，支持用户指定卫星PRN编号以生成对应的正交码；2. 导航电文（D码）生成模块，模拟50bps的卫星导航数据流，构建包含遥测字（TLM）和交接字（HOW）的标准帧结构；3. 信号调制模块，将C/A码与导航电文进行异或运算实现扩频，并采用二进制相移键控（BPSK）方式调制到射频载波上。该项目能够提供可配置的仿真环境，用于验证GPS接收机捕获与跟踪算法的有效性，并直观展示直扩信号的时域波形与频域特性。

详 情 说 明

基于MATLAB的GPS信号与导航电文生成仿真系统

项目介绍

本项目是一个高精度的GPS卫星信号仿真模型，旨在通过MATLAB脚本（.m文件）完整复现GPS L1频点信号的生成、扩频、调制及传输过程。系统以软件定义无线电（SDR）的思路构建，不依赖高层封装工具箱，而是从底层的比特级逻辑开始实现，包括C/A码的LFSR生成、导航电文的构建、DSSS直接序列扩频技术以及BPSK载波调制。该仿真平台支持用户通过修改参数配置不同卫星（PRN）、信噪比环境及采样频率，并提供了时域波形与频域功率谱密度的详细分析功能，用于验证GPS信号的基本特性。

功能特性

• 参数化仿真配置：支持自定义卫星PRN编号（1-32）、载波频率、采样率、仿真时长及信噪比（SNR）。
• 底层C/A码生成器：基于ICD-GPS-200标准，使用双LFSR（线性反馈移位寄存器）手动实现Gold码生成，包含相位抽头选择逻辑。
• 信号扩频与调制：实现了导航电文与C/A码的异或扩频（双极性乘法）以及从基带到射频的BPSK调制。
• AWGN信道模拟：能够根据设定的SNR值计算信号功率与噪声功率，向信号中添加高斯白噪声。
• 多维可视化分析：提供码片序列、基带波形对比、射频时域波形（局部放大）及功率谱密度（PSD）的图表展示。

系统要求

• MATLAB R2016a 或更高版本
• Signal Processing Toolbox（用于信号谱分析函数 pwelch）

使用方法

1. 打开MATLAB软件，定位到项目根目录。
2. 打开主脚本文件，根据需要调整代码顶部的“系统参数设置”部分（如 sat_prn 修改卫星编号，snr_db 修改信噪比）。
3. 运行脚本。
4. 观察控制台输出的参数摘要以及自动弹出的信号分析图表。

---

实现原理与核心算法分析

本项目的核心逻辑在主脚本中顺序执行，主要包含以下五个关键处理模块：

1. C/A码生成模块（LFSR架构）

系统不仅是调用库函数，而是通过模拟硬件寄存器行为生成C/A码。

• 寄存器初始化：初始化两个10级寄存器（G1和G2）为全1状态。
• 多项式反馈逻辑：

* G1生成器：采用多项式 $1 + x^3 + x^{10}$，反馈抽头为第3位和第10位。 * G2生成器：采用多项式 $1 + x^2 + x^3 + x^6 + x^8 + x^9 + x^{10}$，反馈抽头为第2、3、6、8、9、10位。

• 相位选择与Gold码合成：

* 脚本内置了PRN 1-32对应的G2输出相位抽头表（prn_taps）。 * 根据用户设置的 sat_prn，选取G2寄存器的两个特定位置进行模2加，产生延迟的G2序列。 * 将G1的输出与相移后的G2输出进行异或运算（模2加），生成长度为1023的Gold码序列。

• 码型转换：将生成的0/1逻辑序列转换为双极性非归零码（-1/1），以便后续数学运算。

2. 导航电文（数据流）生成模块

该模块模拟卫星下发的低速数据流。

• 数据量计算：根据仿真时长（默认0.040s）和数据速率（50bps）计算所需生成的比特数。
• 随机序列生成：当前实现中采用 randi 函数生成随机的0/1比特流来模拟导航电文，未构建具体的TLM/HOW帧结构，但保留了标准的50bps速率特性。
• 双极性转换：同样将数据比特转换为双极性电平（-1/1）。

3. 信号扩频与调制模块

此部分实现了直扩信号的核心处理。

• 采样点索引映射：

* 创建仿真时间轴 t。 * C/A码过采样：利用 floor(mod(t * code_rate, ca_code_period)) 算法，将1023位的C/A码映射到高采样率的时间轴上，实现码片的周期性重复。 * 数据码过采样：利用 floor(t * data_rate) 算法，将低速导航电文映射到时间轴上，使每个数据比特持续20ms。

• 扩频运算：将过采样后的导航电文序列与C/A码序列进行点对点乘法运算（双极性下的异或），实现频谱扩展。
• BPSK调制：生成频率为 fc（默认4MHz）的余弦载波，并与扩频后的基带信号相乘，完成二进制相移键控调制。

4. 信道模拟（AWGN）

模拟信号在传输过程中的损耗与干扰。

• 功率计算：计算调制信号的平均功率。
• 噪声生成：根据设定的信噪比 snr_db（默认-20dB），反向推算所需的噪声功率。
• 加噪：生成对应功率的高斯白噪声并叠加到发射信号上，形成最终的接收信号 received_signal。

5. 结果可视化

脚本运行结束后会生成一个包含四个子图的分析窗口：

• 子图1：展示生成的C/A码的前50个码片（Stem图），用于验证码序列生成的离散特性。
• 子图2：对比基带导航电文（D码）与扩频码（C/A码）的时域波形，直观展示扩频前后信号频率的巨大差异。
• 子图3：接收信号的局部放大波形图，展示载波细节以及噪声对幅度的影响。
• 子图4：接收信号的功率谱密度（PSD），利用Welch法估算。图中可以清晰观察到以载波频率为中心的主瓣，以及C/A码扩频特有的频谱形状（主瓣宽度约为2MHz），验证了GPS信号的频谱特性。