一种基于不完整时间序列数据集的冶金企业蒸汽流量区间预测方法
专利汇可以提供一种基于不完整时间序列数据集的冶金企业蒸汽流量区间预测方法专利检索，专利查询，专利分析的服务。并且一种基于不完整数据集的 冶金 企业 蒸汽 流量区间预测方法。针对含缺失点的蒸汽流量时间序列数据，利用相空间重构技术构造输入输出训练样本集，并建立一个初始相关向量机回归模型。然后通过贝叶斯方法推导出缺失输出的后验概率分布；并由缺失输出的后验概率的均值来填补与输出相对应的缺失输入点，不断 迭代 上述推导与填补过程来训练 预测模型 。通过单步迭代预测方法来预测未来一段趋势的预测值，其预测区间边界由预测分布的均值和方差决定。这种预测方法可以很好地应用在含有缺失点时间序列数据集中，其预测结果更好地辅助调度人员进行冶金蒸汽 能源 系统的平衡优化工作。，下面是一种基于不完整时间序列数据集的冶金企业蒸汽流量区间预测方法专利的具体信息内容。
1.一种基于不完整时间序列数据集的冶金企业蒸汽流量区间预测方法，其特征在于如下步骤：
(1)数据准备与初始填补
从工业现场实时关系数据库提取一段时间的冶金企业蒸汽流量时间序列数据，检测并记录时间序列中缺失点的位置，并用零值填补所有缺失点；
(2)基于相空间重构的输入输出样本集构建
假设有一条含缺失点的离散时间序列{u(1),u(2),...,u(k),...,u(Nu)}，其中u(k)表示第k个采样时刻的时间序列值，Nu是序列的总长度；由式(1)构造第i个相矢量si(p)：
si(p)＝[u(i),u(i+τ),u(i+2τ),...,u(i+(p-1)τ)]T (1)
其中，i＝1,2,...,N，而N＝Nu-(p-1)τ/Δt，τ表示延迟时间，Δt表示数据的采样间隔，p是相矢量的嵌入维数，也是矢量si(p)的长度；第i个相矢量si(p)同时表示输入样本集中的T
第i个输入样本，即xi＝si(p)，有xi＝[xi1,xi2,...,xic,...,xip] ，其中xic是向量xi中的第c个分量，c＝1,2,…,p；向量xi对应的目标值ti为式(2)所示：
ti＝u(i+p×τ) (2)
此外，N个输入样本的矩阵形式为X＝[x1,x2,...,xi,...,xN]T；N个目标值所组成的向量为t＝[t1,t2,...,ti,...,tN]T；对于由含缺失点的时间序列构造的输入输出样本集，事先记录下输入输出中缺失点的对应关系；
(3)构建初始相关向量机回归模型
假设第i个目标值ti与第i个输入向量xi有如下的映射关系：
其中，j＝0,1,2,…,N，xj是第j个输入向量，且规定K(xi,x0)＝1；另外，第i个基函数向量φi表示为φi＝[1,K(xi,x1),K(xi,x2),...,K(xi,xj),...,K(xi,xN)]T，w＝[w0,w1,...,wj,...,wN]T表示权值向量，K(xi,xj)＝exp{-||xi-xj||2/2l2}是高斯核函数，其中，l≠0是核
2 -2
函数参数；εi表示第i个均值为0，方差为σ的高斯分布噪声项，β≡σ ；对w中的每一个元素分别附加均值为0而方差不同的正态先验分布：
其中，参数α＝[α0,α1,...,αj,...,αN]T，αj是第j个权值wj的先验方差的倒数，且有对角阵A＝diag(α)；并且同时将参数α和β视为随机变量，并假设它们服从伽马分布；(4)基于不完整数据集的相关向量机回归模型构建
构建最终预测模型就是要通过输入输出样本集来迭代计算如下变量：
权值向量w；
缺失输出向量t[m]，即输出集中所有缺失点组成的向量；
缺失输入向量x[m]，即输入集中所有缺失点组成的向量；
噪声方差参数β的点估计值；
具体步骤如下：
2
①初始化：向量α中的元素初始化为值1/N ，β[obs]初始化为10倍的t[obs]的方差，β[obs]是输出集中可观测点的噪声方差的倒数，有β[obs]≡σ-2；其中，t[obs]表示输出集中所有可观测值组成的向量；相关向量集初始化为所有输入样本的集合，此时相关向量集中有N个样本；
②分别使用式(5)和(6)更新w和t[m]的联合后验概率的均值μ[w,m|obs]和协方差矩阵Σ[w,m|obs]；
其中，β＝[β[m],β[m],...,β[m],β[obs],...,β[obs],β[obs]]T，β还可写成 写成对角矩阵的形式是Β[m]＝diag(β[m])，Β[obs]＝diag(β[obs])和Β＝diag(β)；矩阵Β[m]和Β[obs]大小分别为N[yobs]×N[yobs]和N[m]×N[m]，而N[yobs]和N[m]分别表示输出集中可观测点和
8
缺失点的个数，有N＝N[yobs]+N[m]，将β[m]设置为不小于10的值；另外，Φ表示核函数矩阵，有Φ＝[φ1,φ2,...,φi,...,φN]T，并且将Φ表示成分块矩阵的形式为
其中 和 是由缺失输出点和输出可观测点分别在输出样本集中的位置索引所提取的核矩阵Φ中以行为单位的行向量所组成的各自的新矩阵，大小分别为N[m]×(N+1)和N[yobs]×(N+1)；
③采用(7)和(8)式更新αj和β[obs]；
其中， 且 Tr(·)表示矩阵
的求迹运算；
④如果αj＞M，M设置为大于1010的值，将第j个样本从相关向量集中除去；
⑤通过(9)式来计算x[m]中每一个元素的后验概率均值，并填补到输入集中输入缺失点的相应位置，并进一步更新核矩阵Φ；
其中， x[m]集合中共有 个输入缺失点，x[m],v表示x[m]中的第v个元素，
表示 的后验概率均值；而t[m],corr_v表示第v个缺失输入点在缺失输出集t[m]中的对应变量值；t[m],corr_v的后验概率p(t[m],corr_v|t[obs],α,β[obs])由式(5)和(6)得出；
⑥由(10)式计算对数边缘似然度函数的值；
其中，协方差矩阵 矩阵大小为N[yobs]×N[yobs]，I是单位矩阵，
log(·)表示自然对数运算；
⑦如果α中有元素相对上次迭代的变化或边缘似然度函数值的变化小于给定的阈值，终止迭代过程，得出最终预测模型；否则，转到第②步；
(5)单步迭代预测得出未来一段时间的预测区间
给定一个新的输入样本x*，此输入向量中没有缺失点，得出对应的输出t*的预测分布近似为一个正态分布
其中，
μ[w|obs]是权值向量w的后验概率的均值， 由式(8)计算得到；此时，预测上下限分别表示为其中，Ulimit和Llimit分别表示置信区间为(1-χ％)的预测区间上下限，同时 是标准高斯分布中对应概率为 的分位数；未来一段时间的预测均值由单步迭代的方式给出；每一次预测输出的预测均值是一个点，且下一轮的预测会将上一轮的预测结果作为输入，进行循环滚动预测。
一种基于不完整时间序列数据集的冶金企业蒸汽流量区间预
测方法
技术领域
背景技术
发明内容
具体实施方式
此外，所有缺失的和观测到的输入点组成的集合分别表示为向量x[m]和
x[obs](需要注意的是，矩阵 和 与向量x[m]和x[obs]不同，它们中的元素既包含缺失点也有可观测点)。
的联合后验概率分布也是一个多维高斯分布的形式，如下所示。
表示x[m],v的后验概率均值；而t[m],corr_v表示第v个缺失输入点在缺失输出集t[m]中的对应变量；t[m],corr_v的后验概率p(t[m],corr_v|t[obs],α,β[obs])由式(8)和(9)得出。
0；而Ub和Lb分别是第b个点的预测上限和下限。PICP表示被预测区间覆盖的真实目标值的个数，NMPIW是归一化的预测区间平均宽度，其中，G是目标值的取值范围。CWC指标值越小区间预测效果越好。并且当PICP值越大且NMPIW值越小时，区间预测的效果越好。表1三种方法区间预测结果比较(NNI表示最近邻点填补方法；SI表示样条插值填补方法)
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
