电力系统分析潮流计算大作业.docx

电力系统分析潮流计算大作业

第一篇：电力系统分析潮流计算大作业电力系统分析潮流计算大作业（源程序及实验报告）源程序如下：采用直角坐标系的牛顿-拉夫逊迭代functionchaoliujisuan()m=3;%m=PQ节点个数v=1;%v=PV节点个数P=[-0.8055-0.180];%P=PQ节点的P值Q=[-0.5320-0.120];%Q=PQ节点的Q值PP=[0.5];%PP=PV节点的P值V=[1.0];%V=PV节点的U值E=[1111.01.0]';%E=PQ,PV,Vθ节点e的初值F=[00000]';%F=PQ,PV,Vθ节点f的初值G=[6.3110-3.5587-2.752300;-3.55878.5587-500;-2.7523-57.752300;00000;00000];B=[-20.402211.38799.174300;11.3879-31.00937154.98890;9.1743-28.775704.9889;04.988905.24930;004.98890-5.2493];Y=G+j*B;X=[];%X=△Xn=m+v+1;%总的节点数FX=ones(2*n-2,1);%F(x)矩阵F1=zeros(n-1,n-1);%F(x)导数矩阵a=0;%记录迭代次数EF=zeros(n-1,n-1);%最后的节点电压矩阵whilemax(FX)>=10^(-5)fori=1:m%PQ节点FX(i)=P(i);%△PFX(n+i-1)=Q(i);%△Qforw=1:nFX(i)=FX(i)-E(i)*G(i,w)*E(w)+E(i)*B(i,w)*F(w)-F(i)*G(i,w)*F(w)-F(i)*B(i,w)*E(w);%△PFX(n+i-1)=FX(n+i-1)-F(i)*G(i,w)*E(w)+F(i)*B(i,w)*F(w)+E(i)*G(i,w)*F(w)+E(i)*B(i,w)*E(w);%△Qendendfori=m+1:n-1%PV节点FX(i)=PP(i-m);%△PFX(n+i-1)=V(i-m)^2-E(i)^2-F(i)^2;%△Qforw=1:nFX(i)=FX(i)-E(i)*G(i,w)*E(w)+E(i)*B(i,w)*F(w)-F(i)*G(i,w)*F(w)-F(i)*B(i,w)*E(w);%△Pendendfori=1:m%PQ节点forw=1:n-1ifi~=wF1(i,w)=-(G(i,w)*E(i)+B(i,w)*F(i));F1(i,n+w-1)=B(i,w)*E(i)-G(i,w)*F(i);F1(n+i-1,w)=B(i,w)*E(i)-G(i,w)*F(i);F1(n+i-1,n+w-1)=G(i,w)*E(i)+B(i,w)*F(i);elseF1(i,w)=-G(i,i)*E(i)-B(i,i)*F(i);F1(i,n+w-1)=B(i,i)*E(i)-G(i,i)*F(i);F1(n+i-1,w)=B(i,i)*E(i)-G(i,i)*F(i);F1(n+i-1,n+w-1)=G(i,i)*E(i)+B(i,i)*F(i);fork=1:nF1(i,w)=F1(i,w)-G(i,k)*E(k)+B(i,k)*F(k);F1(i,n+w-1)=F1(i,n+w-1)-G(i,k)*F(k)-B(i,k)*E(k);F1(n+i-1,w)=F1(n+i-1,w)+G(i,k)*F(k)+B(i,k)*E(k);F1(n+i-1,n+w-1)=F1(n+i-1,n+w-1)-G(i,k)*E(k)+B(i,k)*F(k);endendendendfori=m+1:n-1%PV节点forw=1:n-1ifi~=wF1(i,w)=-(G(i,w)*E(i)+B(i,w)*F(i));F1(i,n+w-1)=B(i,w)*E(i)-G(i,w)*F(i);F1(n+i-1,w)=0;F1(n+i-1,n+w-1)=0;elseF1(i,w)=-G(i,i)*E(i)-B(i,i)*F(i);F1(i,n+w-1)=B(i,i)*E(i)-G(i,i)*F(i);F1(n+i-1,w)=-2*E(i);F1(n+i-1,n+w-1)=-2*F(i);fork=1:nF1(i,w)=F1(i,w)-G(i,k)*E(k)+B(i,k)*F(k);F1(i,n+w-1)=F1(i,n+w-1)-G(i,k)*F(k)-B(i,k)*E(k);endendendendX=inv(F1)*(-FX);fori=1:n-1E(i)=E(i)+X(i);F(i)=F(i)+X(n+i-1);enda=a+1;fprintf('第%d次迭代后的节点电压分别为:n