用电量的分析与预测

1、用电量的分析与预测对各行业和居民的用电量的分析及今年第二季度用电量的预测第一部分绪论背景：问题的提出：2003年夏季，大面积的拉闸限电波及了全国 21个省（市）;今年一季度，24 个省级电网限电。据有关部门测算，预计2004年全国装机容量缺口在2000万千 瓦以上，除东北、山东电网发电装机容量略有富余外， 其他电网均可能出现不同 程度的缺电。此状况被称为“电荒”。为解“电荒”，去年底至今年初，各地发改委纷纷拿出了庞大的电力建设规 戈上2003年电力建设项目投产、新开工和在建规模均创历史最高水平的基础上， 今年的投资规模和装机总量又创了世界之最。电力投资过热。进入2004年5月份，最新统计数据显

2、示，一季度全国电力供需齐增15%, 共有24个省级电网拉闸限电；国家电网公司预计今年夏季用电高峰供电形势比 去年更加严峻，今年将是近几年电力缺口最大的一年。在这样一个“电荒”年里，国家发展和改革委员会对销售电价水平再次作出 调整，调整幅度按全国平均每千瓦时提高 1.4分钱。面对电力供需矛盾的日益突出，公司积极采取有效措施，加强负荷管理，力 求降低拉闸限电给售电造成损失。一是加强用电市场的调查分析和预测，掌握客 户的需求和供需平衡动态信息，预测电力市场的发展需求，实现电力的增供扩销； 二是加大对缺电的社会宣传力度，赢得用户理解和支持，严格规范自备电厂和地 方电厂管理，加强计划用电，加大对缺电地区

3、的网供能力；三是利用经济杠杆、 负控系统，优化限电方案，引导大用户避峰用电，加大峰谷用电执行力度，合理 错峰、避峰、移峰，最大限度地减少电力缺口的影响；四是用足用好电价政策， 挖潜增收，电价执行到位，提高售电平均单价，增大获利空间等等。通过一系列 的合理措施，公司售电量、平均电价节节攀高，四月份创历史新水平。近年，随三峡竣工发电，大量地方电厂投产，各地区电网改造，我国供电能 力迅速增强，用电量近乎指数趋势增长，但随着经济的快速发展整个社会对电的 需求量也迅速膨胀；气候异常，导致地区水电发电量减少，而（特别是夏季）居 民用电量巨增；国家几次调高电价，又会对我国的用电量有何影响呢？我国的用 电量还

4、会以怎样的趋势增长呢？今年，第二季度的用电量将会达到多少呢？以上信息多数来自国家电力信息网，连接 : HYPERLINK word文档可自由复制编辑研究的意义：电力需求预测是电网规划设计与建设的基础， 预测的准确与否，关系到电源 开发、电网建设、社会安定、居民生活及电力公司本身的发展。其作用与电力行 业的特殊性是密切相关的，主要表现为：1、电力行业的发供用、产供销必须瞬时完成，电量不能储存。即发电厂发 出的电量与当时的负荷水平应保持一致，并随时进行不断的调节以达到供需平 衡。这就要求电力需求的短期预测和中期预测必需精确并具有预见性。2、电力行业是技术资金非常密集的行业，建设周期较长。电力行业是

5、装备型行业，技术和资金非常密集，一般地讲，大型火电厂的建设周期为3-4年，水电厂的建设周期为6 10年。只有准确的预测才能保证如此巨大的资金能够使 用得当，既不会建设滞后造成电量短缺也不会超前建设造成大量的资金积压。3、电力行业是关系到国计民生的公用性事业。随着人民生活水平的日益提高，家庭电气化已是大势所趋，各行业对电力的依赖程度愈来愈高， 电力的短缺 或不能连续稳定的供电不但会严重影响居民的生活甚至会影响安定团结的政治局面。4、准确预测需求，有利于电力企业合理安排经营计划，减少购电成本，提 高经营效益。有助于供电量的预测，电力规划，利于谐调供电与社会需求之间的 关系。5、用电量，变化被人认为

