CN109508818A - 一种基于LSSVM的在线NOx预测方法 - Google Patents
一种基于LSSVM的在线NOx预测方法 Download PDFInfo
- Publication number
- CN109508818A CN109508818A CN201811237111.9A CN201811237111A CN109508818A CN 109508818 A CN109508818 A CN 109508818A CN 201811237111 A CN201811237111 A CN 201811237111A CN 109508818 A CN109508818 A CN 109508818A
- Authority
- CN
- China
- Prior art keywords
- data
- sample
- error
- online
- prediction model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 239000000523 sample Substances 0.000 claims description 66
- 238000005259 measurement Methods 0.000 claims description 26
- 238000012549 training Methods 0.000 claims description 26
- 239000003245 coal Substances 0.000 claims description 21
- 238000012706 support-vector machine Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 6
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 5
- 239000003546 flue gas Substances 0.000 claims description 5
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 4
- 239000001301 oxygen Substances 0.000 claims description 4
- 229910052760 oxygen Inorganic materials 0.000 claims description 4
- 238000012512 characterization method Methods 0.000 claims description 3
- 238000012886 linear function Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 241001269238 Data Species 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 20
- 238000000738 capillary electrophoresis-mass spectrometry Methods 0.000 description 4
- 238000010531 catalytic reduction reaction Methods 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 235000000332 black box Nutrition 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000000714 time series forecasting Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Evolutionary Computation (AREA)
- Game Theory and Decision Science (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Entrepreneurship & Innovation (AREA)
- Bioinformatics & Computational Biology (AREA)
- Operations Research (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于LSSVM的在线NOx预测方法,包括以下步骤:1、使用历史数据建立基于LSSVM的离线NOx预测模型;2、建立一个滑动数据窗口;3、计算在线更新预测模型中样本的个体预测误差;4、获取新的数据并输入在线更新预测模型和滑动窗口中;5、计算滑动窗口中的总体预测误差与个体预测误差;6、当总体预测误差超过阈值ε,转至步骤7;否则转至步骤4;7、使用滑动窗口中个体预测误差超过阈值的数据顶替在线更新预测模型样本中依次最小的拉格朗日值的原数据;8、重新建立在线更新预测模型,转到步骤3。本发明的优越效果是保证了计算机的稳定运行;降低了计算成本且在线更新预测模型能够同时满足实时性和预测精度的要求。
Description
技术领域
本发明属于发电厂技术领域,尤其涉及一种基于LSSVM的在线NOx预测方法。
背景技术
燃煤发电是NOx的主要来源,目前,世界各国严格控制NOx的排放浓度。特别是中国的NOx排放浓度限制在50mg/m3是美国的一半,因此,对燃煤机组进行NOx监测、控制 NOx排放浓度是十分必要的,也是当务之急。
经检索,目前正在广泛使用的测量设备是烟气自动监控系统(CEMS)存在许多缺点,CEMS需要定期的离线维护,成本相对较高,但这仍不能保证在实际运行中其测量值是始终精确有效。CEMS的另一个缺点是它需要较长的测量时间,因此导致控制系统难以使用滞后的CEMS测量值来控制NOx排放。为了克服上述缺点,近年来,提出了基于数据驱动的间接测量模型。这种方法也被称为软测量,该名称来自于与基于硬件的测量相比,测量方法的差异。特别是在软测量中,基于人工智能技术的模型是目前最有效的模型,例如人工神经网络(ANN)和支持向量机(SVM)。这类模型也被称为黑箱模型,它们不需要关于燃煤锅炉机理的基本知识,但是,这些方法需要相当多的数据来建立和训练模型。在过去,由于缺乏数据,黑箱模型是很难完成的,但随着分布式控制系统(DCS)在火电机组的应用,现在大量的运行数据被记录并存储在数据库中,可用于构建数据驱动黑箱模型。
SVM是一种基于统计学习理论的机器学习算法,其思想是通过非线性映射将输入数据映射到高维空间并进行回归。因此,SVM被广泛用于模式分类和回归估计。但是,SVM算法涉及相对复杂的计算。为了解决这个问题,Suykens在1999年提出了基于SVM的最小二乘支持向量机(LSSVM)。其本质为将二次规划问题转化为求解线性方程组,从而加速问题求解,提高了计算的收敛性。现在,LSSVM已被广泛应用于可靠性分析、时间序列预测的相关领域。在实际应用中,虽然LSSVM能够较好地拟合历史数据,但由于缺乏在线学习能力,泛化能力在很大程度上依赖于初始建模样本选择,同时如何提高支持向量的稀疏性并降低复杂度仍然是一个权衡问题。
发明内容
本发明针对上述现有技术的缺点,提出一种基于LSSVM的在线NOx预测方法,本发明所述方法包括以下步骤:
步骤1、使用历史数据建立基于LSSVM的离线NOx预测模型:
步骤1.1、使用电厂的历史数据并对历史数据进行分析处理,SCR入口NOx即为SCR反应器入口NOx的测量值,辅助变量为:烟气含氧量、总风量、总煤量、机组负荷、磨煤机(A~F) 风煤比、炉膛中心温度;磨煤机(A~F)风煤比为进入磨煤机的冷风和热风与煤的比值,其中 (A~F)则是从A到F的5台磨煤机,SCR(Selective Catalytic Reduction)即为选择性催化还原技术;
步骤1.2、样本优化:使用基于相似度函数的方法,当两个样本之间的欧式距离小的情况下,两个样本将会被认定为是重复样本;反之,则会认定为是非重复样本,当产生的重复的样本就舍弃一个重复的样本并对全部的训练样本进行优化,原始数据中的数据量大且存在冗余数据即重复数据,通过样本优化的方法将冗余数据提出并得到的新数据则能代表之前的原始数据;
选择相似度函数:
其中,exp为以e为底的指数,e是自然常数,xi∈Rn,i=1,2,3,…l,n为样本的特征维数, l为样本的个数,归一化参数是δ,向量二范数为||·||2,第i个样本点和第j个样本点之间的相似度值用Fij表征,R为实域;
针对样本的数据信息,为减小数据采集和分析的误差,选择整体的归一化参数δ,计算公式如下:
其中,Ci作为选择样本中的第i个特征的值,n是选择样本特征的维数;
当两个样本的相似度越接近1,则说明两个样本值越接近,则两个样本会含有多个相同的数据,造成数据冗余,加大了计算量,选择舍弃一组样本数据;
存在训练集其中xk∈Rd,yk∈Rd,d是为辅助变量的数量,支持向量机的从k=1到n的训练集;
通过非线性映射函数将输入空间Rd中的输入样本映射到特征空间φ(x)=(φ1(x),φ2(x),…,φn(x)),如下所示是一个引入误差的估计测量模型函数:
其中,ω∈Rdn,b∈R,e∈R,k=1,2,…,n均为常数或常数矩阵;
根据结果风险最小化原则,并引入最小二乘法,式(1)必须满足
其中,e是误差,γ是正则化参数,使用拉格朗日法,式(2)能够改写为:
其中,αk(k=1,2,…,n)是为拉格朗日因子,ei是误差;
寻找满足式(3)的α和b是LSSVM建模的目标;定义核函数,K(Xk,Xi),核函数是一个符合Mercer约束条件的任意对称函数,α和b是测量模型的系数,则测量模型为:
其中,k=1,2,…,n;α和b是测量模型的系数,
核函数包括径向基函数(RBF)、多项式函数、S函数和线性函数。
通过计算式(5)中4个参数的偏导为0,消除参数ω和e后,得到式(7)的线性方程组如下:
其中:
y=[y1,y2,…,yn];
1=[1,1,…,1];
α=[α1,α2,…,αn];
Ωkl=φ(xk)Tφ(xl);
k=1,2,…,n;
上式(7)中,测量模型的项数为训练样本总数加1;
步骤2、建立一个滑动数据窗口,由下式(9)计算的总预测误差来确定在线更新预测模型参数是否更新;
步骤3、计算在线更新预测模型中样本的个体预测误差,定义预测误差为预测值和实际值之间的偏差:
上式中,yi为i时刻的实际值,为i时刻的预测值,ei为i时刻的预测误差;
步骤4、获取新的数据并输入在线更新预测模型和滑动窗口中,旧数据从滑动窗口删除以使数据量稳定,不会因数据量的堆积从而导致运行速度放缓并节约了处理时间;
步骤5、计算滑动窗口中的总体预测误差与个体预测误差,在线更新预测模型总预测误差为:
其中,上式中,m为在线更新预测模型训练样本数量,总体预测误差用于决定在线更新预测模型是否需要更新;
步骤6、当总体预测误差超过阈值ε,转至步骤7;否则转至步骤4;
步骤7、使用滑动窗口中个体预测误差超过阈值的数据顶替在线更新预测模型样本中依次最小的拉格朗日值的原数据;若总预测误差大于阈值ε需要更新时,插入滑动窗口中误差大于阈值的n组训练样本,删除拉格朗日值即公式(5)中的最小的n组训练样本,更新在线更新预测模型参数,生成新的在线更新预测模型以避免删除训练数据中的重要样本;
步骤8、重新建立在线更新预测模型,转到步骤3。
本发明的优越效果是:
1,本发明所述方法用计算机语言去替代实际生产中的硬件设备,通过容易获取的辅助变量去建立在线更新预测模型从而间接的预测主导变量并能够实现对难于测量或无法测量的变量进行精确测量,通过历史数据建立最初的LSSVM模型,然后建立能设定误差精度的滑动窗口,当误差大于设定值时更新在线更新预测模型,提高了测量的准确性。
2,本发明所述方法的在线更新预测模型更新是由总体误差触发并由个体样本的误差指导数据的选取,在线更新预测模型的数据数量始终不变且在线更新预测模型的矩阵规模不变,保证了计算机的稳定运行;在建模的过程不需要对每组输入的新数据进行计算,降低了计算成本。
3,本发明所述方法在保证预测精度的前提下使在线更新预测模型快速对NOx进行在线软测量以符合火电厂NOx软测量的要求并能够准确有效的在线预测NOx的数值。
4,本发明所述方法能够根据预测误差来抛弃冗余数据,增加新的数据来更新训练样本且在线更新预测模型能够同时满足实时性和预测精度的要求。
附图说明
图1是本发明所述方法的LSSVMSCR入口NOx测量的在线更新预测模型结构示意图;
图2是本发明所述方法的在线更新预测模型数据更新示意图;
图3是本发明所述方法的在线更新预测模型流程图。
具体实施方式
下面结合附图对本发明的实施方式进行详细描述。本发明所述方法包括以下步骤:
步骤1、使用历史数据建立基于LSSVM的离线NOx预测模型:
步骤1.1、使用电厂的历史数据并对历史数据进行分析处理,SCR入口NOx即为SCR反应器入口NOx的测量值,辅助变量为:烟气含氧量、总风量、总煤量、机组负荷、磨煤机(A~F) 风煤比、炉膛中心温度;磨煤机(A~F)风煤比为进入磨煤机的冷风和热风与煤的比值,其中 (A~F)则是从A到F的5台磨煤机,SCR(Selective Catalytic Reduction)即为选择性催化还原技术;
步骤1.2、样本优化:使用基于相似度函数的方法,当两个样本之间的欧式距离小的情况下,两个样本将会被认定为是重复样本;反之,则会认定为是非重复样本,当产生的重复的样本就舍弃一个重复的样本并对全部的训练样本进行优化,原始数据中的数据量大且存在冗余数据即重复数据,通过样本优化的方法将冗余数据提出并得到的新数据则能代表之前的原始数据;
选择相似度函数:
其中,exp为以e为底的指数,e是自然常数,xi∈Rn,i=1,2,3,…l,n为样本的特征维数, l为样本的个数,归一化参数是δ,向量二范数为||·||2,第i个样本点和第j个样本点之间的相似度值用Fij表征,R为实域;
针对样本的数据信息,为减小数据采集和分析的误差,选择整体的归一化参数δ,计算公式如下:
其中,Ci作为选择样本中的第i个特征的值,n是选择样本特征的维数;
当两个样本的相似度越接近1,则说明两个样本值越接近,则两个样本会含有多个相同的数据,造成数据冗余,加大了计算量,选择舍弃一组样本数据;
存在训练集其中xk∈Rd,yk∈Rd,d是为辅助变量的数量,支持向量机的从k=1到n的训练集;
通过非线性映射函数将输入空间Rd中的输入样本映射到特征空间φ(x)=(φ1(x),φ2(x),…,φn(x)),如下所示是一个引入误差的估计测量模型函数:
其中,ω∈Rdn,b∈R,e∈R,k=1,2,…,n均为常数或常数矩阵;
根据结果风险最小化原则,并引入最小二乘法,式(1)必须满足
其中,e是误差,γ是正则化参数,使用拉格朗日法,式(2)能够改写为:
其中,αk(k=1,2,…,n)是为拉格朗日因子,ei是误差;
寻找满足式(3)的α和b是LSSVM建模的目标;定义核函数,K(Xk,Xi),核函数是一个符合Mercer约束条件的任意对称函数,α和b是测量模型的系数,则测量模型为:
其中,k=1,2,…,n;α和b是测量模型的系数,
核函数包括径向基函数(RBF)、多项式函数、S函数和线性函数。
通过计算式(5)中4个参数的偏导为0,消除参数ω和e后,得到式(7)的线性方程组如下:
其中:
y=[y1,y2,…,yn];
1=[1,1,…,1];
α=[α1,α2,…,αn];
Ωkl=φ(xk)Tφ(xl);
k=1,2,…,n;
上式(7)中,测量模型的项数为训练样本总数加1;
步骤2、建立一个滑动数据窗口,由下式(9)计算的总预测误差来确定在线更新预测模型参数是否更新;在线更新预测模型数据更新示意图如图2所示;
步骤3、计算在线更新预测模型中样本的个体预测误差,定义预测误差为预测值和实际值之间的偏差:
上式中,yi为i时刻的实际值,为i时刻的预测值,ei为i时刻的预测误差;
步骤4、获取新的数据并输入在线更新预测模型和滑动窗口中,旧数据从滑动窗口删除以使数据量稳定,不会因数据量的堆积从而导致运行速度放缓并节约了处理时间;
步骤5、计算滑动窗口中的总体预测误差与个体预测误差,在线更新预测模型总预测误差为:
其中,上式中,m为在线更新预测模型训练样本数量,总体预测误差用于决定在线更新预测模型是否需要更新;
步骤6、当总体预测误差超过阈值ε,转至步骤7;否则转至步骤4;
步骤7、使用滑动窗口中个体预测误差超过阈值的数据顶替在线更新预测模型样本中依次最小的拉格朗日值的原数据;若总预测误差大于阈值ε需要更新时,插入滑动窗口中误差大于阈值的n组训练样本,删除拉格朗日值即公式(5)中的最小的n组训练样本,更新在线更新预测模型参数,生成新的在线更新预测模型以避免删除训练数据中的重要样本;
步骤8、重新建立在线更新预测模型,转到步骤3。
如图1所示,首先在电厂中采取辅助变量和主导变量的历史数据,将历史数据进行预处理,这样可以减少冗余数据,加快建模速度,节省建模时间。通过最小二乘支持向量机(LSSVM)的方法建立SCR入口NOx测量的在线更新预测模型。图1中的辅助变量为:烟气含氧量,总风量,总煤量,机组负荷,磨煤机(A~F)风煤比,炉膛中心温度,11个辅助变量。图3为在线更新预测模型的流程图。图2是以时间为横坐标轴的在线更新预测模型更新系统,图2表示的是随时间的变化的系统运行流程,当更新开始时,收集的辅助变量以及 SCR入口NOx的实时数据随着时间的积累从而建立出了滑动窗口,当出现总体误差大于阈值时,则更新在线更新预测模型,此时输入当前样本准确的预测SCR入口NOx值。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明公开的范围内,能够轻易想到的变化或替换,都应涵盖在本发明权利要求的保护范围内。
Claims (5)
1.一种基于LSSVM的在线NOx预测方法,其特征在于,包括以下步骤:
步骤1、使用历史数据建立基于LSSVM的离线NOx预测模型:
步骤1.1、使用电厂的历史数据并对历史数据进行分析处理,SCR入口NOx即为SCR反应器入口NOx的测量值,辅助变量为:烟气含氧量、总风量、总煤量、机组负荷、磨煤机风煤比、炉膛中心温度;磨煤机风煤比为是进入磨煤机的冷风和热风与煤的比值;
步骤1.2、样本优化:使用基于相似度函数的方法,当两个样本之间的欧式距离小的情况下,两个样本将会被认定为是重复样本;反之,则会认定为是非重复样本;当产生的重复的样本就舍弃一个重复的样本并对全部的训练样本进行优化,原始数据中的数据量大且存在很多冗余数据即多余的重复数据,通过样本优化的方法将冗余数据提出并得到的新数据则能代表之前的原始数据;
步骤2、建立一个滑动数据窗口,由式(9)计算的总预测误差来确定在线更新预测模型参数是否更新;
步骤3、计算在线更新预测模型中样本的个体预测误差,定义预测误差为预测值和实际值之间的偏差:
上式中,yi为i时刻的实际值,为i时刻的预测值,ei为i时刻的预测误差;
步骤4、获取新的数据并输入在线更新预测模型和滑动窗口中,旧数据从滑动窗口删除;
步骤5、计算滑动窗口中的总体预测误差与个体预测误差,在线更新预测模型总预测误差为:
其中,上式中,m为在线更新预测模型训练样本数量;总体预测误差用于决定在线更新预测模型是否需要更新;
步骤6、当总体预测误差超过阈值ε,转至步骤7;否则转至步骤4;
步骤7、使用滑动窗口中个体预测误差超过阈值的数据顶替在线更新预测模型样本中依次最小的拉格朗日值的原数据;若总预测误差大于阈值ε需要更新时,插入滑动窗口中误差大于阈值的n组训练样本,删除拉格朗日值即公式5中的最小的n组训练样本,更新在线更新预测模型参数,生成新的在线更新预测模型以避免删除训练数据中的重要样本;
步骤8、重新建立在线更新预测模型,转到步骤3。
2.根据权利要求1所述的一种基于LSSVM的在线NOx预测方法,其特征在于,所述步骤1.2中,选择相似度函数:
exp为以e为底的指数,e是自然常数,
其中,xi∈Rn,i=1,2,3,…l,n为样本的特征维数,l为样本的个数,归一化参数是δ,向量二范数为||·||2,第i个样本点和第j个样本点之间的相似度值用Fij表征,R为实域。
3.根据权利要求1所述的一种基于LSSVM的在线NOx预测方法,其特征在于,所述步骤1.2中,针对样本的数据信息,选择整体的归一化参数δ,计算公式如下:
其中,Ci作为选择样本中的第i个特征的值,n是选择样本特征的维数;
当两个样本的相似度越接近1,则说明两个样本值越接近,则两个样本会含有很多相同的数据,这样就会造成数据冗余,从而加大了计算量,使得计算时间变长,所以选择舍弃一组样本数据;
存在训练集其中xk∈Rd,yk∈Rd,d是为辅助变量的数量,支持向量机的从k=1到n的训练集。
4.根据权利要求1所述的一种基于LSSVM的在线NOx预测方法,其特征在于,所述步骤1.2中,通过非线性映射函数将输入空间Rd中的输入样本映射到特征空间φ(x)=(φ1(x),φ2(x),…,φn(x)),如下所示是一个引入误差的估计测量模型函数:
其中,ω∈Rdn,b∈R,e∈R,k=1,2,…,n均为常数或常数矩阵;
根据结果风险最小化原则,并引入最小二乘法,式(1)必须满足
其中,e是误差,γ是正则化参数,使用拉格朗日法,式(2)能够改写为:
其中αk(k=1,2,…,n)是为拉格朗日因子,ei是误差;
寻找满足式(3)的α和b是LSSVM建模的目标;定义核函数,K(Xk,Xi),核函数是一个符合Mercer约束条件的任意对称函数,α和b是测量模型的系数,则测量模型为:
其中,k=1,2,…,n;α和b是测量模型的系数;
通过计算式(5)中4个参数的偏导为0,消除参数ω和e后,得到式(7)的线性方程组
其中:
y=[y1,y2,…,yn];
1=[1,1,…,1];
α=[α1,α2,…,αn];
Ωkl=φ(xk)Tφ(xl);
k=1,2,…,n;
上式(7)中,测量模型的项数为训练样本总数加1。
5.根据权利要求4所述的一种基于LSSVM的在线NOx预测方法,其特征在于,所述核函数包括径向基函数、多项式函数、S函数和线性函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811237111.9A CN109508818B (zh) | 2018-10-23 | 2018-10-23 | 一种基于LSSVM的在线NOx预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811237111.9A CN109508818B (zh) | 2018-10-23 | 2018-10-23 | 一种基于LSSVM的在线NOx预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109508818A true CN109508818A (zh) | 2019-03-22 |
CN109508818B CN109508818B (zh) | 2022-10-25 |
Family
ID=65745934
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811237111.9A Active CN109508818B (zh) | 2018-10-23 | 2018-10-23 | 一种基于LSSVM的在线NOx预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109508818B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110533152A (zh) * | 2019-08-12 | 2019-12-03 | 哈尔滨工程大学 | 一种基于改进烟花算法在cmp架构下的任务调度方法 |
CN111487950A (zh) * | 2020-04-24 | 2020-08-04 | 西安交通大学 | 在线预警和离线诊断的“预测-验证-反馈-优化”闭环系统 |
CN111540412A (zh) * | 2020-04-24 | 2020-08-14 | 大唐环境产业集团股份有限公司 | 一种基于最小二乘法的scr反应器入口烟气软测量方法 |
CN111832171A (zh) * | 2020-07-13 | 2020-10-27 | 沈阳铁路信号有限责任公司 | 一种基于数学模型的铁路信号继电器性能状态预测方法 |
CN112488145A (zh) * | 2019-11-26 | 2021-03-12 | 大唐环境产业集团股份有限公司 | 一种基于智能化方法的NOx在线预测方法及系统 |
CN112541296A (zh) * | 2020-07-22 | 2021-03-23 | 华北电力大学(保定) | 一种基于pso-lssvm的so2预测方法 |
CN113221439A (zh) * | 2021-04-07 | 2021-08-06 | 浙江大学 | 一种基于bp神经网络的排水系统实时率定与动态预测方法 |
CN116341770A (zh) * | 2023-05-30 | 2023-06-27 | 无锡兴达泡塑新材料股份有限公司 | 基于聚苯乙烯生产数据的生产能力预测方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170061305A1 (en) * | 2015-08-28 | 2017-03-02 | Jiangnan University | Fuzzy curve analysis based soft sensor modeling method using time difference Gaussian process regression |
CN106649919A (zh) * | 2016-09-13 | 2017-05-10 | 上海电气电站环保工程有限公司 | 燃煤电站锅炉飞灰含碳量预测模型构建方法及系统 |
CN106680428A (zh) * | 2016-12-19 | 2017-05-17 | 华北电力大学(保定) | 一种用于脱硝控制系统的软测量方法 |
CN106681381A (zh) * | 2017-01-03 | 2017-05-17 | 华北电力大学 | 一种基于智能前馈信号的scr脱硝系统喷氨量优化控制系统及方法 |
CN108647483A (zh) * | 2018-04-26 | 2018-10-12 | 国家电投集团河南电力有限公司技术信息中心 | 一种基于模糊树建模方法的scr入口nox浓度的软测量方法 |
-
2018
- 2018-10-23 CN CN201811237111.9A patent/CN109508818B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170061305A1 (en) * | 2015-08-28 | 2017-03-02 | Jiangnan University | Fuzzy curve analysis based soft sensor modeling method using time difference Gaussian process regression |
CN106649919A (zh) * | 2016-09-13 | 2017-05-10 | 上海电气电站环保工程有限公司 | 燃煤电站锅炉飞灰含碳量预测模型构建方法及系统 |
CN106680428A (zh) * | 2016-12-19 | 2017-05-17 | 华北电力大学(保定) | 一种用于脱硝控制系统的软测量方法 |
CN106681381A (zh) * | 2017-01-03 | 2017-05-17 | 华北电力大学 | 一种基于智能前馈信号的scr脱硝系统喷氨量优化控制系统及方法 |
CN108647483A (zh) * | 2018-04-26 | 2018-10-12 | 国家电投集团河南电力有限公司技术信息中心 | 一种基于模糊树建模方法的scr入口nox浓度的软测量方法 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110533152A (zh) * | 2019-08-12 | 2019-12-03 | 哈尔滨工程大学 | 一种基于改进烟花算法在cmp架构下的任务调度方法 |
CN112488145A (zh) * | 2019-11-26 | 2021-03-12 | 大唐环境产业集团股份有限公司 | 一种基于智能化方法的NOx在线预测方法及系统 |
CN112488145B (zh) * | 2019-11-26 | 2024-05-28 | 大唐环境产业集团股份有限公司 | 一种基于智能化方法的NOx在线预测方法及系统 |
CN111487950A (zh) * | 2020-04-24 | 2020-08-04 | 西安交通大学 | 在线预警和离线诊断的“预测-验证-反馈-优化”闭环系统 |
CN111540412A (zh) * | 2020-04-24 | 2020-08-14 | 大唐环境产业集团股份有限公司 | 一种基于最小二乘法的scr反应器入口烟气软测量方法 |
CN111540412B (zh) * | 2020-04-24 | 2021-03-09 | 大唐环境产业集团股份有限公司 | 一种基于最小二乘法的scr反应器入口烟气软测量方法 |
CN111487950B (zh) * | 2020-04-24 | 2021-11-16 | 西安交通大学 | 在线预警和离线诊断的“预测-验证-反馈-优化”闭环系统 |
CN111832171A (zh) * | 2020-07-13 | 2020-10-27 | 沈阳铁路信号有限责任公司 | 一种基于数学模型的铁路信号继电器性能状态预测方法 |
CN112541296A (zh) * | 2020-07-22 | 2021-03-23 | 华北电力大学(保定) | 一种基于pso-lssvm的so2预测方法 |
CN113221439A (zh) * | 2021-04-07 | 2021-08-06 | 浙江大学 | 一种基于bp神经网络的排水系统实时率定与动态预测方法 |
CN116341770A (zh) * | 2023-05-30 | 2023-06-27 | 无锡兴达泡塑新材料股份有限公司 | 基于聚苯乙烯生产数据的生产能力预测方法及系统 |
CN116341770B (zh) * | 2023-05-30 | 2023-07-21 | 无锡兴达泡塑新材料股份有限公司 | 基于聚苯乙烯生产数据的生产能力预测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN109508818B (zh) | 2022-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109508818A (zh) | 一种基于LSSVM的在线NOx预测方法 | |
CN109146209A (zh) | 基于遗传算法小波神经网络的机床主轴热误差预测方法 | |
Chuentawat et al. | The comparison of PM2. 5 forecasting methods in the form of multivariate and univariate time series based on support vector machine and genetic algorithm | |
CN102053595A (zh) | 一种乙烯装置中控制裂解炉裂解深度的方法 | |
CN113095550A (zh) | 基于变分递归网络和自注意力机制的空气质量预测方法 | |
CN115034129B (zh) | 一种面向火电厂脱硝装置的NOx排放浓度软测量方法 | |
CN112801388B (zh) | 一种基于非线性时间序列算法的电力负荷预测方法及系统 | |
Qiao et al. | Application of PSO LS-SVM forecasting model in oil and gas production forecast | |
Girimurugan et al. | Application of deep learning to the prediction of solar irradiance through missing data | |
CN114418166A (zh) | 一种变压器油中的溶解气体浓度预测方法、装置及介质 | |
CN110738363A (zh) | 一种光伏发电功率预测模型及其构建方法和应用 | |
CN109960146A (zh) | 提高软测量仪表模型预测精度的方法 | |
Zheng et al. | TCN-GAT multivariate load forecasting model based on SHAP value selection strategy in integrated energy system | |
CN117350897A (zh) | 一种考虑空气湿度的转换功率预测的方法和系统 | |
CN105787507B (zh) | 基于预算支持向量集的LS‑SVMs在线学习方法 | |
CN116880201A (zh) | 基于模糊神经网络的水网渠道状态控制系统 | |
CN116861256A (zh) | 一种固废焚烧过程的炉温预测方法、系统、设备及介质 | |
Fuselero et al. | Daily Solar Radiation Forecasting based on a Hybrid NARXGRU Network in Dumaguete, Philippines. | |
CN116757078A (zh) | 一种基于作用力的煤粉流速测量方法及系统 | |
JP2009070200A (ja) | プラント運用適正化装置、プラント運用適正化方法およびプラント運用適正化プログラム | |
Kowalski et al. | The comparison of linear models for PM10 and PM2. 5 forecasting | |
Jianli et al. | Wind power forecasting by using artificial neural networks and Grubbs criterion | |
Zhou et al. | Double-tank liquid level control based on genetic algorithm | |
CN113111588A (zh) | 一种燃气轮机nox排放浓度预测方法及装置 | |
Li et al. | Long short-term memory based on random forest-recursive feature eliminated for hot metal silcion content prediction of blast furnace |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |