GPS信号BPSK与BOC调制功率谱对比仿真系统

资 源 简 介

本项目基于MATLAB环境，开发了一套完整的全球定位系统（GPS）信号调制与频谱分析仿真工具。项目主要功能包括两个核心部分：一是实现传统GPS信号使用的二进制相移键控（BPSK）调制，二是实现现代化GNSS信号（如Galileo和GPS L1C）采用的二进制偏移载波（BOC）调制。程序首先生成标准的伪随机噪声（PRN）序列作为扩频码，然后分别应用BPSK和BOC调制算法生成时域信号。更为关键的是，该项目重点实现了对这两种调制信号功率谱密度（PSD）的计算与仿真，采用周期图法或韦尔奇法（Welch）进行谱估计。系统将BPSK和BOC信号的功率谱绘制在同一坐标系下进行直观对比，清晰展示了BOC信号特有的频谱分裂（Spectrum Splitting）现象，以及其相较于BPSK在频带边缘的能量分布差异。这一比较有助于深入理解BOC调制在提升频谱利用率、增强抗干扰能力以及改善多径抑制性能方面的优势，非常适合用于卫星导航信号体制的研究与教学。

详 情 说 明

GPS信号BPSK与BOC调制功率谱对比仿真系统

项目简介

本项目是一个基于MATLAB环境开发的全球定位系统（GPS）信号仿真工具。它专注于对比分析由于调制方式不同而产生的频域特性差异，具体比较了传统GPS信号使用的二进制相移键控（BPSK）调制与现代化GNSS信号（如Galileo和GPS L1C）采用的二进制偏移载波（BOC）调制。

通过仿真生成时域波形并利用韦尔奇法（Welch）进行功率谱密度估计，本系统直观地展示了BOC调制特有的频谱分裂（Spectrum Splitting）现象，不仅有助于理解卫星导航信号的调制原理，也验证了BOC信号在频带边缘能量分布上的特征。

功能特性

• 参数化仿真设置：基于GPS L1频段基准频率（1.023 MHz）构建仿真环境，支持高采样率（30倍基频）以保证波形保真度。
• 信号生成与调制：

* 生成随机的二进制序列模拟伪随机噪声（PRN）扩频码。 * 实现标准的BPSK调制逻辑。 * 实现BOC(1,1)调制逻辑，生成特定频率的副载波并与扩频码混频。

• 时域波形分析：可视化展示基带PRN码、BOC副载波以及最终调制信号的时域波形细节。
• 功率谱密度（PSD）估计：

* 采用韦尔奇法（Welch's method）进行平滑的谱估计，降低随机噪声影响。 * 支持双边谱计算，将频谱中心移至0Hz，便于观察频谱的对称性。

• 频谱对比可视化：在同一坐标系下绘制BPSK与BOC的归一化功率谱，清晰标注主瓣峰值位置及能量迁移现象。

系统要求

• MATLAB R2016b 或更高版本（推荐使用较新版本以获得更好的图形渲染支持）。
• Signal Processing Toolbox（信号处理工具箱，用于pwelch、rectpulse等函数）。

核心算法与实现逻辑

本项目的主程序脚本（main函数）严谨地按照信号处理流程执行，具体实现细节如下：

1. 系统参数初始化

程序首先定义了基准频率 $f_0 = 1.023 text{ MHz}$。采样频率被设定为 $f_s = 30 times f_0$，以满足奈奎斯特采样定理并确保时域可视化的平滑度。设定了BOC(1,1)调制的关键参数（$alpha=1, beta=1$），确定了副载波频率与码速率一致。

2. 信号生成模型

• PRN序列：利用随机函数生成长度为8192位的双极性（+1, -1）随机序列，用于模拟扩频码。
• 基带成型：使用 rectpulse 函数对离散的码片序列进行上采样（零阶保持），将其转换为具备实际采样率的时间连续信号，即 BPSK 基带信号。
• 副载波生成：根据BOC原理，生成频率为 $f_{sc} = 1.023 text{ MHz}$ 的正弦波，并通过符号函数 sign 将其整形为方波，作为BOC调制的副载波。
• 调制实现：

* BPSK：直接采用基带扩频码流。 * BOC：将基带扩频码流与副载波点对点相乘，实现信号频谱的搬移。

3. 时域可视化

程序截取了信号的前20个码片长度进行绘图，避免因数据过密导致显示不清。绘制了三个子图：

1. 基带PRN数据：展示原始扩频码的跳变。
2. 副载波：展示BOC调制所用的方波时钟信号。
3. BOC调制信号：展示最终生成的调制波形，体现了码片内的电平翻转特性。

4. 功率谱密度（PSD）估计

这是本项目的核心分析部分，采用 韦尔奇法 (Welch's Method) 进行计算：

• 窗口设置：使用长度为1024的汉明窗（Hamming window），重叠率为50%，以平衡频率分辨率和方差。
• FFT参数：FFT点数设为4096，保证了频域显示的细腻度。
• 双边谱处理：通过辅助逻辑 calc_shifted_psd，计算信号的功率谱，并进行 fftshift 操作或使用 centered 参数，将频谱的0频分量移至坐标轴中心。这使得观察正负频率的分布（双边谱）成为可能。
• 归一化：将计算出的功率谱除以最大值并转换为分贝（dB）刻度，使得主瓣峰值归一化为0dB，便于比较波形形状而非绝对功率。

5. 频域对比结果

程序最终生成一张对比图，横轴为相对频率（MHz），纵轴为归一化PSD（dB/Hz）：

• BPSK曲线：显示能量集中在载波中心（0 MHz），呈现典型的 sinc 函数形状，旁瓣随频率增加而衰减。
• BOC(1,1)曲线：直观展示了频谱分裂现象。能量从0 MHz中心向两侧迁移，在 $pm 1.023 text{ MHz}$ 处形成两个主峰。
• 标注：图表中自动标注了BPSK的主瓣峰值和BOC的两个分裂峰值。

使用方法

1. 将代码保存为 MATLAB 脚本文件（例如 main.m）。
2. 在 MATLAB 命令窗口中直接运行该函数 main。
3. 程序运行结束后，将自动弹出两个图形窗口：

* 图1：时域波形对比（PRN码、副载波、BOC信号）。 * 图2：功率谱密度对比（BPSK vs BOC），展示频谱分裂特性。

1. 控制台会输出 "仿真完成" 的提示信息。