13節(jié)點(diǎn)配電網(wǎng)潮流計(jì)算

1、配網(wǎng)前推后代潮流計(jì)算對(duì)給定的13節(jié)點(diǎn)配電網(wǎng)，熟悉MATLAB語言前腿后代潮流計(jì)算程序。13節(jié)點(diǎn)配電網(wǎng)結(jié)構(gòu)如圖1所示，支路參數(shù)見表1;節(jié)點(diǎn)類型及迭代初始電壓及各節(jié)點(diǎn)所接負(fù)荷見表2。要看懂前推回代法計(jì)算程序，報(bào)告敘述計(jì)算原理及計(jì)算流程。繪制計(jì)算流程框圖。確定前推回代支路次序（廣度優(yōu)先，或深度優(yōu)先），編寫前推回代計(jì)算輸入文件。進(jìn)行潮流計(jì)算。整理專利計(jì)算結(jié)果。內(nèi)容包括節(jié)點(diǎn)電壓；線路潮流分布，網(wǎng)損總結(jié)配電網(wǎng)潮流計(jì)算的要點(diǎn)1。O0O9-o1262810311o6O135。圖113節(jié)點(diǎn)配電網(wǎng)結(jié)構(gòu)圖表1系統(tǒng)支路i參數(shù)支路R(Q.)X(Q.)B/2(Tk)(S)123.3673.6850.0232.3562.5

2、410.0341.1451.280.0454.5245.040.0260.8561.140.0672.7452.9650.0283.7434.2510.0892.2372.7560.03104.144.6960.03111.3281.7630.011122.4362.8660.04133.5213.9660.0表2系統(tǒng)負(fù)荷參數(shù)節(jié)點(diǎn)編號(hào)節(jié)點(diǎn)類型節(jié)點(diǎn)初始電壓(kV)P(MVAQ(MVA1根節(jié)點(diǎn)10.4002中間節(jié)點(diǎn)10.00.03420.03013中間節(jié)點(diǎn)10.00.06930.06424中間節(jié)點(diǎn)10.00.08450.07635葉節(jié)點(diǎn)10.00.02950.02616中間節(jié)點(diǎn)10.00.0474

3、0.04097葉節(jié)點(diǎn)10.00.11760.09578中間節(jié)點(diǎn)10.00.09460.08579葉節(jié)點(diǎn)10.00.09160.085910葉節(jié)點(diǎn)10.00.02710.022911中間節(jié)點(diǎn)10.00.06960.064312葉節(jié)點(diǎn)10.00.06760.057913葉節(jié)點(diǎn)10.00.02980.0242主程序活單:PQ,FT,RX=case113();NN=size(PQ,1);NB=size(FT,1);V=PQ(:,1);%調(diào)用數(shù)據(jù)文件%節(jié)點(diǎn)數(shù)%支路數(shù)數(shù)%V初始電壓相量maxd=1k=1whilemaxd0.0001%每一次迭代各節(jié)點(diǎn)的注入有功和無功相同PQ2=PQ;PL=0.0;fori

4、=1:NBkf=FT(i,1);%前推始節(jié)點(diǎn)號(hào)kt=FT(i,2);%前推終節(jié)點(diǎn)號(hào)%計(jì)算沿線電流平方A%計(jì)算支路首端有功/MWRX(i,1)R%計(jì)算沿支路的無功損耗/MvarRX(i,2)X%用PQ1去修正支路末端節(jié)點(diǎn)的有功P單位%用PQ1去修正支路末端節(jié)點(diǎn)的有功Q單位x=(PQ2(kf,2)A2+PQ2(kf,3)A2)/V(kf)/V(kf);PQ1(i,1)=PQ2(kf,2)+RX(i,1)*x;PQ1(i,2)=PQ2(kf,3)+RX(i,2)*x;PQ2(kt,2)=PQ2(kt,2)+PQ1(i,1);MWPQ2(kt,3)=PQ2(kt,3)+PQ1(i,2);MvarPL=

5、PL+RX(i,1)*x;endangle(1)=0.0;fori=NB:-1:1kf=FT(i,2);%回代始節(jié)點(diǎn)號(hào)kt=FT(i,1);%回代終節(jié)點(diǎn)號(hào)%計(jì)算支路電壓損耗的縱分量dv1%計(jì)算支路電壓損耗的橫分量dv2%計(jì)算支路末端電壓/kV%計(jì)算支路dv1=(PQ1(i,1)*RX(i,1)+PQ1(i,2)*RX(i,2)/V(kf);dv2=(PQ1(i,1)*RX(i,2)-PQ1(i,2)*RX(i,1)/V(kf);V2(kt)=sqrt(V(kf)-dv1)A2+dv2A2);angle(kt)=angle(kf)+atand(dv2/(V(kf)-dv1);endmaxd=ab

6、s(V2-V(2);V2(1)=V(1);fori=3:1:NNifabs(V2(i)-V(i)maxd;maxd=abs(V2(i)-V(i);endendmaxdk=k+1PQ1%潮流分布即支路首端潮流MVAV=V2%節(jié)點(diǎn)電壓模計(jì)算結(jié)果kVangle%節(jié)點(diǎn)電壓角度計(jì)算結(jié)果單位度PL%網(wǎng)損單位MWendclear輸入文件清單：functionPQ,FT,RX=case113()PQ=%節(jié)點(diǎn)電壓有功無功10.50010.0.60.45100.40.3100.40.28100.60.4100.40.3100.50.35100.50.4;FT=%首端末端355322221;RX=%RX0.60.351.00.550.650.350.620.50.720.750.900.50.540.65；計(jì)算結(jié)果清單：k=7PQ1=0.40140.28080.40250.30140.50250.35130.80740.58571.52941.07820.50360.40203.53562.6303V=10.500010.