一种钢铁生产过程数据的降噪方法

1.本发明涉及钢铁生产过程的自动控制技术领域，尤其涉及一种钢铁生产过程数据的降噪方法。


背景技术：

2.在现代钢铁生产过程中，钢铁企业每日产生并存储大量生产过程数据，这些数据蕴含了与钢铁生产设备、钢铁生产过程相关的规律性知识。以数据为核心，应用智能分析技术对钢铁生产过程中关键质量指标与工艺参数间进行关联规则挖掘，为钢铁生产过程中质量参数控制提供了新思路。数据是构建模型的基石，只有准确有效的数据才能确保算法的精度，这就要求数据本身具有良好的可用性。随着钢铁生产过程控制技术的发展以及工业大数据应用的普及，往往采用实际生产数据对钢铁生产过程中的质量参数进行分析建模，会使结果更加符合实际情况，避免了传统方法的局限性。但是钢铁生产现场采集到的数据中往往混杂着大量噪声，噪声会很大程度地影响算法的计算效果，对控制量的求解以及质量参数的预测过程造成了困难。因此，结合数据降噪技术，对钢铁生产过程数据进行降噪处理，对提升钢铁生产过程关键质量参数的控制水平具有重要意义。
3.fourier变换是数据降噪中常用的一种频域处理方法，它可以将给定的信号进行展开，分解为多个不同频率的正弦信号的和。19世纪70年代，beltrami提出了一种针对于实正方矩阵的奇异值分解(svd)理论，该理论通过正交变换的手段将信号分解为众多线性分量，以svd 理论为基础的众多分析方法在特征提取和信号降噪领域得到了广泛的应用。1998年，huang 等人以瞬时频率为基础提出了一种经验模态分解(emd)方法，并以此为基础，提出了希尔伯特-黄变换(hht)。同时，为改善emd方法中存在的模态混叠问题，wu等人利用高斯白噪声频率均匀的特点，向原始信号中人为增添了高斯白噪声，提出了集成经验模态分解法 eemd。
4.近些年，各类数据降噪方法也在不断发展，有的方法会针对自身算法缺陷进行改进，也有的方法会对于整体的计算过程进行优化。对于一些实际生产过程而言，现场状况更加复杂多变，数据呈现复杂的非线性，单一的降噪手段不能达到预期效果，因此多种降噪手段联合的方法也逐渐兴起。为了提高钢铁生产过程中质量参数控制水平，解决噪声对控制模型建立的干扰，需要提出一种针对钢铁生产过程的合理、有效的数据降噪手段。


技术实现要素：