6、是衡量一个地区经济发展水平最权威的标志。二.面对现今如此紧张的用电形势，建立一个模型，即能直观地看到用电量的发 展变化情况，从侧面了解社会的发展，居民的生活；与理论值对比，来分析影 响用电量变化的因素；又能对未来的用电量作一个比较精确的预测。因为全国的各季各部门的用电量数据不好找， 选取一个较有代表性的地区来 建立模型，现以缺电最严重的省浙江（绍兴市的用电量）为分析对象，建模型， 分析并对其用电量进行预测。当前的研究状况：由于需电量预测对电源开发、电网建设、社会安定、居民生活及电力公司本 身的发展都有很大的影响，因此世界各地的电力部门都十分重视电力需求预测工 作，设置专门的机构，由经济分析、需

7、电量预测、负荷预测等方面的专业人员来 从事电力需求预测工作，由于起步较早，各自开发出一种或几种适合于本国经济 运行特点的需电量预测方式，而且几种方式可以相互效验。方法有，部门分析法 进行预测，计量模型法，最终需求法，弹性系数法，类比法，积累法，计量经 济模型方法，还有比较复杂的经济模型法如灰色模型法、模糊数学模型法和神 经网络法等。概括介绍论述的内容（摘要）word文档可自由复制编辑用电量同时具有增长性和季节波动性的二重趋势，还受政治，经济，生活水 平等各方面多种因素的影响，这使用电量的变化呈现复杂的非线形组合特征，且 各行业及居民生活用电量的历史数据仍是复杂的非线形组合特征的序列。对各用电群

8、的历史用电数据，先用灰色模型处理，反映用电量的增长性特点；再用季节 变动指数来拟合用电量的季节性趋势；若拟和的精度不够高，再用所得的预测值 建立时间序列的 AR (p)模型，作修整微调，并找历史数据中拟合最好的数据 段，提高精度。影响各行业及居民生活用电量的因素有所不同， 不同用电人群对 应的用电量有不同的增长特点和波动归律。 分别对各用电群的历史用电数据建模 型分析并作预测，对各行业及居民生活用电量的预测值求和作为社会总用电量的 预测。所使用的方法：.灰色预测GM(1, 1)模型：灰色系统是指部分信息已知、部分信息未知的系统。灰色系统理论的实质 是 将无规律的原始数据进行累加生成，得到规律性

9、较强的生成数列再重新建模。由 生成模型得到的数据通过累加生成的逆运算一一累减生成得到还原模型，由还原模型作为预测模型。灰色模型(简称GM真型)的机理是将电力系统的负荷看成 为在一定范围内变化的灰色量，将作用于其中的随机过程看做在一定幅区间和一 定时区间变化的灰色过程。GM(1,1)模型(一次拟合参数模型)，通过对原始数 据进行累加生成，得到规律性较强的序列，用指数曲线去拟合得到预测值。灰色预测具有要求样本数据少、原理简单、运算方便、预测精度高、可检验 等优点，可以较好地对非线性系统进行预测。.季节指数调整的模型：移动平均法：将移动平均值中，周期变动的影响去除。季节变动指数：反映实际值与理论值的

10、差异，是一个比值。对季节指数平均：尽可能消除得出的季节变动指数中存在的不规则变动。.时间序列模型中的自回归序列 AP(p)运用自回归模型中的AP(p)模型，根据时间序列统计规律性的分析，构造拟 合它的最佳数学模型，浓缩时间序列的信息.利用拟合的数学模型给出时间序列的预测值,给出预测结果的精度分析第二部分正文前提假设：.预测的时期没有对用电量影响巨大的突发事件(严重灾情，战争等).用电数据由电力部门给出，电力部门准确测得，能真实反应用电状况.同一用电群中，个体的用电量的增长特点和波动归律有一致性。4,用平均误差来衡量拟合的精度，预测的精度。word文档可自由复制编辑模型组建：1:通过累加技术，使

