【正文】 電力系統(tǒng)潮流計算是研究電力系統(tǒng)穩(wěn)態(tài)運行情況的一種基本計算。6）計算結果輸出：詳細計算結果列表：按節(jié)點、與該節(jié)點相聯(lián)接支路順序，并根據(jù)用戶的要求（通過控制語句控制）可按照字母、分區(qū)或區(qū)域排序輸出潮流計算結果。三、原理和說明PSDBPA電力系統(tǒng)分析程序穩(wěn)態(tài)分析主要是潮流計算，軟件中潮流程序的計算方法有P_Q分解法，牛頓_拉夫遜法，改進的牛頓－拉夫遜算法。e f四、例題及運行結果在上圖所示的簡單電力系統(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ù)課程內容，熟悉MATLAB軟件的使用方法，自行學習MATLAB程序的基礎語法，并根據(jù)所學知識編寫潮流計算牛頓拉夫遜法(或PQ分解法)的計算程序，用相應的算例在MATLAB上進行計算、調試和驗證。232。2SN2246。T2輸入高壓母線的總功率為S39。10kV=232。kV= 231。***SbZII+ZIII+ScZII(+)()+(+)()MVSI==***+ZII+ZIII()=+=*ScZI*+ZIII+SbZI***ZI*+ZII+ZIII()=(+)()+(+)()=++SII=+++=+Sb+Sc=+++=+可見，計算結果誤差很小，無需重算。=(+)180。14)Sc 17=+=+=SAC+SBD=+++=+ SC+SD=++2200+j1200=+ SCD=SBDSD=+=+C點為功率分點，可推算出E點為電壓最低點。試求V功率分布，且找到電壓最低點。==kV=39。=S139。VN248。1102RT==180。母線A的實際運行電壓為117kV，負荷功率：SLDb=30+j12MVA,SLDc=20+j15MVA。P32++=(R23+jX23)=(1+j2)=+++=(R24+jX24)=(+j3)=+則： S23=S3+DS23=+S24=S4+DS24=+39。EF=E+j*F。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。%總的節(jié)點數(shù)FX=ones(2*n2,1)。%PP=PV節(jié)點的P值 V=[]。 0 0。%△P FX(n+i1)=Q(i)。F1(n+i1,w)=F1(n+i1,w)+G(i,k)*F(k)+B(i,k)*E(k)。fprintf(39。%d號節(jié)點電壓的相角度數(shù)為39。(P12R12+Q12X12)180。= 21 XL=180。180。+Sb=+++=+線路中的功率損耗為230。RL+Q139。= (VbDVT)+(dVT)=22()2+=1111Vc=Vc180。(+)=+ 求C點和D點的運算負荷，為 DSCE +=(+)=+ 210=2600+j1600+300+j160++=+SD=600+j200+1600+j1000=2200+j1200kVA循環(huán)功率Sc(V=*A*VBVN()180。+180。104S2DQB213。為了計算線路Ⅰ的電壓損耗，要用A點的電壓和功率SA1。k1248。231。線路和變壓器的參數(shù)如下：線路（每回）：r=, x=, b=180。31500247。1102Mvar=2等值電路中流過線路等值阻抗的功率為S1=S+DST+DS0+jQB2=40+j30++++=++180。)。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ù)已知的電力網的數(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編寫出一個潮流計算的程序。、輸出相關文件 *.dat潮流計算數(shù)據(jù)文件*.bse潮流計算二進制結果文件（可用于潮流計算的輸入或穩(wěn)定計算）*.pfo潮流計算結果文件*.map 供單線圖格式潮流圖及地理接線圖格式潮流圖程序使用的二進制結果文件*.pff，*.pfd 中間文件（正常計算結束后將自動刪除。對應于上述系統(tǒng)及數(shù)據(jù)的潮流計算數(shù)據(jù)（）見例1。從數(shù)學上說，潮流計算是求解一組由潮流方程描述的非線性代數(shù)方程組。* 潮流計算迭代過程和平衡節(jié)點相關信息數(shù)據(jù)計算結果收斂。2）輸入、輸出文件及輸出的內容列表。二、實驗背景隨著科學技術的飛速發(fā)展，電力系統(tǒng)也在不斷地發(fā)展，電網通過互聯(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)。)。8760+252180。=2231。231。VA=，VA=109kV（5）計及從高壓母線輸入變壓器T1和T2的功率。1246。因兩臺變壓器電抗相等，故S1LD=S2LD=11SLD=(16+j12)MVA=8+j6MVA 22（2）求循環(huán)功率。40S=180。106180。7+j1200180。線路LL2和L3導線型號均為LJ185，線路長度分別為10km，4km和3km，線路L4為2km長的LJ70導線；各負荷點負荷amp。180。(RL+jXL)232。202+152247。104S由于線路電壓未知，可用線路額定電壓計算線路產生的充電功率，并將其等分為兩部分，便得112DQB=BcVN=180。+(1+j2)=+ DS23=+(+j3)=+ DS24= S23=S3+DS23=+ S24=S4+DS24=+39。)。,a)。elseF1(i,w)=G(i,i)*E(i)B(i,i)*F(i)。%△Pend endfor i=1:m %PQ節(jié)點for w=1:n1if i~=wF1(i,w)=(G(i,w)*E(i)+B(i,w)*F(i))。00 0。第一篇：電力系統(tǒng)分析潮流計算大作業(yè)電力系統(tǒng)分析潮流計算大作業(yè)（源程序及實驗報告）源程序如下：采用直角坐標系的牛頓拉夫遜迭代 function chaoliujisuan()m=3。 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)。F1(n+i1,n+w1)=0。第%d次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:n39。平衡節(jié)點的功率39。V2DV23=（3）根據(jù)上述求得的線路各點電壓，重新計算各線路的功率損耗和線路始端功率。106S=180。Sc246。247。180。 某一額定電壓為10kV的兩端供電網，如圖所示。3+180。60W=+BI=180。106180。解（1）假定兩臺變壓器變比相同，計算其功率分布。不計電壓降落的橫分量時，按變壓器T1計算可得高壓母線電壓為180。 計及電壓降落的橫分量，按T1和T2計算克分別得。DPs+jsSN247。DP0+j0SN247。3150=172180。請輸入平衡節(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*