(电力系统潮流计算程序设计.docVIP

电力系统潮流计算程序设计 姓名：韦应顺 学号：2011021052 电力工程学院 牛顿—拉夫逊潮流计算方法具有能够将非线性方程线性化的特点，而使用MATLAB语言是由于MATLAB语言的数学逻辑强，易编译。 MATLAB程序【1】【2】 Function tisco %这是一个电力系统潮流计算的程序 n=input（‘\n请输入节点数：n=’）； m=input(‘请输入支路数：m=’); ph=input(‘\n请输入平衡母线的节点号：ph=’)； B1=input（‘\n请输入支路信号：B1=’）; %它以矩阵形式存贮支路的情况，每行存贮一条支路 %第一列存贮支路的一个端点 %第二列存贮支路的另一个端点 %第三列存贮支路阻抗 %第四列存贮支路的对地导纳 %第五列存贮变压器的变比，注意支路为1 %第六列存贮支路的序号 B2=input(‘\n请输入节点信息：B2=’)； %第一列为电源侧的功率 %第二列为负荷侧的功率 %第三列为该点的电压值 %第四列为该点的类型：1为PQ，2为PV节点， 3为平衡节点 A=input（‘\n请输入节点号及对地阻抗：A=’）； ip=input（‘\n请输入修正值：ip=’）； %ip为修正值 Y=zeros（n）； e=zeros（1，n）； f=zeros（1，n）； no=2*ph=1； for i=1：n if A（i2）=0 p=A(i1)； Y(p p)=1./A(i2); End End For i=1:m p=B1(i1); p=B1(i2); Y(p,p)=Y(p,p)+1./(B1(i3)*B1(i5)^2+B1(i4)./2); Y(p,q)=Y(p,q)-1./(B1(i3)*B1(i5); Y(p,q)=Y(p,q); Y(q,q)=Y(q,q)+1./B1(i3)+B1(i4)/2； End G=real(Y); B=imag(Y); for i=1：n e(i)=real(B2(i3)); f(i)=imag(B2(i3)); S(i)=B2(i1)-B2(i2); V(i)=B2(i3); end P=real(S); Q=imag(S); [C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no); J=jacci(Y,G,B,P,Q,e,f,V,C,D,B2,n,ph,no); [De,Di]=hxf(J,D,F,ph,n,no); t=0; while max(abs(De))ipmax(abs(Dfi)ip t=t+1; e=e+De； f=f+Df; [C,D,DF]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)； J=jacci(Y,G,B,P,Q,e,f,V,C,D,B2,n,ph,no) ； [De,Df]=hxf(J,Df,ph,n,no); end v=e+f*j; for i=1:n hh(i)=conj(Y(ph,i)*v(i)); end S(ph)=sum(hh)*v(ph); B2(ph,1)=S(ph); V=abs(v); jd=angle(v)*180/p; resulte1=[A(:,1),real(v),imag(v),V,jd,real(S’), imag(S’),real(B2(:1)),imag(B2(:1)),real(B2(:2)),imag(B2(:,2))]; for i=1:m a(i)=conj((v(B1(i1))/B1(i5)-v(B1(i2))/B1(i3)); b(i)=v(B1(i1))*a(i)-j*B1(i4)*v(B1(i))^2/2; c(i)=-v(B1(i2))*a(i)-j*B1(i4)*v(B1(i2))^2/2; end result2=[B1(:,6),B1(:,1),B1(:,2),real(b’),imag(b’),real(c’),imag(c’), real(b’+c’),imag(b’+c’)]; printcut(result1,S,b,c,result2); type resultm function [C,D,Df]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no) %该子程序是用来求取Df for i=1:n If i=ph C(i)=0; D(i)=0; For j=i:n C(i)=C(i)+G(i,j)*e(j)-B(i,j)*f(j); D(i)=D(i)+G(i,j)*f(j)+B(i,j)*e(j); end P1=C(i)*e(i)+D(i)*f(i); Q1=C(