王　超,李　森

(河南工業(yè)職業(yè)技術(shù)學(xué)院,河南 南陽 473000)

?

變參數(shù)QPSO算法優(yōu)化神經(jīng)網(wǎng)絡(luò)的短期電力負(fù)荷預(yù)測

王超*,李森

(河南工業(yè)職業(yè)技術(shù)學(xué)院,河南 南陽 473000)

摘要:為了提高電力負(fù)荷預(yù)測精度,提出了一種變參數(shù)量子粒子群(VPQPSO)算法優(yōu)化RBF神經(jīng)網(wǎng)絡(luò)的短期負(fù)荷預(yù)測模型(VPQPSO-RBFNN)。首先利用電力負(fù)荷的混沌性,對短期負(fù)荷時間序列進(jìn)行相空間重構(gòu);然后采用變參數(shù)QPSO算法優(yōu)化RBF神經(jīng)網(wǎng)絡(luò)參數(shù)對重構(gòu)后的短期負(fù)荷時間序列進(jìn)行學(xué)習(xí),建立短期電力負(fù)荷最優(yōu)預(yù)測模型;最后采用對某地區(qū)短期電力負(fù)荷進(jìn)行預(yù)測。VPQPSO-RBFNN可以準(zhǔn)確描述復(fù)雜多變的電力負(fù)荷變化趨勢,提高了電力負(fù)荷的預(yù)測精度,仿真結(jié)果驗證了VPQPSO-RBFNN可以用于電力系統(tǒng)負(fù)荷預(yù)測。

關(guān)鍵詞:電力負(fù)荷;RBF神經(jīng)網(wǎng)絡(luò);變參數(shù);量子粒子群算法;相空間重構(gòu)

短期負(fù)荷受到氣象、節(jié)假日等因素影響,規(guī)律難以確切把握,因此短期負(fù)荷預(yù)測一直是國內(nèi)外研究的熱點[1]。

短期負(fù)荷是一種復(fù)雜非線性系統(tǒng),線性方法難以準(zhǔn)確描述其變化趨勢,當(dāng)前主要基于人工神經(jīng)網(wǎng)絡(luò)、支持向量機等機器學(xué)習(xí)方法進(jìn)行建模預(yù)測[2]。其中,BP神經(jīng)網(wǎng)絡(luò)不需要先驗知識,能對系統(tǒng)進(jìn)行非線性、無限的逼近,成為使用最為廣泛的短期負(fù)荷預(yù)測算法[3]。但BP神經(jīng)網(wǎng)絡(luò)預(yù)測性能與初始連接權(quán)值、閾值等參數(shù)密切相關(guān)。目前BP神經(jīng)網(wǎng)絡(luò)參數(shù)優(yōu)化方法主要有:遺傳算法(GA)、粒子群優(yōu)化(PSO)算法、蟻群優(yōu)化算法(ACO)、模擬退火算法(SA)等[4-6],這些算法均存在不足,難以找到最優(yōu)的連接權(quán)值和閾值,導(dǎo)致BP神經(jīng)網(wǎng)絡(luò)的預(yù)測結(jié)果不穩(wěn)定,有時難以令人滿意[7]。2004年,Sun Jun等提出一種新型優(yōu)化算法—量子粒子群QPSO(Quantum Particle Swarm Optimization)算法,具有操作參數(shù)少、易于編程、實現(xiàn)容易和適用性強等優(yōu)點,受到廣泛關(guān)注[8]。

為了提高短期負(fù)荷的預(yù)測精度,提出一種變參數(shù)量子粒子群優(yōu)化VPQPSO(Varying Parameters Quantum Particle Swarm Optimization)算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)的短期負(fù)荷預(yù)測模型(VPQPSO-BPNN),并通過仿真實驗對模型性能進(jìn)行測試。

1　相空間重構(gòu)和BP神經(jīng)網(wǎng)絡(luò)

1.1相空間重構(gòu)

