循环平稳信号处理子程序库

资 源 简 介

该项目旨在开发一套完整的MATLAB子程序集成工具，专门用于分析和处理具有循环平稳特性的非平稳信号。循环平稳性是指信号的统计特性（如均值和自相关函数）随时间呈周期性变化的性质，这种特性广泛存在于通信信号、雷达回波、由于旋转或往复运动产生的机械振动信号中。

详 情 说 明

循环平稳信号处理子程序库

项目介绍

本项目是一个专门针对具有循环平稳特征的非平稳信号进行分析和处理的MATLAB集成工具。循环平稳性是通信信号、雷达回波及旋转机械振动信号中的一种核心属性，其统计特性（如自相关函数）随时间呈周期性变化。本项目通过实现经典的循环自相关函数计算与高效的谱相关密度估计（FAM算法），为研究信号的调制识别、参数估计和特征提取提供了完整的软件基础。

功能特性

• 标准信号仿真：内置BPSK调制信号生成引擎，支持自定义采样频率、载波频率和码速率，并可注入高斯白噪声以模拟真实环境。
• 时延-循环频率域分析：通过二维网格搜索实现循环自相关函数（CAF）的计算，揭示信号在不同时延下的周期相关特性。
• 高效谱相关估计：采用FFT累积法（FAM）实现谱相关密度（SCD）的快速估计，提供了远高于传统时域平滑算法的计算效率。
• 循环平稳特征提取：自动检测Alpha截面上的显著循环频率峰值，并量化计算循环平稳度指标（DCS）。
• 多维度可视化：提供时域波形、CAF三维曲面、SCD二维平视图及三维网格图、Alpha谱线特征图等全方位的视觉分析手段。

实现逻辑说明

项目的主处理流程严格按照循环平稳信号处理的标准理论体系构建，其核心执行逻辑如下：

1. 信号环境构建：

- 首先配置系统参数，包括采样率1000Hz、载波150Hz和码速率50Hz。 - 生成随机比特序列，通过过采样和载波调制合成BPSK信号，最后叠加指定强度的加性高斯白噪声。

1. 循环自相关计算：

- 在预设的循环频率范围（$alpha$）和时延范围（$tau$）内，通过嵌套循环计算离散时间循环自相关函数。 - 实现计算公式为：$R_x^alpha(tau) = E[x(t+tau/2)x^*(t-tau/2)exp(-j2pialpha t)]$。

1. FAM算法执行：

- 第一阶段（信道化）：使用Hamming窗对输入信号进行滑动分段，通过第一级FFT将信号转化为时频分布。 - 第二阶段（复解调与乘积）：对所有频率通道对进行复共轭相乘，提取频率平移之后的成分。 - 第三阶段（循环频率估计）：对乘积序列进行第二级FFT，根据频率差值和采样步进将结果映射到$(alpha, f)$坐标系中。

1. 特征分析与量化：

- 在频谱频率$f=0$处提取SCD的截面，即Alpha谱。 - 使用峰值检测算法识别显著的循环频率分量，这些分量对应于载波频率和码速率的相关特征。 - 计算循环平稳度（DCS）：通过统计非零循环频率区域的能量占总能量的比重，评估信号的周期平稳强度。

1. 图形化输出：

- 自动生成四个图形窗口，涵盖了从原始时域数据到高阶统计特征的全部处理结果。

关键算法与算法细节分析

• FAM算法 (FFT Accumulation Method)：该算法是本项目中最复杂的计算模块。它通过两级FFT结构，将计算复杂度从$O(N^2)$有效降低。算法中引入了滑动长度$L$和通道数$Np$的概念，通过在频域利用相关性来逼近谱相关密度的理论值。代码中精确处理了频率轴和循环频率轴的映射关系，确保了物理单位（Hz）的准确性。
• DCS指标计算：为了量化信号的循环平稳性，代码定义了一个阈值（$alpha$轴中心区域之外），将能量划分为平稳部分和循环平稳部分。这种比例关系对于识别背景噪声中的弱信号具有重要参考价值。
• 峰值检测逻辑：利用findpeaks函数并结合最小峰值高度阈值（最大值的30%），有效过滤了由于有限观测时间带来的估计随机波动，确保检测到的循环频率具有统计显著性。

系统要求

• 环境要求：MATLAB R2016b 或更高版本。
• 所需工具箱：Signal Processing Toolbox（信号处理工具箱）。
• 硬件配置：建议内存4GB以上，以支持FAM算法中大矩阵的存储与运算。

使用方法

1. 确保MATLAB环境中已安装必要的信号处理工具箱。
2. 将主程序代码及相关函数定义保存在同一路径下。
3. 在MATLAB命令行窗口输入该主函数名称并回车。
4. 程序将自动执行信号生成、算法计算并弹出分析图表。
5. 观察命令行输出，获取检测到的循环频率峰值数值。