11、数据形成指数规律，从而建立一阶微分方程，求解微分方 程，再累减还原，得灰色预测值，求得拟合曲线后对对象的将来发展进行预测。 累加生成方法:设原始数据序列：x(1),x(2), -(40)作一次累加生成(记为1-AGQ ,彳4到x1(1),x1(2), x1(4n),t其中 xtX x(i),t=1,2 - 4ni 1一般地，对非负数列，累加生成次数越多，数列的随机性就弱化得越多。 当累加生成次数足够大时，时间序列便由随机转化为非随机了。在GM1型中,一般只对数列作1-AGQ构造一阶线性微分方程：叽 aX (1) = udt其中，a称为模型的发展系数，它反映 X与X的发展趋势，u称为模型的协调系

12、 数，反映了数据间的变化关系。求解微分方程:得X(t+1)=AeH+U (*)a利用最小二乘法求解参数a,u,得：产=(BB),BY一u-x(2)l TOC o 1-5 h z f -2(x1(1)+x1(2)x(3)式中，b=-1(x1(2)+x1(3)12.:_x(4n)_111-2(x (4n-1)+x (4n)将a, u代回微分方程(*),得：x(k +1) =(x(1)(1) u)e7k +u , k =0,1，4n + m ,a a累减还原:得灰色预测值序列：x0 x0(4n+ m)(共4n+m个，m个为预测的将来的用电量)其中：x0(t+1)= (1-ea)x(1)-ueta该部

13、分模型不需要大量样本，也不需要考虑序列的分布规律，能很好地对非 线性系统进行预测，并得到拟合的曲线。根据该模型得到的预测值序列作出的拟 合曲线可以较好地反映原始数据(用电量)的发展趋势。2:由上模型作图可知原始数据序列呈季节性波动，所以再做季节调整的模型。该模型先求季节移动平均值；再用实际值除以季节移动平均值，得季节变动指word文档可自由复制编辑数列；对同一季节的季节变动指数列进行平均，得各季的季节指数；用季节指 数修正对应季节的灰色预测值，得季节指数调整后的预测值序列。I( j)为季节变动指数，j=1,2,3,4.考虑长期趋势条件下I (t) = 出x(t)x(t)为原始数据序列，x(t)

14、为季节移动平均值 1 .x(t) (x(t 2) x(t 1) x(t) x(t -1)4将不同年份的同一季节的季节变动指数进行平均 1I(j)=1(|(j,1)+I(j,2)+. + I( j,n) , n 为历史数据所跨年份。 n再进行平均得季节指数而=而下工、丽j 1最后用季节指数修正，得季节调整后的预测值序列：y(t) = x0(t)父7(下x0(t)为第j季度的预测值3:如果上模型所得的预测值精度仍不令人满意，再对季节指数调整后的预测值 序列建立自回归AR(p)模型。对已得的预测值序列作微小的调整，除去数据中的噪声，得最终预测值序列(z( p + 1),并且该数据段是使预测精度最高的

15、数据段(即与实际曲线拟合得最好的预测值的数据段)。y(t) =%y(t T) +%y(t 2) +. +py(t - p) (t) , t=p+1, ,m+n其中外是参数,i=1,p ; N(t)是白噪声是随机扰动项用最小二乘法估计参数 由(i=1,4n m,p ),(要使Z N(t)2最小)t =p 1N(p+1) IN =:产(4n + m)_一 y(p+1)B =:)(4n + m)-word文档可自由复制编辑-y(p)y(p -1) y(i) Iy(p+i)y(p) y(2)A -.y(4n+m1) y(4n+m 2) y(4n+m p)_ 求得： =(AA),AB,。2是N(t)的方

16、差2 BB-BA(AA)AB cy =4n m p求最终预测值,得：z0(t) = A 二找出使z0(p+1)z0(4n+m)平均误差最小的p,在该p下求出的为最终预测值 z0(p+1)z0(4n+m).4:用平均误差来衡量拟合的精度，预测的精度。以上模型中每做一个模型平均 误差都会减小，也就是说拟合的精度(预测的精度)提高。相对误差=(预测值-实际值)/实际值，平均误差=(三相对误差的绝对值)/相对误差的个数5:该模型不但可以对将来的用电量作出预测，还可以根据拟合出的理论值来分 析当前用电量的情况。模型的推广：模型针对用电量，再较短的时期内，季节波动性较明显，而年与年之间的波 动性不明显，所

17、以模型第2步只作季节调整。可将模型第2步充实为：对数据周 期的调整。对有稳定的发展趋势，有一定波动周期的数据序列都可以用该模型来分析预 测。若所给数据的发展趋势稳定，周期性明显，那么拟合和预测的精度将相当高。 比如，生活中某种电器(空调，冰箱等)，防晒用品的多年的各季的销售量和销 售价格，城市居民用水量，学校食堂每日的消费量等。运用模型,对某地区用电量的历史数据进行分析：该地区用电指标大致在 112万118万千瓦左右，加上地方电厂出力 50万千瓦，电力缺口在 60万千瓦以上。在当前全市工业用电实行“停三开四”，有的企业实行停产的情况下，拉限电的范围还在进一扩大。word文档可自由复制编辑 数据

18、是根据模型用Metlab编程，得到的，程序见：程序文件（论文最后） 从数据表中可以看出，每一步建模，预测精度（平均误差）都有提高。用Metlab作的图可以直观的反映用电量变化的趋势和特征。以下数据表中：di ,d2 ,d3相对误差（）：（实际值-预测值）/实际值*100平均误差=各相对误差的名对值的和/个数da1,da2,da3（其中23=各相对误差的名对值的和/（总个数-p）p表示:p+120,用后20-p个预测数据能使精度更高u,a灰色预测中的参数以下图中：实际用电量曲线（红色），灰色预测曲线（绿色），季节调整后预 测曲线（蓝色），最终预测值曲线（黑色）；横轴：1-21为99年第2季度-2

19、004年 第2季度，21对应的是预测值；纵轴：用电量，单位（亿千瓦时）第一产业，即农业，拟合精度（预测精度）最终达到8.0034%灰色预测季节指数调整yd2预测值相自回归微调x实（亿xo际值预测值 .千瓦（亿千瓦d1相对误差zo对误差预d3测值相对误差时间99年2季度3季度4季度00年1季度2季度3季度4季度01年1季度2季度3季度4季度02年1季度2季度3季度4季度03年1季度2季度3季度4季度04年1季度时）时）0.450.450.52 0.551390.49 0.547080.43 0.542810.59 0.538560.62 0.534350.56 0.530170.52 0.526

20、030.52 0.521920.58 0.517840.53 0.513790.44 0.509770.51 0.505780.6 0.501830.52 0.497910.48 0.494010.56 0.490150.58 0.486320.41 0.482520.31 0.47874（%）0-6.0374-11.65-26.2348.718213.8145.3261-1.1594-0.368510.7183.0591-15.8570.8265916.3624.2488-2.919512.47316.152-17.687-54.434（亿千瓦时）（）（亿千瓦时）（%）p0.46016-2.

21、2570.46016050.62148-19.5150.621480a0.53741-9.67630.5374103.5050.47116-9.57170.471160u0.550726.6580.5507200.557090.602272.85990.602272.8599 da10.52086.99940.52086.999411.4020.456612.1930.456612.193 da20.5337-2.63380.5337-2.63388.38640.58365-0.629830.58365-0.62982 da30.504714.77250.504714.77258.00340.

22、44248-0.564210.44248-0.56420.5172-1.41170.5172-1.41170.565615.73120.565615.73120.489115.94120.489115.94120.4288110.6650.428811.07E+010.5012110.4980.5012110.4980.548135.49470.548135.49470.47399-15.6070.47399-15.6070.41555-34.0490.41555-34.049word文档可自由复制编辑2季度0.4750.485720.48572平均误差11.4028.38648.0034表中

23、y与z0两列是一样的，原因是：自回归模型中随机扰动项W(t及b都很小，在表中看不出差别,在Metlab画的Metlab图中可以看出它们是有细微差别的黑线和蓝线并不完全重合，所以自回归模型对精度的调整仍有微小的作用 局部图由图可以看出：2004年第1季度的用电量明显下降，与今年的气候密切相关，也与电价， 政府限电等因素有关。04年1季度用电量特别低，使得理论值整体都比实际值 低。农业用电年与年之间有很大的波动，原因是农业受气候因素影响很大。因为严重受，去年，今年，特别干热的气候影响，所以对农业数据拟合的不 是太好，精度不是很高。农业的用电量是所有之中拟合得最差的。如果有更多的数据可以在作了季节调

24、整后， 再作一下年度调整，可以提高预 测的精度。该地区的农业用电再逐年下降，说明其农业有衰退的发展趋势，不是该地区 重点发展产业。word文档可自由复制编辑第二产业，工业,拟合精度（预测精度）最终达到1.9969%灰色预测季节指数调整x xo di y d2实际值预测值预测值（亿千瓦（亿千瓦相对误差（亿 千相对误差时间时）时）(%)瓦时）(%)zod3P99年2季度14.8714.87015.401-3.57363季度13.9713.4483.73713.7991.2269a4季度15.1614.097.061214.8272.19683.454200年1季度13.2814.762-11.15

25、713.0771.5313u2季度16.3915.4665.637816.0192.265612.4443季度16.516.2041.79516.626-0.76577da14季度17.7616.9774.409417.865-0.593816.025601年1季度15.4117.787-15.42415.756-2.2482da22季度19.0318.6352.073419.301-1.42611.99693季度19.3919.524-0.6936220.034-3.3192da34季度20.6420.4560.8915721.527-4.295702年1季度18.6121.432-15.1

26、6318.985-2.01732季度23.4622.4544.286523.2570.866023季度24.0123.5262.017224.139-0.537694季度26.2124.6485.959425.9381.037403年1季度23.5225.824-9.795722.8762.73762季度27.8127.0562.711328.023-0.765443季度30.7728.3477.875229.0865.47314季度32.0929.6997.450431.2542.606404年1季度27.6931.116-12.37327.5640.454562季度32.60133.766

27、平均误差6.02561.9969该误差已很小，所以不再做自回归模型word文档可自由复制编辑只做了灰色预测模型，和季节调整模型，因为这样已可以达到1.9969%的精度，即使再做自回归模型也提高不了精度。由图可以看出：该地区的工业在以相当快的速度发展，受，去年，今年，特别干热的气候影 响，水电发不出，是供电量严重不足，使去年和今年的工业用电量偏低（对各产业 都造成用电量偏低的结果），导致拟合的曲线后尾的用电量偏低.但总的来说，工 业的用电量拟合的很好，预测也相应更为可信.在当前全市工业用电实行“停三开四”，有的企业实行停产的情况下，拉限 电的范围还在进一扩大。如果供电不受限制，那工业将以更快的速

28、度发展。第二产业，即服务行业，拟合精度（预测精度）最终达到2.2991%灰色预测1xxod1y1d2zod3实际值预测值预测值预测值（亿千（亿千相对证差（亿千相对误差（亿千相对误差时间瓦时）瓦时）(%)瓦时）(%)瓦时）(%)P99年2季度0.790.7900.711719.9102 0.71171093季度1.07 0.9011315.7821.0935-2.19661.09350a4季度0.94 0.96138-2.2747 0.914232.742 0.9142303.607700年1季度1.011.0257-1.551 0.958665.0834 0.958660u2季度1.031.0

29、942-6.23750.98584.29090.985800.821153季度1.411.167417.2051.4166-0.470611.41660da14季度1.071.2455-16.3991.1844-10.691.184402.299101年1季度1.181.3288-12.6061.2419-5.24931.24190da22季度1.21.4176-18.1331.2771-6.42591.277103.61743季度1.841.512417.8051.83530.257971.83530.25797da34季度1.521.6135-6.1521.5344-0.945141.53

30、44-0.945142.299102年1季度1.581.7214-8.94921.6089-1.83151.6089-1.83152季度1.631.8365-12.6691.6545-1.5031.6545-1.5033季度2.331.959315.912.3776-2.04182.3776-2.04184季度2.052.0903-1.96611.98783.03541.98783.035403年1季度2.162.2301-3.24412.08443.50092.08443.50092季度2.222.3792-7.17062.14343.45022.14343.45023季度3.172.538

31、319.9283.08012.83443.08012.83444季度2.652.708-2.18842.57522.8242.57522.82404年1季度2.622.8891-10.272.7003-3.06562.7003-3.06562季度3.08222.77682.7768平均误差2.29913.61742.2991word文档可自由复制编辑局部图由图可以看出：该地区的服务类行业在以相当快的速度发展 ，说明其服务业发展非常之快， 与第二产业一样，受去年，今年，特别干热的气候影响，水电发不出，是供电量严 重不足，使去年和今年的用电量偏低，导致拟合的曲线后尾的用电量偏低.但拟合 曲线的后半

32、段（即最终预测值的曲线）拟合的相当好 ，预测也相应非常可信.居民生活,拟合精度（预测精度）最终达到4.5979%灰色预测季节指数调整自回归微调xxod1yd2zod3实际值预测值预测值预测值（亿千（亿千相对证差（亿千相对误差（亿千相对误差时间瓦时）瓦时）(%)瓦时）(%)瓦时）(%)P99年2季度1.361.3601.156814.9431.1568053季度1.831.55515.032.039-11.4222.0390a4季度1.531.6096-5.20341.49961.98461.499603.567500年1季度1.571.6662-6.12671.51033.80161.5103

33、0u2季度1.461.7248-18.1341.467-0.48051.46701.48133季度2.341.785423.7022.3412-0.050582.3412-0.05058da14季度1.671.8481-10.6671.7219-3.10531.7219-3.105314.448word文档可自由复制编辑01年1季度1.781.9131-7.4771.73412.57761.73412.5776 da22季度1.671.9803-18.5831.6844-0.86251.6844-0.86255.083季度2.552.049919.612.6881-5.41612.6881-5

34、.4161 da34季度1.82.122-17.8891.977-9.83391.977-9.83394.597902年1季度1.912.1966-15.0041.9911-4.24551.9911-4.24542季度2.022.2738-12.5641.9344.25711.9344.25713季度2.962.353720.4833.0864-4.27183.0864-4.27184季度2.322.4364-5.0192.272.15642.272.156403年1季度2.352.5221-7.32252.28612.71762.28612.71762季度2.22.6107-18.6692.

35、2206-0.936122.2206-0.936113季度4.192.702535.5013.543815.4223.543815.4224季度2.482.7975-12.8022.6063-5.09462.6063-5.094604年1季度2.432.8958-19.1692.6249-8.02082.6249-8.02082季度2.99762.54972.5497平均误差14.4485.084.5979局部图：,说明该地区居民生活水平在提由图可以看出：该地区的居民生活用电量在逐年稳步上升word文档可自由复制编辑 高。去年夏天气特别热，空调等降温设备用电，使用电量迅速上跳。与第二产业 一样

36、，受去年，今年，特别干热的气候影响，水电发不出，是供电量严重不足，使 去年和今年的用电量偏低，导致拟合的曲线后尾的用电量偏低.春季生活对用电没有特别大的需求， 加上政府的一些控制用电的措施，还有 价格的上调，使得今年第1季度得用电量比正常情况下要低，这也导致拟合曲线 后尾的值偏低。观察以上各图，发现2004年1季度的用电量都比预测值还低，今年春雨水非 常少，政府控制用电等原因，今年供电形势比去年更加严峻.最终预测值,及原始数据汇总：1999 年 2000 年2001 年2季3季4季1季2季3季4季1季2季3季第一产业（亿千瓦时）0.450.520.490.430.590.620.560.520

37、.520.58第一产业（亿千瓦时）14.8713.9715.1613.2816.3916.517.7615.4119.0319.39第二产业（亿千瓦时）0.791.070.941.011.031.411.071.181.21.84居民生活（亿千瓦 时）1.361.831.531.571.462.341.671.781.672.55社会用电 量（亿千 瓦时）17.4717.3918.1216.2819.4720.8621.0618.8922.4224.362002 年2003 年2004 年1季2季3季4季1季2季3季4季04年1季04年2季0.440.510.60.520.480.560.58

38、0.410.310.4857218.6123.4624.0126.2123.5227.8130.7732.0927.6933.7661.581.632.332.052.162.223.172.652.622.7768word文档可自由复制编辑1.912.022.962.322.352.24.192.482.432.549739.5782222.5327.6329.8930.128.5132.7838.7137.6333.05最终的预测值（2004年第2季度）：（单位:亿千瓦时）0.48572+33.766+2.7768+2.5497=39.57822的预测值之和，为社会总用电量的预测值精度分析

39、：（单位：亿千瓦时）正常预测偏差0.48572* 8.0034%+33.766* 1.9969%+2.7768* 2.2991%+2.5497* 4.5979%=0.8942运用该模型得到的主要结论：2004年第二季度各产业及居民生活用电量将分别是：0.48572 0.0388 ,33.766 0.6742 ,2.7768 0.06384 ,2.5497 0.01723社会总用电量将是39.57822 土 0.8924（单位：亿千瓦时）改进：若能找到更多的数据，类似模型第2步（季节调整），对数据进行年度调整，可使拟合和预测的精度提高。如某一个用电群受某一因素影响明显，而这一因素又无明显的周期性

40、，寻找更多关于这一因素的信息，寻找一个合理的修正系数对 于预测值进行修正，使预测的精度再一次提高（如农业的用电量，可用气候系数 修正）。还可对第3,第4季度以及以后的进行预测。写了推广程序（见 程序文件（论 文最后）参考文献电力系统负荷灰色预测的新方法。电力系统自动化，周平，周家启灰色预测与决策，邓聚龙，华中理工大学出版社，1992测量数据建模与参数估计，王正明，易东云，国防科技大学出版社，1996计量经济模型与经济预测，Pindyck R.S. , Rubinfeld D.L.机械工业出版社，1999Pardo A,Meneu V,Valor E.Temperture and seasona

41、lity influences onSpanish electricity loadJ.Energy Economics,2002,（24）:55-70数据处理程序及扩展程序:数据处理程序：word文档可自由复制编辑 x=;x1=；%输入原始数据% （灰色预测）x1(1)=x(1);for i=1:1:19x1(i+1)=x1(i)+x(i+1);end;%以上，作一阶累加 1-AGO得x1c=； for i=1:1:19c(i)=-1/2*(x1(i)+x1(i+1);end;B=c,ones(19,1); Y=x(2:20); Y=Y;%以上，构造矩阵 B 和丫b=B*B;d=B;d1=d

42、*B;d2=d*Y;c=d1d2;a=c(1);u=c(2);%以上，利用最小二乘法解参数a和ux0=;x0(1)=x(1); for i=1:1:20 x0(i+1)=(1-exp(a)*(x(1)-u/a)*exp(-a*i);end;%累减还原得到灰色预测模型x0t=1:1:21;plot(t,x0,g,t(1:20),x,r);%做实际售电量曲线（红色），灰色预测曲线（绿色）%分析误差；相对误差（）:（实际值-预测值）/实际值*100%平均误差：各相对误差的绝对值的和/个数xq=x0(1:20);d1=(x-xq)./x.*100;w=d1;W=;for i=1:1:20if w(i)

43、=0W(i)=w(i);elseW(i)=-w(i);end;end;w=;da1=sum(W)/20;xa=;% （季节调整模型）for i=2:1:18xa(i)=0.25*(x(i+2)+x(i+1)+x(i)+x(i-1);I1(i)=x(i)/xa(i);end;%求个季的移动平均值xa,季节变动指数列I1a=0;b=0;word文档可自由复制编辑 c=0;d=0;for i= 0 4 8 12c=c+I1(i+2);d=d+I1(i+3);a=a+I1(i+4);b=b+I1(i+5);end;c=c+I1(18); TOC o 1-5 h z I2=b/4 c/5 d/4 a/4

44、;%以上，对不同年份的变动指数进行平均，得各季的平均季节指数I2e=sum(I2)/4;I=I2./e;%求得各季的季节指数II0=I,I,I,I,I,I(1); HYPERLINK l bookmark105 o Current Document y=x0.*I0;%将灰色预测值用季节指数修正，得季节调整后的预测值yplot(t(1:20),x,r,t,x0,g,t,y,b);yq=y(1:20);d2=(x-yq)./x.*100;w=d2;W=;for i=1:1:20if w(i)=0W(i)=w(i);elseW(i)=-w(i);end;end;w=;da2=sum(W)/20;%

45、作季节调整后预测值的误差分析；相对误差() d2,平均误差da2d3=;da3=0 0;z0=;% (自归模型)for p=3:1:15%循环找使平均误差最小的pA=;A=y(p:20);for i=1:1:p-1B=y(p-i:20-i);A=A,B;end;B=y(p+1:21);e=A;v=e*B;o=e*A;r=B;e=oe;b=ov;d=(r*B-r*A*e*B)/(21-p);Z0=A*b-sqrtm(sqrtm(d*d);% 加入修正项 o-zq=Z0(1:20-p);xq=x(p+1:20);%以上，构造自回归模型，用最小二乘法估计参数值得bw=(xq-zq)./xq.*100

46、;W=;for i=1:1:20-pif w(i)=0 W(i)=w(i);elseword文档可自由复制编辑 end;end;s=sum(W)/(20-p);d3=d3,zeros(1,p),w,0;da3=da3,s;%误差分析；相对误差() d3,平均误差da3,并累记在对应矩阵中z0=z0,zeros(1,p),Z0;end;%结束对p的循环so=da3(3);for i=3:1:15if da3(i)=0 W(i)=w(i);elseW(i)=-w(i);end;end;s=sum(W)/(20-p);d3= zeros(1,p),w,0;da3=s;z0=y(1:p),Z0,;%按

47、模型得出最终预测值z0,相对误差() d3,平均误差da3plot(t(1:20),x,r,t,x0,c,t,y,b,t(p+1:21),z0(p+1:21),k);%作出所有曲线的图：实际售电量曲线(红色)，灰色预测曲线(绿色)，季节调整后预测曲线(蓝色)，最终预测值曲线(黑色)。% (在做自回归模型时由于用最小二乘法，计算的扰动项太小而出现警告)Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 1.794778e-017.word文档可自由复制编辑 (Type wa

48、rning off MATLAB:nearlySingularMatrix to suppress this warning.) Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 1.794778e-017.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning.) Warning: Matrix is close to singular or badly scaled.Results

49、 may be inaccurate. RCOND = 7.015279e-018.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning.) Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 7.015279e-018.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning.) Warning: M

50、atrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 1.489785e-017.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning.) Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 1.489785e-017.(Type warning off MATLAB:nearlySi

51、ngularMatrix to suppress this warning.) Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 8.093471e-018.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning.) Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND

52、 = 8.093471e-018.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning.) Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 3.084048e-018.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning.) Warning: Matrix is close to singula

53、r or badly scaled.Results may be inaccurate. RCOND = 3.084048e-018.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning.) Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 4.001137e-018.(Type warning off MATLAB:nearlySingularMatrix to suppress

54、this warning.) Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 4.001137e-018.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning.) Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 6.361787e-019.(Type wa

55、rning off MATLAB:nearlySingularMatrix to suppress this warning.) Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 6.361787e-019.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning.) Warning: Matrix is close to singular or badly scaled.Results

56、 may be inaccurate. RCOND = 6.365249e-019.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning.) Warning: Matrix is close to singular or badly scaled.word文档可自由复制编辑Results may be inaccurate. RCOND = 6.365249e-019.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning

57、.)Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 8.430305e-019.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning.)Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 8.430305e-019.(Type warning off MATL

58、AB:nearlySingularMatrix to suppress this warning.)Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 2.990102e-019.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning.)Warning: Matrix is close to singular or badly scaled.Results may be inaccura

