快捷信号频谱分析与可视化绘图工具

资 源 简 介

本项目开发了一个高度集成的MATLAB工具函数，旨在彻底解决科研人员在进行信号处理时需要频繁编写冗长的FFT变换代码的痛点。该函数将复杂的频域转换逻辑进行完全封装，核心功能包括：自动计算信号长度并执行高效的快速傅里叶变换（FFT）；根据采样频率自动生成精确的频率轴向量（单位为Hz）；支持自动幅度归一化处理，确保频谱幅度与时域信号能量的一致性。该函数能一键生成标准化的频谱图，支持从单边谱到双边谱的快速切换，并提供对数坐标（dB）与线性坐标的转换选项。其应用场景涵盖了通信仿真数据观测、传感器采集数据快速预览、

详 情 说 明

MATLAB便捷频谱分析与可视化工具

项目介绍

本项目通过封装MATLAB复杂的底层信号处理逻辑，提供了一个高度集成、简洁易用的频谱分析函数。该工具旨在通过极简的调用方式，帮助科研人员和工程师完成从原始时域信号到标准化频域图像的转换。开发者无需手动计算频率分量、处理FFT移位或进行复杂的能量归一化，只需输入原始信号与采样率，即可获得具备出版质量的分析图表。

功能特性

1. 高度封装： 一行代码实现“FFT计算-幅值归一化-频率轴生成-绘图标注”全流程。
2. 多样化显示： 支持单边谱（能量补偿）与双边谱（中心化显示）的无缝切换。
3. 尺度自定义： 提供线性幅值与对数分贝（dB）两种显示尺度，方便观察微弱分量。
4. 自动标注： 内置频域峰值检测逻辑，自动捕捉并标注图中最大幅值点及其频率。
5. 数据交互： 既支持直接绘图可视化，也支持返回频率向量、幅值向量和相位信息进行二次开发。
6. 稳健性： 自动处理输入信号向量方向，内置防止计算对数时产生无穷值的溢出保护。

详细实现逻辑

在演示示例中，系统通过两个主要环节展示了功能：

1. 信号模拟环节 系统设定了固定的采样频率和时长，生成了一个由50Hz和120Hz组成的复合正弦信号。为模拟真实工程环境，在信号中加入了标准正态分布的零均值白噪声，从而建立起一个含有背景干扰的待测信号源。

2. 核心分析环节 核心功能通过解析用户传入的可选参数，动态执行以下逻辑：

• 多参数解析： 采用输入解析器确保参数的灵活性，用户可选择FFT点数、显示模式（单双边）、幅值单位（线性/dB）以及是否静默执行。
• 频域转换与归一化： 执行FFT变换，并严格按照信号长度进行幅值归一化处理（除以L），确保所得幅值与时域原始振幅一致。
• 能量补偿机制： 在单边谱模式下，除直流（DC）和奈奎斯特频率分量外，对其余分量乘以2，以补偿因截取一半频率分量导致的能量折损。
• 中心化处理： 在双边谱模式下，利用位移操作将零频分量移至频谱中心，并对应生成负频率轴。
• 自适应绘图： 根据用户选择的Scale参数动态计算幅值。如果是dB模式，则通过加微小值（eps）防止对零取对数；绘图完成后通过坐标轴微调实现专业视觉效果，并自动在曲线上标记最高峰。

关键函数与算法分析

• 快速傅里叶变换 (FFT)： 算法底层基于FFT实现，支持自定义N点变换。
• 幅值归一化算法： 核心逻辑为 abs(Y/L)，其中L为原始信号长度。这保证了当输入为正弦波时，频谱中对应的峰值点正好等于正弦波的物理幅值。
• 单边谱能量恢复： 实现中对单边谱执行了 P1(2:end-1) = 2*P1(2:end-1)，这是从双边谱导出物理含义明确的单边谱的标准算法。
• 分贝尺度换算： 使用 20*log10 转换，将线性幅值映射为电平强度，通过控制计算精度防止负无穷大情况。
• 自动峰值追踪： 使用索引查找算法实时定位当前显示窗口内的最大值，并通过动态文本框显示其精确频率。

使用方法

在MATLAB环境中，您可以直接调用主演示程序来观察效果。在实际开发中，调用方式如下：

1. 基础调用： 直接传入信号向量和采样频率。
2. 高级配置： 通过 'Mode' 参数切换 'single' 或 'double'；通过 'Scale' 参数切换 'linear' 或 'log'。
3. 静默模式： 设置 'NoPlot' 为 true，函数将不弹出图形窗口，仅返回计算后的频率轴、幅值和相位数组。

系统要求

• 运行环境： MATLAB R2016b 或更高版本。
• 所需工具箱： 仅需MATLAB基础环境（无需专门安装信号处理工具箱，代码采用底层逻辑实现）。