【正文】 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 – 結(jié)果：五、思考討論題？簡述各種算法的優(yōu)缺點。2．了解和掌握PSDBPA電力系統(tǒng)分析程序穩(wěn)態(tài)分析方法（即潮流計算）。采用什么算法以及迭代的最大步數(shù)可以由用戶指定。打開TextEdit應用程序時先讀入該文件。分析報告列表：并根據(jù)用戶的要求（通過控制語句控制），輸出各種潮流分析報告。它的任務是根據(jù)給定的運行條件和網(wǎng)路結(jié)構(gòu)確定整個系統(tǒng)的運行狀態(tài)，如各母線上的電壓（幅值及相角）、網(wǎng)絡中的功率分布以及功率損耗等。PSDBPA仿真軟件中潮流計算模型建模的注意事項?。獲益匪淺。有功出力為負值的發(fā)電機（包括電動機）作為負荷處理，不統(tǒng)計在內(nèi)。5）迭代過程。 定義潮流程序卡片格式文件，用戶可更改及調(diào)整該文件。清晰認知電力系統(tǒng)分析的意義。通過編寫這次程序可以更深一步的理解潮流計算的步驟，也明白了在潮流計算中要注意的一些細節(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)潮流計算的計算程序原代碼由自備移動存儲設備導入計算機。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)試方法。,i)disp(angle(EF(i))*180/pi)。收斂后的節(jié)點電壓用極坐標表示為:39。第%d次迭代后的節(jié)點電壓分別為:n39。F1(n+i1,n+w1)=2*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)=G(i,w)*E(i)+B(i,w)*F(i)。%△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)。X=[]。0000 0。%P=PQ節(jié)點的P值 Q=[ 0]。8760+DPT180。248。I%=2231。247。SDST=DPT+jDQT=2231。=ST39。10248。10248。ST1=S1LD+Sc=8+j6+=8+ST2=S2LDSc=8+=8+（4）計算高壓側(cè)電壓。247。已知低壓母線電壓10kV，負荷功率為16+j12MVA，試求變壓器的功率分布和高壓側(cè)電壓。取SI續(xù)進行計算。40W=18+=180。50W=+B213。線路Ⅰ： ZI=(+)180。進一步可求得E點電壓 +=(+)MVA=+10239。7+2200180。（線路參數(shù)LJ185：z=+；LJ70：z=+）解 線路等值阻抗ZL1=10180。=，Vc= 與計及電壓降落橫分量的計算結(jié)果相比，誤差很小。Pc39。XLQ139。39。V247。=+Sc39。變壓器繞組中的功率損耗為230。W= 2221000SN21000180。180。當變壓器取主軸時，求母線c的電壓。V2DV24=(P23R23+Q23X23)DV23==187。=S23+S24+S2=+S12又39。P32++=(R23+jX23)=(1+j2)=+++=(R24+jX24)=(+j3)=+則： S23=S3+DS23=+S24=S4+DS24=+39。V1DV12=(P24R24+Q24X24)DV24==187。母線A的實際運行電壓為117kV，負荷功率：SLDb=30+j12MVA,SLDc=20+j15MVA。80180。1102RT==180。先按電力網(wǎng)絡的額定電壓計算電力網(wǎng)絡中的功率損耗。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)絡參數(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Ω，其電阻和導納忽略不計。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。%v=PV節(jié)點個數(shù)P=[ 0]。 0 0。Y=G+j*B。%△P FX(n+i1)=Q(i)。F1(n+i1,w)=B(i,w)*E(i)G(i,w)*F(i)。F1(n+i1,w)=F1(n+i1,w)+G(i,k)*F(k)+B(i,k)*E(k)。F1(n+i1,w)=2*E(i)。fprintf(39。enddisp(39。%d號節(jié)點電壓的相角度數(shù)為39。運行結(jié)果：運行結(jié)果復制如下：第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é)點電壓用極坐標表示為: 1號節(jié)點電壓的幅值為:2號節(jié)點電壓的幅值為:3號節(jié)點電壓的幅值為:4號節(jié)點電壓的幅值為:平衡節(jié)點的功率 第三篇：電力系統(tǒng)分析潮流計算的計算機算法潮流計算的計算機算法實驗報告姓名：學號：班級：一、實驗目的掌握潮流計算的計算機算法。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)********************************************************************** w