含有基线漂移和高频噪声的光电传感信号重构方法及系统与流程

本发明涉及信号处理技术领域，特别是涉及一种含有基线漂移和高频噪声的光电传感信号重构方法及系统。

背景技术：

本部分的陈述仅仅是提供了与本发明相关的背景技术信息，不必然构成在先技术。

数据中心机柜等重要小空间设施独特，功能专业，一旦发生火灾，若扑救不及时会造成严重的信息中断和经济社会影响，这类场所亟需在火灾萌芽期就能预警的火灾探测预警系统。国内外均已对吸气式感烟探测技术开展了研究并形成了相关产品，但是在火灾“极早期”探测预警方面基本处于空白。

火灾“极早期”是指物质从被过度加热超过其材质热分解点，到氧化燃烧并开始产生碳烟的阶段，尚无大量的烟雾粒子产生；因此需要对微量的碳烟粒子进行检测，而由此转化的电信号非常微弱，极易掺杂在噪声中无法判断。

理想状况下，在没有烟雾颗粒时，光电接收器的电压信号应该是一直不变的值，但是实际应用中发现，此电压值会随着温度变化、光线衍射、干涉的缓慢变化等产生低频的缓慢漂移，如果漂移过大，会造成误报警，因此需要合适的算法将基线漂移滤除。

此外，硅光电二极管作为光电导器件的原理是当红外入射光照射在pn结上时，光能被吸收然后转化为电能，在入射光传输到硅光电二极管的过程中，不可避免的会伴随着一些随机性的噪声，主要包括器件内部噪声和环境噪声；其中，环境噪声是光电探测器接收信号过程中出现的噪声，其干扰源来自器件外部，例如杂散光的入射、电磁干扰、光路传输介质的湍流等，这些噪声可以通过屏蔽、遮断杂光、选择滤光片等方法进行改善，但是在实际应用中发现，即使应用暗盒，这些干扰仍然无法完全消除；而内部噪声主要有散粒噪声和热噪声，均与频率无关，属于白噪声，上述所有噪声构成原信号的背景噪声。

因此，光电探测器的输出信号并不是一个无掺杂的纯净信号，而是带有低频漂移和白噪声的信号，低频漂移会造成误报警，白噪声会降低系统的检测精度，都会干扰火灾“极早期”阶段的预警，故数据中心机柜“极早期”烟雾探测中存在基线漂移幅值大和背景噪声大的问题。

技术实现要素：

为了解决上述问题，本发明提出了一种含有基线漂移和高频噪声的光电传感信号重构方法及系统，在原始光电传感信号中加入高斯噪声后进行eemd分解，根据imf分量阶数设计自适应低通滤波器，对含有基线漂移的imf分量进行处理，对剩余的imf分量处理高频噪声，得到无漂移、无噪声的重构信号，实现准确的极早期预警。

为了实现上述目的，本发明采用如下技术方案：

第一方面，本发明提供一种含有基线漂移和高频噪声的光电传感信号重构方法，包括：

对光电传感信号进行eemd分解后得到imf分量；

对存在基线漂移的imf分量，根据其阶数选择自适应滤波器以获取基线漂移估计值，对剩余的imf分量获取高频噪声估计值；

对光电传感信号滤除前述两类估计值后，得到重构信号。

第二方面，本发明提供一种含有基线漂移和高频噪声的光电传感信号重构系统，包括：

分解模块，用于对光电传感信号进行eemd分解后得到imf分量；

计算模块，用于对存在基线漂移的imf分量，根据其阶数选择自适应滤波器以获取基线漂移估计值，对剩余的imf分量获取高频噪声估计值；

重构模块，用于对光电传感信号滤除前述两类估计值后，得到重构信号。

第三方面，本发明提供一种电子设备，包括存储器和处理器以及存储在存储器上并在处理器上运行的计算机指令，所述计算机指令被处理器运行时，完成第一方面所述的方法。

第四方面，本发明提供一种计算机可读存储介质，用于存储计算机指令，所述计算机指令被处理器执行时，完成第一方面所述的方法。

与现有技术相比，本发明的有益效果为：

