【文章內(nèi)容簡介】 數(shù)據(jù)結(jié)構(gòu)主要是矩陣，另外功率方程里的求和部分也不像WORD里那樣，可以用公式編輯器寫出，得用循環(huán)實(shí)現(xiàn)。說它不是難得無法下手，是因?yàn)橥ㄟ^我們查閱資料和自己調(diào)試，最后完成了潮流計(jì)算的程序并得到了收斂的結(jié)果。通過這個(gè)任務(wù)，自己在Matlab編程，潮流計(jì)算，WORD文檔的編輯等方面均有提高，但也暴露出了一些問題：理論知識儲(chǔ)備不足，對Matlab的性能和特點(diǎn)還不能有一個(gè)全面的把握，對WORD軟件也不是很熟練，相信通過以后的學(xué)習(xí)能彌補(bǔ)這些不足，達(dá)到一個(gè)新的層次。參考文獻(xiàn)[1]黃平，孫德寶．基于低壓電力載波遠(yuǎn)程集中抄表系統(tǒng)的研制與實(shí)現(xiàn)[J]．交通與計(jì)算機(jī)，2004，22(1)：99101．Huang Ping，Sun Debao．Study of remote collect meter reading system based on power carrier wave of low voltage[J]．Traffic and Computer，2004，22(1)：99101(in Chinese)．[2]謝開貴，周家啟．樹狀網(wǎng)絡(luò)潮流計(jì)算的新算法[J]．中國電機(jī)工程 學(xué)報(bào)，2001，21(9)：116120．Xie Kaigui，Zhou Jiaqi．A new load flow algorithm for radial distribution networks[J]．Proceeding of the CSEE，2001，21(9)： 116120(in Chinese)．[3]張堯，王琴，宋文南，等．樹狀網(wǎng)的潮流算法[J]．中國電機(jī)工程 學(xué)報(bào)，1998，18(3)：217220．Zhang Yao，Wang Qin，SongWennan，et al．A load flow algorithm for radial distribution power networks[J]．Proceeding of the CSEE，1998，18(3)：217220(in Chinese)．[4]車仁飛，李仁?。?一種少環(huán)配電網(wǎng)三相潮流計(jì)算新方法[J]．中國 電機(jī)工程學(xué)報(bào)，2003，23(1)，7479．Che Renfei，Li Renjun．A new threephase power flow method for weakly meshed distribution systems[J]．Proceeding of the CSEE，2003，23(1)，7479(in Chinese)．[5]Cheng Carol S，Shirmohammadi Dariush．A threephase power flow method for realtime distribution system analysis[J]．IEEE Transaction on Power Systems，1995，10(2)：671679．[6]Zimmerman Ray D，Chiang Hsiaodong．Fast decoupled power flow for unbalanced radial distribution systems[J]．IEEE Transaction on Power Systems，1995，10(4)：20452051．附錄計(jì)算潮流的程序：%本程序的功能是用牛頓拉夫遜法進(jìn)行潮流計(jì)算 n=input(39。請輸入節(jié)點(diǎn)數(shù):n=39。)。nl=input(39。請輸入支路數(shù):nl=39。)。isb=input(39。請輸入平衡母線節(jié)電號:isb=39。)。pr=input(39。請輸入誤差精度:pr=39。)。B1=input(39。請輸入由支路參數(shù)形成的矩陣:B1=39。)。%變壓器側(cè)為1，否則為0 B2=input(39。請輸入各節(jié)點(diǎn)參數(shù)形成的矩陣:B2=39。)。X=input(39。請輸入由節(jié)點(diǎn)號及其對地阻抗形成的矩陣:X=39。)。Y=zeros(n)。U=zeros(1,n)。cta=zeros(1,n)。V=zeros(1,n)。O=zeros(1,n)。S1=zeros(nl)。for i=1:n if X(i,2)~=0。p=X(i,1)。Y(p,p)=X(i,2)。end end for i=1:nlif B1(i,6)==0 p=B1(i,1)。q=B1(i,2)。else p=B1(i,2)。q=B1(i,1)。end Y(p,q)=Y(p,q)1./(B1(i,3)*B1(i,5))。Y(q,p)=Y(p,q)。Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2。Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2。end %求導(dǎo)納矩陣G=real(Y)。B=imag(Y)。for i=1:n cta(i)=angle(B2(i,3))。U(i)=abs(B2(i,3))。%V(i)=B2(i,4)。end for i=1:n S(i)=B2(i,1)B2(i,2)。B(i,i)=B(i,i)+B2(i,5)。end P=real(S)。Q=imag(S)。ICT1=0。IT2=1。while IT2~=0 IT2=0。t1=1。t2=1。for i=1:n if i~=isb C(i)=0。D(i)=0。for j1=1:nC(i)=C(i)+U(i)*U(j1)*(G(i,j1)*cos(cta(i)cta(j1))+B(i,j1)*sin(cta(i)cta(j1)))。D(i)=D(i)+U(i)*U(j1)*(G(i,j1)*sin(cta(i)cta(j1))B(i,j1)*cos(cta(i)cta(j1)))。end DP(t1)=P(i)C(i)。t1=t1+1。if B2(i,6)==2 DQ(t2)=Q(i)D(i)。t2=t2+1。end end end t1=t11。t2=t21。DPQ=[DP39。DQ39。]。%求DP,DQ for i=1:t1+t2 if abs(DPQ(i))pr IT2=IT2+1。end end H=zeros(t1,t1)。N=zeros(t1,t2)。K=zeros(t2,t1)。L=zeros(t2,t2)。for i=1:t1 for j1=1:t1 if j1~=isbamp。j1~=iH(i,j1)=0U(i)*U(j1)*(G(i,j1)*sin(cta(i)cta(j1))B(i,j1)*cos(cta(i)cta(j1)))。elseif j1~=isbamp。j1==i H(i,j1)=U(i)^2*B(i,j1)+D(i)。end end end for i=1:t1 for j1=1:t2 if j1~=isbamp。j1~=iN(i,j1)=0U(i)*U(j1)*(G(i,j1)*cos(cta(i)cta(j1))+B(i,j1)*sin(cta(i)cta(j1)))。elseif j1~=isbamp。j1==i N(i,j1)=0U(i)^2*G(i,j1)C(i)。end end end for i=1:t2 for j1=1:t1 if j1~=isbamp。j1~=i K(i,j1)= U(i)*U(j1)*(G(i,j1)*cos(cta(i)cta(j1))+B(i,j1)*sin(cta(i)cta(j1)))。elseif j1~=isbamp。j1==i K(i,j1)=U(i)^2*G(i,j1)C(i)。end end end for i=1:t2 for j1=1:t2 if j1~=isbamp。j1~=iL(i,j1)=0U(i)*U(j1)*(G(i,j1)*sin(cta(i)cta(j1))B(i,j1)*cos(cta(i)cta(j1)))。elseif j1~=isbamp。j1==i L(i,j1)=U(i)^2*B(i,j1)D(i)。end end end J=[H,N。K,L]。%求雅可比矩陣 modify=JDPQ。Dcta=modify([1:t1],:)。t3=U(:,[1:t2])。DU=diag(t3,0)*modify([t1+1:t1+t2],:)。t4=1。for i=1:t1 if B2(i,6)~=1 cta(1,i)=cta(1,i)+Dcta(t4,1)。t4=t4+1。end end t5=1。for i=1:t2 if B2(i,6)==2 U(1,i)=U(1,i)+DU(t5,1)。t5=t5+1。end end ICT1=ICT1+1。end %修正原值 for i=1:n UU(i)=U(i)*cos(cta(i))+1i*U(i)*sin(cta(i))。end for p=1:n c(p)=0。for q=1:n c(p)=c(p)+conj(Y(p,q))*conj(UU(q))。end s(p)=UU(p)*c(p)。end disp(39。39。)。disp(39。各節(jié)點(diǎn)電壓U為（節(jié)點(diǎn)從小到大排列）:39。)。disp(UU)。disp(39。39。)。disp(39。各節(jié)點(diǎn)電壓相角為（節(jié)點(diǎn)從小到大排列）:39。)。disp(180*angle