【正文】 此后，潮流計算曾采用了各種不同的方法，這些方法的發(fā)展主要是圍繞著對潮流計算的一些基本要求進行的。牛頓拉夫遜法迭代次數(shù)為 5次。表1 IEEE 9節(jié)點算例節(jié)點參數(shù)表2 IEEE 9節(jié)點算例線路參數(shù)表3 IEEE 9節(jié)點算例變壓器參數(shù)注：表1表3中功率基準值為100MVA；電阻、電感值為標幺值。3）錯誤信息。除了潮流計算功能外，該軟件還具有自動電壓控制、聯(lián)絡線功率控制、系統(tǒng)事故分析（N1開斷模擬）、網(wǎng)絡等值、靈敏度分析、節(jié)點P－V、Q－V和P－Q曲線、確定系統(tǒng)極限輸送水平、負荷靜特性模型、靈活多樣的分析報告、詳細的檢錯功能等功能。無論是電力系統(tǒng)規(guī)劃、設計還是運行，對其安全穩(wěn)定進行分析都是極其重要的。？為什么？答:不一樣,因為每次迭代的電壓、有功、無功都是與前一次不同的新值，所以每次迭代過程中，雅可比矩陣都是變化的。e=[1 1 ]。for k=1:ls=s(g(i,k)*f(k)+b(i,k)*e(k))。pq(2,2*j)=pq(1,2*j1)。s=0。endend endfunction pq=PQ_out(g,b,e,f,i)%計算pq節(jié)點的雅克比矩陣 l=length(e)。J=zeros(2*l2,2*l2)。for i=1:(l1)s1=0。for i=1:l s1=0。while max(abs(s))10^5 J=J_out(g,b,kind,e,f)。二、實驗內(nèi)容：編制調(diào)試電力系統(tǒng)潮流計算的計算機程序。end end for i=1:nif B2(i,6)==3j=j+1。endif k==(n1)k=0。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))。elseJacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。h=0。DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2real(B2(i,3))^2imag(B2(i,3))^2。B2(i,6)==2h=h+1。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)))。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)。%四：創(chuàng)建PVU 用于存儲PV節(jié)點的初始電壓PVU=zeros(nh1,1)。B2(i,6)==2h=h+1。Y(p,p)=Y(p,p)+1/B1(i,3)。Y(p,q)=Y(p,q)1/B1(i,3)。)。pr=input(39。請輸入節(jié)點數(shù):n=39。3150kWh=2300520kWh線路全年的電能損耗DWL=DPL180。100180。86+j180。2+j180。231。試求線路和變壓器全年的電能損耗。S S39。 VA=230。 VA=230。232。因為阻抗已歸算到低壓側，宜用低壓側的電壓求環(huán)路電勢。由于線路Ⅰ和Ⅲ的功率均流向節(jié)點b為功率分點，且有功功率分點和無公功功率分點都在b點，因此這點的電壓最低。104S =180。50S=180。60S=180。=VC=VADVAC== DVCE=180。3)+Sc17 =++580+j129=+=1(180。(+)=+ ZL4=2180。A=176。==kV=Vc39。==kV=VA117b點電壓為Vb= (VADVL)+(dVL)=22()2+=，橫分量分別為Pc39。線路中電壓降落的縱分量和橫分量分別為P139。N248。=Sc39。(RT+jXT)=DST=231。W= 22100SN2100180。180。80180。=S23+S24+S2=+ 則 S12又+=(+)=+ 從而可得線路始端功率 S12=+ 這個結果與第（1）%，所以第（2）和第（3）的結果可作為最終計算結果；若相差較大，則應返回第（2）步重新計算，直道相差較小為止。故： S12=S12+DS12=+（2）再用已知的線路始端電壓V1=12，求出線路各點電壓。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。%V=PV節(jié)點的U值E=[1 1 1 ]39。 15 0。%F(x)矩陣F1=zeros(n1,n1)。%△P FX(n+i1)=V(im)^2E(i)^2F(i)^2。F1(n+i1,w)=B(i,i)*E(i)G(i,i)*F(i)。F1(n+i1,w)=0。endendend end X=inv(F1)*(FX)。fprintf(39。for i=1:n1 fprintf(39。end fprintf(39。=S23+S24+S2=+S12又39。V2DV24=(P23R23+Q23X23)DV23==187。當變壓器取主軸時，求母線c的電壓。180。W= 2221000SN21000180。變壓器繞組中的功率損耗為230。=+Sc39。V247。39。XLQ139。Pc39。=，Vc= 與計及電壓降落橫分量的計算結果相比，誤差很小。（線路參數(shù)LJ185：z=+；LJ70：z=+）解 線路等值阻抗ZL1=10180。7+2200180。進一步可求得E點電壓 +=(+)MVA=+10239。線路Ⅰ： ZI=(+)180。50W=+B213。40W=18+=180。取SI續(xù)進行計算。已知低壓母線電壓10kV，負荷功率為16+j12MVA，試求變壓器的功率分布和高壓側電壓。247。ST1=S1LD+Sc=8+j6+=8+ST2=S2LDSc=8+=8+（4）計算高壓側電壓。10248。10248。=ST39。SDST=DPT+jDQT=2231。247。I%=2231。248。8760+DPT180。三、實驗要求每人一組，在實驗課時內(nèi)，用MATLAB調(diào)試和修改運行程序，用算例計算輸出潮流結果。isb=input(39。)。Times=1。q=B1(i,2)。%二：創(chuàng)建OrgS，用于存儲初始功率參數(shù)h=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)))。%五：創(chuàng)建DetaS,用于存儲有功功率、無功功率和電壓幅值的不平衡量h=0。t=t+1。end end I。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)+2*real(I(h,1))。for i=1:n%對PV節(jié)點的處理if B2(i,6)==3h=h+1。Jacbi(2*h,2*k1)=0。for i=1:n%對PQ節(jié)點處理if B2(i,6)==2j=j+1。for i=1:nif i~=isbamp。endend end OrgS。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。五、實驗流程六、實驗結果參數(shù)輸入：運行結果：七、實驗體會通過這次實驗，讓我第一次接觸到了MATLAB,并深切體會到了它的強大之處；潮流計算的計算機算法的實現(xiàn)不僅鞏固了我的學過的知識，還讓我學到一些MATLAB的編程，雖然在實驗的過程中出現(xiàn)了很多的錯誤，但在老師的細心指導下，問題都解決啦；計算機為我們省去了大量的人工計算，希望在以后的學習中能接觸到更多的軟件，學習到更多的知識。應用計算例題驗證程序的計算效果。f(i)=f(i)dv(2*i)。q(i)=f(i)*s1e(i)*s2。if kind(1,i)==1s(2*i)=kind(3,i)f(i)*s1+e(i)*s2。endelses=PV_out(g,b,e,f,i)。s=0。endpq(2,2*i)=s+g(i,i)*e(i)+b(i,i)*f(i)。for k=1:ls=s(g(i,k)*e(k)b(i,k)*f(k))。pv(1,2*j)=b(i,j)*e(i)g(i,j)*f(i)。 容許誤差ε=105節(jié)點導納矩陣：各節(jié)點電壓：節(jié)點efvζ 各節(jié)點功率：節(jié)點PQ 2 – 結果：五、思考討論題？簡述各種算法的優(yōu)缺點。2．了解和掌握PSDBPA電力系統(tǒng)分析程序穩(wěn)態(tài)分析方法（即潮流計算）。采用什么算法以及迭代的最大步數(shù)可以由用戶指定。打開TextEdit應用程序時先讀入該文件