【正文】 電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種基本計算。6）計算結(jié)果輸出：詳細計算結(jié)果列表：按節(jié)點、與該節(jié)點相聯(lián)接支路順序，并根據(jù)用戶的要求（通過控制語句控制）可按照字母、分區(qū)或區(qū)域排序輸出潮流計算結(jié)果。三、原理和說明PSDBPA電力系統(tǒng)分析程序穩(wěn)態(tài)分析主要是潮流計算，軟件中潮流程序的計算方法有P_Q分解法，牛頓_拉夫遜法，改進的牛頓－拉夫遜算法。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。for j=1:(l1)if j==i s=0。endpq(1,2*i1)=sg(i,i)*e(i)b(i,i)*f(i)。ends(2*i1)=kind(2,i)e(i)*s1f(i)*s2。for i=1:le(i)=e(i)dv(2*i1)。%迭代次數(shù)加1 end Times。elseJacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。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))。for i=1:nif i~=isbamp。j=0。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(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))。B2(i,6)==3h=h+1。B2(i,6)==3h=h+1。else%含有變壓器的支路p=B1(i,1)。請輸入支路參數(shù):B1=39。二、實驗準備根據(jù)課程內(nèi)容，熟悉MATLAB軟件的使用方法，自行學習MATLAB程序的基礎語法，并根據(jù)所學知識編寫潮流計算牛頓拉夫遜法(或PQ分解法)的計算程序，用相應的算例在MATLAB上進行計算、調(diào)試和驗證。)。for k=1:nF1(i,w)=F1(i,w)G(i,k)*E(k)+B(i,k)*F(k)。elseF1(i,w)=G(i,i)*E(i)B(i,i)*F(i)。%X=△X n=m+v+1。%Q=PQ節(jié)點的Q值 PP=[]。線路末端充電功率QB2=2blV2=180。247。1+ST22=8++8+=16+ 計算所得功率分布，如圖所示。按變壓器T2計算可得180。232。SIII=+=SbSI=+=+由此得到功率初分布，如圖所示。=180。 SAC=+++=+DVAC=180。(+)=+ZL2=4180。180。+DSL=+++=+由母線A輸出的功率為SA=S139。=Sc+DPT+jDQT=20+j15++=+39。151Vs%180。解（1）計算參數(shù)并作出等值電路。39。=== V2187。= Bc=2180。(+j)MVA=+變壓器的勵磁功率也作為接于節(jié)點b的負荷，于是節(jié)點b的負荷Sb=SLDb+jDQB+(DP0+jDQ0)=30+j12++=+ Sc=20+j15MVA 這樣就得到圖所示的等值電路節(jié)（2）計算母線A輸出的功率。DSL=231。+180。kV=11011039。()Z229。網(wǎng)絡各組件參數(shù)為：變電所b SNXT==20MVA，DS0=+，RT=，變電所c SNXT=127W=10MVA，DS0=+，RT=，負荷功率 SLDb=24+j18MVA，SLDc=12+j9MVA試求電力網(wǎng)絡的功率分布及最大電壓損耗。104180。+180。VBDE10180。k2=(10+)180。246。kVA247。180。 0 0。%最后的節(jié)點電壓矩陣 while max(FX)=10^(5)for i=1:m %PQ節(jié)點FX(i)=P(i)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。end a=a+1。fprintf(39。)。X=input(39。Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3))。t=0。%七：創(chuàng)建I，用于存儲節(jié)點電流參數(shù)i=zeros(n1,1)。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)。%九：求解修正方程，獲取節(jié)點電壓的不平衡量DetaU=zeros(2*n2,1)。endend end for i=1:nif i~=isbamp。end end DetaS。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。endendendend end Jacbi。程序要求根據(jù)已知的電力網(wǎng)的數(shù)學模型（節(jié)點導納矩陣）及各節(jié)點參數(shù)，完成該電力系統(tǒng)的潮流計算，要求計算出節(jié)點電壓、功率等參數(shù)。s2=0。for i=1:(l1)。for k=1:ls=s+(g(i,k)*f(k)+b(i,k)*e(k))。endpv(1,2*i)=s+b(i,i)*e(i)g(i,i)*f(i)。六、實驗心得這次實驗是通過matlab編寫出一個潮流計算的程序。、輸出相關(guān)文件 *.dat潮流計算數(shù)據(jù)文件*.bse潮流計算二進制結(jié)果文件（可用于潮流計算的輸入或穩(wěn)定計算）*.pfo潮流計算結(jié)果文件*.map 供單線圖格式潮流圖及地理接線圖格式潮流圖程序使用的二進制結(jié)果文件*.pff，*.pfd 中間文件（正常計算結(jié)束后將自動刪除。對應于上述系統(tǒng)及數(shù)據(jù)的潮流計算數(shù)據(jù)（）見例1。從數(shù)學上說，潮流計算是求解一組由潮流方程描述的非線性代數(shù)方程組。* 潮流計算迭代過程和平衡節(jié)點相關(guān)信息數(shù)據(jù)計算結(jié)果收斂。2）輸入、輸出文件及輸出的內(nèi)容列表。二、實驗背景隨著科學技術(shù)的飛速發(fā)展，電力系統(tǒng)也在不斷地發(fā)展，電網(wǎng)通過互聯(lián)變得越來越復雜，同時也使系統(tǒng)穩(wěn)定問題越來越突出。b=[ 0 0 0 ]。pq(1,2*j)=b(i,j)*e(i)g(i,j)*f(i)。J(2*i,j)=s(2,j)。s=zeros(2*l2,1)。k=0。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。for j=1:nif j~=isbk=k+1。%創(chuàng)建Jacbi Jacbi=zeros(2*n2)。for i=1:nif i~=isbamp。end end for i=1:n%對PV節(jié)點的處理if B2(i,6)==3j=j+1。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))。h=0。B2(i,6)==2h=h+1。for i=1:n%對PQ節(jié)點的處理if i~=isbamp。q=B1(i,2)。)。平衡節(jié)點的功率39。第%d次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:n39。F1(n+i1,n+w1)=0。%△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)。 0。3150=180。31500kVA 247。100232。T182+=8++180。231。若取其假定正方向為順時針方向，則可得k2 DE187。104180。104S2DQBI=180。10+2200180。kV、Vamp。RT+Qc39。+=(+)MVA=+39。(+)MVA2231。104180。 如圖所示簡單系統(tǒng)，額定電壓為110KV 雙回輸電線路，長度為80km，采用LGJ150導線，其單位長度的參數(shù)為：r=，x=,b=180。第一篇：電力系統(tǒng)分析潮流計算例題3．1 電網(wǎng)結(jié)構(gòu)如圖3—11所示，其額定電壓為10KV。=S23+S24+S2=+ 則 S12又+=(+)=+ 從而可得線路始端功率 S12=+ 這個結(jié)果與第（1）%，所以第（2）和第（3）的結(jié)果可作為最終計算結(jié)果；若相差較大，則應返回第（2）步重新計算，直道相差較小為止。180。(RT+jXT)=DST=231。N248。==kV=VA117b點電壓為Vb= (VADVL)+(dVL)=22()2+=，橫分量分別為Pc39。A=176。3)+Sc17 =++580+j129=+=1(180。60S=180。104S =180。因為阻抗已歸算到低壓側(cè)，宜用低壓側(cè)的電壓求環(huán)路電勢。 VA=230。S S39。231。86+j180。3150kWh=2300520kWh線路全年的電能損耗DWL=DPL180。 15 0。%△P FX(n+i1)=V(im)^2E(i)^2F(i)^2。F1(n+i1,w)=0。fprintf(39。end fprintf(39。請輸入平衡節(jié)點號:isb=39。%一：創(chuàng)建節(jié)點導納矩陣for i=1:n1if B1(i,6)==0%不含變壓器的支路p=B1(i,1)。j=0。for i=1:n%對PQ節(jié)點的處理if i~=isbamp。%八：創(chuàng)建Jacbi(雅可比矩陣)Jacbi=zeros(2*n2)。for j=1:nif j~=isbk=k+1。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。%創(chuàng)建DetaS h=0。end end I。for i=1:nif B2(i,6)==3h=h+1。for i=1:nif B2(i,6)==2j=j+1。三、實驗程序：function [e,f,p,q]=flow_out(g,b,kind,e,f)%計算潮流后efpq的終值 s=flow(g,b,kind,e,f)。endfunction s=flow(g,b,kind,e,f)%計算當前ef與規(guī)定的pqv的差值 l=length(e)。for j=1:(2*l2)J(2*i1,j)=s(1,j)。elsepq(1,2*j1)=(g(i,j)*e(i)+b(i,j)*f(i))。end end%數(shù)據(jù)輸入g=[ 0 0 0 0 0 ]。3．了解并掌握PSDBPA電力系統(tǒng)分析程序單線圖和地理接線圖的使用。常用的控制語句主要包括：(1)