相空間重構(gòu)的目的在于從實測時間序列中的某一分量了解非線性動力系統(tǒng)相空間的幾何特性,并在高維相空間中恢復(fù)混沌吸引子,Packard和Takens等[9]提出了采用延遲坐標(biāo)法對一維混沌時間序列x(i),i=1,2,…,n進(jìn)行相空間重構(gòu),對于短期負(fù)荷時間序列:x(t),t=1,2,…,N,通過選擇合適的嵌入維數(shù)m和延遲時間τ,可以產(chǎn)生一組新的多維時間序列:

X(i)=(x(i),x(i+τ),…,x(i+(m-1)τ)T

i=1,2,…,N

(1)

式中,m為嵌入維數(shù);τ為延遲時間。

從式(1)可知,在相空間重構(gòu)中,m和τ選取十分重要,本文采用互信息法計算τ,采用CAO方法選擇m。

1.2BP神經(jīng)網(wǎng)絡(luò)

BP神經(jīng)網(wǎng)絡(luò)是一種多層前饋神經(jīng)網(wǎng)絡(luò)。網(wǎng)絡(luò)結(jié)構(gòu)如圖1所示,由輸入層、輸出層、一個或多個隱含層組成。每層由許多并行計算的簡單神經(jīng)元組成,層與層之間的神經(jīng)元采用全互連的方式,兩層之間的每個神經(jīng)元都關(guān)聯(lián)一個權(quán)值,同層神經(jīng)元之間無連接。各神經(jīng)元的傳遞函數(shù)通常為連續(xù)可微的S型函數(shù),它可以實現(xiàn)從輸入到輸出的任意非線性映射[10]。

圖1　BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)

BP神經(jīng)網(wǎng)絡(luò)模型的主要思想就是使網(wǎng)絡(luò)輸出層的誤差平方和達(dá)到最小。在網(wǎng)絡(luò)的學(xué)習(xí)訓(xùn)練過程中,網(wǎng)絡(luò)首先隨機產(chǎn)生一組神經(jīng)元之間的連接權(quán)值和閾值,訓(xùn)練樣本輸入網(wǎng)絡(luò)通過正向傳播產(chǎn)生一個輸出結(jié)果,將它與期望值比較,若兩者誤差未達(dá)到設(shè)定的誤差,誤差信號沿原路反傳回來,同時修正各層神經(jīng)元的權(quán)值和閾值,以減小誤差,正向傳播的輸出計算與反向傳播的權(quán)值修正反復(fù)交替進(jìn)行,直到網(wǎng)絡(luò)的實際輸出與期望輸出之間的誤差值達(dá)到要求,從而得到滿意的連接權(quán)值和閥值。

設(shè)輸出層有m神經(jīng)元,BP網(wǎng)絡(luò)的實際輸出是y,期望輸出是y′,函數(shù)E為:

(2)

每個權(quán)值的修正值為:

(3)

式中,wij是從神經(jīng)元i到神經(jīng)元j的連接權(quán)值,η是學(xué)習(xí)速率,Ij是中間第j個隱藏層的傳輸函數(shù)。

在BP神經(jīng)網(wǎng)絡(luò)開始訓(xùn)練前,需要選擇最合適的連接初始權(quán)值,通常情況下,采用隨機方式初始化,易使BP神經(jīng)網(wǎng)絡(luò)出現(xiàn)收斂速度慢、易陷入局部最優(yōu)解等缺陷,影響B(tài)P神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)和泛化能力,為此,本文采用量子粒子群算法(QPSO)優(yōu)化BP神經(jīng)網(wǎng)絡(luò)參數(shù)。

2　VPQPSO算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)模型

2.1VPQPSO算法

雖然PSO算法擁有控制參數(shù)少、易于實現(xiàn)等優(yōu)點,但不能以概率100%收斂于全局最優(yōu)解,為此,SunJun等從量子力學(xué)的角度出發(fā)提出了一種新的PSO算法—量子粒子群(QPSO)算法,以保證算法收斂到全局最優(yōu)解。量子粒子群中所有粒子的平均最好位置點為pmbest(MeanBestPosition),pmbest定義為:

(4)

在QPSO算法中,每一個粒子必須收斂于各自的隨機點P,第i個粒子第d維的位置為:

(5)

式中,uid是均勻分布在(0,1)上的隨機數(shù);Lid(t)的值由下式確定:

