[发明专利]激光诱导击穿光谱的小波降噪的最佳分解层数选择方法

说明书

技术领域

本发明涉及激光诱导击穿光谱的信号预处理分析领域，具体是一种基于噪声能量分布的数学模型，通过小波能量熵计算得到噪声能量值的预测区间并用于确定激光诱导击穿光谱的小波降噪的最佳分解层数。

背景技术

激光诱导击穿光谱(LIBS)是一种常用的激光烧蚀光谱分析技术。激光经过透镜聚焦到样品上就会在样品局部产生等离子体。用光谱仪直接收集样品表面等离子体产生的发射谱线信号，根据发射光谱的强度对样品进行元素分析。LIBS发射光谱的强度受到多种噪声源因素的影响，比如光电耦合器的暗电流、电子电路的热噪声以及光谱仪的杂散光和等离子体的连续辐射等。这些噪声往往与LIBS发射光谱信号混杂在一起，不利于后续的样品元素分析。

为了解决LIBS光谱信号的降噪问题，通常采用的方法有两类：硬件优化方法和软件方法。硬件优化方法通常包括设计滤波器以及隔离器来降噪，但是受到开发周期较长和耗费多的影响，造成了硬件优化方法使用起来受到了较大的局限性。软件方法不但具备硬件优化方法不可比拟的可靠性和精度，而且还具有硬件优化方法无法达到的性能。软件滤波降噪的方法有很多，其中以具有出色的时频局部特性的小波降噪方法为最常用的方法。小波降噪方法主要包括三种：系数相关性方法、模极大值方法和阈值方法。这三种小波滤波方法的降噪效果均受到小波分解层数的影响。分解层数过多会造成信号有用信息的丢失，信噪比反而下降且增加运算量；分解层数过少则降噪效果不理想，信噪比提高有限。常用的确定小波最佳分解层数的方法是白噪声检测方法，它基于小波系数的相关性系数白化检验来确定小波降噪的最佳分解层数。但是，这种方法在较大的分解层上会产生自相关性较好的子信号，与真实信号混淆，造成了小波最佳分解层数的选择出现偏差影响了小波降噪效果。

发明内容

为了解决现有小波降噪中选择最佳分解层数的困难所引起的降噪效果不理想的不足，本发明的目的在于提出了一种基于噪声能量分布的数学模型，通过小波能量熵计算得到噪声能量值的预测区间并用于确定激光诱导击穿光谱的小波降噪的最佳分解层数。

本发明为实现上述目的所采用的技术方案是：一种激光诱导击穿光谱的小波降噪的最佳分解层数选择方法，包括以下步骤：

步骤1)输入整个激光诱导击穿光谱信号Y，并对激光诱导击穿光谱信号Y进行标准化计算得到标准化后的激光诱导击穿光谱信号y；

步骤2)利用蒙特卡洛方法产生与步骤1)中激光诱导击穿光谱信号Y长度一致、统计特性一致的高斯白噪声N；

步骤3)利用Shannon熵来选择小波基函数并计算得到最大的分解层数M；

步骤4)根据小波基函数对标准化后的激光诱导击穿光谱信号y进行离散小波变换，计算得到标准化后的激光诱导击穿光谱信号y在每一个分解层上的细节系数cD_j；根据小波基函数对高斯白噪声N进行离散小波变换，计算得到高斯白噪声N在每一个分解层上的细节系数cD’_j，j＝1,…,M；

步骤5)根据标准化后的激光诱导击穿光谱信号y在每一个分解层上的细节系数cD_j重构标准化后的激光诱导击穿光谱信号y在每一层的高频子信号S_j；根据高斯白噪声N在每一个分解层上的细节系数cD'_j重构高斯白噪声N在每一层的高频子信号S'_j，j＝1,…,M；

步骤6)根据标准化后的激光诱导击穿光谱信号y在每一层的高频子信号S_j，j＝1,…,M，计算小波能量熵WEE_j及其差分导数D(WEE_j)，j＝2,…,M；根据高斯白噪声N在每一层的高频子信号S'_j，j＝1,…,M，计算小波能量熵WEE′_j，j＝1,…,M，及其差分导数D(WEE′_j)，j＝2,…,M；

步骤7)根据D(WEE_j)在置信水平为0.95下，计算预测区间；

步骤8)如果某一D(WEE′_j)超出所述预测区间，就得到了最佳分解层数j-1。

所述的标准化后的激光诱导击穿光谱信号y的计算公式是其中，μ是Y的平均值，σ是Y的标准方差。

所述统计特性是激光诱导击穿光谱信号Y的平均值和标准方差。

所述的最大分解层数M的计算公式是M＝[log₂(L)]，其中，L是激光诱导击穿光谱信号y的长度，[.]是取整函数。