【文章內(nèi)容簡介】 39。 Total loss 39。)fprintf(39。%39。, real(SLT)), fprintf(39。%\n39。, imag(SLT))clear Ik In SL SLT Skn Snk% Power flow solution by NewtonRaphson methodns=0。 ng=0。 Vm=0。 delta=0。 yload=0。 deltad=0。nbus = length(busdata(:,1))。kb=[]。Vm=[]。 delta=[]。 Pd=[]。 Qd=[]。 Pg=[]。 Qg=[]。 Qmin=[]。 Qmax=[]。 Pk=[]。 P=[]。 Qk=[]。 Q=[]。 S=[]。 V=[]。for k=1:nbusn=busdata(k,1)。kb(n)=busdata(k,2)。 Vm(n)=busdata(k,3)。 delta(n)=busdata(k, 4)。Pd(n)=busdata(k,5)。 Qd(n)=busdata(k,6)。 Pg(n)=busdata(k,7)。 Qg(n) = busdata(k,8)。Qmin(n)=busdata(k, 9)。 Qmax(n)=busdata(k, 10)。Qsh(n)=busdata(k, 11)。 if Vm(n) = 0 Vm(n) = 。 V(n) = 1 + j*0。 else delta(n) = pi/180*delta(n)。 V(n) = Vm(n)*(cos(delta(n)) + j*sin(delta(n)))。 P(n)=(Pg(n)Pd(n))/basemva。 Q(n)=(Qg(n)Qd(n)+ Qsh(n))/basemva。 S(n) = P(n) + j*Q(n)。 endendfor k=1:nbusif kb(k) == 1, ns = ns+1。 else, endif kb(k) == 2 ng = ng+1。 else, endngs(k) = ng。nss(k) = ns。endYm=abs(Ybus)。 t = angle(Ybus)。m=2*nbusng2*ns。maxerror = 1。 converge=1。iter = 0。 mline=ones(nbr,1)。for k=1:nbr for m=k+1:nbr if((nl(k)==nl(m)) amp。 (nr(k)==nr(m)))。 mline(m)=2。 elseif ((nl(k)==nr(m)) amp。 (nr(k)==nl(m)))。 mline(m)=2。 else, end end end %雅可比矩陣clear A DC J DXwhile maxerror = accuracy amp。 iter = maxiter for ii=1:mfor k=1:m A(ii,k)=0。 %初始化雅可比矩陣end, enditer = iter+1。for n=1:nbusnn=nnss(n)。lm=nbus+nngs(n)nss(n)ns。J11=0。 J22=0。 J33=0。 J44=0。 for ii=1:nbr if mline(ii)==1 if nl(ii) == n | nr(ii) == n if nl(ii) == n , l = nr(ii)。 end if nr(ii) == n , l = nl(ii)。 end J11=J11+ Vm(n)*Vm(l)*Ym(n,l)*sin(t(n,l) delta(n) + delta(l))。 J33=J33+ Vm(n)*Vm(l)*Ym(n,l)*cos(t(n,l) delta(n) + delta(l))。 if kb(n)~=1 J22=J22+ Vm(l)*Ym(n,l)*cos(t(n,l) delta(n) + delta(l))。 J44=J44+ Vm(l)*Ym(n,l)*sin(t(n,l) delta(n) + delta(l))。 else, end if kb(n) ~= 1 amp。 kb(l) ~=1 lk = nbus+lngs(l)nss(l)ns。 ll = l nss(l)。 % J1的非對角元素 A(nn, ll) =Vm(n)*Vm(l)*Ym(n,l)*sin(t(n,l) delta(n) + delta(l))。 if kb(l) == 0 % J2的非對角元素 A(nn, lk) =Vm(n)*Ym(n,l)*cos(t(n,l) delta(n) + delta(l))。end if kb(n) == 0 % J3的非對角元素 A(lm, ll) =Vm(n)*Vm(l)*Ym(n,l)*cos(t(n,l) delta(n)+delta(l))。 end if kb(n) == 0 amp。 kb(l) == 0 % J4的非對角元素 A(lm, lk) =Vm(n)*Ym(n,l)*sin(t(n,l) delta(n) + delta(l))。end else end else , end else, end end Pk = Vm(n)^2*Ym(n,n)*cos(t(n,n))+J33。 Qk = Vm(n)^2*Ym(n,n)*sin(t(n,n))J11。 if kb(n) == 1 P(n)=Pk。 Q(n) = Qk。 end % Swing bus P if kb(n) == 2 Q(n)=Qk。 if Qmax(n) ~= 0 Qgc = Q(n)*basemva + Qd(n) Qsh(n)。 if iter = 7 if iter 2 if Qgc Qmin(n), Vm(n) = Vm(n) + 。 elseif Qgc Qmax(n), Vm(n) = Vm(n) 。end else, end else,end else,end end if kb(n) ~= 1 A(nn,nn) = J11。 % J1對角元素 DC(nn) = P(n)Pk。 end if kb(n) == 0 A(nn,lm) = 2*Vm(n)*Ym(n,n)*cos(t(n,n))+J22。 % J2對角元素 A(lm,nn)= J33。 % J3對角元素 A(lm,lm) =2*Vm(n)*Ym(n,n)*sin(t(n,n))J44。 % J4對角元素 DC(lm) = Q(n)Qk。 endendDX=A\DC39。for n=1:nbus nn=nnss(n)。 lm=nbus+nngs(n)nss(n)ns。 if kb(n) ~= 1 delta(n) = delta(n)+DX(nn)。 end if kb(n) == 0 Vm(n)=Vm(n)+DX(lm)。 end end maxerror=max(abs(DC))。 if iter == maxiter amp。 maxerror accuracy fprintf(39。\nWARNING: Iterative solution did not converged after 39。) fprintf(39。%g39。, iter), fprintf(39。 iterations.\n\n39。) fprintf(39。Press Enter to terminate the iterations and print the results \n39。) converge = 0。 pause, else, end end if converge ~= 1 tech= (39。 ITERATIVE SOLUTION DID NOT CONVERGE39。)。 else, tech=(39。 Power Flow Solution by NewtonRaphson Method39。)。end V = Vm.*cos(delta)+j*Vm.*sin(delta)。deltad=180/pi*delta。i=sqrt(1)。k=0。for n = 1:nbus if kb(n) == 1 k=k+1。 S(n)= P(n)+j*Q(n)。 Pg(n) = P(n)*basemva + Pd(n)。 Qg(n) = Q(n)*basemva + Qd(n) Qsh(n)。 Pgg(k)=Pg(n)。 Qgg(k)=Qg(n)。