【正文】 80。31500247。231。kVA247。kVA=252+j41661002180。I%=2231。DP0+j0SN247。=2231。86+j180。31500kVA 247。kVA=172+j1701100100232。248。232。248。線路末端充電功率QB2=2blV2=180。106180。100180。1102Mvar=2等值電路中流過線路等值阻抗的功率為S1=S+DST+DS0+jQB2=40+j30++++=++180。180。180。100MW= DPL=2RL=2V1102已知cosj=，Tmax=4500h，從表中查得t=3150h，假定變壓器全年投入運(yùn)行，則變壓器全年的電能損耗 DWT=2DP0180。8760+DPT180。3150=172180。8760+252180。3150kWh=2300520kWh線路全年的電能損耗DWL=DPL180。3150=180。3150kWh=5946885kWh 輸電系統(tǒng)全年的總電能損耗DWT+DWL=2300520+5946885kWh=8247405kWh第二篇：電力系統(tǒng)分析潮流計算大作業(yè)電力系統(tǒng)分析潮流計算大作業(yè)（源程序及實(shí)驗(yàn)報告）源程序如下：采用直角坐標(biāo)系的牛頓拉夫遜迭代 function chaoliujisuan()m=3。%m=PQ節(jié)點(diǎn)個數(shù) v=1。%v=PV節(jié)點(diǎn)個數(shù)P=[ 0]。%P=PQ節(jié)點(diǎn)的P值 Q=[ 0]。%Q=PQ節(jié)點(diǎn)的Q值 PP=[]。%PP=PV節(jié)點(diǎn)的P值 V=[]。%V=PV節(jié)點(diǎn)的U值E=[1 1 1 ]39。%E=PQ,PV,Vθ節(jié)點(diǎn)e的初值F=[0 0 0 0 0]39。%F=PQ,PV,Vθ節(jié)點(diǎn)f的初值 G=[ 0 0。 0 0。 0 0。0000 0。0000 0]。B=[ 00。 15 0。 0。00 0。000]。Y=G+j*B。X=[]。%X=△X n=m+v+1。%總的節(jié)點(diǎn)數(shù)FX=ones(2*n2,1)。%F(x)矩陣F1=zeros(n1,n1)。%F(x)導(dǎo)數(shù)矩陣 a=0。%記錄迭代次數(shù)EF=zeros(n1,n1)。%最后的節(jié)點(diǎn)電壓矩陣 while max(FX)=10^(5)for i=1:m %PQ節(jié)點(diǎn)FX(i)=P(i)。%△P FX(n+i1)=Q(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)。%△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)。%△Qend end for i=m+1:n1 %PV節(jié)點(diǎn) FX(i)=PP(im)。%△P FX(n+i1)=V(im)^2E(i)^2F(i)^2。%△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)。%△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))。F1(i,n+w1)=B(i,w)*E(i)G(i,w)*F(i)。F1(n+i1,w)=B(i,w)*E(i)G(i,w)*F(i)。F1(n+i1,n+w1)=G(i,w)*E(i)+B(i,w)*F(i)。elseF1(i,w)=G(i,i)*E(i)B(i,i)*F(i)。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。F1(n+i1,w)=B(i,i)*E(i)G(i,i)*F(i)。F1(n+i1,n+w1)=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)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。F1(n+i1,w)=F1(n+i1,w)+G(i,k)*F(k)+B(i,k)*E(k)。F1(n+i1,n+w1)=F1(n+i1,n+w1)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))。F1(i,n+w1)=B(i,w)*E(i)G(i,w)*F(i)。F1(n+i1,w)=0。F1(n+i1,n+w1)=0。elseF1(i,w)=G(i,i)*E(i)B(i,i)*F(i)。F1(i,n+w1)=B(i,i)*E(i)G(i,i)*F(i)。F1(n+i1,w)=2*E(i)。F1(n+i1,n+w1)=2*F(i)。for k=1:nF1(i,w)=F1(i,w)G(i,k)*E(k)+B(i,k)*F(k)。F1(i,n+w1)= F1(i,n+w1)G(i,k)*F(k)B(i,k)*E(k)。endendend end X=inv(F1)*(FX)。for i=1:n1E(i)=E(i)+X(i)。F(i)=F(i)+X(n+i1)。end a=a+1。fprintf(39。第%d次迭代后的節(jié)點(diǎn)電壓分別為:n39。,a)。disp(E+j*F)。fprintf(39。第%d次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:n39。,a)。disp(FX)。enddisp(39。收斂后的節(jié)點(diǎn)電壓用極坐標(biāo)表示為:39。)。EF=E+j*F。for i=1:n1 fprintf(39。%d號節(jié)點(diǎn)電壓的幅值為:39。,i)disp(abs(EF(i)))。fprintf(39。%d號節(jié)點(diǎn)電壓的相角度數(shù)為39。,i)disp(angle(EF(i))*180/pi)。end PPH=0。for i=1:nPPH=PPH+EF(n)*conj(Y(n,i))*conj(EF(i))。end fprintf(39。平衡節(jié)點(diǎn)的功率39。)。disp(PPH)。運(yùn)行結(jié)果：運(yùn)行結(jié)果復(fù)制如下：第1次迭代后的節(jié)點(diǎn)電壓分別為:第1次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:00第2次迭代后的節(jié)點(diǎn)電壓分別為:第2次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:第3次迭代后的節(jié)點(diǎn)電壓分別為:第3次迭代后功率偏差△P △Q電壓偏差△V的平方分別為:第4次迭代后的節(jié)點(diǎn)電壓分別為:第4次迭代后功率偏差△P △Q電壓偏差△V的平方分別為: *收斂后的節(jié)點(diǎn)電壓用極坐標(biāo)表示為: 1號節(jié)點(diǎn)電壓的幅值為:2號節(jié)點(diǎn)電壓的幅值為:3號節(jié)點(diǎn)電壓的幅值為:4號節(jié)點(diǎn)電壓的幅值為:平衡節(jié)點(diǎn)的功率 第三篇：電力系統(tǒng)分析潮流計算的計算機(jī)算法潮流計算的計算機(jī)算法實(shí)驗(yàn)報告姓名：學(xué)號：班級：一、實(shí)驗(yàn)?zāi)康恼莆粘绷饔嬎愕挠嬎銠C(jī)算法。熟悉MATLAB，并掌握MATLAB程序的基本調(diào)試方法。二、實(shí)驗(yàn)準(zhǔn)備根據(jù)課程內(nèi)容，熟悉MATLAB軟件的使用方法，自行學(xué)習(xí)MATLAB程序的基礎(chǔ)語法，并根據(jù)所學(xué)知識編寫潮流計算牛頓拉夫遜法(或PQ分解法)的計算程序，用相應(yīng)的算例在MATLAB上進(jìn)行計算、調(diào)試和驗(yàn)證。三、實(shí)驗(yàn)要求每人一組，在實(shí)驗(yàn)課時內(nèi)，用MATLAB調(diào)試和修改運(yùn)行程序，用算例計算輸出潮流結(jié)果。四、實(shí)驗(yàn)程序clear。%清空內(nèi)存n=input(39。請輸入節(jié)點(diǎn)數(shù):n=39。)。n1=input(39。請輸入支路數(shù):n1=39。)。isb=input(39。請輸入平衡節(jié)點(diǎn)號:isb=39。)。pr=input(39。請輸入誤差精度:pr=39。)。B1=input(39。請輸入支路參數(shù):B1=39。)。B2=input(39。請輸入節(jié)點(diǎn)參數(shù):B2=39。)。X=input(39。節(jié)點(diǎn)號和對地參數(shù):X=39。)。Y=zeros(n)。Times=1。%一：創(chuàng)建節(jié)點(diǎn)導(dǎo)納矩陣for i=1:n1if B1(i,6)==0%不含變壓器的支路p=B1(i,1)。q=B1(i,2)。Y(p,q)=Y(p,q)1/B1(i,3)。Y(q,p)=Y(p,q)。Y(p,p)=Y(p,p)+1/B1(i,3)+*B1(i,4)。Y(q,q)=Y(q,q)+1/B1(i,3)+*B1(i,4)。else%含有變壓器的支路p=B1(i,1)。q=B1(i,2)。Y(p,q)=Y(p,q)1/(B1(i,3)*B1(i,5))。Y(q,p)=Y(p,q)。Y(p,p)=Y(p,p)+1/B1(i,3)。Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3))。end end Y。%將OrgS、DetaS初始化OrgS=zeros(2*n2,1)。DetaS=zeros(2*n2,1)。%二：創(chuàng)建OrgS，用于存儲初始功率參數(shù)h=0。j=0。for i=1:n%對PQ節(jié)點(diǎn)的處理if i~=isbamp。B2(i,6)==2h=h+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)))。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)))。endend end %三：對PV節(jié)點(diǎn)的處理，注意這時不可再將h初始化為0 for i=1:nif i~=isbamp。B2(i,6)==3h=h+1。for j=1:nO