Lid(t)=2α(t)|pmbestid(t)-xid(t)|

(6)

那么,QPSO算法的粒子位置更新公式為:

xid(t+1)=pid(t)±α(t)×|pmbestid(t)-

Xid(t)|×ln(1/uid(t))

(7)

式中,α稱為收縮擴張系數(shù)。

粒子的當(dāng)前最優(yōu)位置Pi與全局最優(yōu)位置Pg的更新方式為:

(8)

式中,f()為目標(biāo)函數(shù)。

在QPSO算法中,唯一控制參數(shù)為式(9)中收縮擴張系數(shù)α,它在算法中起到的作用是能夠調(diào)整算法的收斂過程?？刂茀?shù)α最簡單的方法是固定α的值,但會降低算法性能。當(dāng)前較好的方法是線性變換α,如式(9)所示,但性能提升有限。

α(t)=m-(m-n)*t/Maxtimes

(9)

式中,m、n為常數(shù),Maxtimes為最大迭代數(shù)。

為了使QPSO算法具有更好的適應(yīng)性,本文提出一種變參數(shù)量子粒子群優(yōu)化(VPQPSO)算法。VPQPSO算法的參數(shù)α定義為:

(10)

式中,使z=lg(ΔF),ΔF為誤差函數(shù),其定義如下。

(11)

式中,f(pi)為Pi的適應(yīng)度值,f(pg)為pg的適應(yīng)度值,Min()為取最小值函數(shù)。

VPQPSO算法的偽代碼為:

初始化種群的每個粒子的位置向量;

do

fori=1to種群規(guī)模m

iff(xi)

Pg=min(Pi);

計算式(8)計算pmbest的值;

ford=1to維數(shù)d

u=rand(0,1);

ifrand(0,1)<0.5

xid(t+1)=Pid-a(t)*|pmbestid-

xid(t)|*ln(1/uid(t));

else

xid(t+1)=Pid+a(t)*|pmbestid(t)-

xid(t)|*ln(1/uid(t));

直到終止條件滿足。

選用3個常用基準(zhǔn)函數(shù)對VPQPSO算法和QPSO算法性能進(jìn)行對比,具體如下:

(1)Griewank函數(shù):

(12)

(2)Rastrigin函數(shù):

(13)

(3)Schaffer函數(shù)

(14)

圖1為3個測試函數(shù)適應(yīng)度對數(shù)值進(jìn)化曲線,從圖1可知,對所有函數(shù),VPQPSO的收斂速度明顯優(yōu)于QPSO。對多峰值Rastrigin和強烈振蕩的多峰值Schaffer,VPQPSO能很快達(dá)到理論極小點0和-1,且避免了QPSO陷入局部最優(yōu)的缺點,這表明了通過自適應(yīng)的改變參數(shù)α值,可以使VPQPSO的全局搜索能力、收斂精度和收斂速度均優(yōu)于傳統(tǒng)的QPSO算法。

圖2　QPSO和VPQPSO算法的收斂性能比較

3　仿真實驗

3.1數(shù)據(jù)來源

硬件驗環(huán)境為:Pentium E5300 2.8 GHz,4GRAM平臺上采用MATLAB 2012編程實現(xiàn)仿真實驗。數(shù)據(jù)來自于河南鄭州2013年7月1日至2013年7月30日的小時短期負(fù)荷數(shù)據(jù),共采集到720個短期負(fù)荷數(shù)據(jù),具體如圖3所示,前600個數(shù)據(jù)作為訓(xùn)練集,最后120個數(shù)據(jù)作為測試集。

圖3　原始短期負(fù)荷數(shù)據(jù)

3.2對比模型和評價標(biāo)準(zhǔn)

為了使VPQPSO-BPNN的預(yù)測結(jié)果具有可比性,采用量子粒子群算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)(QPSO-BPNN)作為對比模型。評價標(biāo)準(zhǔn)為平均相對誤差(MAPE)和均方根誤差(RMSE)。它們具體定義為:

(15)

(16)

3.3數(shù)據(jù)預(yù)處理

