数字调制信号循环谱相关分析与可视化系统

资 源 简 介

本项目旨在构建一个基于MATLAB的数字调制信号循环平稳特性分析平台。项目主要包含三个核心模块：信号生成、核心算法计算及结果可视化。首先，信号生成模块能够模拟产生BPSK、QPSK、AM、FM等多种常见的数字调制信号，支持用户自定义采样率、载波频率、符号速率及信噪比（SNR），以模拟真实的通信环境。其次，核心计算模块基于循环平稳信号处理理论，利用时间平滑或者频率平滑算法（如FFT累加法FAM或条带谱相关算法SSCA）计算信号的谱相关密度（Spectral Correlation Density, SCD），该算法能有效抑制平稳噪声并提取信号的周期性特征。最后，可视化模块将计算得到的四维数据（循环频率、谱频率、幅值、相位）进行图形化展示，重点功能是绘制以循环频率（alpha）和常规频率（f）为平面的三维循环谱幅值图，直观呈现信号的特定循环平稳特征峰。同时，系统支持交互式或自动切割功能，提取并画出特定循环频率处的二维截面图（如循环特征切片）以及特定谱频率下的切片，用于辅助信号的参数估计、盲识别及特征检测。

详 情 说 明

数字调制信号循环谱相关分析与可视化系统

1. 项目介绍

本项目是一个基于MATLAB构建的数字调制信号循环平稳特性分析平台。该系统集成了信号模拟生成、信道噪声叠加、核心循环谱计算（基于FAM算法）以及结果可视化功能。项目旨在通过计算谱相关密度（Spectral Correlation Density, SCD），提取信号在循环频率域和常规频率域的独特的二维特征，从而在低信噪比环境下有效区分和识别BPSK、QPSK、AM、FM等不同调制类型的信号。

2. 功能特性

• 多类型信号生成：支持产生BPSK、QPSK、AM（常规调幅）、FM（调频）四种典型的通信信号。
• 参数灵活配置：允许自定义采样频率、信号时长、载波频率、符号/码元速率以及信噪比（SNR）。
• 模拟真实信道：内置高斯白噪声（AWGN）生成模块，可根据指定的SNR对信号进行加噪处理。
• 高效核心算法：实现了基于FFT累加法（FAM）的循环谱计算核心逻辑。采用时间抽取（Channelization）与频率平滑相结合的方式，通过汉明窗加窗处理和通道化FFT，计算信号的循环谱密度。
• 独立运行设计：代码内部实现了基础信号处理所需的辅助函数（如矩形脉冲成型），减少了对MATLAB特定工具箱（如通信工具箱）的依赖。

3. 系统要求

• 软件环境：MATLAB
• 工具箱依赖：核心算法主要依赖基础MATLAB函数。代码中已手动实现 rectpulse 等函数以替代工具箱依赖，增强了兼容性。

4. 使用方法

直接运行主脚本函数即可启动系统。程序将按照以下流程自动执行：

1. 初始化系统参数（采样率、载波、调制类型等）。
2. 生成对应的数字调制信号并添加高斯白噪声。
3. 利用FAM算法计算信号的循环谱矩阵。
4. 统计并输出核心算法的计算耗时。
5. 调用可视化模块绘制分析结果（注：可视化函数的具体实现在提供的代码片段之外，但在主流程中被调用）。

5. 核心代码逻辑与实现细节分析

以下是对系统主要模块及其在代码中实际实现逻辑的详细分析：

5.1 系统主控模块

主函数作为程序的入口，负责全局资源的调度。它首先定义了关键的系统参数结构体 sys_params，包括1000Hz的采样率、250Hz的载波频率、50 Baud的码元速率以及10dB的信噪比。同时，设定了FAM算法的配置参数，如窗口点数（Np=64）和滑动步长（L=16）。主流程清晰地展示了"参数配置 -> 信号生成 -> 噪声叠加 -> FAM计算 -> 结果显示"的线性处理逻辑，并包含了一个计时器用于监测核心算法的性能。

5.2 信号生成模块

该模块根据用户指定的调制类型，通过数学公式直接合成基带和频带信号：

• BPSK：生成随机比特流，变换为双极性码（+1/-1），通过自定义的矩形脉冲成型函数扩展采样点，最后乘以余弦载波。
• QPSK：生成0-3的随机整数，映射到复数星座图（$pi/4$ 相位偏移），分离实部（I路）和虚部（Q路）分别进行脉冲成型，最后执行正交调制（I路乘余弦，Q路乘负正弦）。
• AM：基于单频正弦波（20Hz）构造基带信号，并添加直流偏置（调制度0.8），实现常规双边带调幅。
• FM：基于单频正弦波（15Hz）作为消息信号，对其进行累积求和（模拟积分）以获得相位偏移，结合调频灵敏度生成调频信号。

5.3 噪声处理与辅助功能

• 高斯白噪声添加：代码包含一个独立的噪声函数，参考信号的平均功率，根据设定的SNR分贝值反推噪声功率，利用 randn 生成标准正态分布序列并按比例叠加到原始信号上。
• 自定义脉冲成型：为了避免依赖信号处理工具箱中的 rectpulse 函数，代码内部通过矩阵重排（repmat 和 reshape 操作）实现了一个等效的矩形脉冲扩展函数，确保代码在基础版MATLAB中也能运行。

5.4 核心FAM算法实现

这是系统的核心计算模块，代码实现了FFT累加法（FAM）的主要步骤，用于估计循环谱密度：

1. 数据加窗与通道化：

* 首先计算分块数量，初始化通道矩阵。 * 使用汉明窗（Hamming window）对信号进行加窗处理。 * 执行短时傅里叶变换（STFT）风格的处理，将长信号分割、加窗并在时间轴上进行FFT，将一维时域信号转换为“频率-时间”的二维矩阵 $X(f, t)$。

1. 谱相关计算：

* 代码采用了基于频率差遍历的优化逻辑。外层循环遍历频率通道差值 $p$（对应粗略的循环频率）。 * 根据差值 $p$，筛选出对应的通道对索引 $(k, l)$，使得 $k - l = p$。 * 时域相乘：对选定的通道行向量进行共轭相乘，计算 $X_k(t) cdot X_l^*(t)$。 * FFT变换：对上述乘积结果沿时间轴执行FFT，从而获得特定循环频率下的谱特征。 * 映射与累加：根据通道索引 $(k+l)$ 确定结果在常规谱频率轴（$f$）上的位置，将计算出的幅度谱累加到输出矩阵 S_out 中。该步骤实现了将计算出的相关值正确映射到 $(f, alpha)$ 平面的逻辑。

5.5 可视化接口

代码末尾调用了一个可视化函数，传入了时间轴、含噪信号、计算出的循环谱矩阵（Sx）、循环频率向量（alpha）和谱频率向量（f）。虽然该函数的具体绘图代码未包含在主要逻辑块中，但接口的存在表明系统支持将计算得到的二维/三维数据进行图形化展示。