【正文】 7。=Sc+DPT+jDQT=20+j15++=+39。Sc246。151Vs%180。106S=180。解（1）計(jì)算參數(shù)并作出等值電路。V2DV23=（3）根據(jù)上述求得的線路各點(diǎn)電壓，重新計(jì)算各線路的功率損耗和線路始端功率。39。平衡節(jié)點(diǎn)的功率39。%d號(hào)節(jié)點(diǎn)電壓的幅值為:39。第%d次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:n39。for i=1:n1E(i)=E(i)+X(i)。F1(n+i1,n+w1)=0。F1(n+i1,n+w1)=G(i,i)*E(i)+B(i,i)*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)。%F(x)導(dǎo)數(shù)矩陣 a=0。 0。%E=PQ,PV,Vθ節(jié)點(diǎn)e的初值F=[0 0 0 0 0]39。第一篇：電力系統(tǒng)分析潮流計(jì)算大作業(yè)電力系統(tǒng)分析潮流計(jì)算大作業(yè)（源程序及實(shí)驗(yàn)報(bào)告）源程序如下：采用直角坐標(biāo)系的牛頓拉夫遜迭代 function chaoliujisuan()m=3。%F=PQ,PV,Vθ節(jié)點(diǎn)f的初值 G=[ 0 0。00 0。%記錄迭代次數(shù)EF=zeros(n1,n1)。%△Pend endfor i=1:m %PQ節(jié)點(diǎn)for w=1:n1if i~=wF1(i,w)=(G(i,w)*E(i)+B(i,w)*F(i))。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)。F(i)=F(i)+X(n+i1)。,a)。,i)disp(abs(EF(i)))。)。P12++=(R12+jX12)=(+)22VN10=+39。+(1+j2)=+ DS23=+(+j3)=+ DS24= S23=S3+DS23=+ S24=S4+DS24=+39。輸電線路的等值電阻、電抗和電納分別為 RL=1180。104S由于線路電壓未知，可用線路額定電壓計(jì)算線路產(chǎn)生的充電功率，并將其等分為兩部分，便得112DQB=BcVN=180。1102RT==180。202+152247。39。(RL+jXL)232。+jDQB=+=+（3）計(jì)算各節(jié)點(diǎn)電壓。180。180。線路LL2和L3導(dǎo)線型號(hào)均為L(zhǎng)J185，線路長(zhǎng)度分別為10km，4km和3km，線路L4為2km長(zhǎng)的LJ70導(dǎo)線；各負(fù)荷點(diǎn)負(fù)荷amp。(+)=+ZL3=3180。7+j1200180。+180。106180。106180。40S=180。（4）計(jì)算電壓損耗。因兩臺(tái)變壓器電抗相等，故S1LD=S2LD=11SLD=(16+j12)MVA=8+j6MVA 22（2）求循環(huán)功率。10248。1246。1246。VA=，VA=109kV（5）計(jì)及從高壓母線輸入變壓器T1和T2的功率。3．6 如圖所示網(wǎng)絡(luò)，變電所低壓母線上的最大負(fù)荷為40MW，cosj=，Tmax=4500h。231。248。=2231。106180。8760+252180。%清空內(nèi)存n=input(39。)。請(qǐng)輸入節(jié)點(diǎn)參數(shù):B2=39。q=B1(i,2)。Y(q,p)=Y(p,q)。for i=1:n%對(duì)PQ節(jié)點(diǎn)的處理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%非對(duì)角元素的處理Jacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。if i==j%對(duì)角元素的處理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%對(duì)PV節(jié)點(diǎn)的處理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)。通過(guò)實(shí)驗(yàn)教學(xué)加深學(xué)生對(duì)復(fù)雜電力系統(tǒng)潮流計(jì)算計(jì)算方法的理解，學(xué)會(huì)運(yùn)用電力系統(tǒng)的數(shù)學(xué)模型，掌握潮流計(jì)算的過(guò)程及其特點(diǎn)，熟悉各種常用應(yīng)用軟件，熟悉硬件設(shè)備的使用方法，加強(qiáng)編制調(diào)試計(jì)算機(jī)程序的能力，提高工程計(jì)算的能力，學(xué)習(xí)如何將理論知識(shí)和實(shí)際工程問(wèn)題結(jié)合起來(lái)。k=0。endl=length(e)。s=zeros(2*l2,1)。end endfunction J=J_out(g,b,kind,e,f)%計(jì)算節(jié)點(diǎn)的雅克比矩陣 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 ]。高斯迭代法是直接迭代，對(duì)初值要求比較低，程序簡(jiǎn)單，內(nèi)存小，但收斂性差，速度慢，多用于配電網(wǎng)或輻射式網(wǎng)絡(luò)中；牛頓拉夫遜法是將非線性方程線性化之后再迭代的，對(duì)初值要求比較高，收斂性好，速度快，迭代次數(shù)少，運(yùn)行時(shí)間短，被廣泛使用；PQ分解法是在極坐標(biāo)牛頓法的基礎(chǔ)上進(jìn)行三個(gè)簡(jiǎn)化所得，有功、無(wú)功分開迭代，迭代次數(shù)比牛頓多一倍但運(yùn)算量小，整體速度更快，運(yùn)行時(shí)間更短，多用于110KV以上的高壓電網(wǎng)中，電力網(wǎng)絡(luò)的節(jié)點(diǎn)分幾類？各類節(jié)點(diǎn)的已知量和待求量是什么？答: PQ節(jié)點(diǎn)：P、Q為已知量，V、d為待求量；PV節(jié)點(diǎn)：給定P、V，求Q、d；平衡節(jié)點(diǎn)：給定V、d，求P、Q。二、實(shí)驗(yàn)背景隨著科學(xué)技術(shù)的飛速發(fā)展，電力系統(tǒng)也在不斷地發(fā)展，電網(wǎng)通過(guò)互聯(lián)變得越來(lái)越復(fù)雜，同時(shí)也使系統(tǒng)穩(wěn)定問(wèn)題越來(lái)越突出?？蛇M(jìn)行交流系統(tǒng)潮流計(jì)算，也可進(jìn)行包括雙端和多端直流系統(tǒng)的交直流混合潮流計(jì)算。2）輸入、輸出文件及輸出的內(nèi)容列表。IEEE 9節(jié)點(diǎn)例題：圖1 IEEE9節(jié)點(diǎn)系統(tǒng)接線圖節(jié)點(diǎn)參數(shù)、線路參數(shù)及變壓器參數(shù)分別見表1～表3。* 潮流計(jì)算迭代過(guò)程和平衡節(jié)點(diǎn)相關(guān)信息數(shù)據(jù)計(jì)算結(jié)果收斂。潮流計(jì)算的方法有最基本的手算迭代方法，而利用電子計(jì)算機(jī)進(jìn)行潮流計(jì)算從20世紀(jì)50年代中期就已經(jīng)開始。從數(shù)學(xué)上說(shuō)，潮流計(jì)算是求解一組由潮流方程描述的非線性代數(shù)方程組。各區(qū)域平衡機(jī)出力數(shù)據(jù)列表區(qū)域平衡機(jī)電壓額定有功有功出力無(wú)功出力有功負(fù)荷無(wú)功負(fù)荷所屬分區(qū)SYSTEM發(fā)電機(jī)1 01* 沒有遇到錯(cuò)誤信息 23:03:48 3）單線圖：（二）課本習(xí)題：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é)的變壓器等。對(duì)應(yīng)于上述系統(tǒng)及數(shù)據(jù)的潮流計(jì)算數(shù)據(jù)（）見例1。如為致命性錯(cuò)誤，則中斷計(jì)算。、輸出相關(guān)文件 *.dat潮流計(jì)算數(shù)據(jù)文件*.bse潮流計(jì)算二進(jìn)制結(jié)果文件（可用于潮流計(jì)算的輸入或穩(wěn)定計(jì)算）*.pfo潮流計(jì)算結(jié)果文件*.map 供單線圖格式潮流圖及地理接線圖格式潮流圖程序使用的二進(jìn)制結(jié)果文件*.pff，*.pfd 中間文件（正常計(jì)算結(jié)束后將自動(dòng)刪除。PSDBPA軟件包主要由潮流和暫穩(wěn)程序構(gòu)成，具有計(jì)算規(guī)模大、計(jì)算速度快、數(shù)值穩(wěn)定性好、功能強(qiáng)等特點(diǎn)，已在我國(guó)電力系統(tǒng)規(guī)劃、調(diào)度、生產(chǎn)運(yùn)行及科研部門得到了廣泛應(yīng)用。六、實(shí)驗(yàn)心得這次實(shí)驗(yàn)是通過(guò)matlab編寫出一個(gè)潮流計(jì)算的程序。f=[0 0 0 0]。endpv(1,2*i)=s+b(i,i)*e(i)g(i,i)*f(i)。pq(2,2*j1)=pq(1,2*j)。for k=1:ls=s+(g(i,k)*f(k)+b(i,k)*e(k))。pq=zeros(2,2*l2)。for i=1:(l1)。s2=0。s2=0。J_ni=inv(J)。程序要求根據(jù)已知的電力網(wǎng)的數(shù)學(xué)模型（節(jié)點(diǎn)導(dǎo)納矩陣）及各節(jié)點(diǎn)參數(shù)，完成該電力系統(tǒng)的潮流計(jì)算，要求計(jì)算出節(jié)點(diǎn)電壓、功率等參數(shù)。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。endendendend end Jacbi。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))。k=0。end end DetaS。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。endend end for i=1:nif i~=isbamp。end end B2。%九：求解修正方程，獲取節(jié)點(diǎn)電壓的不平衡量