【正文】 +247。=232。 計(jì)及電壓降落的橫分量，按T1和T2計(jì)算克分別得。S S39。j1MVA=8+10282+=8++180。T2輸入高壓母線的總功率為S39。1+ST22=8++8+=16+ 計(jì)算所得功率分布，如圖所示。試求線路和變壓器全年的電能損耗。106S/km變壓器（每臺(tái)）:DP0=86kW,DPs=200kW,I0%=,Vs%=解 最大負(fù)荷時(shí)變壓器的繞組功率損耗為V%230。230。DPs+jsSN247。231。248。2SN2246。247。2+j180。231。kVA=252+j41661002180。DP0+j0SN247。86+j180。kVA=172+j1701100100232。232。線路末端充電功率QB2=2blV2=180。100180。180。100MW= DPL=2RL=2V1102已知cosj=，Tmax=4500h，從表中查得t=3150h，假定變壓器全年投入運(yùn)行，則變壓器全年的電能損耗 DWT=2DP0180。3150=172180。3150kWh=2300520kWh線路全年的電能損耗DWL=DPL180。3150kWh=5946885kWh 輸電系統(tǒng)全年的總電能損耗DWT+DWL=2300520+5946885kWh=8247405kWh第三篇：電力系統(tǒng)分析潮流計(jì)算的計(jì)算機(jī)算法潮流計(jì)算的計(jì)算機(jī)算法實(shí)驗(yàn)報(bào)告姓名：學(xué)號(hào)：班級(jí)：一、實(shí)驗(yàn)?zāi)康恼莆粘绷饔?jì)算的計(jì)算機(jī)算法。二、實(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)證。四、實(shí)驗(yàn)程序clear。請(qǐng)輸入節(jié)點(diǎn)數(shù):n=39。n1=input(39。)。請(qǐng)輸入平衡節(jié)點(diǎn)號(hào):isb=39。pr=input(39。)。請(qǐng)輸入支路參數(shù):B1=39。B2=input(39。)。節(jié)點(diǎn)號(hào)和對(duì)地參數(shù):X=39。Y=zeros(n)。%一：創(chuàng)建節(jié)點(diǎn)導(dǎo)納矩陣for i=1:n1if B1(i,6)==0%不含變壓器的支路p=B1(i,1)。Y(p,q)=Y(p,q)1/B1(i,3)。Y(p,p)=Y(p,p)+1/B1(i,3)+*B1(i,4)。else%含有變壓器的支路p=B1(i,1)。Y(p,q)=Y(p,q)1/(B1(i,3)*B1(i,5))。Y(p,p)=Y(p,p)+1/B1(i,3)。end end Y。DetaS=zeros(2*n2,1)。j=0。B2(i,6)==2h=h+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)))。B2(i,6)==3h=h+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)))。%四：創(chuàng)建PVU 用于存儲(chǔ)PV節(jié)點(diǎn)的初始電壓PVU=zeros(nh1,1)。for i=1:nif B2(i,6)==3t=t+1。end end PVU。for i=1:n%對(duì)PQ節(jié)點(diǎn)的處理if i~=isbamp。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。end end t=0。B2(i,6)==3h=h+1。DetaS(2*h1,1)=real(B2(i,2))OrgS(2*h1,1)。end end DetaS。h=0。I(h,1)=(OrgS(2*h1,1)OrgS(2*h,1)*sqrt(1))/conj(B2(i,3))。%八：創(chuàng)建Jacbi(雅可比矩陣)Jacbi=zeros(2*n2)。k=0。for j=1:nif j~=isbk=k+1。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*h,2*k)=Jacbi(2*h1,2*k1)2*imag(I(h,1))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k)=Jacbi(2*h1,2*k1)。endendendend end k=0。for j=1:nif j~=isbk=k+1。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*h,2*k)=2*real(B2(i,3))。Jacbi(2*h1,2*k)=real(Y(i,j))*real(B2(i,3))+imag(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k)=0。endendendend end Jacbi。DetaU=inv(Jacbi)*DetaS。%修正節(jié)點(diǎn)電壓 j=0。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j1,1)*sqrt(1)。%十：開始循環(huán)********************************************************************** while abs(max(DetaU))pr OrgS=zeros(2*n2,1)。j=0。B2(i,6)==2h=h+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)))。B2(i,6)==3h=h+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)))。%創(chuàng)建DetaS h=0。B2(i,6)==2h=h+1。DetaS(2*h,1)=imag(B2(i,2))OrgS(2*h,1)。for i=1:nif i~=isbamp。t=t+1。DetaS(2*h,1)=real(PVU(t,1))^2+imag(PVU(t,1))^2real(B2(i,3))^2imag(B2(i,3))^2。%創(chuàng)建I i=zeros(n1,1)。for i=1:nif i~=isbh=h+1。end end I。h=0。for i=1:nif B2(i,6)==2h=h+1。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))。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)+2*real(I(h,1))。elseJacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k1)=Jacbi(2*h1,2*k)。endif k==(n1)k=0。for i=1:nif B2(i,6)==3h=h+1。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))。Jacbi(2*h,2*k1)=2*imag(B2(i,3))。elseJacbi(2*h1,2*k1)=imag(Y(i,j))*real(B2(i,3))+real(Y(i,j))*imag(B2(i,3))。Jacbi(2*h,2*k1)=0。endif k==(n1)k=0。DetaU=zeros(2*n2,1)。DetaU。for i=1:nif B2(i,6)==2j=j+1。end end for i=1:nif B2(i,6)==3j=j+1。end end B2。%迭代次數(shù)加1 end Times。第四篇：電力系統(tǒng)分析潮流實(shí)驗(yàn)報(bào)告南昌大學(xué)實(shí)驗(yàn)報(bào)告電力系統(tǒng)潮流計(jì)算實(shí)驗(yàn) 學(xué)生姓名： 學(xué) 號(hào)： 專業(yè)班級(jí)： 實(shí)驗(yàn)類型：□ 驗(yàn)證 □ 綜合 ■ 設(shè)計(jì) □ 創(chuàng)新 實(shí)驗(yàn)日期： 實(shí)驗(yàn)成績(jī)：一、實(shí)驗(yàn)?zāi)康模罕緦?shí)驗(yàn)通過對(duì)電力系統(tǒng)潮流計(jì)算的計(jì)算機(jī)程序的編制與調(diào)試，獲得對(duì)復(fù)雜電力系統(tǒng)進(jìn)行潮流計(jì)算的計(jì)算機(jī)程序，使系統(tǒng)潮流計(jì)算能夠由計(jì)算機(jī)自行完成，即根據(jù)已知的電力網(wǎng)的數(shù)學(xué)模型（節(jié)點(diǎn)導(dǎo)納矩陣）及各節(jié)點(diǎn)參數(shù)，由計(jì)算程序運(yùn)行完成該電力系統(tǒng)的潮流計(jì)算。二、實(shí)驗(yàn)內(nèi)容：編制調(diào)試電力系統(tǒng)潮流計(jì)算的計(jì)算機(jī)程序。在各種潮流計(jì)算的算法中選擇一種，按照計(jì)算方法編制程序。在相應(yīng)的編程環(huán)境下對(duì)程序進(jìn)行組織調(diào)試。三、實(shí)驗(yàn)程序：function [e,f,p,q]=flow_out(g,b,kind,e,f)%計(jì)算潮流后efpq的終值 s=flow(g,b,kind,e,f)。while max(abs(s))10^5 J=J_out(g,b,kind,e,f)。dv=J_ni*s。for i=1:le(i)=e(i)dv(2*i1)。ends=flow(g,b,kind,e,f)。for i=1:l s1=0。for j=1:ls1=s1+g(i,j)*e(j)b(i,j)*f(j)。endp(i)=e(i)*s1+f(i)*s2。endfunction s=flow(g,b,kind,e,f)%計(jì)算當(dāng)前ef與規(guī)定的pqv的