【正文】 。分析報告列表：并根據(jù)用戶的要求（通過控制語句控制），輸出各種潮流分析報告。它的任務(wù)是根據(jù)給定的運行條件和網(wǎng)路結(jié)構(gòu)確定整個系統(tǒng)的運行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡(luò)中的功率分布以及功率損耗等。PSDBPA仿真軟件中潮流計算模型建模的注意事項?。獲益匪淺。有功出力為負(fù)值的發(fā)電機(jī)（包括電動機(jī)）作為負(fù)荷處理，不統(tǒng)計在內(nèi)。5）迭代過程。 定義潮流程序卡片格式文件，用戶可更改及調(diào)整該文件。清晰認(rèn)知電力系統(tǒng)分析的意義。通過編寫這次程序可以更深一步的理解潮流計算的步驟，也明白了在潮流計算中要注意的一些細(xì)節(jié)。[e,f,p,q]=flow_out(g,b,kind,e,f)。pv(2,2*i)=2*f(i)。pv=zeros(2,2*l2)。s=0。for k=1:ls=s(g(i,k)*e(k)b(i,k)*f(k))。for j=1:(2*l2)J(2*i1,j)=s(1,j)。s2=s2+g(i,j)*f(j)+b(i,j)*e(j)。s2=s2+g(i,j)*f(j)+b(i,j)*e(j)。l=length(dv)/2。將事先編制好的電力系統(tǒng)潮流計算的計算程序原代碼由自備移動存儲設(shè)備導(dǎo)入計算機(jī)。Times=Times+1。DetaU=inv(Jacbi)*DetaS。Jacbi(2*h,2*k)=2*real(B2(i,3))。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)。for j=1:nif j~=isbk=k+1。h=0。end end 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)))。%！!初始功率參數(shù)在迭代過程中是不累加的，所以在這里必須將其初始化為零矩陣h=0。DetaU。else%非對角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。endif k==(n1)%將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行k=0。if i==j%對角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))+imag(I(h,1))。for i=1:nif i~=isbh=h+1。for i=1:n%六：對PV節(jié)點的處理，注意這時不可再將h初始化為0if i~=isbamp。PVU(t,1)=B2(i,3)。endend end %三：對PV節(jié)點的處理，注意這時不可再將h初始化為0 for i=1:nif i~=isbamp。%將OrgS、DetaS初始化OrgS=zeros(2*n2,1)。Y(q,q)=Y(q,q)+1/B1(i,3)+*B1(i,4)。)。B1=input(39。請輸入支路數(shù):n1=39。熟悉MATLAB，并掌握MATLAB程序的基本調(diào)試方法。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é)點電壓均為額定電壓，求線路始端功率。end PPH=0。)。,a)。for k=1:nF1(i,w)=F1(i,w)G(i,k)*E(k)+B(i,k)*F(k)。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))。elseF1(i,w)=G(i,i)*E(i)B(i,i)*F(i)。%△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)。%X=△X n=m+v+1。0000 0]。%Q=PQ節(jié)點的Q值 PP=[]。%P=PQ節(jié)點的P值 Q=[ 0]。0000 0。X=[]。%△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)。F1(n+i1,n+w1)=G(i,w)*E(i)+B(i,w)*F(i)。F1(n+i1,n+w1)=F1(n+i1,n+w1)G(i,k)*E(k)+B(i,k)*F(k)。F1(n+i1,n+w1)=2*F(i)。第%d次迭代后的節(jié)點電壓分別為:n39。收斂后的節(jié)點電壓用極坐標(biāo)表示為:39。,i)disp(angle(EF(i))*180/pi)。已知各節(jié)點的負(fù)荷功率及參數(shù)：S2=(+)MVA，S3=(+)MVA，S4=(+)MVAZ12=(+)W，Z23=(+)W，Z24=(+)W試求電壓和功率分布。+180。106S/km。80180。1102Mvar= 和b，作為節(jié)點負(fù)荷的一部分。15DPo+jDQo=2180。247。S1246。39。180。180。=180。B=176。10=(+)kVA=580+j129kVA =*17180。14+180。=VE=VCDVCE== 圖所示110kV閉式電網(wǎng)，A點為某發(fā)電廠的高壓母線，其運行電壓為117kV。104180。=180。1102Mvar=1(+)W=+ 2 2DQBIII 變電所b：ZTb= DS0b=2(+)MVA=+變電所b：ZTc=1(+j127)W=+ 2 DS0c=2(+)MVA=+ 等值電路如圖所示（2）計算節(jié)點b和c的運算負(fù)荷。SA1=SI++(+)MVA=+=++1102DVI=PA1RI+180。247。故循環(huán)功率為 Sc187。10+247。10+247。j1MVA=8+10282+=8++180。106S/km變壓器（每臺）:DP0=86kW,DPs=200kW,I0%=,Vs%=解 最大負(fù)荷時變壓器的繞組功率損耗為V%230。248。231。kVA=172+j1701100100232。180。3150kWh=5946885kWh 輸電系統(tǒng)全年的總電能損耗DWT+DWL=2300520+5946885kWh=8247405kWh第三篇：電力系統(tǒng)分析潮流計算的計算機(jī)算法潮流計算的計算機(jī)算法實驗報告姓名：學(xué)號：班級：一、實驗?zāi)康恼莆粘绷饔嬎愕挠嬎銠C(jī)算法。n1=input(39。)。節(jié)點號和對地參數(shù):X=39。Y(p,p)=Y(p,p)+1/B1(i,3)+*B1(i,4)。end end Y。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 B2(i,6)==3t=t+1。end end t=0。h=0。for j=1:nif j~=isbk=k+1。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)。Jacbi(2*h,2*k)=2*real(B2(i,3))。DetaU=inv(Jacbi)*DetaS。%十：開始循環(huán)********************************************************************** while abs(max(DetaU))pr OrgS=zeros(2*n2,1)。B2(i,6)==3h=h+1。DetaS(2*h,1)=imag(B2(i,2))OrgS(2*h,1)。%創(chuàng)建I i=zeros(n1,1)。for i=1:nif B2(i,6)==2h=h+1。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)。Jacbi(2*h,2*k1)=2*imag(B2(i,3))。DetaU=zeros(2*n2,1)。end end B2。在各種潮流計算的算法中選擇一種，按照計算方法編制程序。dv=J_ni*s。for j=1:ls1=s1+g(i,j)*e(j)b(i,j)*f(j)。for j=1:ls1=s1+g(i,j)*e(j)b(i,j)*f(j)。if kind(1,i)==1s=PQ_out(g,b,e,f,i)。for j=1:(l1)if j==i s=0。endpq(2,2*i1)=s+b(i,i)*e(i)g(i,i)*f(i)。end endfunction pv=PV_out(g,b,e,f,i)%計算pv節(jié)點的雅克比矩陣 l=length(e)。pv(2,2*i1)=2*e(i)。kind=[