【正文】 第一篇：電力系統(tǒng)分析潮流計(jì)算大作業(yè)電力系統(tǒng)分析潮流計(jì)算大作業(yè)（源程序及實(shí)驗(yàn)報(bào)告）源程序如下：采用直角坐標(biāo)系的牛頓拉夫遜迭代 function chaoliujisuan()m=3。%m=PQ節(jié)點(diǎn)個(gè)數(shù) v=1。%v=PV節(jié)點(diǎn)個(gè)數(shù)P=[ 0]。%P=PQ節(jié)點(diǎn)的P值 Q=[ 0]。%Q=PQ節(jié)點(diǎn)的Q值 PP=[]。%PP=PV節(jié)點(diǎn)的P值 V=[]。%V=PV節(jié)點(diǎn)的U值E=[1 1 1 ]39。%E=PQ,PV,Vθ節(jié)點(diǎn)e的初值F=[0 0 0 0 0]39。%F=PQ,PV,Vθ節(jié)點(diǎn)f的初值 G=[ 0 0。 0 0。 0 0。0000 0。0000 0]。B=[ 00。 15 0。 0。00 0。000]。Y=G+j*B。X=[]。%X=△X n=m+v+1。%總的節(jié)點(diǎn)數(shù)FX=ones(2*n2,1)。%F(x)矩陣F1=zeros(n1,n1)。%F(x)導(dǎo)數(shù)矩陣 a=0。%記錄迭代次數(shù)EF=zeros(n1,n1)。%最后的節(jié)點(diǎn)電壓矩陣 while max(FX)=10^(5)for i=1:m %PQ節(jié)點(diǎn)FX(i)=P(i)。%△P FX(n+i1)=Q(i)。%△Qfor w=1:n FX(i)= FX(i)E(i)*G(i,w)*E(w)+E(i)*B(i,w)*F(w)F(i)*G(i,w)*F(w)F(i)*B(i,w)*E(w)。%△PFX(n+i1)=FX(n+i1)F(i)*G(i,w)*E(w)+F(i)*B(i,w)*F(w)+E(i)*G(i,w)*F(w)+E(i)*B(i,w)*E(w)。%△Qend end for i=m+1:n1 %PV節(jié)點(diǎn) FX(i)=PP(im)。%△P FX(n+i1)=V(im)^2E(i)^2F(i)^2。%△Qfor w=1:n FX(i)= FX(i)E(i)*G(i,w)*E(w)+E(i)*B(i,w)*F(w)F(i)*G(i,w)*F(w)F(i)*B(i,w)*E(w)。%△Pend endfor i=1:m %PQ節(jié)點(diǎn)for w=1:n1if i~=wF1(i,w)=(G(i,w)*E(i)+B(i,w)*F(i))。F1(i,n+w1)=B(i,w)*E(i)G(i,w)*F(i)。F1(n+i1,w)=B(i,w)*E(i)G(i,w)*F(i)。F1(n+i1,n+w1)=G(i,w)*E(i)+B(i,w)*F(i)。elseF1(i,w)=G(i,i)*E(i)B(i,i)*F(i)。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。F1(n+i1,w)=B(i,i)*E(i)G(i,i)*F(i)。F1(n+i1,n+w1)=G(i,i)*E(i)+B(i,i)*F(i)。for k=1:nF1(i,w)=F1(i,w)G(i,k)*E(k)+B(i,k)*F(k)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。F1(n+i1,w)=F1(n+i1,w)+G(i,k)*F(k)+B(i,k)*E(k)。F1(n+i1,n+w1)=F1(n+i1,n+w1)G(i,k)*E(k)+B(i,k)*F(k)。endendend end for i=m+1:n1 %PV節(jié)點(diǎn)for w=1:n1if i~=wF1(i,w)=(G(i,w)*E(i)+B(i,w)*F(i))。F1(i,n+w1)=B(i,w)*E(i)G(i,w)*F(i)。F1(n+i1,w)=0。F1(n+i1,n+w1)=0。elseF1(i,w)=G(i,i)*E(i)B(i,i)*F(i)。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。F1(n+i1,w)=2*E(i)。F1(n+i1,n+w1)=2*F(i)。for k=1:nF1(i,w)=F1(i,w)G(i,k)*E(k)+B(i,k)*F(k)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。endendend end X=inv(F1)*(FX)。for i=1:n1E(i)=E(i)+X(i)。F(i)=F(i)+X(n+i1)。end a=a+1。fprintf(39。第%d次迭代后的節(jié)點(diǎn)電壓分別為:n39。,a)。disp(E+j*F)。fprintf(39。第%d次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:n39。,a)。disp(FX)。enddisp(39。收斂后的節(jié)點(diǎn)電壓用極坐標(biāo)表示為:39。)。EF=E+j*F。for i=1:n1 fprintf(39。%d號(hào)節(jié)點(diǎn)電壓的幅值為:39。,i)disp(abs(EF(i)))。fprintf(39。%d號(hào)節(jié)點(diǎn)電壓的相角度數(shù)為39。,i)disp(angle(EF(i))*180/pi)。end PPH=0。for i=1:nPPH=PPH+EF(n)*conj(Y(n,i))*conj(EF(i))。end fprintf(39。平衡節(jié)點(diǎn)的功率39。)。disp(PPH)。運(yùn)行結(jié)果：運(yùn)行結(jié)果復(fù)制如下：第1次迭代后的節(jié)點(diǎn)電壓分別為:第1次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:00第2次迭代后的節(jié)點(diǎn)電壓分別為:第2次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:第3次迭代后的節(jié)點(diǎn)電壓分別為:第3次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:第4次迭代后的節(jié)點(diǎn)電壓分別為:第4次迭代后功率偏差△P △Q電壓偏差△V的平方分別為: *收斂后的節(jié)點(diǎn)電壓用極坐標(biāo)表示為: 1號(hào)節(jié)點(diǎn)電壓的幅值為:2號(hào)節(jié)點(diǎn)電壓的幅值為:3號(hào)節(jié)點(diǎn)電壓的幅值為:4號(hào)節(jié)點(diǎn)電壓的幅值為:平衡節(jié)點(diǎn)的功率 第二篇：電力系統(tǒng)分析潮流計(jì)算例題3．1 電網(wǎng)結(jié)構(gòu)如圖3—11所示，其額定電壓為10KV。已知各節(jié)點(diǎn)的負(fù)荷功率及參數(shù)：S2=(+)MVA，S3=(+)MVA，S4=(+)MVAZ12=(+)W，Z23=(+)W，Z24=(+)W試求電壓和功率分布。解：（1）先假設(shè)各節(jié)點(diǎn)電壓均為額定電壓，求線路始端功率。P32++=(R23+jX23)=(1+j2)=+++=(R24+jX24)=(+j3)=+則： S23=S3+DS23=+S24=S4+DS24=+39。=S23+S24+S2=+S12又39。39。P12++=(R12+jX12)=(+)22VN10=+39。故： S12=S12+DS12=+（2）再用已知的線路始端電壓V1=12，求出線路各點(diǎn)電壓。(P12R12+Q12X12)180。+180。=== V2187。V1DV12=(P24R24+Q24X24)DV24==187。V2DV24=(P23R23+Q23X23)DV23==187。V2DV23=（3）根據(jù)上述求得的線路各點(diǎn)電壓，重新計(jì)算各線路的功率損耗和線路始端功率。+(1+j2)=+ DS23=+(+j3)=+ DS24= S23=S3+DS23=+ S24=S4+DS24=+39。=S23+S24+S2=+ 則 S12又+=(+)=+ 從而可得線路始端功率 S12=+ 這個(gè)結(jié)果與第（1）%，所以第（2）和第（3）的結(jié)果可作為最終計(jì)算結(jié)果；若相差較大，則應(yīng)返回第（2）步重新計(jì)算，直道相差較小為止。 如圖所示簡(jiǎn)單系統(tǒng)，額定電壓為110KV 雙回輸電線路，長度為80km，采用LGJ150導(dǎo)線，其單位長度的參數(shù)為：r=，x=,b=180。106S/km。變電所中裝有兩臺(tái)三相110/11kV的變壓器，每臺(tái)的容量為15MVA,其參數(shù)為：DP0=，DPs=128kW,Vs%=,Io%=。母線A的實(shí)際運(yùn)行電壓為117kV，負(fù)荷功率：SLDb=30+j12MVA,SLDc=20+j15MVA。當(dāng)變壓器取主軸時(shí)，求母線c的電壓。解（1）計(jì)算參數(shù)并作出等值電路。輸電線路的等值電阻、電抗和電納分別為 RL=1180。80180。= 21 XL=180。80180。= Bc=2180。80180。180。106S=180。104S由于線路電壓未知，可用線路額定電壓計(jì)算線路產(chǎn)生的充電功率，并將其等分為兩部分，便得112DQB=BcVN=180。180。104180。1102Mvar= 和b，作為節(jié)點(diǎn)負(fù)荷的一部分。兩臺(tái)變壓器并聯(lián)運(yùn)行時(shí)，它們的等值電阻、電抗及勵(lì)磁功率分別為1DPsVN21128180。1102RT==180。W= 2221000SN21000180。151Vs%180。1102RT==180。W= 22100SN2100180。180。15DPo+jDQo=2180。(+j)MVA=+變壓器的勵(lì)磁功率也作為接于節(jié)點(diǎn)b的負(fù)荷，于是節(jié)點(diǎn)b的負(fù)荷Sb=SLDb+jDQB+(DP0+jDQ0)=30+j12++=+ Sc=20+j15MVA 這樣就得到圖所示的等值電路節(jié)（2）計(jì)算母線A輸出的功率。先按電力網(wǎng)絡(luò)的額定電壓計(jì)算電力網(wǎng)絡(luò)中的功率損耗。變壓器繞組中的功率損耗為230。Sc246。202+152247。(RT+jXT)=DST=231。(+)MVA2231。247。110232。VN248。=+Sc39。=Sc+DPT+jDQT=20+j15++=+39。39。=Sc39。+Sb=+++=+線路中的功率損耗為230。S1246。DSL=231。231。V247。247。(RL+jXL)232。N248。+=(+)MVA=+39。39。2 于是可得 S139。=S139。39。+DSL=+++=+由母線A輸出的功率為SA=S139。+jDQB=+=+（3）計(jì)算各節(jié)點(diǎn)電壓。線路中電壓降落的縱分量和橫分量分別為P139。RL+Q139。180。+180。==kV=P139。XLQ139。180。180。==kV=VA117b點(diǎn)電壓為Vb= (VADVL)+(dVL)=22()2+=，橫分量分別為Pc39。RT+Qc39。180。+180。==kV=39。Pc39。180。180。==kV=Vc39。= (VbDVT)+(dVT)=22()2+=1111Vc=Vc180。=180。kV=11011039。如果在上述計(jì)算中都不計(jì)電