【正文】 ile) if dfile(I)==39。 break。 end outfile=strcat(outfile,39。)。w39。 %打開輸出文件 %**********************給節(jié)點(diǎn)編號 ********************* function change global nb。 global nPQ。 global line。 nSW = 0。 nPQ = 0。 if type == 3, %求平衡節(jié)點(diǎn)數(shù) nSW = nSW + 1。 elseif type == 2, %求 PV節(jié)點(diǎn)數(shù) nPV = nPV +1。 else %求 PQ節(jié)點(diǎn)數(shù) nPQ = nPQ + 1。 end end bus=[PQ。SW]。 nodenum=[newbus bus(:,1)]。 for I=1:nl %根據(jù)節(jié)點(diǎn)號對應(yīng)表重新排列線路矩陣節(jié)點(diǎn)編號 for J=1:2 for k=1:nb if line(I,J)==nodenum(k,2) line(I,J)=nodenum(k,1)。 global nb。 global line。 global myf。 %對導(dǎo)納矩陣賦初值 0 for k=1:nl I=line(k,1)。 Zt=line(k,3)+j*line(k,4)。 %非接地支路為阻抗的倒數(shù) end Ym=line(k,5)+j*line(k,6)。 %求 K if (K==0)amp。 Y(I,I)=Y(I,I)+Yt+Ym。 Y(I,J)=Y(I,J)Yt。 end if (K==0)amp。 Y(I,I)=Y(I,I)+Ym。 Y(J,J)=Y(J,J)+Yt/K/K。 Y(J,I)=Y(I,J)。 Y(J,J)=Y(J,J)+K*K*Yt。 Y(J,I)=Y(I,J)。 global nPQ。 global line。 global myf。 %最大迭代次數(shù) eps1=。 fprintf(myf, 39。)。%10f+j*(%10f) 39。 end fprintf(myf, 39。)。 u(:,1)=bus(1:nPQ,2)。u]。 %得到雅可比矩陣 del=dPQ(Y,bus)。 %求得 dx fprintf(myf,39。,i)。第 %d次迭代的雅比矩陣 J\n39。 for I=1:nb+nPQ1 for J=1:nb+nPQ1 fprintf(myf, 39。, Jac(I,J))。\n39。 end fprintf(myf,39。,i)。%10e 39。 %輸出每次迭代的功率偏差 fprintf(myf, 39。)。 第 %d 次 迭 代 的 節(jié) 點(diǎn) 相 角 和 電 壓 的 偏 差dx\n39。 for I=1:nb+nPQ1 fprintf(myf, 39。, dx(I,1))。\n39。 end for I=nb:nb+nPQ1 dx(I,1)=dx(I,1)*x(I,1)。 %修正 fprintf(myf,39。,i)。 %將合一起的 x，拆成相角、電壓兩個(gè)列向量 u=x(nb:nb+nPQ1,1)。 fprintf(myf, 39。, I,angl(I,1))。 fprintf(myf, 39。, I,u(I,1))。(max(abs(dPQ(Y,bus)))eps2) %判斷是否達(dá)到所指定的精度要求 break end end if i==max1 %超出最大迭代次數(shù)。超過最大迭代次數(shù)，不收斂停機(jī)！ 39。 end %*******************計(jì)算節(jié)點(diǎn)注入功率 ************ function PQ 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)。 %求取的平衡節(jié)點(diǎn)的 P存回 bus矩陣 end sum=0。 B=real(Y(I,J))*sin(ang)imag(Y(I,J))*cos(ang)。 end bus(I,5)=sum。 global lPQ。 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 15 global line。 %讀入線路參數(shù) J=line(k,2)。 %用 lPQ矩陣第一、二列保存線路兩端節(jié)點(diǎn)編號 lPQ(k,2)=J。 Yt=1/Zt。 K=line(k,7)。 %用復(fù)數(shù)表示 Ui if J~=0 Uj=bus(J,2)*(cos(bus(J,3))+j*sin(bus(J,3)))。(J~=0) %普通線路非接地支路 。 Iji=Uj*(Yt+Ym)Ui*Yt。(J==0) %接地支路 Iij=Ui*Ym。 end if K0 %變壓器線路 : Zt和 Ym為折算到 i側(cè)的值 ,K在 j側(cè) Iij=(UiUj)*Yt/K+Ui*(Ym+Yt*(K1)/K)。 end if K0 %變壓器線路 : Zt和 Ym為折算到 K側(cè)的值 ,K在 i側(cè) K=K。 Iji=(UjUi)*Yt*K+Uj*Yt*K*(K1)。 %求取節(jié)點(diǎn)注入功率 Sji=Uj*conj(Iji)。 %求線路損耗 lPQ(k,3:5)=[Sij Sji delS]。 global nb。 global line。 global nodenum。 for I=1:nb %bus矩陣中的節(jié)點(diǎn) 號恢復(fù)原來的編號 for J=1:nb %查找其原編號 if nodenum(J,2)==I break。%bus矩陣中的節(jié)點(diǎn)號恢復(fù)原來的編號 tem(I,1)=I。 for I=1:nl %line矩陣中的節(jié)點(diǎn)號恢復(fù)原來的 編號 for J=1:2 for k=1:nb if lPQ(I,J)==nodenum(k,1) lPQ(I,J)=nodenum(k,2)。 global nl。 global lPQ。 fprintf(myf, 39。)。 節(jié)點(diǎn)計(jì)算結(jié)果： \n39。 %輸出節(jié)點(diǎn)計(jì)算結(jié)果到文件 fprintf(myf, 39。)。 %2d %10f %10f %10f+j%10f \n39。 end fprintf(myf, 39。)。節(jié)點(diǎn) I 節(jié)點(diǎn) J 線路功率 S(I,J) 線路功率 S(J,I) 線路損耗 dS(I,J)\n39。 for I=1:nl, fprintf(myf, 39。, lPQ(I,1),lPQ(I,2),real(lPQ(I,3)),imag(lPQ(I,3)),real(lPQ(I,4)),imag(lPQ(I,4)),real(lPQ(I,5)),imag(lPQ(I,5)))。 例題數(shù)據(jù)代碼 % (bus)(volt) (ang) (p) (q) (bus type) bus=[ 1 1。 3 1。 5 3]。 1 3 0。 5 3 。 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 18 第五章 運(yùn)行結(jié)果分析 運(yùn)用 MATLAB 運(yùn)行結(jié)果 節(jié)點(diǎn)導(dǎo)納矩陣 Y +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*() +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*() +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*() +j*( ) +j*( ) +j*( ) +j*( ) +j*( ) +j*() 第 1 次迭代結(jié)果 第 1 次迭代的雅比矩陣 J 第 1 次迭代的功率偏差 dP 和 dQ +000 +000 +000 +000 +000 +000 第 1 次迭代的節(jié)點(diǎn)相角和電壓的偏差 dx 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 19 第 1 次迭代的節(jié)點(diǎn)相角 delta 和電壓 U ang1 ang2 ang3 ang4 U1 U2 U3 第 2 次迭代結(jié)果 第 2 次迭代的雅比矩陣 J