【正文】 Jacbi=zeros(2*n2)。h=0。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。end end t=0。for i=1:n%對PQ節(jié)點的處理if i~=isbamp。for i=1:nif B2(i,6)==3t=t+1。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)))。j=0。end end Y。Y(p,q)=Y(p,q)1/(B1(i,3)*B1(i,5))。Y(p,p)=Y(p,p)+1/B1(i,3)+*B1(i,4)。%一：創(chuàng)建節(jié)點導(dǎo)納矩陣for i=1:n1if B1(i,6)==0%不含變壓器的支路p=B1(i,1)。節(jié)點號和對地參數(shù):X=39。B2=input(39。)。請輸入平衡節(jié)點號:isb=39。n1=input(39。四、實驗程序clear。3150kWh=5946885kWh 輸電系統(tǒng)全年的總電能損耗DWT+DWL=2300520+5946885kWh=8247405kWh第三篇：電力系統(tǒng)分析潮流計算的計算機算法潮流計算的計算機算法實驗報告姓名：學(xué)號：班級：一、實驗?zāi)康恼莆粘绷饔嬎愕挠嬎銠C算法。3150=172180。180。線路末端充電功率QB2=2blV2=180。kVA=172+j1701100100232。DP0+j0SN247。231。247。248。DPs+jsSN247。106S/km變壓器（每臺）:DP0=86kW,DPs=200kW,I0%=,Vs%=解 最大負(fù)荷時變壓器的繞組功率損耗為V%230。1+ST22=8++8+=16+ 計算所得功率分布，如圖所示。j1MVA=8+10282+=8++180。 計及電壓降落的橫分量，按T1和T2計算克分別得。10+247。按變壓器T2計算可得180。10+247。不計電壓降落的橫分量時，按變壓器T1計算可得高壓母線電壓為180。故循環(huán)功率為 Sc187。232。247。解（1）假定兩臺變壓器變比相同，計算其功率分布。SA1=SI++(+)MVA=+=++1102DVI=PA1RI+180。SIII=+=SbSI=+=+由此得到功率初分布，如圖所示。1102Mvar=1(+)W=+ 2 2DQBIII 變電所b：ZTb= DS0b=2(+)MVA=+變電所b：ZTc=1(+j127)W=+ 2 DS0c=2(+)MVA=+ 等值電路如圖所示（2）計算節(jié)點b和c的運算負(fù)荷。106180。=180。=180。104180。60W=+BI=180。=VE=VCDVCE== 圖所示110kV閉式電網(wǎng)，A點為某發(fā)電廠的高壓母線，其運行電壓為117kV。 SAC=+++=+DVAC=180。14+180。3+180。10=(+)kVA=580+j129kVA =*17180。(+)=+ZL2=4180。B=176。 某一額定電壓為10kV的兩端供電網(wǎng)，如圖所示。=180。180。180。180。180。+DSL=+++=+由母線A輸出的功率為SA=S139。39。247。S1246。=Sc+DPT+jDQT=20+j15++=+39。247。Sc246。15DPo+jDQo=2180。151Vs%180。1102Mvar= 和b，作為節(jié)點負(fù)荷的一部分。106S=180。80180。解（1）計算參數(shù)并作出等值電路。106S/km。V2DV23=（3）根據(jù)上述求得的線路各點電壓，重新計算各線路的功率損耗和線路始端功率。+180。39。已知各節(jié)點的負(fù)荷功率及參數(shù)：S2=(+)MVA，S3=(+)MVA，S4=(+)MVAZ12=(+)W，Z23=(+)W，Z24=(+)W試求電壓和功率分布。平衡節(jié)點的功率39。,i)disp(angle(EF(i))*180/pi)。%d號節(jié)點電壓的幅值為:39。收斂后的節(jié)點電壓用極坐標(biāo)表示為:39。第%d次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:n39。第%d次迭代后的節(jié)點電壓分別為:n39。for i=1:n1E(i)=E(i)+X(i)。F1(n+i1,n+w1)=2*F(i)。F1(n+i1,n+w1)=0。F1(n+i1,n+w1)=F1(n+i1,n+w1)G(i,k)*E(k)+B(i,k)*F(k)。F1(n+i1,n+w1)=G(i,i)*E(i)+B(i,i)*F(i)。F1(n+i1,n+w1)=G(i,w)*E(i)+B(i,w)*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)。%△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。X=[]。 0。0000 0。%E=PQ,PV,Vθ節(jié)點e的初值F=[0 0 0 0 0]39。%P=PQ節(jié)點的P值 Q=[ 0]。第一篇：電力系統(tǒng)分析潮流計算大作業(yè)電力系統(tǒng)分析潮流計算大作業(yè)（源程序及實驗報告）源程序如下：采用直角坐標(biāo)系的牛頓拉夫遜迭代 function chaoliujisuan()m=3。%Q=PQ節(jié)點的Q值 PP=[]。%F=PQ,PV,Vθ節(jié)點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é)點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é)點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。)。解：（1）先假設(shè)各節(jié)點電壓均為額定電壓，求線路始端功率。P12++=(R12+jX12)=(+)22VN10=+39。=== V2187。+(1+j2)=+ DS23=+(+j3)=+ DS24= S23=S3+DS23=+ S24=S4+DS24=+39。變電所中裝有兩臺三相110/11kV的變壓器，每臺的容量為15MVA,其參數(shù)為：DP0=，DPs=128kW,Vs%=,Io%=。輸電線路的等值電阻、電抗和電納分別為 RL=1180。= Bc=2180。104S由于線路電壓未知，可用線路額定電壓計算線路產(chǎn)生的充電功率，并將其等分為兩部分，便得112DQB=BcVN=180。兩臺變壓器并聯(lián)運行時，它們的等值電阻、電抗及勵磁功率分別為1DPsVN21128180。1102RT==180。(+j)MVA=+變壓器的勵磁功率也作為接于節(jié)點b的負(fù)荷，于是節(jié)點b的負(fù)荷Sb=SLDb+jDQB+(DP0+jDQ0)=30+j12++=+ Sc=20+j15MVA 這樣就得到圖所示的等值電路節(jié)（2）計算母線A輸出的功率。202+152247。110232。39。DSL=231。(RL+jXL)232。2 于是可得 S139。+jDQB=+=+（3）計算各節(jié)點電壓。+180。180。+180。180。kV=11011039。線路LL2和L3導(dǎo)線型號均為LJ185，線路長度分別為10km，4km和3km，線路L4為2km長的LJ70導(dǎo)線；各負(fù)荷點負(fù)荷amp。kV時的初始如圖所示。(+)=+ZL3=3180。()Z229。7+j1200180。10+j1200180。+180。網(wǎng)絡(luò)各組件參數(shù)為：變電所b SNXT==20MVA，DS0=+，RT=，變電所c SNXT=127W=10MVA，DS0=+，RT=，負(fù)荷功率 SLDb=24+j18MVA，SLDc=12+j9MVA試求電力網(wǎng)絡(luò)的功率分布及最大電壓損耗。106180。1102Mvar=線路Ⅱ： Z213。106180。104180。40S=180。DSTb242+182(+)MVA=+ =2110Sb=SLDb+DSTb+DSob+jDQBI+jDQBIII=24+j18++++=++92(+)MVA=+ =1102Sc=SLDc+DSTc+DSoc+jDQBIII+jDQBII=12+j9++++=+（3）計算閉式網(wǎng)絡(luò)的功率分布。（4）計算電壓損耗。+180。因兩臺變壓器電抗相等，故S1LD=S2LD=11SLD=(16+j12)MVA=8+j6MVA 22（2）求循環(huán)功率。=101247。10248。VBDE10180。1246。k1=(10+)180。1246。k2=(10+)180。VA=，VA=109kV（5）計及從高壓母線輸入變壓器T1和T2的功率。j1MVA=8+ 21039。3．6 如圖所示網(wǎng)絡(luò)，變電所低壓母線上的最大負(fù)荷為40MW，cosj=，Tmax=4500h。246。231。232。248。kVA247。=2231。248。106180。180。8760+252180。熟悉MATLAB，并掌握MATLAB程序的基本調(diào)試方法。%清空內(nèi)存n=input