集合经验模态分解eemd算法是处理信号噪声比较常用的方法，目前较多的应用于图像分析处理、心电信号ecg去噪等方面，本发明将其应用于烟雾探测方面，本发明提出一种改进的eemd分解方法，在原始信号中加入高斯噪声后进行分解，同时引入一组自适应低通滤波器，根据imf分量阶数设计自适应低通滤波器，依次对含有基线漂移的imf分量进行处理，对剩余的imf分量处理高频噪声的影响，实现滤除基线漂移、降低背景噪声、提高信噪比的目的，实现准确的极早期预警。

本发明附加方面的优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。

附图说明

构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。

图1为本发明实施例1提供的含有基线漂移和高频噪声的光电传感信号重构方法流程图；

图2(a)-2(n)为本发明实施例1提供的eemd分解结果示意图；

图3(a)-3(b)为本发明实施例1提供的原始信号与分解后的信号对比示意图。

具体实施方式：

下面结合附图与实施例对本发明做进一步说明。

应该指出，以下详细说明都是示例性的，旨在对本发明提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。

需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。

在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。

实施例1

如图1所示，本实施例提供一种含有基线漂移和高频噪声的光电传感信号重构方法，包括：

s1：对光电传感信号进行eemd分解后得到imf分量；

s2：对存在基线漂移的imf分量，根据其阶数选择自适应滤波器以获取基线漂移估计值，对剩余的imf分量获取高频噪声估计值；

s3：对光电传感信号滤除前述两类估计值后，得到重构信号。

在本实施例中，为了方便且准确的计算信噪比，在原始光电传感信号中加入高斯噪声，对光电传感信号进行eemd分解，分解成n级imf分量和残留分量。

在本实施例中，根据光电传感信号基线漂移和高频白噪声的特点，即基线漂移分布在最后几阶imf分量中，噪声信号主要分布在前几阶imf分量中，而删除最后几阶imf分量可能会造成重构信号的失真，故本实施例提出eemd自适应滤波，引入一组自适应低通滤波器依次处理最后几阶imf分量，获取基线漂移信号估计值，以此滤除基线漂移信号。

在本实施例中，所述步骤s2中，对于imf分量，选择与imf阶数自适应的低通滤波器处理基线漂移信号，将处理过后的各阶imf分量进行重构，获取基线漂移估计值和高频噪声估计值；具体包括：

对原始光电传感信号加入白噪声后得到信号x(t)，进行eemd分解处理，x(t)分解为多个imf分量，设为xk(t)：

其中，n是imf分量的数量；

对于存在基线漂移的m阶imf分量，通过低通滤波器hk(t)将xk(t)由高阶处理为低阶，该低通滤波器hk(t)的输出为：

yk(t)＝hk(t)*xk(t)

其中，*表示卷积，hk(t)是低通滤波器。

随着imf阶数的降低，每个imf分量中包含的慢振荡分量越来越少，而信号分量越来越多；因此，设计一组截止频率ωk随imf阶数的减小而减小的自适应低通滤波器，具体为：首先从最后一阶imf分量开始，即xn(t)分量的自适应滤波器截止频率为ωn(t)，因此，第k个滤波器的截止频率为：

ωk-1＝ωk*α

其中，α是0.1到0.99之间的数，k＝n,…,2,1。

根据imf各阶数分量的特性，基线幅值随着阶数的减少而逐渐减小，因此，从每个imf分量中提取出包含低频成分的滤波器输出yk(t)，计算yk的标准差std(yk)和评价系数pk作为停止准则，以确定m值：

其中，k＝1,2,…n，flip代表反转函数，设置基线漂移的阈值系数ε，范围是0-0.1，如果pk＜ε，则截止阶数m＝n+1-k；

在本实施例中，设置ωn为10，α为0.9，ε为0.01，将m+1～n阶imf滤波器的和作为重构基线估计。

采用阈值函数法，设置高频噪声的阈值系数为δ，采用δ与剩余imf分量中的表示高频噪声的imf分量相乘得到高频噪声信号估计，因此基线漂移和高频噪声的估计值总和为：

最后，从原始光电传感信号中减去基线漂移和高频噪声估计值总和，得到无漂移、无噪声的重构信号。

