[发明专利]基于历史数据稳态值的动态系统增益估计方法有效
申请号: | 201811638956.9 | 申请日: | 2018-12-29 |
公开(公告)号: | CN109753634B | 公开(公告)日: | 2020-08-11 |
发明(设计)人: | 王建东;杨子江;曹鹏飞;周东华 | 申请(专利权)人: | 山东科技大学 |
主分类号: | G06F17/18 | 分类号: | G06F17/18 |
代理公司: | 济南圣达知识产权代理有限公司 37221 | 代理人: | 李琳 |
地址: | 266590 山东*** | 国省代码: | 山东;37 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 基于 历史数据 稳态 动态 系统 增益 估计 方法 | ||
1.一种基于历史数据稳态值的动态系统增益估计方法,其特征是:包括以下步骤:
(1)采集机组的有功功率、主蒸汽流量的控制器输出和主蒸汽压力的历史数据样本;采用线性分段表示方法找到机组处于稳定状态的有功功率、主蒸汽流量的控制器输出和主蒸汽压力的数据段:
具体地,对历史样本数据中输入和输出的长连续时间序列进行分割,形成数据段,引入最小振幅变化和最小样本数,作为阈值参数表征各数据段的稳态条件,保留处于输入和输出均在稳态下的数据段;
(2)根据数据样本的平均值来计算保留数据段的稳态值;
(3)在稳态值中,确定在统计意义上与相同稳态增益相关联的元素,将这些元素分成一组并分别估计每组的稳态增益;
(4)计算所述稳态增益估计参数在某个置信度下的区间估计;
(5)得到不同的稳态增益对应不同工况,估计出在不同工况下的多个稳态增益,确定动态系统机组此时可能处于的运行点;
所述步骤(3)中,具体过程包括:
步骤3.1:通过求解多元线性方程组,从稳态值中估计参数;
步骤3.2:找到稳态值中与其估计值最大偏差的稳态值;
步骤3.3:若该稳态值与对应的估计值的差值大于等于最小振幅变化,将对应的稳态值删除;
步骤3.4:重复步骤3.1-3.3,直到稳态值与对应的估计值的差值不大于最小振幅变化为止,表示选择出来的稳态值和相应的估计参数;
步骤3.5:对剩余的稳态值进行重复步骤3.1-3.4的操作,直到剩余的稳态值的个数小于设定数值。
2.如权利要求1所述的一种基于历史数据稳态值的动态系统增益估计方法,其特征是:所述步骤(1)中,采用自底向上的分段线性表示方法从历史数据样本中分别将输入和输出的长连续时间序列分割成多个短数据段,每个短数据段用直线表示。
3.如权利要求1所述的一种基于历史数据稳态值的动态系统增益估计方法,其特征是:所述步骤(1)中,当某数据段同时满足其振幅变化值小于等于设定的最小振幅变化,且该数据段的样本总数大于等于设置的最小样本数时,认为所述数据段处于稳定状态。
4.如权利要求1所述的一种基于历史数据稳态值的动态系统增益估计方法,其特征是:所述步骤(1)中,构造整体指示序列,整体指示序列为各输入数据段稳态指示值与输出数据段稳态指示值乘积,计算得到整体指示序列值,其中,当某数据段处于稳定状态,对应的数据段稳态指示值为1,否则为0。
5.如权利要求4所述的一种基于历史数据稳态值的动态系统增益估计方法,其特征是:所述步骤(2)中,确定稳态条件下的整体指示序列的起止位置,确定起止位置之间的数据段的样本平均值确定相应的稳态值。
6.如权利要求1所述的一种基于历史数据稳态值的动态系统增益估计方法,其特征是:所述步骤(4)中,根据中心极限定理,相应的稳态值的估计参数收敛于高斯分布。
7.如权利要求1所述的一种基于历史数据稳态值的动态系统增益估计方法,其特征是:所述步骤(4)中,计算估计参数的均值和协方差矩阵,根据所述均值和协方差矩阵,给出设定置信度值的所有估计参数的区间估计。
8.如权利要求1或7所述的一种基于历史数据稳态值的动态系统增益估计方法,其特征是:所述步骤(4)中,设定置信度值为95%。
9.一种动态系统稳态增益估计系统,其特征是:运行于处理器或存储器上,被配置为执行以下指令:
(1)接收采集机组的有功功率、主蒸汽流量的控制器输出和主蒸汽压力的历史数据样本,采用线性分段表示方法找到机组处于稳定状态的有功功率、主蒸汽流量的控制器输出和主蒸汽压力的数据段;
对历史样本数据中输入和输出的长连续时间序列进行分割,形成数据段,引入最小振幅变化和最小样本数,作为阈值参数表征各数据段的稳态条件,保留处于输入和输出均在稳态下的数据段;
(2)根据数据样本的平均值来计算保留数据段的稳态值;
(3)在稳态值中,确定在统计意义上与相同稳态增益相关联的元素,将这些元素分成一组并分别估计每组的稳态增益;
(4)计算所述稳态增益估计参数在某个置信度下的区间估计;
(5)得到不同的稳态增益对应不同工况,估计出在不同工况下的多个稳态增益,确定动态系统机组此时可能处于的运行点;
所述步骤(3)中,具体过程包括:
步骤3.1:通过求解多元线性方程组,从稳态值中估计参数;
步骤3.2:找到稳态值中与其估计值最大偏差的稳态值;
步骤3.3:若该稳态值与对应的估计值的差值大于等于最小振幅变化,将对应的稳态值删除;
步骤3.4:重复步骤3.1-3.3,直到稳态值与对应的估计值的差值不大于最小振幅变化为止,表示选择出来的稳态值和相应的估计参数;
步骤3.5:对剩余的稳态值进行重复步骤3.1-3.4的操作,直到剩余的稳态值的个数小于设定数值。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于山东科技大学,未经山东科技大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201811638956.9/1.html,转载请声明来源钻瓜专利网。