基于迭代更新的ukfnn铝电解功耗模型构建方法
【专利摘要】本发明公开了一种基于迭代更新的UKFNN铝电解功耗模型构建方法，包括如下步骤：对铝电解功耗模型的结构参数进行初始化；计算所述功耗模型的一组采样点；根据所述功耗模型上一时刻的状态估计值，由状态方程F对所述功耗模型下一时刻的状态和协方差进行预测；重新采样以预测状态值为中心、以预测方差为协方差所产生的样点，用观测函数对采样点和协方差进行更新，然后进行状态变量和观测值的协方差更新；用更新后的下一时刻的状态和下一时刻的协方差矩阵重新进行状态估计更新。所述方法将迭代方法与UKFNN相结合，从而得到了更加准确的估计值。
【专利说明】
基于迭代更新的UKFNN错电解功耗模型构建方法
技术领域
[0001] 本发明设及特别适用于特定功能的数字计算设备或数据处理设备或数据处理方 法技术领域，尤其设及一种基于迭代更新的UKFNN侣电解功耗模型构建方法。
【背景技术】
[0002] 经过多年的快速发展，中国已逐步发展成为一个侣冶金大国。根据侣冶金能源成 本的统计，侣±矿的开采超过原侣成本的=分之一，生产吨侣消费量是生产吨钢能耗的4.5 倍，平均能耗约为182-212mj/kg，电解过程大约占到64%，节能的要求日渐强烈，侣电解的 节能技术的应用在我国侣冶金工业的可持续发展中占据着重要的地位。
[0003] 为了解决侣冶金工业的能耗问题，实现侣工业的可持续发展，W及节能减排。仍然 需要采取一系列的措施，可W从两方面进行改进：
[0004] a.硬件方面的改进，如加大阳极的长度、增强电解槽的电流密度或者加大阴极炭 块的尺寸等方法。
[0005] b.软件方面的改进，如采用优化建模、对决策参数的调整等，对侣电解能耗的核屯、 算法进行精确的改造优化，使侣电解的能耗逐步的减少。
[0006] 我国侣冶金行业能耗比较大，针对如何选择生产过程中的工艺参数、实现工艺过 程的最优控制，建立精确的侣电解能耗模型是实现侣电解最低能耗的依据，也同时是国内 外正在积极探索的问题。传统的非线性系统建模是利用BP神经网络和无迹卡尔曼神经网络 法建模。
[0007] BP神经网络模型具有较强的非线性映射能力、高度自学习和自适应的能力，鉴于 BP神经网络的运些优点，B巧巾经网络模型应用比较广泛，与此同时也暴露出了越来越多的 不足，如局部极小化、收敛速度慢、神经网络结构选择不一等。
[000引无迹卡尔曼神经网络建模法具有应用范围广和数学建模简单易行，不需计算雅克 比矩阵，它可估计信号的过去和当前状态，甚至能估计将来的状态，当新的数据被观测到 后，只要根据新的数据和前一时刻的估计量，即可算出新的估计量。运种方法在非线性估计 工程中应用比较广泛，运种方法主要是对神经网络的权值和阔值进行动态调整，修改，它可 W根据生产条件的实时变化而建立动态演化模型，。但UKRW仍然存在卡尔曼滤波器自身存 在的缺陷，UKFr^N在协方差矩阵的转移过程中，因为不稳定的计算而造成状态协方差矩阵失 去对称性，使UKFr^N算法失效;使用线性最小均方估计测量更新的精度不高。
【发明内容】
[0009]本发明所要解决的技术问题是提供一种基于迭代更新的UKF順侣电解功耗模型构 建方法，所述方法将迭代方法与UKFNN相结合，利用迭代更新的方法提高估计精度，即求得 UKF順算法中状态估计值后，再将估计值返回测量更新阶段进行重新采样计算，从而得到更 加准确的估计值。
[0010]为解决上述技术问题，本发明所采取的技术方案是:一种基于迭代更新的UKF順侣 电解功耗模型构建方法，其特征在于包括如下步骤：
[0011] 1)对侣电解功耗模型的结构参数进行初始化
其中，XO表示侣电解功耗模型权值和阔值的初始状态，Il是对XO的估计，C。表示XO的方差； [001^ 2)计算所述功耗模型的一组采样点，求取2L+1个点，L为状态X日的维数；
[0013] 3)根据所述功耗模型上一时刻的状态估计值，由状态方程F对所述功耗模型下一 时刻的状态和协方差进行预测；
[0014] 4)重新采样W预测值为中屯、、W预测方差为协方差所产生的样点，用观测函数对 采样点和协方差进行更新，然后进行状态变量和观测值的协方差更新，求卡尔曼增益，状态 更新得到下一时刻状态，协方差矩阵更新得到下一时刻协方差矩阵；用更新后的下一时刻 的状态和下一时刻的协方差矩阵重新进行状态估计更新。
[0015] 进一步的技术方案在于:所述步骤1)中初始化状态变量XO包括电解侣槽电压、系 列电流、电解质水平、分子比、侣水平、出侣量、槽溫、效应间隔和下料间隔。
[0016] 进一步的技术方案在于:所述步骤2)中：
[0017] 计算所述功耗模型的一组采样点，求取2L+1个点，L为状态XO的维数；
[001引其4
[0019] 站-1是所述功耗模型K-I时刻的采样值，Illl是所述功耗模型K-I时刻的估计值，L 是状态变量的维数，A是缩放比例参数，Pk-I是所述功耗模型K-I时刻的协方差矩阵。
[0020] 进一步的技术方案在于:所述步骤3)中：
[0021 ]根据上一时刻状态由状态方程F对下一时刻进行状态预测
[0022]
[0023] 其中：xk|k-i表示根据K-I时刻的状态对下一时刻的状态进行预测，f是状态转移函 数，W是各个采样点对应的权重，Il表示各个采样点加权后的更新值；
[0024]
[00巧]&表示即寸刻的协方差矩阵，Q是过程噪声的协方差。
[0026] 进一步的技术方案在于:所述的步骤4)中
[0027] 测量更新：
[002引 当j = 0时
[0029]
[0030] 其中Xj表示重新采样后的采样矩阵，Il表示第j次各个采样点加权后的更新值， 重新采样W预测值为中屯、、W预测方差为协方差所产生的样点对观测值和协方差进行更 新；
[0031]
[0032] ' 个测量值
[0033]
[0034]
[0035]
[0036]
[0037]
[00；3 引
[0039]
[0040]
[0041]
[0042]
[0043] 协方差矩阵更新，得到下一时刻协方差矩阵，Pk是更新后的协方差矩阵。
[0044] 采用上述技术方案所产生的有益效果在于:所述方法将迭代与UKF順相结合，利用 迭代更新的方法提高估计精度，即求得UKF順算法中状态估计值后，再将估计值返回测量更 新阶段进行重新采样计算，从而得到更加准确的估计值。
【附图说明】
[0045] 下面结合附图和【具体实施方式】对本发明作进一步详细的说明。
[0046] 图1是本发明所述功耗模型构建方法的流程图；
[0047] 图2是50组基于迭代更新的UKF順神经网络预测图；
[004引图3是50组基于迭代更新的UKF順神经网络预测误差图；
[0049]图4是80组基于迭代更新的UKF順神经网络预测图；
[0化0]图5是80组基于迭代更新的UKRW神经网络预测误差图；
[0化1 ]图6是110组基于迭代更新的UKF順神经网络预测图；
[0化2]图7是110组基于迭代更新的UKRW神经网络预测误差图；
[0化3]图8是200组基于迭代更新的UKRW神经网络预测图；
[0化4] 图9是200组基于迭代更新的UKRW神经网络预测误差图。
【具体实施方式】
[0055] 下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整 地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0056] 在下面的描述中阐述了很多具体细节W便于充分理解本发明，但是本发明还可W 采用其他不同于在此描述的其它方式来实施，本领域技术人员可W在不违背本发明内涵的 情况下做类似推广，因此本发明不受下面公开的具体实施例的限制。
[0057] 如图1所示，本发明公开了一种基于迭代更新的UKF順侣电解功耗模型构建方法， 其特征在于包括如下步骤：
[0化引1)对侣电解功耗模型的结构参数进行初始本
其中，X质示侣电解功耗模型权值和阔值的初始状态，；6是对XO的估计，巧。表示XO的方差；
[0059] 其中初始化状态变量XO可W包括电解侣槽电压、系列电流、电解质水平、分子比、 侣水平、出侣量、槽溫、效应间隔和下料间隔，当然还可W为其它本领域技术人员需要的参 数。
[0060] 2)计算所述功耗模型的一组采样点，求取2L+1个点，L为状态XO的维数；
[OOW] 其中
[006^ 站-1是所述功耗模型K-I时刻的采样值III是所述功耗模型K-I时刻的估计值，L是 状态变量的维数，A是缩放比例参数，Pk-I是所述功耗模型K-I时刻的协方差矩阵。
[0063] 3)根据所述功耗模型上一时刻的状态估计值，由状态方程F对所述功耗模型下一 时刻的状态和协方差进行预测；
[0064] 根据上一时刻状态由状态方程F对下一时刻进行状态预测
[00 化]
[0066] 其中：xk|k-i表示根据K-I时刻的状态对下一时刻的状态进行预测，f是状态转移函 数，W是各个采样点对应的权重，.?表示各个采样点加权后的更新值；
[0067]
[006引 Il表示K时刻的协方差矩阵，Q是过程噪声的协方差。
[0069] 4)重新采样W预测值为中屯、、W预测方差为协方差所产生的样点，用观测函数对 采样点和协方差进行更新，然后进行状态变量和观测值的协方差更新，求卡尔曼增益，状态 更新得到下一时刻状态，协方差矩阵更新得到下一时刻协方差矩阵；用更新后的下一时刻 的状态和下一时刻的协方差矩阵重新进行状态估计更新。
[0070] 所述的步骤4)中测量更新的过程如下：
[0071] 当 J = O时
[0072]
[0073] 新采样I
[0074]
[0075] 个测量'
[0076]
[0077]
[007引
[0079]
[0080]
[0081]
[0082]
[0083]
[0084]
[0085]
[0086] 协方差矩阵史新，得到下一时刻协方差矩阵，Pk是史新后的协方差矩阵。仿具与分 析：
[0087] 经过实验分析，本方法隐层节点采用8个，根据IUKF順中的状态维数计算可知，维 数为8*9+8+8+1 = 89。当训练样本数据为50组、80组、110组、200组时IUKFr^N模型预测效果如 下分析：
[0088] a.当训练样本数据为50JUKFNN模型预测输出和误差如图2-3所示。
[0089] b.当训练样本数据为80JUKFN饰莫型预测输出和误差如图4-5所示。
[0090] C.当训练样本数据为110，IUKFN饰莫型预测输出和误差如图6-7所示。
[0091] d.当训练样本数据为200，IUKFN饰莫型预测输出和误差如图8-9所示。
[0092] 实验结果表明，随着训练样本的增加，IUKF順神经网络的误差逐渐减小。在实时生 产过程中，由于其能进行动态调整，测量更新部分能不断的更新新的状态估计值和方差重 新开始采样，解决了 BP順、UKFN饰莫型精度差的问题，几乎能够达到精准的预测效果。
[0093] 所述方法将迭代与UKF顺相结合，利用迭代更新的方法提高估计精度，即求得 UKF順算法中状态估计值后，再将估计值返回测量更新阶段进行重新采样计算，从而得到更 加准确的估计值。
【主权项】
1. 一种基于迭代更新的UKFNN铝电解功耗模型构建方法，其特征在于包括如下步骤： 1) 对铝电解功耗模型的结构参数进行初始化：^ = -元]h - ，其 中，χο表示铝电解功耗模型权值和阈值的初始状态，Il是对χο的估计，圍表示χο的方差； 2) 计算所述功耗模型的一组采样点，求取2L+1个点，L为状态XO的维数； 3) 根据所述功耗模型上一时刻的状态估计值，由状态方程F对所述功耗模型下一时刻 的状态和协方差进行预测； 4) 重新采样以预测状态值为中心、以预测方差为协方差所产生的样点，用观测函数对 采样点和协方差进行更新，然后进行状态变量和观测值的协方差更新，求卡尔曼增益，状态 更新得到下一时刻状态，协方差矩阵更新得到下一时刻协方差矩阵；用更新后的下一时刻 的状态和下一时刻的协方差矩阵重新进行状态估计更新。2. 如权利要求1所述的基于迭代更新的UKFNN铝电解功耗模型构建方法，其特征在于： 所述步骤1)中初始化状态变量XQ包括电解铝槽电压、系列电流、电解质水平、分子比、铝水 平、出铝量、槽温、效应间隔和下料间隔。3. 如权利要求1所述的基于迭代更新的UKFNN铝电解功耗模型构建方法，其特征在于， 所述步骤2)中： 计算所述功耗模型的一组采样点，求取2L+1个点，L为状态XO的维数； 其中Xk-^所述功耗模型κ-1时刻的采样值HI是所述功耗模型K-I时刻的估计值，L是状态 变量的维数，λ是缩放比例参数，？!^是所述功耗模型K-I时刻的协方差矩阵。4. 如权利要求1所述的基于迭代更新的UKFNN铝电解功耗模型构建方法，其特征在于， 所述步骤3)中： 根据上一时刻状态由状态方程F对下一时刻进行状态预测其中：处^:表示根据K-I时刻的状态对下一时刻的状态进行预测，f是状态转移函数，W 是各个米样点对应的权重，Jp表不各个米样点加权后的更新值；III表不K时刻的协方差矩阵，Q是过程噪声的协方差。5. 如权利要求1所述的基于迭代更新的UKFNN铝电解功耗模型构建方法，其特征在于， 所述的步骤4)中 测量更新： 当j = 〇时其中^表示重新采样后的采样矩阵，表示第j次各个采样点加权后的更新值，重新 采样以预测值为中心、以预测方差为协方差所产生的样点对观测值和协方差进行更新；H是观测函数，^表示对各个采样点进行计算的测量值构成的采样矩阵，Il表示各个测 量值加权后的Ir新倌。协方差矩阵更新，得到下一时刻协方差矩阵，Pk是更新后的协方差矩阵。
【文档编号】G06F17/50GK106021698SQ201610325327
【公开日】2016年10月12日
【申请日】2016年5月17日
【发明人】姚立忠, 张恒健, 李太福, 苏盈盈, 易军, 黄迪, 曹旭鹏
【申请人】重庆科技学院