本实施例模拟实际火灾过程中“极早期”烟雾状态，将设备开机后预热足够时间使其充分稳定，在正常运行环境中不间断的取6000个数据点；同时，为了验证所用算法的有效性，计算信噪比，在获取6000个光电信号数据的过程中，通过标准烟箱三次加入等量的烟雾颗粒模拟实际火灾过程中“极早期”的烟雾状态。eemd的处理结果如图2(a)-2(n)所示，原信号与处理过后的信号对比如图3(a)-3(b)所示。

采用“相关系数”计算原始信号与能反映低频漂移趋势的分解量的关系，相关系数体现变量之间线性相关程度的量，相关系数的绝对值一般在0.8以上认为a和b有强相关性；在0.3到0.8之间认为有弱相关性；在0.3以下认为没有相关性。

为了计算高频分量与原信号基线漂移趋势的贴近程度，去除有效烟雾信号的干扰，首先在原始数据中手动去除有效烟雾信号数据，只留下无烟雾信号的背景数据，采用matlab计算改进的eemd方法的低频漂移分量和原始信号的关系，结果为：

r_eemd＝0.9925

采用信噪比的方法进行对比，信噪比反映信号与噪声能量的关系，一般来说，信噪比高且绝对平均误差低的方法比信噪比低且绝对平均误差高的方法好；信噪比与绝对平均误差的公式为：

snr＝10lg{∑x(t)²/∑[x(t)-x(t)']²}

采用matlab计算可得利用eemd方法处理数据后的信噪比为：

snr_eemd＝22.3986

在本实施例中，采用改进的eemd方法，分解计入噪声信号的原始信号，信号被分解为12级imf分量，基线漂移分布主要在最后一级imf分量中，与原始信号的低频漂移趋势非常近似；选取合适的自适应低通滤波器对最后一级imf进行处理，滤波器输出的结果视为基线漂移估计；高频噪声信号主要分布在前几级imf分量中，选取合适的阈值对前几级imf分量进行处理，获取预估的纯噪声信号；最后通过从信号中减去基线漂移估计值和纯噪声信号估计值，获得无漂移、无噪声的重构信号，更好的实现了滤除基线漂移、降低背景噪声、提高信噪比的目的。

实施例2

本实施例提供一种含有基线漂移和高频噪声的光电传感信号重构系统，包括：

分解模块，用于对光电传感信号进行eemd分解后得到imf分量；

计算模块，用于对存在基线漂移的imf分量，根据其阶数选择自适应滤波器以获取基线漂移估计值，对剩余的imf分量获取高频噪声估计值；

重构模块，用于对光电传感信号滤除前述两类估计值后，得到重构信号。

此处需要说明的是，上述模块对应于实施例1中的步骤s1至s3，上述模块与对应的步骤所实现的示例和应用场景相同，但不限于上述实施例1所公开的内容。需要说明的是，上述模块作为系统的一部分可以在诸如一组计算机可执行指令的计算机系统中执行。

在更多实施例中，还提供：

一种电子设备，包括存储器和处理器以及存储在存储器上并在处理器上运行的计算机指令，所述计算机指令被处理器运行时，完成实施例1中所述的方法。为了简洁，在此不再赘述。

应理解，本实施例中，处理器可以是中央处理单元cpu，处理器还可以是其他通用处理器、数字信号处理器dsp、专用集成电路asic，现成可编程门阵列fpga或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。

存储器可以包括只读存储器和随机存取存储器，并向处理器提供指令和数据、存储器的一部分还可以包括非易失性随机存储器。例如，存储器还可以存储设备类型的信息。

一种计算机可读存储介质，用于存储计算机指令，所述计算机指令被处理器执行时，完成实施例1中所述的方法。

实施例1中的方法可以直接体现为硬件处理器执行完成，或者用处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器、闪存、只读存储器、可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器，处理器读取存储器中的信息，结合其硬件完成上述方法的步骤。为避免重复，这里不再详细描述。

本领域普通技术人员可以意识到，结合本实施例描述的各示例的单元即算法步骤，能够以电子硬件或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本申请的范围。

以上仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。