5.针对现有技术的不足，本发明提出了一种钢铁生产过程数据的降噪方法。
6.一种钢铁生产过程数据的降噪方法，具体为以下步骤：
7.步骤1：采集含噪声的钢铁生产过程原始数据，建立eemd-wt降噪模型；
8.步骤1.1：对含噪声的钢铁生产过程原始数据进行集成经验模态分解eemd，得到频率从高至低的若干组噪声主导和信号主导的各本征模态函数imf分量；
9.10.式中wi(t)为第i次分解选取的高斯白噪声，xi(t)为第i次加入高斯白噪声后的数据，x(t)为含噪声的钢铁生产过程原始数据，ri为余项，n为eemd分解的imf数目，imf
ij
为第i次加入高斯白噪声时分解出的第j个imf分量；
11.步骤1.1.1：向钢铁生产过程原始数据中加入正态分布的高斯白噪声；
12.步骤1.1.2：提取出加入高斯白噪声后数据xi(t)的所有局部极值，包括局部最大值以及局部最小值；
13.步骤1.1.3：采用三次样条插值的方法，分别对局部最大值和局部最小值进行拟合，得到上、下包络线；
14.步骤1.1.4：求上、下包络线的平均值，并将其定义为均值包络线m1(t)；
15.步骤1.1.5：定义h1(t)＝xi(t)-m1(t)，判断h1(t)是否为imf；如果h1(t)是imf，则一阶imf 分量为c1(t)＝h1(t)；如果h1(t)不是imf，则将其视为加入高斯白噪声的数据，并重步骤 1.1.2-1.1.4，直到标准差sd小于设定值停止；
[0016][0017]
式中t为每个时刻点，t为总时刻，hk(t)为k-1阶含噪声的数据与k阶imf包络线的差值；
[0018]
步骤1.1.6：剔除一阶imf分量，得到剩余的新数据r1(t)＝xi(t)-c1(t)，对r1(t)重复进行步骤1.1.2-1.1.5，得到二阶imf分量c2(t)，以此类推，直到第n阶分量小于预设值；
[0019]
步骤1.1.7：重复m次操作步骤1.1.1到步骤1.1.6，每次重新选取新的高斯白噪声序列，得到m组imf分量；
[0020]
步骤1.1.8：将m组imf分量进行集成均值处理，得到eemd分解后的最终结果：
[0021][0022]
式中m为总体平均次数；
[0023]
得到频率从高至低的若干组噪声主导和信号主导的各本征模态函数imf分量；
[0024]
步骤1.2：采用连续均方误差的方法cmse计算噪声能量突变点b；
[0025]
信号主导的分量和噪声信号主导的分量的分界点b为噪声能量突变点：
[0026][0027][0028]
式中n为数据样本量，c为本征模态函数总数，cmse(imfk，imf
k+1
)为k阶imf与k+1 阶imf的连续均方误差，imfk(ti)为第k阶本征模态分量中第i个本征模态分量， argfirstlocalmin为cmse(imfk，imf
k+1
)最小时，求imf阶数k的值；
[0029]
步骤1.3：计算出噪声能量突变点b值后，将噪声主导的高频imf分量进行小波变换wt 降噪处理，除去高频中的噪声，保留高频数据段的剩余信息；
[0030]
步骤1.4：最后高频数据段的剩余信息与信号主导的低频imf分量重构，完成eemd-wt 降噪模型的建立；
[0031]
步骤2：把含噪声的钢铁生产过程数据导入eemd-wt降噪模型，得到利用eemd-wt 降噪模型降噪后的数据。
[0032]
本发明的有益效果：本发明提供一种钢铁生产过程数据的降噪方法，涉及钢铁生产过程的自动控制技术领域。
[0033]
本发明提出的集成经验模态分解结合小波变换eemd-wt的降噪方法降噪效果好，相比于 svd降噪方法和均值降噪方法能更好地还原出原始数据，可以广泛地投入到钢铁生产过程当中。
附图说明
[0034]
图1为本发明具体实施时含噪原始数据经eemd分解后的imfs图。
[0035]
图2为本发明具体实施时利用小波变换对imf1～imf4进行降噪的结果图。
[0036]
图3为本发明集成经验模态分解结合小波变换(eemd-wt)降噪模型的算法流程图。
具体实施方式
[0037]
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。
[0038]
为验证降噪模型有效性，以钢铁生产的轧制过程为例，将轧制现场中带有噪声的弯辊力数据作为对象进行降噪性能测试。
[0039]
步骤1：采集含噪声的钢铁生产过程原始数据，建立eemd-wt降噪模型；
[0040]
步骤1.1：对含噪声的钢铁生产过程原始数据进行集成经验模态分解eemd，得到频率从高至低的若干组各本征模态函数imf分量；分解后结果如图1所示，本实施例把数据分解为 10个频率从高至低排序的本征模态函数；
[0041][0042]
式中wi(t)为第i次分解选取的高斯白噪声，xi(t)为第i次加入高斯白噪声后的数据，x(t)为含噪声的钢铁生产过程原始数据，ri为余项，n为eemd分解的imf数目，imf
ij
为第i次加入高斯白噪声时分解出的第j个imf分量；
[0043]
步骤1.1.1：向钢铁生产过程原始数据中加入正态分布的高斯白噪声；
[0044]
步骤1.1.2：提取出加入高斯白噪声后数据xi(t)的所有局部极值，包括局部最大值以及局部最小值；
[0045]
步骤1.1.3：采用三次样条插值的方法，分别对局部最大值和局部最小值进行拟合，得到上、下包络线；
[0046]
步骤1.1.4：求上、下包络线的平均值，并将其定义为均值包络线m1(t)；
[0047]
步骤1.1.5：定义h1(t)＝xi(t)-m1(t)，判断h1(t)是否为imf；如果h1(t)是imf，则一阶imf 分量为c1(t)＝h1(t)；如果h1(t)不是imf，则将其视为加入高斯白噪声的数据，并重步骤 1.1.2-1.1.4，直到标准差sd小于设定值停止；
[0048][0049]
式中t为每个时刻点，t为总时刻，hk(t)为k-1阶含噪声的数据与k阶imf包络线的差值；
[0050]
步骤1.1.6：剔除一阶imf分量，得到剩余的新数据r1(t)＝xi(t)-c1(t)，对r1(t)重
复进行步骤1.1.2-1.1.5，得到二阶imf分量c2(t)，以此类推，直到第n阶分量小于预设值；
[0051]
步骤1.1.7：重复m次操作步骤1.1.1到步骤1.1.6，每次重新选取新的高斯白噪声序列，得到m组imf分量；
[0052]
步骤1.1.8：将m组imf分量进行集成均值处理，得到eemd分解后的最终结果：
[0053][0054]
式中m为总体平均次数；
[0055]
得到频率从高至低的若干组噪声主导主导和信号主导的各本征模态函数imf分量；
[0056]
步骤1.2：采用连续均方误差的方法cmse计算噪声能量突变点b；
[0057]
信号主导的分量和噪声信号主导的分量的分界点b为噪声能量突变点：
[0058][0059][0060]
式中n为数据样本量，c为本征模态函数总数，cmse(imfk，imf
k+1
)为k阶imf与k+1 阶imf的连续均方误差，imfk(ti)为第k阶本征模态分量中第i个本征模态分量， argfirstlocalmin为cmse(imfk，imf
k+1
)最小时，求imf阶数k的值；计算各相邻imf 分量间的连续均方差cmse结果如表1所示，cmse(imf3，imf4)为出现的第一个极小值，因此判定b＝4，需要对前四个本征模态分量采取进一步的降噪处理；
[0061]
表1相邻imf的连续均方误差；
[0062][0063]
步骤1.3：计算出噪声能量突变点b值后，将噪声主导的高频imf分量进行小波变换wt 降噪处理，除去高频中的噪声，保留高频数据段的剩余信息；
[0064]
将imf1～imf4进行小波变换降噪处理，得到降噪后的高频数据a1～a4，如图2所示；
[0065]
步骤1.4：最后高频数据段的剩余信息与信号主导的低频imf分量重构，完成eemd-wt 降噪模型的建立；
[0066]
将小波变换处理后的数据a1～a4与imf5～imf10进行重构，得到eemd-wt降噪后的数据；eemd-wt降噪模型的算法流程如图3所示；
[0067]
步骤2：把含噪声的钢铁生产过程数据导入eemd-wt降噪模型，得到利用eemd-wt 降噪模型降噪后的数据。
[0068]
采用常见的svd降噪方法和工程领域常用的均值降噪方法与本发明的eemd-wt降噪模型进行比较。并根据信噪比(snr)来衡量三种方法的降噪性能，其计算公式为：
[0069]
snr(db)＝10log
10
(p
signal
/p
noise
)＝20log
10
(a
signal
/a
noise
)
[0070]
式中，psignal为信号功率；pnoise为噪声功率；asignal为信号幅度；anoise为噪声幅度；
[0071]
表2三种降噪方法信噪比对比结果。
[0072][0073]
eemd-wt降噪方法的信噪比为24.94db，优于svd降噪方法的22.19db和均值降噪方法20.10db，更能说明eemd-wt降噪方法的准确性。

