【正文】 fprintf(fid,39。迭代次數(shù)k為:%d\n39。)。fprintf(fid,39。wt39。39。,num2str(etime(clock,t1))])。pause(3*rand)disp([39。 endendend%以下是輸出數(shù)據(jù)，注意涉及k，k=1:n時，把fprintf放在循環(huán)中，否則只會讀出k=n時的數(shù)據(jù)tic。else kq=0。 %修正電壓模值 end end kp=1。 for i=1:n if B2(i,5)==2 %i節(jié)點為PQ節(jié)點 L=L+1。 for k=IC1:na DQ(k)=DQ(k)A(k,i)*DQ(i)。 DQ(I)=DQ(I)A(I,L)*DQ(L)。 IC4=L1。 for i=1:na DQ(i)=A(i,i)*DQ(i)。 end %=============== 用因子表求解修正方程的修正量電壓模 =============== Nq(K)=ICT3。 if DET1(K)=pr %節(jié)點功率誤差是否滿足要求 ICT3=ICT3+1。 end end % DET1(K)=max(DE1)。 %節(jié)點功率Q計算,ΔQ DQ(L)=DQ1(i)./V(i)。 for k=1:n %所有列（包括PV節(jié)點及平衡節(jié)點列） C(i)=C(i)+V(k)*(G(i,k)*sin(O(i)O(k))BI(i,k)*cos(O(i)O(k)))。 for i=1:n if B2(i,5)==2 %i節(jié)點為PQ節(jié)點 C(i)=0。 end %================= 對無功Q及電壓U的關(guān)系進行解算 =========== kq=1。 %按列回代。 for MZ=2:IC4 I=IC4+2MZ。這里的DP為Δθi end else for LZ=3:i L=i+3LZ。 for k=IC1:n DP(k)=DP(k)B(k,i)*DP(i)。 if ICT2~=0 for i=2:n DP(i)=B(i,i)*DP(i)。 if DET(K)=pr %節(jié)點功率誤差是否滿足要求 ICT2=ICT2+1。 end end %DET(K)=max(DE)。 %節(jié)點功率P計算 DP(i)=DP1(i)./V(i)。 for k=1:n %所有列（包括PV節(jié)點及平衡節(jié)點列） C(i)=C(i)+V(k)*(G(i,k)*cos(O(i)O(k))+BI(i,k)*sin(O(i)O(k)))。ICT3=0。ICT3為不滿足無功功率精度值的個數(shù)。ICT3=1。DET=0。kq=1。ICT1=0。for k=i+1:nafor j1=i+1:naA(k,j1)=A(k,j1)A(k,i)*A(i,j1)。for j1=k:naA(i,j1)=A(i,j1)./A(i,i)。%這里的A為系數(shù)矩陣B″end endendend%*********************以上是將Y矩陣中的元素存在A矩陣中**************for i=1:naif i==naA(i,i)=1./A(i,i)。%J1為PQ節(jié)點時k=k+1。k=0。for i=1:nif B2(i,5)==2。%消去運算endendendend%**********以上是用YI矩陣形成因子表，并求逆陣（計算有功功率）**********p=0。%B（i,j1）是從第二行起的上三角矩陣，j1iendB(i,i)=1./B(i,i)。%去倒數(shù)并賦值else IC1=i+1。%該節(jié)點的幅值=節(jié)點電壓初始值的實部和虛部的算術(shù)平方根；因為此節(jié)點不是PV節(jié)點而是PQ節(jié)點，所以第75行求Vi的式子求不出電壓來，所以需此行求電壓O(i)=atan(f(i)./e(i))。%節(jié)點電壓初始值的虛部V(i)=B2(i,4)。for i=1:ne(i)=real(B2(i,3))。%發(fā)電機功率節(jié)點負荷功率=線路功率endP=real(S)。BI=imag(Y)。G=real(Y)。%Y的實部虛部，YI的虛部。B=imag(YI)。%自導(dǎo)納=節(jié)點對地自導(dǎo)納+支路導(dǎo)納+線路對地導(dǎo)納YI(p,p)=YI(p,p)+1./B1(i,3)。%對變壓器支路求自導(dǎo)納YI(q,q)=YI(q,q)+1./B1(i,3)。%求Y矩陣互導(dǎo)納YI(q,p)=YI(p,q)。%求Y矩陣互導(dǎo)納YI(p,q)=YI(p,q)1./B1(i,3)。q=B1(i,1)。q=B1(i,2)。%Y(p,p)=1./X(i,2)。%for i=1:n%if X(i,2)~=0。V=zeros(1,n)。e=zeros(1,n)。 %*********************求節(jié)點導(dǎo)納矩陣Y**********************Y=zeros(n)。 tic。B1(:,3)=B1(:,3)/ZB。SB=10。)。,39。%由節(jié)點號及其對地阻抗形成的矩陣[na]=xlsread(39。)。,39。%節(jié)點編號%z(:,2)=xlsread(39。A9:A100039。39。%節(jié)點類型（1平衡節(jié)點，2PQ節(jié)點，） B2=[(y(:,1)+y(:,2)*i),(y(:,3)+y(:,4)*i),y(:,5),y(:,6),y(:,7)]。T3:T100039。39。)。,39。%節(jié)點電壓的初始值(Vm)y(:,6)=xlsread(39。R3:R100039。39。)。,39。%節(jié)電負荷的有功功率(Pd)y(:,4)=xlsread(39。P3:P100039。39。)。,39。%節(jié)點所接發(fā)電機的有功功率(Pg)y(:,2)=xlsread(39。N3:N100039。39。 B1=[x(:,1),x(:,2),(x(:,3)+x(:,4)*i),x(:,5)*i,x(:,6),x(:,7)]。J3:J100039。39。)。,39。%支路對地導(dǎo)納(b)x(:,6)=xlsread(39。H3:H100039。39。)。,39。%支路電阻(r)x(:,4)=xlsread(39。F3:F100039。39。)。,39。%起始節(jié)點編號(fbus)x(:,2)=xlsread(39。D3:D100039。39。 %把矩陣數(shù)值輸?shù)紼xcel里，然后從Matlab中調(diào)出該Excel里的矩陣數(shù)值file=[pname fname]。39。)。,39。%isb為平衡節(jié)點編號[pr]=xlsread(39。A4:A439。39。)。,39。%n為節(jié)點數(shù) [nl]=xlsread(39。A2:A239。39?？梢詮囊韵聨讉€方面入手：①對于BX法可以在形成系數(shù)矩陣的時候，對考慮電阻和忽略電阻進行不同的組合，找出收斂情況最好的組合：②可以使用支路分解法來解決PQ算法的收斂性問題【15】；③可以更多的參考國內(nèi)外的參考文獻，研究出新的算法來改善Q算法的不足。 在本次的畢業(yè)設(shè)計中，改進的PQ算法雖然在五節(jié)點和九節(jié)點的配電網(wǎng)絡(luò)中收斂的效果比較好，但是在十四節(jié)點以及更大的配電網(wǎng)絡(luò)中仍然會出現(xiàn)不收斂的問題。但是。本次設(shè)計的主要目的就上改善PQ算法，使該算法能夠應(yīng)用于配電網(wǎng)中，由于該算法對于電阻與電抗比值的敏感性，所以它并不適用于配電網(wǎng)的潮流計算中，但由于該算法的快速性以及減少容量要求，所以就意圖改善它的不足，使其能夠應(yīng)用于配電網(wǎng)的潮流計算中。第五章 總結(jié)與展望 總結(jié)隨著智能電網(wǎng)的發(fā)展，分布式電源必然會在配電網(wǎng)的應(yīng)用中越來越廣泛，它的應(yīng)用有利于改善環(huán)境、節(jié)約能源，符合可持續(xù)發(fā)展的國策。 將本章中所有表格之間的對比，可看出如改進PQ算法，配網(wǎng)中所得潮流計算結(jié)果不收斂。由于時間的原因，在程序中也就沒能夠加入新的算法，所以十四節(jié)點結(jié)果的網(wǎng)絡(luò)未能夠得到改善。 十四節(jié)點算例結(jié)果分析 在十四節(jié)點的網(wǎng)絡(luò)中，其系統(tǒng)的基準電壓，系統(tǒng)的基準容量為,。發(fā)現(xiàn)在用BX法進行改進的算法中結(jié)果收斂，在此基礎(chǔ)上加入二范數(shù)時計算的速度和迭代次數(shù)都有所改善。下面是對其電網(wǎng)的線路損耗和電壓以及運算結(jié)果的情況。 九節(jié)點算例結(jié)果分析在9節(jié)點的配電網(wǎng)絡(luò)系統(tǒng)中，節(jié)點1為平衡節(jié)點，其余的節(jié)點都是PQ節(jié)點，其中分布式電源加在節(jié)點2上面。而且在二范數(shù)的加入PQ算法的時候，比較兩組數(shù)據(jù)明顯可以看出所得出的結(jié)果幾近相似，而且所得的總網(wǎng)損的結(jié)果是完全相同，因此二范數(shù)的加入并沒有結(jié)果的精度。而從表43中可以看出當使用BX法對PQ算法進行改進的時候，算法的結(jié)果收斂，然后在此基礎(chǔ)上加入二范數(shù)的時候，迭代的次數(shù)有所減少，并且計算的速度也加快了。使用一般的PQ法進行潮流計算，結(jié)果不收斂。五節(jié)點網(wǎng)絡(luò)拓撲結(jié)構(gòu)見附錄。第四章 潮流計算結(jié)果的分析 開始 讀取數(shù)據(jù) 形成系數(shù)矩陣 數(shù)據(jù)標幺化形成因子表 計算有有功不平衡量Y 有功滿足精度要求N 修正電壓相角 計算無功不平衡量N 無功滿足要求Y 計算電壓、網(wǎng)損等 結(jié)束 圖41 PQ算法的總流程圖 五節(jié)點的算例結(jié)果分析在5節(jié)點配電網(wǎng)的計算中，節(jié)點1位平衡節(jié)點，其余的節(jié)點都是PQ節(jié)點。通過以上的對系數(shù)矩陣在形成的過程中的簡化方法，來改變PQ算法比值的敏感性問題。而在本次的設(shè)計中分別對這個系數(shù)矩陣進行了不同的處理。在此設(shè)計中采用的BX法其實是指在PQ算法運算的最初階段形成系數(shù)矩陣和的時候?qū)λ木仃噮?shù)進行改進【8】。 BX法的介紹及應(yīng)用 在本次的畢業(yè)設(shè)計中是利用PQ算法來對配電網(wǎng)進行計算，但是一般說來PQ算法并不適應(yīng)于配電網(wǎng)的潮流計算，故在此不再進行說明。 二范數(shù)的應(yīng)用 本次算法中是討論利用PQ算法進行配電網(wǎng)的潮流計算的時候，所得到的結(jié)果是否收斂，判斷的依據(jù)是,原有的計算中是在一組解中找出最大的然后進行判斷，這樣做固然可以使所有的滿足精度要求，但是誤差太大。以下是三者的具體運算公式： 1范數(shù)：設(shè)有向量X=（）,；2范數(shù)：設(shè)有向量X=（）,； 范數(shù)：設(shè)有向量X=（）。在不同的情況下可以根據(jù)定義設(shè)置不同的范數(shù)。由于本次的設(shè)計正是有關(guān)算法的收斂性問題，因此才加入了二范數(shù)的概念，在本文的算法驗證中也體現(xiàn)了二范數(shù)對于結(jié)果的收斂性有所改善。 二范數(shù)的介紹及應(yīng) 二范數(shù)的介紹 范數(shù)為在向量空間內(nèi)對所有的向量賦一個非零的大小。在無功功率不平衡量的計算中與其原理是相同的。相當于等式右邊的常數(shù)項，相當于等式左邊的系數(shù)矩陣，()過程如下（為了方便起見用來表示）：設(shè)系數(shù)矩陣 ， ⑴對系數(shù)矩陣中的元素進行規(guī)格化: (32) 其中j=(i+1),(i+2),…,n (33) ⑵對矩陣中的元素進行按列消去運算：（p=k+1，k+2,…n）；這兩步使矩陣,其中 （34） （35） （36） ⑶用變化后的系數(shù)矩陣的下三角元素對常數(shù)項矩陣中的元素進行消去運算：如對第i個元素進行消去運算，則 (37) (38)根據(jù)式（37）、（3—8）對矩陣中的每一個元素進行消去運算； ⑷對變化后的矩陣用矩陣中的元素進行回代運算， (39) (310)根據(jù)以上（39）、（310）則可求出未知數(shù),即。而線路電阻、電抗等參數(shù)所形成的的系數(shù)矩陣，幾位每次迭代時所求得的電壓相角。下面利用求取有功的不平衡量來進行舉例說明。因此在解方程的時候首先對這兩個矩陣進行因子表化。而本次畢業(yè)設(shè)計的主要任務(wù)就是改善PQ算法的收斂性問題，因此采用了因子表法。 因子表法的應(yīng)用在本次的畢業(yè)設(shè)計中最主要的程序就是因子表的建立，然后利用因子表來解方程組。設(shè)方程組為AX=B，其中矩陣為，消去步驟為： ⑴首先對第i行的元素來進行規(guī)格化：, （其中