【正文】 III 變電所b：ZTb= DS0b=2(+)MVA=+變電所b：ZTc=1(+j127)W=+ 2 DS0c=2(+)MVA=+ 等值電路如圖所示（2）計(jì)算節(jié)點(diǎn)b和c的運(yùn)算負(fù)荷。SIII=+=SbSI=+=+由此得到功率初分布，如圖所示。SA1=SI++(+)MVA=+=++1102DVI=PA1RI+180。解（1）假定兩臺(tái)變壓器變比相同，計(jì)算其功率分布。247。232。故循環(huán)功率為 Sc187。不計(jì)電壓降落的橫分量時(shí)，按變壓器T1計(jì)算可得高壓母線電壓為180。10+247。按變壓器T2計(jì)算可得180。10+247。 計(jì)及電壓降落的橫分量，按T1和T2計(jì)算克分別得。j1MVA=8+10282+=8++180。1+ST22=8++8+=16+ 計(jì)算所得功率分布，如圖所示。106S/km變壓器（每臺(tái)）:DP0=86kW,DPs=200kW,I0%=,Vs%=解 最大負(fù)荷時(shí)變壓器的繞組功率損耗為V%230。DPs+jsSN247。248。247。231。DP0+j0SN247。kVA=172+j1701100100232。線路末端充電功率QB2=2blV2=180。180。3150=172180。3150kWh=5946885kWh 輸電系統(tǒng)全年的總電能損耗DWT+DWL=2300520+5946885kWh=8247405kWh第二篇：電力系統(tǒng)分析潮流計(jì)算大作業(yè)電力系統(tǒng)分析潮流計(jì)算大作業(yè)（源程序及實(shí)驗(yàn)報(bào)告）源程序如下：采用直角坐標(biāo)系的牛頓拉夫遜迭代 function chaoliujisuan()m=3。%Q=PQ節(jié)點(diǎn)的Q值 PP=[]。%F=PQ,PV,Vθ節(jié)點(diǎn)f的初值 G=[ 0 0。0000 0]。00 0。%X=△X n=m+v+1。%記錄迭代次數(shù)EF=zeros(n1,n1)。%△PFX(n+i1)=FX(n+i1)F(i)*G(i,w)*E(w)+F(i)*B(i,w)*F(w)+E(i)*G(i,w)*F(w)+E(i)*B(i,w)*E(w)。%△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))。elseF1(i,w)=G(i,i)*E(i)B(i,i)*F(i)。for k=1:nF1(i,w)=F1(i,w)G(i,k)*E(k)+B(i,k)*F(k)。endendend end for i=m+1:n1 %PV節(jié)點(diǎn)for w=1:n1if i~=wF1(i,w)=(G(i,w)*E(i)+B(i,w)*F(i))。elseF1(i,w)=G(i,i)*E(i)B(i,i)*F(i)。for k=1:nF1(i,w)=F1(i,w)G(i,k)*E(k)+B(i,k)*F(k)。F(i)=F(i)+X(n+i1)。,a)。,a)。)。,i)disp(abs(EF(i)))。end PPH=0。)。二、實(shí)驗(yàn)準(zhǔn)備根據(jù)課程內(nèi)容，熟悉MATLAB軟件的使用方法，自行學(xué)習(xí)MATLAB程序的基礎(chǔ)語法，并根據(jù)所學(xué)知識(shí)編寫潮流計(jì)算牛頓拉夫遜法(或PQ分解法)的計(jì)算程序，用相應(yīng)的算例在MATLAB上進(jìn)行計(jì)算、調(diào)試和驗(yàn)證。請(qǐng)輸入節(jié)點(diǎn)數(shù):n=39。)。pr=input(39。請(qǐng)輸入支路參數(shù):B1=39。)。Y=zeros(n)。Y(p,q)=Y(p,q)1/B1(i,3)。else%含有變壓器的支路p=B1(i,1)。Y(p,p)=Y(p,p)+1/B1(i,3)。DetaS=zeros(2*n2,1)。B2(i,6)==2h=h+1。B2(i,6)==3h=h+1。%四：創(chuàng)建PVU 用于存儲(chǔ)PV節(jié)點(diǎn)的初始電壓PVU=zeros(nh1,1)。end end PVU。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。B2(i,6)==3h=h+1。end end DetaS。I(h,1)=(OrgS(2*h1,1)OrgS(2*h,1)*sqrt(1))/conj(B2(i,3))。k=0。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))。endendendend end k=0。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))。endendendend end Jacbi。%修正節(jié)點(diǎn)電壓 j=0。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。j=0。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)))。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)))。B2(i,6)==2h=h+1。for i=1:nif i~=isbamp。DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2real(B2(i,3))^2imag(B2(i,3))^2。for i=1:nif i~=isbh=h+1。h=0。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))。elseJacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。endif k==(n1)k=0。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))。elseJacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。endif k==(n1)k=0。DetaU。end end for i=1:nif B2(i,6)==3j=j+1。%迭代次數(shù)加1 end Times。二、實(shí)驗(yàn)內(nèi)容：編制調(diào)試電力系統(tǒng)潮流計(jì)算的計(jì)算機(jī)程序。在相應(yīng)的編程環(huán)境下對(duì)程序進(jìn)行組織調(diào)試。while max(abs(s))10^5 J=J_out(g,b,kind,e,f)。for i=1:le(i)=e(i)dv(2*i1)。for i=1:l s1=0。endp(i)=e(i)*s1+f(i)*s2。for i=1:(l1)s1=0。ends(2*i1)=kind(2,i)e(i)*s1f(i)*s2。J=zeros(2*l2,2*l2)。J(2*i,j)=s(2,j)。endend endfunction pq=PQ_out(g,b,e,f,i)%計(jì)算pq節(jié)點(diǎn)的雅克比矩陣 l=length(e)。endpq(1,2*i1)=sg(i,i)*e(i)b(i,i)*f(i)。s=0。for k=1:ls=s(g(i,k)*e(k)b(i,k)*f(k))。pq(2,2*j)=pq(1,2*j1)。for j=1:(l1)if j==i s=0。for k=1:ls=s(g(i,k)*f(k)+b(i,k)*e(k))。elsepv(1,2*j1)=(g(i,j)*e(i)+b(i,j)*f(i))。e=[1 1 ]。e f四、例題及運(yùn)行結(jié)果在上圖所示的簡(jiǎn)單電力系統(tǒng)中，系統(tǒng)中節(jié)點(diǎn)2為PQ節(jié)點(diǎn)，節(jié)點(diǎn)3為PV節(jié)點(diǎn)，節(jié)點(diǎn)4為平衡節(jié)點(diǎn)，已給定 P1s+jQ1s= P2s+jQ2s= P3s= V3s= V4s=∠0176。？為什么？答:不一樣,因?yàn)槊看蔚碾妷?、有功、無功都是與前一次不同的新值，所以每次迭代過程中，雅可比矩陣都是變化的。第五篇：電力系統(tǒng)分析計(jì)算機(jī)算法PSDBPA潮流計(jì)算實(shí)驗(yàn)報(bào)告電力系統(tǒng)分析的計(jì)算機(jī)算法實(shí)驗(yàn)報(bào)告學(xué)生姓名 課 程 電力系統(tǒng)分析的計(jì)算機(jī)算法 學(xué) 號(hào)專 業(yè) 電氣工程及其自動(dòng)化 指導(dǎo)教師 邱曉燕二Ο一四 年 六 月 二日實(shí)驗(yàn)一潮流計(jì)算一、實(shí)驗(yàn)?zāi)康?．了解并掌握電力系統(tǒng)計(jì)算機(jī)算法的相關(guān)原理。無論是電力系統(tǒng)規(guī)劃、設(shè)計(jì)還是運(yùn)行，對(duì)其安全穩(wěn)定進(jìn)行分析都是極其重要的。三、原理和說明PSDBPA電力系統(tǒng)分析程序穩(wěn)態(tài)分析主要是潮流計(jì)算，軟件中潮流程序的計(jì)算方法有P_Q分解法，牛頓_拉夫遜法，改進(jìn)的牛頓－拉夫遜算法。除了潮流計(jì)算功能外，該軟件還具有自動(dòng)電壓控制、聯(lián)絡(luò)線功率控制、系統(tǒng)事故分析（N1開斷模擬）、網(wǎng)絡(luò)等值、靈敏度分析、節(jié)點(diǎn)P－V、Q－V和P－Q曲線、確定系統(tǒng)極限輸送水平、負(fù)荷靜特性模型、靈活多樣的分析報(bào)告、詳細(xì)的檢錯(cuò)功能等功能。該文件安裝時(shí)放在與潮流程序相同的目錄中。3）錯(cuò)誤信息。6）計(jì)算結(jié)果輸出：詳細(xì)計(jì)算結(jié)果列表：按節(jié)點(diǎn)、與該節(jié)點(diǎn)相聯(lián)接支路順序，并根據(jù)用戶的要求（通過控制語句控制）可按照字母、分區(qū)或區(qū)域排序輸出潮流計(jì)算結(jié)果。表1 IEEE 9節(jié)點(diǎn)算例節(jié)點(diǎn)參數(shù)表2 IEEE 9節(jié)點(diǎn)算例線路參數(shù)表3 IEEE 9節(jié)點(diǎn)算例變壓器參數(shù)注：表1表3中功率基準(zhǔn)值為100MVA；電阻、電感值為標(biāo)幺值。當(dāng)最大出力值小于實(shí)際出力時(shí)，統(tǒng)計(jì)時(shí)最大出力值用實(shí)際出力值代替。牛頓拉夫遜法迭代次數(shù)為 5次。電力系統(tǒng)潮流計(jì)算是研究電力系統(tǒng)穩(wěn)態(tài)運(yùn)行情況的一種基本計(jì)算。此后，潮流計(jì)算曾采用了各種不同的方法，這些方法的發(fā)展主要是圍繞著對(duì)潮流計(jì)算的一些基本要求進(jìn)