电力负荷预测第八灰色系统预测方法.pptx

教学要求：清楚灰色系统、生成数等基本概念；掌握灰色关联度的计算；掌握GM(1,1)的建模及预测；了解GM(1,1)残差模型的作用及其实现；了解其它灰色预测模型的用途；教学重点：关联度计算，GM(1,1)预测。教学难点：GM(1,1)残差模型。一.灰色系统理论概述 80年代，中国华中理工大学，邓聚龙教授创建的；概率统计：要求大样本，事先需知道分布规律；时间序列：数据的拟合；灰色理论：少量数据（4个以上）发现规律；是系统理论的新分支 系统理论工程与实践（一级刊物）；The Journal of Grey System（SCI）；作为横断学科，广泛应用于社会经济系统，控制系 统等的预测、决策、控制；1.灰色系统的概念部分信息已知、部分信息未知的系统，。相对于白色系统，其系统内部特征完全已知或系统信息是充足的。相对于黑色系统，其系统内部信息一无所知，只能从它与外部的联系来观测。举例：电力供求系统 计划体制下白色系统.电价确定；电量需求由用电计划指标决定；市场条件下灰色系统.电价不确定；电量需求受生产经营状况等因素影响等；2.2.灰色预测的分类灰色预测的分类时间序列预测（重点介绍）用等时距观测到的，反映对象特征的一系列数据，构造出灰色预测模型，并预测未来某一时刻的特征量，或者是达到某一特征量所需要的时间，。畸变预测 通过灰色模型预测异常值出现的时刻，预测异常值什么时候出现在特定的时区，。波形预测（拓扑预测）通过灰色模型，预测对于未来变动的轨迹，。系统预测 对系统行为特征指标，建立一组相互关联的灰色预测模型，在预测系统整体变化的同时，预测系统各个环节的变化，。3.生成数特点：不找概率分布；不寻求统计特征；处理方法：累加生成；累减生成；均值生成；级比生成；与概率与数理统计中各随机变量不同。减少波动性。还原用。序列中空穴数的插值。序列两头空穴数的补齐。计算过程：原始序列 生成序列 其中：原始序列生成序列序号例题 例1:令 累加生成2.28+2.98=5.26 例2：令 累减生成(1IAGO)5.26-2.28=2.98二.灰色关联度的分析与计算 分析系统中各因素关联程度。分析系统中各因素关联程度。计算过程：step1:计算关联系数;设 多个参考序列被比较序列则关联系数 注意：若单位不一，初值不同的序列，需事先进行初始化，即将序列中所有的数据分别除以第一个数据。Step2:计算关联度 表示被比较数列与参考数列间的关联度；为各关联系数的平均值。灰色关联分析的主要优势另一类量化分析方法数理统计类 （回归分析，方差分析，主成分分析）样本量大；样本具有较好的分布规律和确定的发展趋势；计算量大；可能出现量化结果与定性分析结果不符的现象；灰色关联分析灰色关联分析 依据因素间发展态势的几何相似或相异程度，来衡量因素间关联程度的。样本多少没有过多要求；不需要典型的分布规律；计算量少；不会出现关联量化结果与定性分析不一致的情况；算例 已知：参考序列 被比较数列 求：关联度。yo解:Step 1:初始化xoStep2:求绝对差序列Step3 求关联系数与k有关,k=16.Step4:求关联度 结论：y1和y0的关联程度大于y2和y0的关联程度。三.GM(1,1)模型及预测 GM(1,1)（Grey Model）1阶1元微分方程；不是差分方程；是一个线性动态模型，常用于时间序列。1.理论依据把随机量当作是在一定范围内变化的灰色量；把随机过程当作是在一定范围、一定时区内变化的灰色过程；灰色理论将无规律的历史数据，经累加生成后，使其变为具有指数增长规律的上升形状数列，由于一阶微分方程解的形式，即是指数增长形式，所以可对生成后的数列，建立微分方程模型。故灰色模型实际上是生成数列所建的模型；灰色理论通过灰数的不同生成方式，数据的不同取舍，不同级别的残差GM模型，来调整，修正，提高精度；对高阶系统建模，灰色理论是通过GM(1,n)模型群解决的。GM(1,n)模型群也即一阶微分方程组组成的灰色模型；GM模型所得到的数据，必须经过逆生成，即累减生成还原后才能应用；2.GM(1,1)模型的建立 原始时间序列 累加生成序列 生成序列X1的微分方程为：其中，a发展灰数；u内生控制灰数；（灰色作用量）待估参数-a反映了 的发展态势。一般情况下，系统作用量应是外生的或前定的，而GM(1,1)是单序列建模，只用到系统行为序列（或称输出序列，背景值），而无外作用序列（或称输入序列，驱动值）。b为从背景值挖掘出来的数据，反映了数据变化的关系，其确切内涵是灰的。灰色作用量是内涵外延化的具体体现，它的存在，是区别灰色建模于一般输入输出建模（黑箱建模）的分水岭。运用最小二乘法（推导参看牛东晓p172)，估计参数其中：生成序列的值原始序列的值原微分方程的解为：变量GM(1,1)模型的预测方程。原微分方程的还原解为：GM(1,1)模型的预测方程的还原解。几点说明给定原始序列X0中的数据，不一定要全部用来建模，对原始序列的取舍不同，模型会不同，即a，u会不同。建模的数据取舍，应保证建模序列等时距、相连的，不得有跳跃出现。一般建模数据序列，应当由最新数据及其相邻数 据构成。当出现新数据时，处理方法：将新信息加入原始序列，重估参数。去掉原始序列中最老的一个数据，再加工最新数据 （等维信息），所形成序列和原序列维数相同，重估 参数。3.GM(1,1)模型的检验 残差检验 关联度检验 后验差检验 残差检验按预测模型，计算 ;将 累减生成 ;计算原始序列 与 的绝对残差 及相对误差 关联度检验：按关联度计算方法，算出 与原始序列 的关联 系数，再算出关联度。依经验，当 时，若关联度大于0.6,模型满意。后验差检验计算原始序列的平均值计算原始序列的均方差计算残差 的均值 计算残差的均方差计算方差比计算小误差概率p 若记 查表若相对误差、关联度、后验差检验均在允许范围内，则可用所建模型进行预测，否则应进行残差修正。P C 模型评价 0.95 0.80 0.70 0.65 勉强合格 0.70 0.65 不合格4.例题已知：试用GM(1,1)建模，并预测第8期的值。拐点解:构造累加生成序列 构造数据矩阵B和数据向量Yn 计算 ，和 。即 预测模型 指数函数形式残差检验因为 故i1234562.673.133.253.263.563.722.675.809.0512.4115.9719.692.675.789.0312.43015.9719.682.673.113.253.403.543.71GM(1,1)模型计算一览表故绝对误差 相对误差小于4%关联度检验：计算 与 的绝对误差 。注：由于只有2个序列，故不必寻求二级最大或最小。关联系数为 故 关联度为 ，模型满足要求。后验差检验X0的均值X0的均方差残差的均值残差的均方差 计算方差比计算小概率 因为所有 都小于 ，故P=1预测作业8-1.对以下数据进行GM(1,1)建模，并进行模型检验。X0=62.9,58.8,61.4,87.2,104.9,124.8,110.7,129.0,155.3,219.038-2.采用编制的GM(1,1)程序，探讨采用不同长度序列进行GM(1,1)建模的差异性。四.GM(1,1)残差模型1.问题的提出：当用X0建立的GM(1,1)模型检验不合格时；当用X0建立的GM(1,1)模型精度不理想时；2.残差模型的建立 X0构建的GM(1,1)模型为 为预测序列；为生成序列；对e1建立GM(1,1)模型表示 的偏差GM(1,1)修正模型的还原解为导数型残差修正注意：一般不使用全部残差建模。修正前的值？参见例题4.算例已知：试：作下一期的预测。解：step1：普通GM(1,1)建模；序号i1043.4543.4543.4500.00%2147.0590.589.930.570.63%3252.75143.25140.792.461.72%4357.14200.39196.443.951.97%5462.64263.03257.325.712.17%6568.52331.55323.957.62.29%未修正前的建模计算偏差较大且有逐步增大的趋势；问：为什么与 X1 比较？step2：基于残差的GM(1,1)建模；取n=5step3：GM(1,1)残差修正模型；残差未发生变化序号i1043.4543.4543.4500.00%2147.0590.589.930.570.63%3252.75143.25140.790.590.41%4357.14200.39196.440.330.17%5462.64263.03257.320.310.12%6568.52331.55323.950.850.26%修正后的建模计算结论：偏差大大降低；step4：采用残差GM(1,1)模型预测；注意：本残差修正GM(1,1)模型，采用的是残差求导 方式；残差修正是基于 预测模型之上进行的。本残差修正模型，适用于偏差均为正的情况。五.其它灰色预测模型 1.GM(1,N)预测模型 2.GM(2,1)预测模型 3.灰色系统预测模型1.GM(1,N)模型多个影响因素的建模问题系统有n个影响变量每个变量有m个数据，则有：对每一个原始序列进行累加形成生成数列由这n个序列组成的系统，围绕着主要因素X1建立微分方程：一阶n个变量的微分方程。由最小二乘法知：GM(1,N)的模型为注：例：对以下两组数据建立GM(1,2)模型，且 进行预测。i123456X102.2782.0373.2613.6793.851X207.047.6458.0758.538.7749.54解：step1:累加生成新序列。i12345X112.2784.9158.17611.85515.706X217.0414.68522.7631.2940.064step2:计算B和Yn阵。step3:计算a和b2。step4:写微分方程且写出解的形式。step5:模型的检验。当k=0时，GM(1,2)检验一览表i12345X12.2782.0373.2613.6793.851X1预测2.2782.1523.573.8483.871绝对偏差00.4830.3090.170.034相对偏差(%)018.399.484.620.88step5:模型的预测2.GM(2,1)预测模型问题的提出 GM(1.1):适用于具有较强指数规律的序列，只能描述 单调的变化过程。GM(2.1):解决非单调的摆动发展序列或有饱和的S性序 列问题。推导例题b/a2证明（略）3.灰色系统预测模型对于含有多个相互关联的因素与多个自主控制变量 的复杂系统，当任何单个模型都不能反映系统的发展 变化时。求出系统发展的数值解；弄清楚各影响因素的制约关系与程度；协同关系预测；灰色系统预测模型的建模步骤：Step1:确定系统的主导因素和关联因素；Step2:建立GM模型群，一般情况下，对主导因素建立 GM(1,1)模型，对影响因素建立GM(1,h)模型；Step3:根据GM模型群，将灰参数列成系统状态方程的矩 阵；Step4:用龙格-库塔法，求解状态方程；Step5:对状态方程的解作累减还原，击球的系统各因素的 拟合值和预测值；常微分方程的数值求解方法例题：农村电力需求预测90年补充：龙格-库塔法作业8-2.阅读文献，清楚灰色预测理论在电力负荷预测中的实际应用。小 结1.掌握灰色系统的基本概念。2.掌握灰色关联度计算，GM(1,1)预测模型。3.清楚GM(1,1)残差预测模型及GM(1,N)预测模型。3.了解GM(2,1)和灰色系统预测模型。