基于EMD和SVM的短期负荷预测

  • 格式:pdf
  • 大小:1.62 MB
  • 文档页数:5

基于EMD和SVM的短期负荷预测

祝志慧,孙云莲,季󰀁宇(武汉大学电气工程学院,武汉430072)摘󰀁要:为了有效预测具有一定周期性和随机性的非平稳时间序列性的电力负荷,提高预测精度,提出了结合经验模式分解(EMD)和支持向量机(SVM)的短期负荷预测法。该法运用EMD将负荷序列分解成若干个不同频率的平稳分量,分解后的分量突出了原负荷的局部特征,能更明显的看出原负荷序列的周期项、随机项和趋势项;根据各个分量的变化规律,选择合适的SVM参数和核函数构造不同的支持向量机分别预测;由SVM对各分量的预测值组合得到最终预测值。仿真试验表明,此方法与单一的SVM预测法及BP神经网络预测法相比,具有较高的精度和较强的推广能力。关键词:短期负荷;经验模式分解;本征模式分量;支持向量机;核函数;组合预测中图分类号:TM715文献标志码:A文章编号:1003󰀁6520(2007)05󰀁0118󰀁05Short󰀁termLoadForecastingBasedonEMDandSVMZHUZhi󰀁hui,SUNYun󰀁lian,JIYu(SchoolofElectricalEngineering,WuhanUniversity,Wuhan430072,China)Abstract:Thepowerloadisinherentlynon󰀁stationarytimeseries,andhasacertainperiodicityandrandomnessbyitself,soitisdifficulttoconstructtheforecastingmodel.Thetraditionalmodelsareconstructedonthebasisofthesuppositionthattheloadseriesarelinearandstationary,whichcannotpredictaccuratelytherealnon󰀁stationaryloadseries.Inordertopredicttheshort󰀁termpowerloadeffectivelyandleveluptheforecastprecision,ahybridforecas󰀁tingmethodbasedonEmpiricalModeDecomposition(EMD)andSupportVectorMachine(SVM)ispresentedinthispaper.EMDcandecomposenon󰀁stationarysignalsintosomesmoothandstationaryintrinsicmodefunctions(1MF)withdifferentfrequencyinthedifferentscalespacebythesiftingprocess.Whichisregardedasnewadaptivewaveletdecompositionmethod.SVM,anovelmachinelearningmethodbasedonthestructuralriskminimization(SRM)principle,ispowerfulfortheproblemwithsmallsample,nonlinearity,highdimensionandlocalminima.AccordingtotheoutstandingfeatureofEMDalgorithm,firstly,thepowerloadtimeseriesisdecomposedintoase󰀁riesofstationaryintrinsicmodefunctionsindifferentscalespaceviaEMDsiftingprocedure.Thelocalfeaturesofo󰀁riginalloadseriesareprominentintheintrinsicmodefunctionssothatitismoreobvioustoobservethecycle,ran󰀁domandtrendpartsoftheoriginalloadsequence.Secondly,accordingtothechangeregulationofeachofallresul󰀁tedintrinsicmodefunctions,therightparameterandkernelfunctionsarechosentobuilddifferentSVMrespectivelytoforecasteachintrinsicmodefunctions.Atlast,theseforecastingresultsofeachIMFarecombinedwithSVMtoobtainfinalforecastingresult.ThesimulationresultsshowthatthehybridmethodbasedonEMDandSVMhasfas󰀁terspeed,higherprecisionandgreatergeneralizationabilitythanthatofthesingleSVMmethodandthatoftheBPneuralnetworkmethod,whichprovesthatitisaneffectivemethod.Keywords:short󰀁termload;empiricalmodedecomposition;intrinsicmodefunctions;supportvectormachine;ker󰀁nelfunctions;hybridforecasting

