【正文】
⑵對第i列的元素進行消去運算:(p=k+1,k+2,…n)然后按照步驟⑴、 ⑵對其進行矩陣的消去運算,直到對所有的列進行消去運算為止。 二范數(shù)的介紹及應(yīng) 二范數(shù)的介紹 范數(shù)為在向量空間內(nèi)對所有的向量賦一個非零的大小。通過以上的對系數(shù)矩陣在形成的過程中的簡化方法,來改變PQ算法比值的敏感性問題。發(fā)現(xiàn)在用BX法進行改進的算法中結(jié)果收斂,在此基礎(chǔ)上加入二范數(shù)時計算的速度和迭代次數(shù)都有所改善??梢詮囊韵聨讉€方面入手:①對于BX法可以在形成系數(shù)矩陣的時候,對考慮電阻和忽略電阻進行不同的組合,找出收斂情況最好的組合:②可以使用支路分解法來解決PQ算法的收斂性問題【15】;③可以更多的參考國內(nèi)外的參考文獻,研究出新的算法來改善Q算法的不足。%isb為平衡節(jié)點編號[pr]=xlsread(39。,39。H3:H100039。N3:N100039。)。%節(jié)點類型(1平衡節(jié)點,2PQ節(jié)點,) B2=[(y(:,1)+y(:,2)*i),(y(:,3)+y(:,4)*i),y(:,5),y(:,6),y(:,7)]。)。%Y(p,p)=1./X(i,2)。%Y的實部虛部,YI的虛部。%B(i,j1)是從第二行起的上三角矩陣,j1iendB(i,i)=1./B(i,i)。ICT1=0。 end end %DET(K)=max(DE)。 for i=1:n if B2(i,5)==2 %i節(jié)點為PQ節(jié)點 C(i)=0。 DQ(I)=DQ(I)A(I,L)*DQ(L)。39。)。 %修正電壓模值 end end kp=1。 end end % DET1(K)=max(DE1)。 for k=IC1:n DP(k)=DP(k)B(k,i)*DP(i)。ICT3=1。k=0。%發(fā)電機功率節(jié)點負荷功率=線路功率endP=real(S)。%求Y矩陣互導(dǎo)納YI(p,q)=YI(p,q)1./B1(i,3)。 tic。%節(jié)點編號%z(:,2)=xlsread(39。%節(jié)點電壓的初始值(Vm)y(:,6)=xlsread(39。)。)。F3:F100039。39。%n為節(jié)點數(shù) [nl]=xlsread(39。 將本章中所有表格之間的對比,可看出如改進PQ算法,配網(wǎng)中所得潮流計算結(jié)果不收斂。使用一般的PQ法進行潮流計算,結(jié)果不收斂。以下是三者的具體運算公式: 1范數(shù):設(shè)有向量X=(),;2范數(shù):設(shè)有向量X=(),; 范數(shù):設(shè)有向量X=()。因此在解方程的時候首先對這兩個矩陣進行因子表化。直接法經(jīng)過有限次的算術(shù)運算就能得到答案,運算的次數(shù)與采用的計算方法和方程組的階數(shù)有關(guān)。對于狀態(tài)量的約束條件,也就是說的數(shù)值不宜過大,而且電納遠遠大于電導(dǎo)所以可以認為: (26) (27)最終可以將公式將化為: (28) (29)從上式可以看出在雅可比矩陣中的兩個子陣H、L中的元素有相同的表達式,但是二者的階數(shù)不同,前者為(n1)階,后者是(m1)階。因此對該算法需要進行改進。不同容量的分布式電源提出了一種基于網(wǎng)絡(luò)結(jié)構(gòu)的分布式松弛母線結(jié)構(gòu)的潮流算法,這種算法適用于含有多個分布式電源的網(wǎng)絡(luò),并且在IEEE14節(jié)點的網(wǎng)絡(luò)中得以驗證其有效性。傳統(tǒng)的發(fā)電機的模型由PV、PQ、和平衡節(jié)點,但是分布式電源的節(jié)點類型需要取何種類型需要視具體情況而定。⑶供電質(zhì)量提高。針對PQ算法用于配電網(wǎng)的不足提出了改進的方法(BX法),并且在此基礎(chǔ)上加入了二范數(shù)。但整體上來說,分布式電源是規(guī)模比較小的(一般的發(fā)電功率為數(shù)千瓦到50MW),能夠單獨的輸出電能且分布位置一般在用戶附近的電能系統(tǒng)。但是當分布式電源在并入電網(wǎng)的時候也會給電網(wǎng)帶來影響,因此研究分布式電源的潮流計算是很有必要的。文獻【3】對PQ算法進行了改進,并且指出了在形成系數(shù)矩陣的時候考慮支路電阻在形成系數(shù)矩陣的時候忽略支路電阻這樣的組合使得PQ算法的收斂性得到了很大的提高,在本次畢業(yè)設(shè)計中就用的該理論。之所以會出現(xiàn)不收斂的情況是因為當比值很大時,分解算法的假設(shè)條件均不成立,因此分解法失去了簡化條件的基礎(chǔ),導(dǎo)致出現(xiàn)了對比值的敏感性。為了建立修正方程必須對其節(jié)點進行約定:①假設(shè)網(wǎng)絡(luò)中有n個節(jié)點,它們的編號為1,2,3,…,n,在這些節(jié)點中包含一個平衡節(jié)點為d;②在這網(wǎng)絡(luò)中有(a1)個PV節(jié)點,編號為1,2,3,…,a,在這中間包含平衡節(jié)點d;③在此網(wǎng)絡(luò)中有(na)個PQ節(jié)點,編號為a+1,a+2, …,n.由于平衡節(jié)點的注入功率不能事先給定,因此也就不能列出的表達式;平衡節(jié)點的電壓;也不需要求取。這樣就相當于在用戶側(cè)有了電源,潮流的方向也就不想以前那樣單向的由電源側(cè)流向負荷側(cè),有可能就會出現(xiàn)電流回流的情況,因此需要對配電網(wǎng)的潮流在進行重新分布。其運算過程可以分為常數(shù)項的消去和系數(shù)矩陣的消去,而這兩者的消去過程是分開進行的。然后進行有功功率不平衡量的判斷,看其是否滿足精度要求,滿足則進行其他配電網(wǎng)中的其他的參數(shù)的運算,不滿足則形成新的電壓相角,并利用新的電壓相角參數(shù)來進行下一次的迭代計算。原本的兩個系數(shù)矩陣的參數(shù)是電網(wǎng)中的導(dǎo)納矩陣的虛部。九節(jié)點網(wǎng)絡(luò)拓撲結(jié)構(gòu)見附錄。在十四節(jié)點的網(wǎng)絡(luò)中這種方法并沒有使結(jié)果收斂,因此本次設(shè)計還有不足,有待提高。,39。)。%支路電抗(x)x(:,5)=xlsread(39。%由支路參數(shù)形成的矩陣y(:,1)=xlsread(39。39。,39。39。O=zeros(1,n)。%自導(dǎo)納=節(jié)點對地自導(dǎo)納+支路導(dǎo)納+線路對地導(dǎo)納 end %求導(dǎo)納矩陣%*********************求節(jié)點導(dǎo)納矩陣Y**********************% G=real(Y)。%相角endend%*****************以上是求PQ節(jié)點電壓初值,各節(jié)點電壓相角初值********for i=2:nif i==nB(i,i)=1./B(i,i)。endA(i,i)=1./A(i,i)。 end DP1(i)=P(i)V(i)*C(i)。這里的DP為Δθi end end end end %=========================== 對角度進行修正 ================= for i=2:n O(i)=O(i)DP(i)。 %規(guī)格化,求ΔVi if i==na for LZ=2:i L=i+2LZ。etime程序總運行時間:39。,K)。 %無功功率都滿足要求 if kp~=0 %有功功率存在不滿足的節(jié)點 K=K+1。 %本次迭代無功不滿足的節(jié)點個數(shù) end else kp=0。 IC4=L1。%*********************迭代次數(shù)K置零,Kp=1,Kq=1開始***************** while (ICT2~=0|ICT3~=0) ICT2=0。%PQ節(jié)點個數(shù)A(p,k)=BI(i,j1)。%節(jié)點電壓初始值的實部f(i)=imag(B2(i,3))。%求YI矩陣互導(dǎo)納 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。YI=zeros(n)。B9:B100039。S3:S100039。,39。,39。39。x(:,1)=xlsread(39。B2:B239。最重要的就是其在電網(wǎng)發(fā)生故障的時候可以為那些重要的負荷供電,減少電力事故帶來的損失。說明在五節(jié)點的配電網(wǎng)中,BX法對于算法的改進是很適用的。因此利用二范數(shù)來處理所得到的,利用二范數(shù)小于無窮范數(shù)來減小迭代次數(shù),加快潮流計算的時間并且取到了較好的效果。有功計算的公式: (31)其中利用式(33)可以求得,在最初的時候是一直或者自己設(shè)定的,設(shè)定的一般都為1。迭代法主要用于解非線性方程組。 必須指出雖然PQ分解法做了很多的簡化步驟,但是并沒有對系統(tǒng)最后的算法的精確度有所影響,因為其收斂判據(jù)依然是無功功率的不平衡量和有功功率的不平衡量都要滿足題目所要求的精度。在最開始的階段我們假設(shè)全網(wǎng)的電壓為額定的電壓,由于已知線路末端的負荷功率,可以根據(jù)公式(21)由末端向始端進行推算,在此過程中只計算各元件的功率損耗而不計算節(jié)點的電壓,可以求出各個支路的電流和功率損耗,最終求得線路的始端功率,這是其回代過程:然后再根據(jù)求得的始端功率和已知的電壓,由始端想末端根據(jù)公式(22)進行電壓降落的推算,從而求得各個節(jié)點的電壓,這是前推的過程。 文獻【14】介紹了分布式電源的處理方法以及分布式電源加在不同的位置上對電網(wǎng)的影響 本文的主要內(nèi)容 在摘要中已經(jīng)指出本次畢業(yè)設(shè)計采用的潮流計算方法是PQ分解法,它是基于牛頓—拉夫遜法所提出來的。此外,王成山、鄭海峰等人對于風力發(fā)電、太陽能發(fā)電等隨機性的發(fā)電模型提出了相應(yīng)的潮流計算方法,運用了概率統(tǒng)計方法處理系統(tǒng)中的隨機變化的因素,給出系統(tǒng)運行電壓、支路潮流等概率分布情況,能夠更清楚地表明系統(tǒng)的運行狀況,為系統(tǒng)的安全運行提供更加完整的數(shù)據(jù)信息。在分布式電源中一般都會采用較為先進的控制設(shè)備,性能比較好,因此其在開機或停機的時候易于操作,在負荷調(diào)節(jié)的方面也比較靈活。關(guān)鍵詞:分布式電源、PQ算法、配電網(wǎng)、BX法AbstractWith rapid development of the power system, the application of distributed generator supply bee ordinary. Power grid with distributed generator will influence the voltage ,network loss of the power distribution the research of power grid with distributed generator ‘s power flow calculation is particularly important.In this article, the definition of distributed generator supply and the situation of both at home and abroad is produced firstly. Used for mathematical algorithm are introduced, and then summarizes the mon trend of the algorithm Aiming at the shortings of the PQ algorithm used in distribution network is put forward to improve the method, the method (BX). For 5, 9, 14 and 33 nodes of network has carried on the example analysis, in 5 and 9 nodes in the network to improve the algorithm can make the PQ good convergence, but at 14, 33 nodes in the network effect is not obvious, so the algorithm needs to be impr