首先根據(jù)互信息法得到時間延遲(τ)的值在[1,20]區(qū)間內(nèi)的互信息函數(shù),該互信息函數(shù)的第1極小值為最佳時間延遲τ,互信息函數(shù)圖如圖4所示。從圖4可以看出當(dāng)時間延遲為6時,互信息函數(shù)達(dá)到第1極小值,所以短期負(fù)荷時間序列的最佳時間延遲為6。

圖4　互信息函數(shù)圖

然后采用CAO算法計算在不同嵌入維數(shù)下的關(guān)聯(lián)維數(shù),如圖5所示。從圖5可知,隨著嵌入維數(shù)的增加,關(guān)聯(lián)維數(shù)不斷增加,當(dāng)嵌入維數(shù)m=7時,關(guān)聯(lián)維數(shù)達(dá)到飽和狀態(tài),說明短期負(fù)荷時間序列的最佳嵌入維數(shù)為7。

圖5　不同嵌入維數(shù)下的關(guān)聯(lián)維

最后采用τ=6,m=7,根據(jù)式(1)對圖2中的短期電力負(fù)荷數(shù)據(jù)進(jìn)行相空間重構(gòu),得一個多維的短期電力負(fù)荷數(shù)據(jù)。BP神經(jīng)網(wǎng)絡(luò)訓(xùn)練過程中,對0～1間的數(shù)據(jù)最為敏感,為了提高訓(xùn)練效率,對收集短期負(fù)荷數(shù)據(jù)進(jìn)行歸一化處理,具體為:

(17)

式中,x(i)和x′(i)分別表示原始和歸一化后值,Ex和σx分別表示原始時間序列均值和標(biāo)準(zhǔn)差。

3.4結(jié)果與分析

最后50個樣本作為測試樣本,其余的作為訓(xùn)練樣本。VPSO-BPNN和QPSO-BPNN的擬合結(jié)果和預(yù)測結(jié)果分別如圖6和7所示,它們擬合結(jié)果和預(yù)測結(jié)果的誤差見表1。

圖6　PSO-BPNN和CPSO-BPNN的擬合結(jié)果對比

圖7　PSO-BPNN和CPSO-BPNN的預(yù)測結(jié)果對比

表1QPSO-BPNN和VPQPSO-BPNN的性能對比

評價指標(biāo)QPSO?BPNNVPQPSO?BPNN擬合RMSE4．412．01擬合MAPE6．54%3．48%預(yù)測RMSE5．602．52預(yù)測MAPE8．30%4．12%

對圖3、圖4和表2兩種模型對短期負(fù)荷擬合和預(yù)測結(jié)果進(jìn)行分析,可以得到如下結(jié)論:

(1)VPQPSO-BPNN對訓(xùn)練樣本擬合的精度相當(dāng)高,擬合誤差要小于QPSO-BPNN,這表明VPQPSO算法較好的克服了QPSO算法存在的缺陷,提高了短期負(fù)荷預(yù)測模型的擬合精度。

(2)VPQPSO-BPNN對測試樣本集的預(yù)測結(jié)果要優(yōu)于QPSO-BPNN,這說明表明采用VPQPSO-BPNN建立的短期負(fù)荷預(yù)測模型更加全面、準(zhǔn)確的刻畫了短期電力負(fù)荷的時變性、周期性等變化趨勢,泛化能力更優(yōu),從而提高了短期電力負(fù)荷預(yù)測精度。

4　結(jié)語

為了提高短期電力負(fù)荷預(yù)測精度,提出一種基于量子粒子群優(yōu)化算法優(yōu)化BP網(wǎng)絡(luò)參數(shù)的短期電力負(fù)荷預(yù)測模型。仿真結(jié)果表明,相對于其他預(yù)測模型,VPQPSO-BPNN更能映電力負(fù)荷復(fù)雜的變化特點,獲得了更優(yōu)的預(yù)測效果,在非線性電力負(fù)荷預(yù)測有著廣泛的應(yīng)用前景。

參考文獻(xiàn):

[1]康重慶,夏清,張伯明.電力系統(tǒng)負(fù)荷預(yù)測研究綜述與發(fā)展方向的探討[J].電力系統(tǒng)自動化,2004,28(17):1-9.

