【正文】 PQ。 ret。 打開bus line的數(shù)據(jù)文件的子程序：function openfile global nb nl。global bus。 [dfile,pathname]=uigetfile(39。,39。)。 you must select a valid data file39。 % strip off .m eval(dfile(1:lfile2))。 [nb,mb]=size(bus)。39。.39。 end outfile=strcat(outfile,dfile(I))。.txt39。 myf=fopen(outfile,39。)。global nl。global bus。global nodenum。 nPV = 0。 for I = 1:nb, type= bus(I,6)。 SW(nSW,:)=bus(I,:)。 PV(nPV,:)=bus(I,:)。 PQ(nPQ,:)=bus(I,:)。PV。 newbus=[1:nb]39。 bus(:,1)=newbus。 break end end endend 建立節(jié)點(diǎn)導(dǎo)納矩陣：function Y = ybus(bus,line) global nl。global bus。global Y。 Y=zeros(nb,nb)。 J=line(k,2)。if J~=0 Yt=1/Zt。 K=line(k,7)。(J~=0) Y(I,I)=Y(I,I)+Yt+Ym。 Y(I,J)=Y(I,J)Yt。 endif (K==0)amp。 endif K0 Y(I,I)=Y(I,I)+Yt+Ym。 Y(I,J)=Y(I,J)Yt/K。 endif K0 Y(I,I)=Y(I,I)+Yt+Ym。 Y(I,J)=Y(I,J)+K*Yt。 endend 牛拉法解方程組：function NR global nb。global bus。global Y。 max1=100。 eps2=。節(jié)點(diǎn)導(dǎo)納矩陣Y\n39。 for I=1:nb for J=1:nb fprintf(myf, 39。, real(Y(I,J)),imag(Y(I,J)))。\n39。end for i=1:max1 angl(:,1)=bus(1:nb1,3)。 x=[angl。 Jac=form_jac(bus,Y)。 dx=Jac\del。第%d次迭代結(jié)果\n39。 fprintf(myf,39。,i)。%10f 39。 end fprintf(myf, 39。)。第%d次迭代的功率偏差dP和dQ\n39。 for I=1:nb+nPQ1 fprintf(myf, 39。, del(I,1))。\n39。 end fprintf(myf,39。,i)。%10e 39。 fprintf(myf, 39。)。 end x=xdx。第%d次迭代的節(jié)點(diǎn)相角delta和電壓U\n39。 angl=x(1:nb1,1)。 for I=1:nb1 bus(I,3)=angl(I,1)。ang%d %10f \n39。 end for I=1:nPQ % bus(I,2)=u(I,1)。U%d %10f \n39。 end if(max(abs(dx))eps1)amp。超過最大迭代次數(shù)，不收斂停機(jī)！39。end計(jì)算節(jié)點(diǎn)注入功率：global nb。global bus。global Y。 for J=1:nb ang=bus(I,3)bus(J,3)。 sum=sum+bus(I,2)*bus(J,2)*(A)。 end sum=0。 B=real(Y(I,J))*sin(ang)imag(Y(I,J))*cos(ang)。 endbus(I,5)=sum。global lPQ。global line。 J=line(k,2)。 lPQ(k,2)=J。 Yt=1/Zt。 K=line(k,7)。 if J~=0 Uj=bus(J,2)*(cos(bus(J,3))+j*sin(bus(J,3)))。(J~=0) Iij=Ui*(Yt+Ym)Uj*Yt。 end if (K==0)amp。 Iji=0。 Iji=(UjUi)*Yt/K+Uj*Yt*(1K)/K^2。 Iij=(UiUj)*Yt*K+Ui*(Ym+Yt*(1K))。 end Sij=Ui*conj(Iij)。delS=Sij+Sji。 end恢復(fù)原節(jié)點(diǎn)編號：global nl。global bus。global Y。global lPQ。 end end tem(I,:)=bus(nodenum(J,1),:)。endbus=tem。 break end end endend 輸出：function output global nb。global bus。global myf。\n牛頓－拉夫遜法潮流計(jì)算結(jié)果\n39。fprintf(myf, 39。)。節(jié)點(diǎn) 節(jié)點(diǎn)電壓 節(jié)點(diǎn)相角（角度） 節(jié)點(diǎn)注入功率\n39。for I=1:nb, fprintf(myf, 39。, bus(I,1),bus(I,2),bus(I,3)*180/pi,bus(I,4),bus(I,5))。\n 線路計(jì)算結(jié)果：\n39。 fprintf(myf, 39。)。 %2d %2d %10f+j%10f %10f+j%10f %10f+j%10f\n39。endfclose(myf)。 2 1。 4 2。 % b1 b2 ( R ) ( X ) (G) ( B ) ( K )line = [ 1 2 0。 4 5 0。 5 1 ]。 各節(jié)點(diǎn)，支路的運(yùn)行狀態(tài) 非正常運(yùn)行結(jié)果分析 若某發(fā)電站發(fā)電量減半時對系統(tǒng)斷電影響假設(shè)當(dāng)發(fā)電機(jī)1（節(jié)點(diǎn)1上的發(fā)電機(jī)，屬于PQ節(jié)點(diǎn)）發(fā)電量減半時，即有功和無功都減半，此時電路的潮流分布圖如下：各節(jié)點(diǎn)，支路的運(yùn)行狀態(tài)：節(jié)點(diǎn)支路由上面三幅圖結(jié)合正常運(yùn)行時的數(shù)據(jù)可以看出，當(dāng)發(fā)電機(jī)1的發(fā)電量減半時，PQ節(jié)點(diǎn)的電壓會降低，而PV，平衡節(jié)點(diǎn)的電壓基本保持不變。但是由于系統(tǒng)的波動不是很大，尚在系統(tǒng)的調(diào)節(jié)能力范圍內(nèi)，故系統(tǒng)不會出現(xiàn)電壓崩潰，不過2節(jié)點(diǎn)的電壓已經(jīng)很低了，如果短時間內(nèi)不采取相應(yīng)措施升高電壓將會對用戶造成影響。第七章總結(jié)通過本次課程設(shè)計(jì)讓我有復(fù)習(xí)了一次潮流計(jì)算的相關(guān)知識，更加家清晰了什么是潮流計(jì)算以及潮流計(jì)算的在電力系統(tǒng)的重要性。電力系統(tǒng)分析的潮流計(jì)算是電力系統(tǒng)分析的一個重要的部分。同時潮流分布也是電力系統(tǒng)規(guī)劃設(shè)計(jì)的一項(xiàng)基礎(chǔ)工作