技术特征：
1.一种钢铁生产过程数据的降噪方法，其特征在于，包括以下步骤：步骤1：采集含噪声的钢铁生产过程原始数据，建立eemd-wt降噪模型；步骤1.1：对含噪声的钢铁生产过程原始数据进行集成经验模态分解eemd，得到频率从高至低的若干组噪声主导和信号主导的各本征模态函数imf分量；步骤1.2：采用连续均方误差的方法cmse计算噪声能量突变点b；步骤1.3：计算出噪声能量突变点b值后，将噪声主导的高频imf分量进行小波变换wt降噪处理，除去高频中的噪声，保留高频信号段的剩余信息；步骤1.4：最后高频信号段的剩余信息与信号主导的低频imf分量重构，完成eemd-wt降噪模型的建立；步骤2：把含噪声的钢铁生产过程数据导入eemd-wt降噪模型，得到利用eemd-wt降噪模型降噪后的数据。2.根据权利要求1所述的一种钢铁生产过程数据的降噪方法，其特征在于，所述的步骤1.1对含噪声的钢铁生产过程原始数据进行集成经验模态分解eemd，得到频率从高至低的若干组噪声主导和信号主导的各本征模态函数imf分量；具体为：式中w
i
(t)为第i次分解选取的高斯白噪声，x
i
(t)为第i次加入高斯白噪声后的数据，x(t)为含噪声的钢铁生产过程原始数据，r
i
为余项，n为eemd分解的imf数目，imf
ij
为第i次加入高斯白噪声时分解出的第j个imf分量。3.根据权利要求1所述的一种钢铁生产过程数据的降噪方法，其特征在于，所述的步骤1.1具体为：步骤1.1.1：向钢铁生产过程原始数据中加入正态分布的高斯白噪声；步骤1.1.2：提取出加入高斯白噪声后数据x
i
(t)的所有局部极值；步骤1.1.3：采用三次样条插值的方法，分别对局部最大值和局部最小值进行拟合，得到上、下包络线；步骤1.1.4：求上、下包络线的平均值，并将其定义为均值包络线m1(t)；步骤1.1.5：定义h1(t)＝x
i
(t)-m1(t)，判断h1(t)是否为imf；步骤1.1.6：剔除一阶imf分量，得到剩余的新数据r1(t)＝x
i
(t)-c1(t)，对r1(t)重复进行步骤1.1.2-1.1.5，得到二阶imf分量c2(t)以此类推，直到第n阶分量小于预设值；步骤1.1.7：重复m次操作步骤1.1.1到步骤1.1.6；步骤1.1.8：将m组imf分量进行集成均值处理，得到eemd分解后的最终结果：式中m为总体平均次数。4.根据权利要求3所述的一种钢铁生产过程数据的降噪方法，其特征在于，所述的步骤1.1.2中的局部极值具体为：所述局部极值包括局部最大值以及局部最小值。5.根据权利要求3所述的一种钢铁生产过程数据的降噪方法，其特征在于，所述的步骤1.1.5中的判断h1(t)是否为imf具体为：
如果h1(t)是imf，则一阶imf分量为c1(t)＝h1(t)；如果h1(t)不是imf，则将其视为加入高斯白噪声的信号，并重步骤1.1.2-1.1.4，直到标准差sd小于设定值停止；式中t为每个时刻点，t为总时刻，h
k
(t)为k-1阶含噪声的信号与k阶imf包络线的差值。6.根据权利要求3所述的一种钢铁生产过程数据的降噪方法，其特征在于，所述的步骤1.1.7重复m次操作步骤1.1.1到步骤1.1.6；具体为：每次重新选取新的高斯白噪声序列，得到m组imf分量。7.根据权利要求1所述的一种钢铁生产过程数据的降噪方法，其特征在于，所述的步骤1.2采用连续均方误差的方法cmse计算噪声能量突变点b；具体为：信号主导的分量和噪声信号主导的分量的分界点b为噪声能量突变点：信号主导的分量和噪声信号主导的分量的分界点b为噪声能量突变点：式中n为数据样本量，c为本征模态函数总数，cmse(imf
k
，imf
k+1
)为k阶imf与k+1阶imf的连续均方误差，imf
k
(t
i
)为第k阶本征模态分量中第i个本征模态分量，arg firstlocalmin为cmse(imf
k
，imf
k+1
)最小时，求imf阶数k的值。

技术总结
本发明提供一种钢铁生产过程数据的降噪方法，涉及钢铁生产过程的自动控制技术领域。首先建立EEMD-WT降噪模型，对含有噪声的钢铁生产过程数据进行EEMD分解，然后采用连续均方误差的方法(CMSE)计算噪声能量突变点b，计算出b值后，将前面相应的高频IMF分量进行小波变换降噪处理，除去高频中的噪声，保留高频信号段的剩余信息，最后与低频其余的IMF分量重构，完成EEMD-WT降噪模型的建立；将含噪的钢铁生产过程数据导入到EEMD-WT降噪模型中，得到最后的降噪后数据。本发明提出的EEMD-WT降噪方法降噪效果好，相比于SVD降噪方法和均值降噪方法能更好地还原出原始数据，可以广泛地投入到钢铁生产过程当中。到钢铁生产过程当中。


技术研发人员：孙杰 单鹏飞 丁肇印 李梦琴 乔继柱 李树 刘云霄 李霄剑 彭文 张殿华
受保护的技术使用者：东北大学
技术研发日：2021.12.08
技术公布日：2022/3/8