[2]Shi B,Li Y X,Yu X H,et al.A Modified Particle Swarm Optimization and Radial Basis Function Neural Network Hybrid Algorithm Model and Its Application[C]//2009 WRI Global Congress on Intelligent Systems,2009(1):134-138.

[3]Wang H,Li B S,Han X Y,et al.Study of Neural Networks for Electric Power Load Forecasting[C]//The Third International Symposium on Neural Networks Proceedings,2010:1277-1283.

[4]王濤,彭云華,吳小亞.基于相關(guān)因素修正的短期電力負(fù)荷回歸預(yù)測模型[J].電氣自動化,2013,35(1):65-54.

[5]蔣剛.基于模糊支持向量核回歸方法的短期峰值負(fù)荷預(yù)測[J].控制理論與應(yīng)用,2007,24(6):986-990.

[6]陳國初,劉軍.基于神經(jīng)網(wǎng)絡(luò)的中長期用電量預(yù)測模型[J].上海電機學(xué)院學(xué)報,2009,12(1):20-24.

[7]俞阿龍.基于遺傳算法的RBF神經(jīng)網(wǎng)絡(luò)在渦流傳感器非線性補償中應(yīng)用[J].電子器件,2008,31(3):1039-1041.

[8]師彪,李郁俠,于新花,等.改進(jìn)粒子群-BP神經(jīng)網(wǎng)絡(luò)模型的短期電力負(fù)荷預(yù)測[J].計算機應(yīng)用,2009,29(4):1036-1039.

[9]肖文平,葉家瑋.混沌粒子群優(yōu)化神經(jīng)網(wǎng)絡(luò)算法應(yīng)用于SRG建模[J].計算機工程與應(yīng)用,2010,46(27):238-241.

[10]Coelho L S.Gaussian Quantum Behaved Particle Swarm Optimization Approaches for Constrained Engineering Design Problems[J].Expert Systems with Applications,2010,37(2):1676-1683

[11]魏秀業(yè),潘宏霞.粒子群優(yōu)化及智能故障診斷[M].北京:國防工業(yè)出版社,2010:17-18.

[12]Sun J,Fang W,Xu X J,et al.Quantum-Behaved Particle Swarm Optimization:Analysis of the Individual Particle’s Behavior and Parameter Selection[J].Evolutionary Computation,2012,20(3):349-393.

[13]許少華,王皓,王穎,告示.一種改進(jìn)的量子粒子群優(yōu)化算法及其應(yīng)用[J].計算機工程與應(yīng)用,2011,47(20):34-36.

王超(1975-),男,碩士,講師,研究方向計算機網(wǎng)絡(luò)與應(yīng)用技術(shù),wangchao770212@sohu.com;

李森(1980-),男,碩士,講師,研究方向計算機網(wǎng)絡(luò)與軟件。

Short-TermLoadForecastingBasedonNeuralNetworkOptimizedbyVaryingParametersQuantumParticleSwarmOptimizationAlgorithm

WANGChao*,LISen

(Henan Polytechnic Institute,Nanyang He’nan 473000,China)

Abstract:In order to forecast short-term load accurately and quickly,a short load forecasting model is proposed,based on BP neural network optimized by quantum particle swarm optimization algorithm.Firstly,the data of short load are reconstructed by chaotic theory,and then the parameters of BPNN were considered the position vector of quantum particle,the optimal parameters of BPNN are found by quantum particle optimization algorithm to reduce blindness and inefficiency.Lastly,the optimal model for network traffic is built and the performance of mode are tested by short load data.The simulation results show that VPQPSO algorithm has solved the problems of the BP neural network,and the proposed model improved the power load forecasting accuracy.

Key words:power load;forecasting accuracy;quantum particle swarm optimization algorithm;BP neural network;varying parameters

doi:EEACC:815010.3969/j.issn.1005-9490.2014.04.043

中圖分類號:TP391

文獻(xiàn)標(biāo)識碼:A

文章編號:1005-9490(2014)04-0782-05

收稿日期:2013-08-18修改日期:2013-09-22