【正文】 2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)))。j=0。Y(p,q)=Y(p,q)1/(B1(i,3)*B1(i,5))。%一：創(chuàng)建節(jié)點導納矩陣for i=1:n1if B1(i,6)==0%不含變壓器的支路p=B1(i,1)。B2=input(39。請輸入平衡節(jié)點號:isb=39。四、實驗程序clear。end fprintf(39。for i=1:n1 fprintf(39。fprintf(39。endendend end X=inv(F1)*(FX)。F1(n+i1,w)=0。F1(n+i1,w)=B(i,i)*E(i)G(i,i)*F(i)。%△P FX(n+i1)=V(im)^2E(i)^2F(i)^2。%F(x)矩陣F1=zeros(n1,n1)。 15 0。%V=PV節(jié)點的U值E=[1 1 1 ]39。3150kWh=2300520kWh線路全年的電能損耗DWL=DPL180。100180。86+j180。2+j180。231。試求線路和變壓器全年的電能損耗。S S39。 VA=230。 VA=230。232。因為阻抗已歸算到低壓側(cè)，宜用低壓側(cè)的電壓求環(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=+ 這個結(jié)果與第（1）%，所以第（2）和第（3）的結(jié)果可作為最終計算結(jié)果；若相差較大，則應返回第（2）步重新計算，直道相差較小為止。故： S12=S12+DS12=+（2）再用已知的線路始端電壓V1=12，求出線路各點電壓。第一篇：電力系統(tǒng)分析潮流計算例題3．1 電網(wǎng)結(jié)構(gòu)如圖3—11所示，其額定電壓為10KV。(P12R12+Q12X12)180。 如圖所示簡單系統(tǒng)，額定電壓為110KV 雙回輸電線路，長度為80km，采用LGJ150導線，其單位長度的參數(shù)為：r=，x=,b=180。= 21 XL=180。104180。180。(+)MVA2231。+Sb=+++=+線路中的功率損耗為230。+=(+)MVA=+39。RL+Q139。RT+Qc39。= (VbDVT)+(dVT)=22()2+=1111Vc=Vc180。kV、Vamp。(+)=+ 求C點和D點的運算負荷，為 DSCE +=(+)=+ 210=2600+j1600+300+j160++=+SD=600+j200+1600+j1000=2200+j1200kVA循環(huán)功率Sc(V=*A*VBVN()180。10+2200180。+180。104S2DQBI=180。104S2DQB213。104180。為了計算線路Ⅰ的電壓損耗，要用A點的電壓和功率SA1。若取其假定正方向為順時針方向，則可得k2 DE187。k1248。231。231。T182+=8++180。線路和變壓器的參數(shù)如下：線路（每回）：r=, x=, b=180。100232。31500247。31500kVA 247。1102Mvar=2等值電路中流過線路等值阻抗的功率為S1=S+DST+DS0+jQB2=40+j30++++=++180。3150=180。%E=PQ,PV,Vθ節(jié)點e的初值F=[0 0 0 0 0]39。 0。%F(x)導數(shù)矩陣 a=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)=G(i,i)*E(i)+B(i,i)*F(i)。F1(n+i1,n+w1)=0。for i=1:n1E(i)=E(i)+X(i)。第%d次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:n39。%d號節(jié)點電壓的幅值為:39。平衡節(jié)點的功率39。%清空內(nèi)存n=input(39。)。請輸入節(jié)點參數(shù):B2=39。q=B1(i,2)。Y(q,p)=Y(p,q)。for i=1:n%對PQ節(jié)點的處理if i~=isbamp。endend end OrgS。B2(i,6)==2h=h+1。DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2real(B2(i,3))^2imag(B2(i,3))^2。h=0。else%非對角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。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))。endif k==(n1)%將用于內(nèi)循環(huán)的指針置于初始值，以確保雅可比矩陣換行k=0。end end for i=1:n%對PV節(jié)點的處理if B2(i,6)==3j=j+1。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)))。for i=1:nif i~=isbamp。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。%創(chuàng)建Jacbi Jacbi=zeros(2*n2)。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)2*imag(I(h,1))。for j=1:nif j~=isbk=k+1。Jacbi(2*h,2*k)=0。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。通過實驗教學加深學生對復雜電力系統(tǒng)潮流計算計算方法的理解，學會運用電力系統(tǒng)的數(shù)學模型，掌握潮流計算的過程及其特點，熟悉各種常用應用軟件，熟悉硬件設(shè)備的使用方法，加強編制調(diào)試計算機程序的能力，提高工程計算的能力，學習如何將理論知識和實際工程問題結(jié)合起來。k=0。endl=length(e)。s=zeros(2*l2,1)。end endfunction J=J_out(g,b,kind,e,f)%計算節(jié)點的雅克比矩陣 l=length(e)。J(2*i,j)=s(2,j)。endpq(1,2*i)=s+b(i,i)*e(i)g(i,i)*f(i)。pq(1,2*j)=b(i,j)*e(i)g(i,j)*f(i)。s=0。b=[ 0 0 0 ]。高斯迭代法是直接迭代，對初值要求比較低，程序簡單，內(nèi)存小，但收斂性差，速度慢，多用于配電網(wǎng)或輻射式網(wǎng)絡中；牛頓拉夫遜法是將非線性方程線性化之后再迭代的，對初值要求比較高，收斂性好，速度快，迭代次數(shù)少，運行時間短，被廣泛使用；PQ分解法是在極坐標牛頓法的基礎(chǔ)上進行三個簡化所得，有功、無功分開迭代，迭代次數(shù)比牛頓多一倍但運算量小，整體速度更快，運行時間更短，多用于110KV以上的高壓電網(wǎng)中，電力網(wǎng)絡的節(jié)點分幾類？各類節(jié)點的已知量和待求量是什么？答: PQ節(jié)點：P、Q為已知量，V、d為待求量；PV節(jié)點：給定P、V，求Q、d；平衡節(jié)點：給定V、d，求P、Q。二、實驗背景隨著科學技術(shù)的飛速發(fā)展，電力系統(tǒng)也在不斷地發(fā)展，電網(wǎng)通過互聯(lián)變得越來越復雜，同時也使系統(tǒng)穩(wěn)定問題越來越突出?？蛇M行交流系統(tǒng)潮流計算，也可進行包括雙端和多端直流系統(tǒng)的交直流混合潮流計算。2）輸入、輸出文件及輸出的內(nèi)容列表。IEEE 9節(jié)點例題：圖1 IEEE9節(jié)點系統(tǒng)接線圖節(jié)點參數(shù)、線路參數(shù)及變壓器參數(shù)分別見表1～表3。* 潮流計算迭代過程和平衡節(jié)點相關(guān)信息數(shù)據(jù)計算結(jié)果收斂。潮流計算的方法有最基本的手算迭代方法，而利用電子計算機進行潮流計算從20世紀50年代中期就已經(jīng)開始。從數(shù)學上說，潮流計算是求解一組由潮流方程描述的非線性代數(shù)方程組。各區(qū)域平衡機出力數(shù)據(jù)列表區(qū)域平衡機電壓額定有功有功出力無功出力有功負荷無功負荷所屬分區(qū)SYSTEM發(fā)電機1 01* 沒有遇到錯誤信息 23:03:48 3）單線圖：（二）課本習題：E25 ：：(POWERFLOW,CASEID=IEEE9,PROJECT=IEEE_9BUS_TEST_SYSTEM)/SOL_ITER,DECOUPLED=2,NEWTON=15,OPITM=0 /P_OUTPUT_LIST,ZONES=ALL /RPT_SORT=ZONE /NEW_BASE,FILE= /PF_MAP,FILE = /母線4999999B母線1 B母線2 BE母線3 999.L支路數(shù)據(jù)L母線1母線2 L母線2母線4 L母線4母線1 .T變壓器數(shù)據(jù)，包括普通變壓器、移相器、帶調(diào)節(jié)的變壓器等。對應于上述系統(tǒng)及數(shù)據(jù)的潮流計算數(shù)據(jù)（）見例1。如為致命性錯誤，則中斷計算。、輸出相關(guān)文件 *.dat潮流計算數(shù)據(jù)文件*.bse潮流計算二進制結(jié)果文件（可用于潮流計算的輸入或穩(wěn)定計算）*.pfo潮流計算結(jié)果文件*.map 供單線圖格式潮流圖及地理接線圖格式潮流圖程序使用的二進制結(jié)果文件*.pff，*.pfd 中間文件（正常計算結(jié)束后將自動刪除。PSDBPA軟件包主要由潮流和暫穩(wěn)程序構(gòu)成，具有計算規(guī)模大、計算速度快、數(shù)值穩(wěn)定性好、功能強等特點，已在我國電力系統(tǒng)規(guī)劃、調(diào)度、生產(chǎn)運行及科研部門得到了廣泛應用。六、實驗心得這次實驗是通過matlab編寫出一個潮流