小波变换与高斯拟合在光谱重叠峰解析中的应用
文件信息
- 云盘路径：/工作相关/Personnal/2017工作/质谱仪软件架构资料调研/相关论文/小波变换与高斯拟合在光谱重叠峰解析中的应用_冯飞.pdf
- 文件大小：1021 KB
- 页数：6页
- 提取方式：PDF文本提取(pdftotext)
- 日期：2026-09-20
文献信息
- 期刊：光子学报 (ACTA PHOTONICA SINICA)
- 卷期：第44卷第6期，2015年6月
- DOI：10.3788/gzxb20154406.0630001
- 作者：冯飞，王府北，谢非，吴琼水，曾立波
- 单位：武汉大学电子信息学院，武汉430072
- 基金：国家科技支撑计划课题（No.2011BAF02B02）
核心观点
本文提出了一种基于小波变换与高斯拟合的光谱重叠峰解析算法，能够同时对含噪声和基线漂移的光谱信号进行去噪、基线校正和峰位检测，并通过高斯拟合求解峰强和峰宽信息。
方法步骤
第一步：小波变换预处理
- 小波函数选择：sym8小波（综合去噪和基线校正效果最佳）
- 分解层数：8层（6-9层对比实验确定）
- 阈值方法：Heursure软阈值法
- 处理流程：
- 对低频系数置零 → 基线校正
- 对高频系数阈值量化 → 谱线去噪
- 小波重构 → 得到预处理后的光谱信号
第二步：多尺度小波变换寻峰
- 母小波函数：Mexican Hat (Mexh) 函数
- 尺度参量：1至20，间隔为1
- 脊线搜索：
- 搜寻各尺度小波系数模极大值
- 连接相邻尺度的模极大值形成脊线
- 复合脊线校正：剔除偏离过大的极值点
- 峰位确定：脊线中小波尺度参量最大时对应的位置即为峰位
第三步：高斯函数谱峰拟合
- 单峰高斯函数：g(u,σ)(t) = A·exp(-(t-u)²/2σ²)
- 多峰拟合：f’(t) = Σg(uᵢ,σᵢ)(t)
- 优化目标：在最小二乘意义下寻找最佳拟合光谱
- 输出参数：峰位(u)、峰强(A)、峰宽(σ)
实验结果
仿真光谱参数
- 6个高斯峰，包含：
- 2个完全独立峰（峰位50、100）
- 2个完全重叠峰（峰位150、160）
- 2个部分重叠峰（峰位230、250）
- 叠加随机噪声（信噪比20dB）和基线漂移
峰位检测精度
| 峰号 | 真实值 | 拟合值 | 偏差 |
|---|---|---|---|
| 1 | 50 | 50.014 | +0.014 |
| 2 | 100 | 100.031 | +0.031 |
| 3 | 150 | 149.067 | -0.933 |
| 4 | 160 | 158.977 | -1.023 |
| 5 | 230 | 230.342 | +0.342 |
| 6 | 250 | 249.617 | -0.383 |
- 峰位偏差：±1.3以内
- 峰强偏差：不超过8.5%
与OrIGIN软件对比
- 峰位拟合：本文算法除峰5外均优于OrIGIN
- 峰强拟合：本文算法除峰1和峰6外均优于OrIGIN
- 峰宽拟合：本文算法在峰1、2、4、5上优于OrIGIN
关键创新点
- 同时处理噪声和基线：通过小波变换低频系数置零实现基线校正，高频系数阈值化实现去噪
- 多尺度脊线搜索：利用小波系数模极大值的多尺度特性，精确定位峰位
- 复合脊线校正：提出复合脊线的判定和校正方法，提高重叠峰检测准确性
- 高斯拟合参数优化：基于小波变换得到的峰位数和峰位值作为高斯拟合的初始参数
应用价值
- 通用性强：适用于多种光谱类型（原子发射光谱、红外光谱、拉曼光谱等）
- 非平稳信号：不仅适用于平稳信号，还可推广到非平稳信号处理
- 不依赖先验知识：无需了解待测样品的成分和含量信息
- 可推广性：可应用于其他谱图重叠峰解析场景
研究局限
- 对于完全重叠峰（峰3、4），峰位拟合精度相对较低
- 需要人工设定小波函数类型和分解层数等参数
- 计算复杂度较高，实时性受限
参考文献
- 方勇华, 孔超, 兰天鸽, 等. 应用小波变换实现光谱的噪声去除和基线校正[J]. 光学精密工程, 2006, 14(6): 1089-1091.
- 蔡涛, 王先培, 杜双育, 等. 基于多尺度小波变换的红外光谱谱峰识别算法[J]. 分析化学研究简报, 2011, 39(6): 911-914.
- 陈鹏飞, 田地, 乔淑君, 等. 一种基于连续小波变换的LIBS自动寻峰方法[J]. 光谱学与光谱分析, 2014, 34(7): 1969-1972.