资源描述:
《中国石油大学电力系统分析matlab编程.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、电力系统分析大作业m=2;%m=PQ节点个数v=1;%v=PV节点个数P=[-0.3-0.55];%P=PQ节点的P值Q=[-0.18-0.13];%Q=PQ节点的Q值PP=[0.5];%PP=PV节点的P值V=[1.1];%V=PV节点的U值E=[111.11.05]';%E=PQ,PV,Vθ节点e的初值F=[0000]';%F=PQ,PV,Vθ节点f的初值G=[1.-0.0-0.;-0.1.0-0.;0000;-0.-0.00.];B=[-8.248342.3.1.;2.-4.02.;3.0-3.0;1.2.0-4.26142];Y=G+j*B;X=[];%X=△Xn=m+v+1
2、;%总的节点数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(2*i-1)=P(i);%△PFX(2*i)=Q(i);%△Qforw=1:nFX(2*i-1)=FX(2*i-1)-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(2*i)=FX(2*i)-F(i)*G(i,w)*E(w
3、)+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(2*i-1)=PP(i-m);%△PFX(2*i)=V(i-m)^2-E(i)^2-F(i)^2;%△Qforw=1:nFX(2*i-1)=FX(2*i-1)-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(2*i-1,2*w-1)=-(G(i,w)*E(
4、i)+B(i,w)*F(i));F1(2*i-1,2*w)=B(i,w)*E(i)-G(i,w)*F(i);F1(2*i,2*w-1)=B(i,w)*E(i)-G(i,w)*F(i);F1(2*i,2*w)=G(i,w)*E(i)+B(i,w)*F(i);elseF1(2*i-1,2*w-1)=-G(i,i)*E(i)-B(i,i)*F(i);F1(2*i-1,2*w)=B(i,i)*E(i)-G(i,i)*F(i);F1(2*i,2*w-1)=B(i,i)*E(i)-G(i,i)*F(i);F1(2*i,2*w)=G(i,i)*E(i)+B(i,i)*F(i);fork=1:nF1
5、(2*i-1,2*w-1)=F1(2*i-1,2*w-1)-G(i,k)*E(k)+B(i,k)*F(k);F1(2*i-1,2*w)=F1(2*i-1,2*w)-G(i,k)*F(k)-B(i,k)*E(k);F1(2*i,2*w-1)=F1(2*i,2*w-1)+G(i,k)*F(k)+B(i,k)*E(k);F1(2*i,2*w)=F1(2*i,2*w)-G(i,k)*E(k)+B(i,k)*F(k);endendendendfori=m+1:n-1%PV节点forw=1:n-1ifi~=wF1(2*i-1,2*w-1)=-(G(i,w)*E(i)+B(i,w)*F(i));F
6、1(2*i-1,2*w)=B(i,w)*E(i)-G(i,w)*F(i);F1(2*i,2*w-1)=0;F1(2*i,2*w)=0;elseF1(2*i-1,2*w-1)=-G(i,i)*E(i)-B(i,i)*F(i);F1(2*i-1,2*w)=B(i,i)*E(i)-G(i,i)*F(i);F1(2*i,2*w-1)=-2*E(i);F1(2*i,2*w)=-2*F(i);fork=1:nF1(2*i-1,2*w-1)=F1(2*i-1,2*w-1)-G(i,k)*E(k)+B(i,k)*F(k);F1(2*i-1,2*w)=F1(2*i-1,2*w)-G(i,k)*F(k)
7、-B(i,k)*E(k);endendendendX=inv(F1)*(-FX);fori=1:n-1E(i)=E(i)+X(2*i-1);F(i)=F(i)+X(2*i);enda=a+1;fprintf('第%d次迭代后的节点电压分别为:',a);disp(E+j*F);fprintf('第%d次迭代后功率偏差△P△Q电压偏差△V的平方分别为:',a);disp(FX);enddisp('收敛后的节点电压用极坐标表示为:');EF=E+j*F;