信号稀疏表示与原子分解分析系统

资 源 简 介

本项目旨在通过MATLAB平台实现一个功能完善的信号原子分解软件，系统通过将复杂信号分解为过完备字典中若干原子的线性组合，从而实现信号的高度稀疏表示。该软件集成了多种经典的变换字典，包括但不限于离散余弦变换（DCT）字典、各类常用小波（Wavelet）字典以及Gabor字典，能够根据信号特征灵活选择最优的稀疏基。在算法实现层面，软件内置了标准的匹配追踪（Matching Pursuit, MP）算法和经过优化的正交匹配追踪（Orthogonal Matching Pursuit, OMP）算法，通过迭代迭代过程自适应地寻找与原始信号投影最大的原子项。该系统广泛应用于非平稳信号处理、图像信号压缩、噪声滤除以及复杂特征提取等科研与工程场景。软件不仅提供了高效的计算内核，还通过可视化的方式向用户展示分解过程中的残差收敛曲线、稀疏系数分布以及重构信号与原始信号的对比误差分析，帮助用户直观评估原子分解的效果与精度。

详 情 说 明

信号稀疏表示与原子分解系统

项目介绍

本项目是一个基于 MATLAB 开发的信号稀疏表示与原子分解系统。其核心目标是将复杂的时域信号分解为过完备字典中少量原子的线性组合，从而在保证重构质量的前提下实现信号的高度压缩和特征提取。系统通过迭代搜索的方式，自适应地从包含多种数学特性的原子库中匹配最能代表原始信号特征的分量，广泛适用于非平稳信号处理、噪声滤除和数据压缩等领域。

功能特性

1. 复合信号生成：系统能够生成包含简谐正弦波、余弦成分、瞬态突变脉冲以及局部高频振荡的复杂测试信号，用于模拟真实的工程信号环境。
2. 混合过完备字典构建：集成离散余弦变换（DCT）原子库与多尺度 Gabor 原子库，涵盖了全局频域特征和局部时频特征。
3. 双算法对比：内置正交匹配追踪（OMP）与标准匹配追踪（MP）两种经典的原子分解算法，支持收敛性能的对比分析。
4. 自动化重构与评估：自动计算稀疏度、压缩比、均方误差（MSE）以及峰值信噪比（PSNR），量化评估分解效果。
5. 综合可视化分析：提供多维度的图形界面，包括时域重构对比、残差收敛曲线、稀疏系数分布、误差时域分布以及主要原子贡献度分析。

系统要求

1. 软件环境：MATLAB R2016b 或更高版本。
2. 基础工具箱：MATLAB 基准功能（核心算法基于矩阵运算实现，不依赖特定工具箱）。

实现逻辑详解

系统的运行流程遵循信号处理的标准生命周期，具体步骤如下：

1. 环境初始化与信号合成

系统首先定义采样长度为 256 点，通过叠加不同频率的正弦波和余弦波构建基础信号。随后，在特定时间点植入瞬态脉冲，并在特定区间叠加局部高频信号。最后引入高斯白噪声，以测试算法在含噪环境下的鲁棒性。

1. 过完备字典 D 的构建

字典由两部分拼接而成，确保了对不同性质信号的稀疏表达能力：

• DCT 字典：通过离散余弦变换公式生成 N 个基向量，并进行过归一化处理，主要负责捕捉信号中的平稳频率成分。
• Gabor 字典：采用多尺度策略（尺度因子包含 8, 16, 32）和多频率策略（0.05 至 0.3 范围），利用高斯窗函数对正弦/余弦波进行调制，并通过时移参数在全时域内滑动。这种结构能够极好地匹配信号中的瞬态成分和局部振荡。

1. 稀疏分解计算

系统并行运行两种分解模式：

• 正交匹配追踪模式：在每轮迭代中寻找与当前残差相关性最高的原子，并利用最小二乘法将信号投影到已选原子集合构成的子空间。这种“正交化”处理确保了残差与已选原子始终保持正交，从而加快收敛速度。
• 标准匹配追踪模式：仅通过简单的投影更新系数和残差，不进行正交化处理。通过对比两者的残差能量下降轨迹，验证 OMP 的高效性。

1. 性能指标采集

系统记录每次迭代后的残差 L2 范数，并最终提取非零系数的个数作为稀疏度。根据重构信号与含噪原始信号的差异，计算 MSE 和 PSNR，评估系统对原始信息的还原程度。

1. 结果可视化

系统生成一个包含六个子图的看板：

• 左上：对比原始含噪信号与 OMP 重构信号，展示拟合精度。
• 右上：绘制 OMP 与 MP 的残差随迭代次数下降的对数曲线，展示收敛速度。
• 中左：展示全字典索引下的稀疏系数，反映系数的稀疏分布特性。
• 中右：展示重构后的点对点误差波动。
• 下左：柱状图展示贡献度最大的前 10 个原子及其权重。
• 下右：以文本形式输出详细的性能分析报告。

核心算法与实现细节

1. 字典归一化逻辑

系统在构建 DCT 和 Gabor 原子时，均强制执行归一化操作（L2 范数为 1）。这是确保匹配追踪算法中相关性计算（内积运算）能够准确反映原子能量贡献的关键前提。

1. OMP 中的最小二乘更新

在 OMP 的实现中，系统并非简单地累加系数，而是通过求解 (D_active' * D_active) (D_active' * y) 来重新评估所有已选原子的最优系数。这一步骤有效消除了原子间非正交带来的重叠干扰。

1. 残差迭代终止条件

系统设置了双重退出机制：一是最大迭代次数（默认 50 次，即期望稀疏度），二是目标误差阈值（1e-4）。只要满足其一，计算即停止，从而平衡了计算效率与分解精度。

1. Gabor 原子的时频构造

Gabor 原子通过高斯包络与正余弦调制的乘积实现。通过 floor(N/8) 的步长进行时移，既保证了时域覆盖的完备性，又控制了字典的总体规模，避免了过高的计算复杂度。

使用方法

1. 启动 MATLAB 软件。
2. 将项目相关的脚本文件放置在 MATLAB 的当前工作路径下。
3. 运行主处理脚本。
4. 观察自动生成的图形界面。
5. 在图形界面的右下角区域查看系统自动生成的“信号稀疏表示分析报告”，获取具体的稀疏度、压缩比和 PSNR 指标。