【正文】 S23+S24+S2=+ 則 S12又+=(+)=+ 從而可得線路始端功率 S12=+ 這個結果與第（1）%，所以第（2）和第（3）的結果可作為最終計算結果；若相差較大，則應返回第（2）步重新計算，直道相差較小為止。母線A的實際運行電壓為117kV，負荷功率：SLDb=30+j12MVA,SLDc=20+j15MVA。80180。80180。180。1102RT==180。W= 22100SN2100180。先按電力網絡的額定電壓計算電力網絡中的功率損耗。(RT+jXT)=DST=231。VN248。=Sc39。231。N248。=S139。線路中電壓降落的縱分量和橫分量分別為P139。==kV=P139。==kV=VA117b點電壓為Vb= (VADVL)+(dVL)=22()2+=，橫分量分別為Pc39。==kV=39。==kV=Vc39。如果在上述計算中都不計電壓降落的橫分量，所得結果為Vb=，Vc39。A=176。試求V功率分布，且找到電壓最低點。(+)=+ ZL4=2180。)1(180。3)+Sc17 =++580+j129=+=1(180。14)Sc 17=+=+=SAC+SBD=+++=+ SC+SD=++2200+j1200=+ SCD=SBDSD=+=+C點為功率分點，可推算出E點為電壓最低點。=VC=VADVAC== DVCE=180。解（1）計算網絡參數及制定等值電路。60S=180。=(+)180。50S=180。1102Mvar=線路Ⅱ： ZIII BIII=(+)180。104S =180。***SbZII+ZIII+ScZII(+)()+(+)()MVSI==***+ZII+ZIII()=+=*ScZI*+ZIII+SbZI***ZI*+ZII+ZIII()=(+)()+(+)()=++SII=+++=+Sb+Sc=+++=+可見，計算結果誤差很小，無需重算。由于線路Ⅰ和Ⅲ的功率均流向節(jié)點b為功率分點，且有功功率分點和無公功功率分點都在b點，因此這點的電壓最低。== VA117變電所b高壓母線的實際電壓為 Vb=VADVI== 變比分別為k1=110/11和k2=，如圖所示，兩臺變壓器歸算到低壓側的電抗均為1Ω，其電阻和導納忽略不計。因為阻抗已歸算到低壓側，宜用低壓側的電壓求環(huán)路電勢。kV= 231。232。=MVA= **ZT1+ZT2j1j1（3）計算兩臺變壓器的實際功率分布。 VA=230。10kV=232。 VA=230。=232。S S39。T2輸入高壓母線的總功率為S39。試求線路和變壓器全年的電能損耗。230。231。2SN2246。2+j180。kVA=252+j41661002180。86+j180。232。100180。100MW= DPL=2RL=2V1102已知cosj=，Tmax=4500h，從表中查得t=3150h，假定變壓器全年投入運行，則變壓器全年的電能損耗 DWT=2DP0180。3150kWh=2300520kWh線路全年的電能損耗DWL=DPL180。二、實驗準備根據課程內容，熟悉MATLAB軟件的使用方法，自行學習MATLAB程序的基礎語法，并根據所學知識編寫潮流計算牛頓拉夫遜法(或PQ分解法)的計算程序，用相應的算例在MATLAB上進行計算、調試和驗證。請輸入節(jié)點數:n=39。)。pr=input(39。請輸入支路參數:B1=39。)。Y=zeros(n)。Y(p,q)=Y(p,q)1/B1(i,3)。else%含有變壓器的支路p=B1(i,1)。Y(p,p)=Y(p,p)+1/B1(i,3)。DetaS=zeros(2*n2,1)。B2(i,6)==2h=h+1。B2(i,6)==3h=h+1。%四：創(chuàng)建PVU 用于存儲PV節(jié)點的初始電壓PVU=zeros(nh1,1)。end end PVU。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。B2(i,6)==3h=h+1。end end DetaS。I(h,1)=(OrgS(2*h1,1)OrgS(2*h,1)*sqrt(1))/conj(B2(i,3))。k=0。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))。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))+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))。endendendend end Jacbi。%修正節(jié)點電壓 j=0。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。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)))。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,6)==2h=h+1。for i=1:nif i~=isbamp。DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2real(B2(i,3))^2imag(B2(i,3))^2。for i=1:nif i~=isbh=h+1。h=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))。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))。endif k==(n1)k=0。DetaU。end end for i=1:nif B2(i,6)==3j=j+1。%迭代次數加1 end Times。二、實驗內容：編制調試電力系統潮流計算的計算機程序。在相應的編程環(huán)境下對程序進行組織調試。while max(abs(s))10^5 J=J_out(g,b,kind,e,f)。for i=1:le(i)=e(i)dv(2*i1)。for i=1:l s1=0。endp(i)=e(i)*s1+f(i)*s2。for i=1:(l1)s1=0。ends(2*i1)=kind(2,i)e(i)*s1f(i)*s2。J=zeros(2*l2,2*l2)。J(2*i,j)=s(2,j)。endend endfunction pq=PQ_out(g,b,e,f,i)%計算pq節(jié)點的雅克比矩陣 l=length(e)。endpq(1,2*i1)=sg(i,i)*e(i)b(i,i)*f(i)。s=0。for k=1:ls=s(g(i,k)*e(k)b(i,k)*f(k))。pq(2,2*j)=pq(1,2*j1)。for j=1:(l1)if j==i s=0。for k=1:ls=s(g(i,k)*f(k)+b(i,k)*e(k))。elsepv(1,2*j1)=(g(i,j)*e(i)+b(i,j)*f(i))。e=[1 1 ]。e f四、例題及運行結果在上圖所示的簡單電力系統中，系統中節(jié)點2為PQ節(jié)點，節(jié)點3為PV節(jié)點，節(jié)點4為平衡節(jié)點，已給定 P1s+jQ1s= P2s+jQ2s= P3s= V3s= V4s=∠0176。？為什么？答:不一樣,因為每次迭代的電壓、有功、無功都是與前一次不同的新值，所以每次迭代過程中，雅可比矩陣都是變化的。第五篇：電力系統分析計算機算法PSDBPA潮流計算實驗報告電力系統分析的計算機算法實驗報告學生姓名 課 程 電力系統分析的計算機算法 學 號專 業(yè) 電氣工程及其自動化 指導教師 邱曉燕二Ο一四 年 六 月 二日實驗一潮流計算一、實驗目的1．了解并掌握電力系統計算機算法的相關原理。無論是電力系統規(guī)劃、設計還是運行，對其安全穩(wěn)定進行分析都是極其重要的。三、原理和說明PSDBPA電力系統分析程序穩(wěn)態(tài)分析主要是潮流計算，軟件中潮流程序的計算方法有P_Q分解法，牛頓_拉夫遜法，改進的牛頓－拉夫遜算法。除了潮流計算功能外，該軟件還具有自動電壓控制、聯絡線功率控制、系統事故分析（N1開斷模擬）、網絡等值、靈敏度分析、節(jié)點P－V、Q－V和P－Q曲線、確定系統極限輸送水平、負荷靜特性模型、靈活多樣的分析報告、詳細的檢錯功能等功能。該文件安裝時放在與潮流程序相同的目錄中。3）錯誤信息。6）計算結果輸出：詳細計算結果列表：按節(jié)點、與該節(jié)點相聯接支路順序，并根據用戶的要求（通過控制語句控制）可按照字母、分區(qū)或區(qū)域排序輸出潮流計算結果。表1 IEEE 9節(jié)點算例節(jié)點參數表2 IEEE 9節(jié)點算例線路參數表3 IEEE 9節(jié)點算例變壓器參數注：表1表3中功率基準值為100MVA；電阻、電感值為標幺值。當最大出力值小于實際出力時，統計時最大出力值用實際出力值代替。牛頓拉夫遜法迭代次數為 5次。電力系統潮流計算是研究電力系統穩(wěn)態(tài)運行情況的一種基本計算。此后，潮流計算曾采用了各種不同的方法，這些方法的發(fā)展主要是圍繞著對潮流計算的一些基本要求進