0󰀁引󰀁言电力短期负荷预测是能量管理系统(EMS)中一个重要组成部分,随着电力市场改革的深入,负荷预测的作用愈来愈重要。负荷预测精度越高,越有利于提高发电设备的利用率和经济调度的有效性[1]。现有的负荷预测方法主要有回归分析法[2]、时间序列法[3]、人工神经网络法[4󰀁7]及小波󰀁神经网络组合法[8,9]。其中神经网络预测法因固有的非线性,在电力负荷预测中被广泛应用。但其存在网络结构难以确定、易陷入局部极小值、收敛速度慢等缺点。Vapnik等提出的支持向量机(SVM)是建立在VC维理论和结构化风险最小原则基础上,结构简单、学习速度快、全局最优、泛化性好,能较好的解决小样本、非线性、高维数和局部极小点等问题,被认为是神经网络的替代方法,已在模式识别、函数估计和信号处理领域广泛应用[10,11]。电力负荷为非平稳时间序列,但现有的预测方法是假设它为平稳信号,预测精度难以提高,造成了负荷预测的难度。为了更有效的把握负荷序列变化的信息对其进行分解,根据各分量的特点建模提高预测精度。经验模式分解(EMD)是一种新的信号󰀁118󰀁第33卷第5期2007年󰀁󰀁5月高󰀁电󰀁压󰀁技󰀁术HighVoltageEngineeringVol.33No.5May󰀁2007󰀁󰀁处理方法,它将非平稳信号按不同尺度的波动或趋势逐级分解成若干个本征模式分量(IMF),对信号作了平稳化处理,减少了信号间的特征信息的干涉或耦合。该方法吸取了小波变换的多分辨优势,同时克服小波变换中需选取小波基的困难,从信号本身的尺度特征出发对信号分解[12],具有自适应、正交性、完备性和瞬时频率等特点,属自适应小波分解方法。本文尝试将EMD和SVM相结合,对短期负荷进行预测。首先运用EMD将非平稳的负荷时间序列分解成具有不同特征尺度的平稳分量,然后用SVM对每个IMF分量进行负荷预测。最后,所有分量的预测值通过SVM组合得到最终的预测值。1󰀁原󰀁理1.1󰀁经验模式分解EMD是1998年由Huang提出,它通过对信号的󰀁筛选󰀁将信号分解成不同频率的IMF,IMF具如下特点:󰀁极值(极大值和极小值)数与过零点的数目相等或最多相差一个;󰀁在任意时刻,其上、下包络线的均值必须是零[13]。具体的分解过程是[14]:1)根据信号x(t)的局部极值求出其上、下包络线的平均值m1;2)将原数据序列减去平均包络后即可得一个去掉低频的新数据序列h1=x(t)-m1;判断h1是否为IMF,若不满足IMF条件,将h1看作新x(t),重复上述处理过程,直到h1满足IMF条件时,记c1=h1,视为IMF1;3)将r=x(t)-h1看作新的x(t),重复以上1)和2)步骤,即可依次得到IMF2,IMF3,󰀁,直到cn或r满足给定的终止条件时筛选结束;最后,原始的数据序列x(t)可表示为:x(t)=󰀁ni=1ci+r。(1)󰀁󰀁式(1)说明,可以将信号x(t)分解成频率从大到小排列的n个IMF分量ci和一个趋势项r之和。因每个IMF分量代表一个特征尺度的数据序列,故󰀁筛󰀁过程实际上将原始数据序列分解为各种不同特征波动序列的叠加。1.2󰀁支持向量机回归算法SVM回归算法的基本思想是[15,16]:通过用内积函数定义的非线性变换将输入空间变换到高维空间,使样本线性可分,在样本线性可分的情况下在高维空间作线性回归。对于给定的训练数据集s=((x1,y1),(x2,y2),󰀁,(xl,yl)),其SVM估计:f(x)=w󰀁󰀁(x)+b,(2)式中,󰀁(x)为从输入空间到高维特征空间的非线性映射,w为权系数向量,b为偏置。采用󰀁不敏感损失函数进行回归,则w和b由最小化式󰀁来估计:minimise12󰀁w󰀁2+C󰀁1l󰀁li=1L󰀁(yi,f(xi)),(3)

L󰀁(y,f(x))=|y-f(x)|-󰀁0,|y-f(x)|󰀁󰀁|y-f(x)|<󰀁。为寻找系数,引入松弛变量󰀁i,󰀁*i,使式(4)最小化minimise12󰀁w󰀁2+C󰀁li=1(󰀁i+󰀁*i),(4)

subjecttow󰀁󰀁(xi)+bi-yi󰀁󰀁+󰀁*i-w󰀁󰀁(xi)-bi+yi󰀁󰀁+󰀁i󰀁i,󰀁*i󰀁0󰀁󰀁i=1,2,󰀁,l,

式中C为惩罚因子,实现在经验风险和置信范围的折中。式󰀁所确定的约束优化问题是一个凸最优化问题,由拉格朗日理论可知,关于w求鞍点有:w=󰀁li=1(󰀁i-󰀁*i)󰀁󰀁(xi),(5)式中,󰀁i,󰀁*i为拉格朗日乘子,将式(5)代入式(2):f(x)=󰀁li=1(󰀁i-󰀁*i)󰀁k(xi,x)+b,(6)其中,k(xi,x)=(󰀁(xi)󰀁󰀁(x))称为核函数,任何函数只要满足Mercer条件,均可作为核函数。常用核函数有三:󰀁多项式核函数k(x,y)=((x󰀁y)+1)d(d为多项式阶数);󰀁径向基(RBF)核函数k(x,y)=exp(-󰀁x-y󰀁2/󰀁2)(󰀁为核函数宽度);󰀁Sig󰀁moid核函数k(x,y)=tanh(󰀁(x,y)+c)。为了求解式(4)的解,通常采用对偶定理,得到式󰀁优化问题的对偶式:maximise-12󰀁li,j=1(󰀁*i-󰀁i)(󰀁*j-󰀁j)k(xi,xj)+