59、te. RCOND = 2.990102e-019.(Type warning off MATLAB:nearlySingularMatrix to suppress this warning.) da2 da3n=?;%n输入数据的年数（本文n=5）m=?;%mW4往后预测的个数（本文 m=1）x=;%输入连续n年的各季度的4n个售电量的数据（实际值 x）x1=;% （灰色预测）word文档可自由复制编辑 x1(1)=x(1);for i=1:1:4*n-1x1(i+1)=x1(i)+x(i+1);end;%以上，作一阶累加1-AGO得x1c=；for i=1:1:4*n-1c(i)=-1/

60、2*(x1(i)+x1(i+1);end;B=c,ones(4*n-1),1); Y=x(2:4*n); Y=Y;%以上，构造矩阵B 和丫b=B*B;d=B;d1=d*B;d2=d*Y;c=d1d2;a=c(1);u=c(2);%以上，利用最小二乘法解参数a和ux0=;x0(1)=x(1); for i=1:1:4*n+m-1x0(i+1)=(1-exp(a)*(x(1)-u/a)*exp(-a*i);end;%累减还原得到灰色预测模型x0t=1:1:4*n+m;plot(t(1:4*n+m),x0,g,t(1:4*n),x,r);%做实际售电量曲线(红色)，灰色预测曲线(绿色)xq=x0(1