【正文】 DetaU=zeros(2*n2,1)。Jacbi(2*h,2*k1)=2*imag(B2(i,3))。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)。for i=1:n%對PQ節(jié)點的處理if B2(i,6)==2h=h+1。%七：創(chuàng)建I，用于存儲節(jié)點電流參數(shù)i=zeros(n1,1)。DetaS(2*h,1)=imag(B2(i,2))OrgS(2*h,1)。t=0。for j=1:nOrgS(2*h1,1)=OrgS(2*h1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*Imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3))。Y(q,p)=Y(p,q)。X=input(39。請輸入誤差精度:pr=39。)。3150=180。1102Mvar=2等值電路中流過線路等值阻抗的功率為S1=S+DST+DS0+jQB2=40+j30++++=++180。31500kVA 247。31500247。100232。線路和變壓器的參數(shù)如下：線路（每回）：r=, x=, b=180。T182+=8++180。231。231。k1248。若取其假定正方向為順時針方向，則可得k2 DE187。為了計算線路Ⅰ的電壓損耗，要用A點的電壓和功率SA1。104180。104S2DQB213。104S2DQBI=180。+180。10+2200180。(+)=+ 求C點和D點的運算負(fù)荷，為 DSCE +=(+)=+ 210=2600+j1600+300+j160++=+SD=600+j200+1600+j1000=2200+j1200kVA循環(huán)功率Sc(V=*A*VBVN()180。kV、Vamp。= (VbDVT)+(dVT)=22()2+=1111Vc=Vc180。RT+Qc39。RL+Q139。+=(+)MVA=+39。+Sb=+++=+線路中的功率損耗為230。(+)MVA2231。180。104180。= 21 XL=180。 如圖所示簡單系統(tǒng)，額定電壓為110KV 雙回輸電線路，長度為80km，采用LGJ150導(dǎo)線，其單位長度的參數(shù)為：r=，x=,b=180。(P12R12+Q12X12)180。運行結(jié)果：運行結(jié)果復(fù)制如下：第1次迭代后的節(jié)點電壓分別為:第1次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:00第2次迭代后的節(jié)點電壓分別為:第2次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:第3次迭代后的節(jié)點電壓分別為:第3次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:第4次迭代后的節(jié)點電壓分別為:第4次迭代后功率偏差△P △Q電壓偏差△V的平方分別為: *收斂后的節(jié)點電壓用極坐標(biāo)表示為: 1號節(jié)點電壓的幅值為:2號節(jié)點電壓的幅值為:3號節(jié)點電壓的幅值為:4號節(jié)點電壓的幅值為:平衡節(jié)點的功率 第二篇：電力系統(tǒng)分析潮流計算例題3．1 電網(wǎng)結(jié)構(gòu)如圖3—11所示，其額定電壓為10KV。%d號節(jié)點電壓的相角度數(shù)為39。enddisp(39。fprintf(39。F1(n+i1,w)=2*E(i)。F1(n+i1,w)=F1(n+i1,w)+G(i,k)*F(k)+B(i,k)*E(k)。F1(n+i1,w)=B(i,w)*E(i)G(i,w)*F(i)。%△P FX(n+i1)=Q(i)。Y=G+j*B。 0 0。%v=PV節(jié)點個數(shù)P=[ 0]。%PP=PV節(jié)點的P值 V=[]。B=[ 00。%總的節(jié)點數(shù)FX=ones(2*n2,1)。%△Qend end for i=m+1:n1 %PV節(jié)點 FX(i)=PP(im)。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。F1(i,n+w1)=B(i,w)*E(i)G(i,w)*F(i)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。disp(E+j*F)。EF=E+j*F。for i=1:nPPH=PPH+EF(n)*conj(Y(n,i))*conj(EF(i))。P32++=(R23+jX23)=(1+j2)=+++=(R24+jX24)=(+j3)=+則： S23=S3+DS23=+S24=S4+DS24=+39。V1DV12=(P24R24+Q24X24)DV24==187。母線A的實際運行電壓為117kV，負(fù)荷功率：SLDb=30+j12MVA,SLDc=20+j15MVA。80180。1102RT==180。先按電力網(wǎng)絡(luò)的額定電壓計算電力網(wǎng)絡(luò)中的功率損耗。VN248。231。=S139。==kV=P139。==kV=39。如果在上述計算中都不計電壓降落的橫分量，所得結(jié)果為Vb=，Vc39。試求V功率分布，且找到電壓最低點。)1(180。14)Sc 17=+=+=SAC+SBD=+++=+ SC+SD=++2200+j1200=+ SCD=SBDSD=+=+C點為功率分點，可推算出E點為電壓最低點。解（1）計算網(wǎng)絡(luò)參數(shù)及制定等值電路。=(+)180。1102Mvar=線路Ⅱ： ZIII BIII=(+)180。***SbZII+ZIII+ScZII(+)()+(+)()MVSI==***+ZII+ZIII()=+=*ScZI*+ZIII+SbZI***ZI*+ZII+ZIII()=(+)()+(+)()=++SII=+++=+Sb+Sc=+++=+可見，計算結(jié)果誤差很小，無需重算。== VA117變電所b高壓母線的實際電壓為 Vb=VADVI== 變比分別為k1=110/11和k2=，如圖所示，兩臺變壓器歸算到低壓側(cè)的電抗均為1Ω，其電阻和導(dǎo)納忽略不計。kV= 231。=MVA= **ZT1+ZT2j1j1（3）計算兩臺變壓器的實際功率分布。10kV=232。=232。T2輸入高壓母線的總功率為S39。230。2SN2246。kVA=252+j41661002180。232。100MW= DPL=2RL=2V1102已知cosj=，Tmax=4500h，從表中查得t=3150h，假定變壓器全年投入運行，則變壓器全年的電能損耗 DWT=2DP0180。二、實驗準(zhǔn)備根據(jù)課程內(nèi)容，熟悉MATLAB軟件的使用方法，自行學(xué)習(xí)MATLAB程序的基礎(chǔ)語法，并根據(jù)所學(xué)知識編寫潮流計算牛頓拉夫遜法(或PQ分解法)的計算程序，用相應(yīng)的算例在MATLAB上進行計算、調(diào)試和驗證。)。請輸入支路參數(shù):B1=39。Y=zeros(n)。else%含有變壓器的支路p=B1(i,1)。DetaS=zeros(2*n2,1)。B2(i,6)==3h=h+1。end end PVU。B2(i,6)==3h=h+1。I(h,1)=(OrgS(2*h1,1)OrgS(2*h,1)*sqrt(1))/conj(B2(i,3))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))+real(I(h,1))。endendendend end k=0。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。%修正節(jié)點電壓 j=0。j=0。OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))imag(Y(i,j))*imag(B2(j,3)))real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。for i=1:nif i~=isbamp。for i=1:nif i~=isbh=h+1。if i==jJacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1))。endif k==(n1)k=0。elseJacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。DetaU。%迭代次數(shù)加1 end Times。在相應(yīng)的編程環(huán)境下對程序進行組織調(diào)試。for i=1:le(i)=e(i)dv(2*i1)。endp(i)=e(i)*s1+f(i)*s2。ends(2*i1)=kind(2,i)e(i)*s1f(i)*s2。J(2*i,j)=s(2,j)。endpq(1,2*i1)=sg(i,i)*e(i)b(i,i)*f(i)。for k=1:ls=s(g(i,k)*e(k)b(i,k)*f(k))。for j=1:(l1)if j==i s=0。elsepv(1,2*j1)=(g(i,j)*e(i)+b(i,j)*f(i))。e f四、例題及運行結(jié)果在上圖所示的簡單電力系統(tǒng)中，系統(tǒng)中節(jié)點2為PQ節(jié)點，節(jié)點3為PV節(jié)點，節(jié)點4為平衡節(jié)點，已給定 P1s+jQ1s= P2s+jQ2s= P3s= V3s= V4s=∠0176。第五篇：電力系統(tǒng)分析計算機算法PSDBPA潮流計算實驗報告電力系統(tǒng)分析的計算機算法實驗報告學(xué)生姓名 課 程 電力系統(tǒng)分析的計算機算法 學(xué) 號專 業(yè) 電氣工程及其自動化 指導(dǎo)教師 邱曉燕二Ο一四 年 六 月 二日實驗一潮流計算一、實驗?zāi)康?．了解并掌握電力系統(tǒng)計算機算法的相關(guān)原理。三、原理和說明PSDBPA電力系統(tǒng)分析程序穩(wěn)態(tài)分析主要是潮流計算，軟件中潮流程序的計算方法有P_Q分解法，牛頓_拉夫遜法，改進的牛頓－拉夫遜算法。該文件安裝時放在與潮流程序相同的目錄中。6）計算結(jié)果輸出：詳細(xì)計算結(jié)果列表：按節(jié)點、與該節(jié)點相聯(lián)接支路順序，并根據(jù)用戶的要求（通過控制語句控制）可按照字母、分區(qū)或區(qū)域排序輸出潮流計算結(jié)果。當(dāng)最大出力值小于實際出力時，統(tǒng)計時最大出力值用實際出力值代替。電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種基本