【正文】 計算的程序。f=[0 0 0 0]。endpv(1,2*i)=s+b(i,i)*e(i)g(i,i)*f(i)。pq(2,2*j1)=pq(1,2*j)。for k=1:ls=s+(g(i,k)*f(k)+b(i,k)*e(k))。pq=zeros(2,2*l2)。for i=1:(l1)。s2=0。s2=0。J_ni=inv(J)。程序要求根據(jù)已知的電力網(wǎng)的數(shù)學(xué)模型（節(jié)點導(dǎo)納矩陣）及各節(jié)點參數(shù)，完成該電力系統(tǒng)的潮流計算，要求計算出節(jié)點電壓、功率等參數(shù)。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。endendendend end Jacbi。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))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。k=0。end end DetaS。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。endend end for i=1:nif i~=isbamp。end end B2。%九：求解修正方程，獲取節(jié)點電壓的不平衡量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。)。disp(PPH)。fprintf(39。disp(FX)。end a=a+1。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。F1(i,n+w1)=B(i,w)*E(i)G(i,w)*F(i)。%最后的節(jié)點電壓矩陣 while max(FX)=10^(5)for i=1:m %PQ節(jié)點FX(i)=P(i)。000]。 0 0。%m=PQ節(jié)點個數(shù) v=1。180。248。kVA247。232。246。j1MVA=8+ 21039。k2=(10+)180。k1=(10+)180。VBDE10180。=101247。+180。DSTb242+182(+)MVA=+ =2110Sb=SLDb+DSTb+DSob+jDQBI+jDQBIII=24+j18++++=++92(+)MVA=+ =1102Sc=SLDc+DSTc+DSoc+jDQBIII+jDQBII=12+j9++++=+（3）計算閉式網(wǎng)絡(luò)的功率分布。104180。1102Mvar=線路Ⅱ： Z213。網(wǎng)絡(luò)各組件參數(shù)為：變電所b SNXT==20MVA，DS0=+，RT=，變電所c SNXT=127W=10MVA，DS0=+，RT=，負(fù)荷功率 SLDb=24+j18MVA，SLDc=12+j9MVA試求電力網(wǎng)絡(luò)的功率分布及最大電壓損耗。10+j1200180。()Z229。kV時的初始如圖所示。kV=11011039。+180。+180。2 于是可得 S139。DSL=231。110232。(+j)MVA=+變壓器的勵磁功率也作為接于節(jié)點b的負(fù)荷，于是節(jié)點b的負(fù)荷Sb=SLDb+jDQB+(DP0+jDQ0)=30+j12++=+ Sc=20+j15MVA 這樣就得到圖所示的等值電路節(jié)（2）計算母線A輸出的功率。兩臺變壓器并聯(lián)運行時，它們的等值電阻、電抗及勵磁功率分別為1DPsVN21128180。= Bc=2180。變電所中裝有兩臺三相110/11kV的變壓器，每臺的容量為15MVA,其參數(shù)為：DP0=，DPs=128kW,Vs%=,Io%=。=== V2187。解：（1）先假設(shè)各節(jié)點電壓均為額定電壓，求線路始端功率。39。V2DV23=（3）根據(jù)上述求得的線路各點電壓，重新計算各線路的功率損耗和線路始端功率。解（1）計算參數(shù)并作出等值電路。106S=180。151Vs%180。Sc246。=Sc+DPT+jDQT=20+j15++=+39。247。+DSL=+++=+由母線A輸出的功率為SA=S139。180。180。 某一額定電壓為10kV的兩端供電網(wǎng)，如圖所示。(+)=+ZL2=4180。3+180。 SAC=+++=+DVAC=180。60W=+BI=180。=180。106180。SIII=+=SbSI=+=+由此得到功率初分布，如圖所示。解（1）假定兩臺變壓器變比相同，計算其功率分布。232。不計電壓降落的橫分量時，按變壓器T1計算可得高壓母線電壓為180。按變壓器T2計算可得180。 計及電壓降落的橫分量，按T1和T2計算克分別得。1+ST22=8++8+=16+ 計算所得功率分布，如圖所示。DPs+jsSN247。247。DP0+j0SN247。線路末端充電功率QB2=2blV2=180。3150=172180。%Q=PQ節(jié)點的Q值 PP=[]。0000 0]。%X=△X n=m+v+1。%△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)。elseF1(i,w)=G(i,i)*E(i)B(i,i)*F(i)。endendend end for i=m+1:n1 %PV節(jié)點for w=1:n1if i~=wF1(i,w)=(G(i,w)*E(i)+B(i,w)*F(i))。for k=1:nF1(i,w)=F1(i,w)G(i,k)*E(k)+B(i,k)*F(k)。,a)。)。end PPH=0。二、實驗準(zhǔn)備根據(jù)課程內(nèi)容，熟悉MATLAB軟件的使用方法，自行學(xué)習(xí)MATLAB程序的基礎(chǔ)語法，并根據(jù)所學(xué)知識編寫潮流計算牛頓拉夫遜法(或PQ分解法)的計算程序，用相應(yīng)的算例在MATLAB上進(jìn)行計算、調(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)境下對程序進(jì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分解法，牛頓_拉夫遜法，改進(jìn)的牛頓－拉夫遜算法。該文件安裝時放在與潮流程序相同的目錄中。6）計算結(jié)果輸出：詳細(xì)計算結(jié)果列表：按節(jié)點、與該節(jié)點相聯(lián)接支路順序，并根據(jù)用戶的要求（通過控制語句控制）可按照字母、分區(qū)或區(qū)域排序輸出潮流計算結(jié)果。當(dāng)最大出力值小于實際出力時，統(tǒng)計時最大出力值用實際出力值代替。電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種基本