【文章內(nèi)容簡介】 nPQ = 0。 for I = 1:nb, % nb為總節(jié)點(diǎn)數(shù) type= bus(I,6)。 if type == 3, %求平衡節(jié)點(diǎn)數(shù) nSW = nSW + 1。 SW(nSW,:)=bus(I,:)。 elseif type == 2, %求 PV節(jié)點(diǎn)數(shù) nPV = nPV +1。 PV(nPV,:)=bus(I,:)。 else %求 PQ節(jié)點(diǎn)數(shù) nPQ = nPQ + 1。 PQ(nPQ,:)=bus(I,:)。 end end bus=[PQ。PV。SW]。 %按 PQ PV SW順序排列的節(jié)點(diǎn)數(shù)據(jù)矩陣 newbus=[1:nb]39。 nodenum=[newbus bus(:,1)]。 %形成節(jié)點(diǎn)號對應(yīng)表 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 11 bus(:,1)=newbus。 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)。 break end end end end %******************建立節(jié)點(diǎn)導(dǎo)納矩陣 ***************** function Y = ybus(bus,line) global nl。 global nb。 global bus。 global line。 global Y。 global myf。 Y=zeros(nb,nb)。 %對導(dǎo)納矩陣賦初值 0 for k=1:nl I=line(k,1)。 %讀入線路參數(shù) J=line(k,2)。 Zt=line(k,3)+j*line(k,4)。 if J~=0 %接地支路分母為零 Yt=1/Zt。 %非接地支路為阻抗的倒數(shù) end Ym=line(k,5)+j*line(k,6)。 %求 Ym K=line(k,7)。 %求 K if (K==0)amp。(J~=0) % 普通線路 : K=0。 Y(I,I)=Y(I,I)+Yt+Ym。 Y(J,J)=Y(J,J)+Yt+Ym。 Y(I,J)=Y(I,J)Yt。 Y(J,I)=Y(I,J)。 end if (K==0)amp。(J==0) % 對地支路 : K=0,J=0,R=X=0。 Y(I,I)=Y(I,I)+Ym。 end 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 12 if K0 % 變壓器線路 :Zt和 Ym為折算到 i側(cè)的值 ,K在 j側(cè) Y(I,I)=Y(I,I)+Yt+Ym。 Y(J,J)=Y(J,J)+Yt/K/K。 Y(I,J)=Y(I,J)Yt/K。 Y(J,I)=Y(I,J)。 end if K0 % 變壓器線路 :Zt和 Ym為折算到 K所在側(cè)的值 ,K在 i側(cè) Y(I,I)=Y(I,I)+Yt+Ym。 Y(J,J)=Y(J,J)+K*K*Yt。 Y(I,J)=Y(I,J)+K*Yt。 Y(J,I)=Y(I,J)。 end end %***************牛頓拉夫遜法解非線性方程序組 *************** function NR global nb。 global nPQ。 global bus。 global line。 global Y。 global myf。 max1=100。 %最大迭代次數(shù) eps1=。 %收斂精度 eps2=。 fprintf(myf, 39。節(jié)點(diǎn)導(dǎo)納矩陣 Y\n39。)。 %輸出節(jié)點(diǎn)導(dǎo)納矩陣 for I=1:nb for J=1:nb fprintf(myf, 39。%10f+j*(%10f) 39。, real(Y(I,J)),imag(Y(I,J)))。 end fprintf(myf, 39。\n39。)。 end for i=1:max1 angl(:,1)=bus(1:nb1,3)。 u(:,1)=bus(1:nPQ,2)。 x=[angl。u]。 %從 bus矩陣中得到 PQ PV節(jié)點(diǎn)的相角、 PQ節(jié)點(diǎn)的電壓形成的列向量作為 x Jac=form_jac(bus,Y)。 %得到雅可比矩陣 del=dPQ(Y,bus)。 %得到功率偏差列向量 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 13 dx=Jac\del。 %求得 dx fprintf(myf,39。第 %d次迭代結(jié)果 \n39。,i)。 fprintf(myf,39。第 %d次迭代的雅比矩陣 J\n39。,i)。 for I=1:nb+nPQ1 for J=1:nb+nPQ1 fprintf(myf, 39。%10f 39。, Jac(I,J))。 %輸出每次迭代的雅可比矩陣 end fprintf(myf, 39。\n39。)。 end fprintf(myf,39。第 %d次迭代的功率偏差 dP和 dQ\n39。,i)。 for I=1:nb+nPQ1 fprintf(myf, 39。%10e 39。, del(I,1))。 %輸出每次迭代的功率偏差 fprintf(myf, 39。\n39。)。 end fprintf(myf,39。 第 %d 次 迭 代 的 節(jié) 點(diǎn) 相 角 和 電 壓 的 偏 差dx\n39。,i)。 for I=1:nb+nPQ1 fprintf(myf, 39。%10e 39。, dx(I,1))。 %輸出每次迭代的節(jié)點(diǎn)相角和電壓的偏差 fprintf(myf, 39。\n39。)。 end for I=nb:nb+nPQ1 dx(I,1)=dx(I,1)*x(I,1)。 %由于求得的是 delU/U，故還需作相應(yīng)處理 end x=xdx。 %修正 fprintf(myf,39。第 %d次迭代的節(jié)點(diǎn)相角 delta和電壓 U\n39。,i)。 angl=x(1:nb1,1)。 %將合一起的 x，拆成相角、電壓兩個(gè)列向量 u=x(nb:nb+nPQ1,1)。 for I=1:nb1 %輸出修正的 bus矩陣中 PQ、 PV結(jié)點(diǎn)的相角 bus(I,3)=angl(I,1)。 fprintf(myf, 39。ang%d %10f \n39。, I,angl(I,1))。 end for I=1:nPQ %輸出修正的 bus矩陣中 PQ結(jié)點(diǎn)的電壓 bus(I,2)=u(I,1)。 fprintf(myf, 39。U%d %10f \n39。, I,u(I,1))。 end 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 14 if (max(abs(dx))eps1)amp。(max(abs(dPQ(Y,bus)))eps2) %判斷是否達(dá)到所指定的精度要求 break end end if i==max1 %超出最大迭代次數(shù)。出錯(cuò)指示 error(39。超過最大迭代次數(shù)，不收斂停機(jī)！ 39。)。 end %*******************計(jì)算節(jié)點(diǎn)注入功率 ************ function PQ global nb。 global nPQ。 global bus。 global line。 global Y。 for I=nPQ+1:nb %求 PV結(jié)點(diǎn)的 Q，平衡節(jié)點(diǎn)的 P,Q if bus(I,6)==3 %對于平衡節(jié)點(diǎn)，求其 P sum=0。 for J=1:nb ang=bus(I,3)bus(J,3)。 A=real(Y(I,J))*cos(ang)+imag(Y(I,J))*sin(ang)。 sum=sum+bus(I,2)*bus(J,2)*(A)。 end bus(I,4)=sum。 %求取的平衡節(jié)點(diǎn)的 P存回 bus矩陣 end sum=0。 for J=1:nb ang=bus(I,3)bus(J,3)。 B=real(Y(I,J))*sin(ang)imag(Y(I,J))*cos(ang)。 sum=sum+bus(I,2)*bus(J,2)*B。 end bus(I,5)=sum。 %求取的 PV、平衡節(jié)點(diǎn)的 Q存回 bus矩陣 end %***************計(jì)算線路功率及損耗 ****************** function flow global nl。 global lPQ。 global bus。 貴州大學(xué)電氣工程學(xué)院《電力系統(tǒng)分析》課程設(shè)計(jì) 15 global line。 for k=1:nl I=line(k,1)。 %讀入線路參數(shù) J=line(k,2)。 lPQ(k,1)=I。 %用 lPQ矩陣第一、二列保存線路