资源描述：
《电力系统潮流计算程序设计》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、电力系统潮流计算程序设计姓名：韦应顺学号：2011021052电力工程学院牛顿—拉夫逊潮流计算方法具有能够将非线性方程线性化的特点，而使用MATLAB语言是由于MATLAB语言的数学逻辑强，易编译。1.MATLAB程序【1】【2】Functiontisco%这是一个电力系统潮流计算的程序n=input（‘请输入节点数：n=’）；m=input(‘请输入支路数：m=’);ph=input(‘请输入平衡母线的节点号：ph=’)；B1=input（‘请输入支路信号：B1=’）;%它以矩阵形式存贮支路的情况，每行存贮一条支路%第一列存贮支路的一个端点%第二列存贮支路的另一个端点%
2、第三列存贮支路阻抗%第四列存贮支路的对地导纳%第五列存贮变压器的变比，注意支路为1%第六列存贮支路的序号B2=input(‘请输入节点信息：B2=’)；%第一列为电源侧的功率%第二列为负荷侧的功率%第三列为该点的电压值%第四列为该点的类型：1为PQ，2为PV节点，3为平衡节点A=input（‘请输入节点号及对地阻抗：A=’）；ip=input（‘请输入修正值：ip=’）；%ip为修正值Y=zeros（n）；e=zeros（1，n）；f=zeros（1，n）；no=2*ph=1；fori=1：nifA（i2）=0p=A(i1)；Y(pp)=1./A(i2);EndEndFo
3、ri=1:mp=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；EndG=real(Y);B=imag(Y);fori=1：ne(i)=real(B2(i3));f(i)=imag(B2(i3));S(i)=B2(i1)-B2(i2);V(i)=B2(i3);endP=real(S);Q=imag(S);[C,D,DF]=xxf(G,B,e,f,P,Q,n
4、,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;whilemax(abs(De))>ip&max(abs(Dfi)>ipt=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);endv=e+f*j;fori=1:nhh(i)=conj(Y(ph,i)*v(i));endS(p
5、h)=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))];fori=1:ma(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)*
6、v(B1(i2))^2/2;endresult2=[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);typeresultmfunction[C,D,Df]=xxf(G,B,e,f,P,Q,n,B2,ph,V,no)%该子程序是用来求取Dffori=1:nIfi=phC(i)=0;D(i)=0;Forj=i:nC(i)=C(i)+G(i,j)*e(j)-B(i,j)*f(j);D(i)=D(i)+G(i
7、,j)*f(j)+B(i,j)*e(j);endP1=C(i)*e(i)+D(i)*f(i);Q1=C(i)*f(i)-D(i)*e(i);V1=e(i)^2+f(i)^2;IfB2(i4)=2p=2*i-1;Df(p)=P(i)-P1;p=p+1;elsep=2*i-1;Df(p)=P(i)-P1;p=p+1;Df(p)=Q(i)-Q1;endendendDf=Df’;Ifph=nDf(noJ=[];endfunction[De,Df]=hxf(J,Df
