【正文】 TLAB運(yùn)行結(jié)果，文件內(nèi)容如下：% (bus)(volt) (ang) (p) (q) (bus type) bus=[ 1 1。, 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)))。for I=1:nl, fprintf(myf, 39。節(jié)點(diǎn)I 節(jié)點(diǎn)J 線路功率S(I,J) 線路功率S(J,I) 線路損耗dS(I,J)\n39。)。endfprintf(myf, 39。 %2d %10f %10f %10f+j%10f \n39。)。 fprintf(myf, 39。 節(jié)點(diǎn)計(jì)算結(jié)果：\n39。)。 fprintf(myf, 39。global lPQ。global nl。 for I=1:nl for J=1:2 for k=1:nb if lPQ(I,J)==nodenum(k,1) lPQ(I,J)=nodenum(k,2)。 tem(I,1)=I。for I=1:nb for J=1:nb if nodenum(J,2)==I break。global nodenum。global line。global nb。 lPQ(k,3:5)=[Sij Sji delS]。 Sji=Uj*conj(Iji)。 Iji=(UjUi)*Yt*K+Uj*Yt*K*(K1)。 end if K0 K=K。 end if K0 Iij=(UiUj)*Yt/K+Ui*(Ym+Yt*(K1)/K)。(J==0) Iij=Ui*Ym。 Iji=Uj*(Yt+Ym)Ui*Yt。 end if (K==0)amp。 Ui=bus(I,2)*(cos(bus(I,3))+j*sin(bus(I,3)))。 end Ym=line(k,5)+j*line(k,6)。if J~=0 Zt=line(k,3)+j*line(k,4)。 lPQ(k,1)=I。 for k=1:nlI=line(k,1)。global bus。 end計(jì)算線路功率及損耗：global nl。 sum=sum+bus(I,2)*bus(J,2)*B。 for J=1:nb ang=bus(I,3)bus(J,3)。 end bus(I,4)=sum。 A=real(Y(I,J))*cos(ang)+imag(Y(I,J))*sin(ang)。 for I=nPQ+1:nb if bus(I,6)==3 sum=0。global line。global nPQ。)。(max(abs(dPQ(Y,bus)))eps2) break endendif i==max1 error(39。, I,u(I,1))。 fprintf(myf, 39。, I,angl(I,1))。 fprintf(myf, 39。 u=x(nb:nb+nPQ1,1)。,i)。 fprintf(myf,39。 end for I=nb:nb+nPQ1 dx(I,1)=dx(I,1)*x(I,1)。\n39。, dx(I,1))。 for I=1:nb+nPQ1 fprintf(myf, 39。第%d次迭代的節(jié)點(diǎn)相角和電壓的偏差dx\n39。)。 fprintf(myf, 39。%10e 39。,i)。 end fprintf(myf,39。\n39。, Jac(I,J))。 for I=1:nb+nPQ1 for J=1:nb+nPQ1 fprintf(myf, 39。第%d次迭代的雅比矩陣J\n39。,i)。 fprintf(myf,39。 del=dPQ(Y,bus)。u]。 u(:,1)=bus(1:nPQ,2)。)。 end fprintf(myf, 39。%10f+j*(%10f) 39。)。 fprintf(myf, 39。 eps1=。global myf。global line。global nPQ。 Y(J,I)=Y(I,J)。 Y(J,J)=Y(J,J)+K*K*Yt。 Y(J,I)=Y(I,J)。 Y(J,J)=Y(J,J)+Yt/K/K。(J==0) Y(I,I)=Y(I,I)+Ym。 Y(J,I)=Y(I,J)。 Y(J,J)=Y(J,J)+Yt+Ym。 if (K==0)amp。 end%Ym=line(k,5)+j*line(k,6)。 Zt=line(k,3)+j*line(k,4)。 for k=1:nlI=line(k,1)。global myf。global line。global nb。 for I=1:nl for J=1:2 for k=1:nb if line(I,J)==nodenum(k,2) line(I,J)=nodenum(k,1)。nodenum=[newbus bus(:,1)]。SW]。 endend bus=[PQ。else nPQ = nPQ + 1。elseif type == 2, nPV = nPV +1。if type == 3, nSW = nSW + 1。 nPQ = 0。 nSW = 0。global line。global nPQ。 節(jié)點(diǎn)重新編號(hào)：function change global nb。w39。)。 endoutfile=strcat(outfile,39。 break。for I=1:length(dfile) if dfile(I)==39。 outfile=39。 end [nl,ml]=size(line)。)else lfile =length(dfile)。 if pathname == 0 error(39。Select Data File39。*.m39。global myf。global line。 output。 flow。 NR。 change。 global myf。 global nodenum。 global line。