基于墨西哥草帽小波的连续小波变换系统

资 源 简 介

该项目实现了基于墨西哥草帽小波（Mexican Hat Wavelet）函数的连续小波变换（CWT）算法。 墨西哥草帽小波是高斯函数的二阶导数，具有极佳的时频局部化特性，被广泛用于捕捉信号特征。 系统核心功能是对比输入信号与不同尺度的小波基函数，通过卷积运算提取信号在不同频率窗口下的能量分布。 它可以有效地应用于非平稳信号的奇异性检测、峰值识别以及复杂动态系统的瞬态特征提取。 用户能够通过调节尺度参数，实现对信号从微观细节到宏观趋势的全方位多尺度透视，为信号处理提供可靠依据。

详 情 说 明

基于墨西哥草帽小波的连续小波变换系统

项目介绍

本项目实现了一个基于墨西哥草帽（Mexican Hat）小波函数的连续小波变换（CWT）处理系统。该系统通过数学建模和卷积运算，能够将一维时域信号分解至二阶导数高斯基函数构建的尺度空间中。这种变换方式特别适用于检测信号中的突变点、瞬态冲击以及分析非平稳信号在高低频率下的能量分布特征，为信号处理提供从微观细节到宏观趋势的全方位多尺度视角。

功能特性

1. 多分量信号模拟：系统内置了一个包含平稳简谐波、局部高频正弦波以及高斯脉冲冲击的复合测试信号，模拟复杂的实际工业或物理传感数据。
2. 自定义小波核构建：严格按照墨西哥草帽小波的物理公式生成基函数，并根据不同尺度参数动态调整支撑范围，确保时频分析的准确性。
3. 时频能量映射：通过计算变换系数的模方，生成能量分布矩阵，直观展示信号能量在时间-尺度（频率）平面上的演化。
4. 多尺度特征提取：支持提取特定尺度下的变换系数。小尺度用于捕捉信号的极值点和奇异性（高频细节），大尺度用于提取信号的包络和慢变趋势（低频成分）。
5. 多架构数值可视化：通过四分图布局，同步对比原始信号、时频能量热力图以及多尺度分量提取结果。

系统要求

• 软件环境：MATLAB R2016b 或更高版本。
• 工具箱需求：基础 MATLAB 核心模块（主要利用其矩阵运算与绘图功能，无需额外的小波工具箱，算法为原生实现）。

核心实现逻辑

系统运行遵循以下标准流程：

1. 参数定义与环境初始化：设置 1000Hz 采样频率，生成 1 秒长度的时间序列。
2. 信号合成：

* 叠加基础 50Hz 余弦波。 * 在 0.5 秒后的区间加入 120Hz 正弦信号，模拟频率突变。 * 在 0.3 秒处注入高斯脉冲，模拟瞬态干扰或系统奇异性。

1. 尺度空间定义：设定 1 至 100 的连续尺度向量，用于控制小波函数的伸缩。
2. 循环卷积计算：

* 针对每一个尺度，确定小波的支撑范围（通常为尺度的 5 倍），保证核函数的完整性。 * 计算归一化时间变量，并应用墨西哥草帽小波公式生成核函数。 * 使用 L2 归一化对核函数进行能量补偿。 * 利用卷积运算将核函数与原信号匹配，获取该尺度下的变换系数。

1. 能量谱转化：对系数矩阵进行平方运算，获取解析能量谱。
2. 可视化输出：利用颜色映射（Jet Colormap）和图像缩放技术，将系数矩阵转化为易于观察的时频图像。

算法与实现细节

• 墨西哥草帽小波公式：算法采用了高斯函数的二阶导数形式，数学表达为 $psi(t) = frac{2}{sqrt{3}pi^{1/4}}(1-t^2)e^{-t^2/2}$。该函数具备良好的局部化特性，且其均值为零。
• 尺度与频率的关系：在代码实现中，尺度 $s$ 直接影响小波核的宽度。较小的尺度对应高频，能够精准定位脉冲位置；较大的尺度对应低频，能够平滑噪声并显示信号的大致轮廓。
• 卷积模式：采用 'same' 模式卷积，确保变换后的每一行系数与原始输入信号的时间点完全对齐。
• 特征对比机制：系统专门选取了尺度 5（细节）和尺度 40（趋势）进行同框对比，展示了 CWT 在信号分解和特征分离方面的强大能力。

使用方法

1. 启动 MATLAB 软件。
2. 将系统程序脚本置于当前工作路径。
3. 在命令行窗口输入主函数名称并回车运行。
4. 程序将自动计算并弹出波形分析窗口，同时在命令行输出信号采样点数、尺度范围及能量分布矩阵的维度信息。
5. 用户可以通过观察“时频能量分布谱”来识别信号中 0.3s 处的脉冲能量集中区域，以及 0.5s 后高频分量的加入。