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