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

电力系统分析潮流计算大作业（源程序及实验报告）源程序如下：采用直角坐标系的牛顿-拉夫逊迭代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.174315-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',a);disp(E+j*F);fprintf('第%d次迭代后功率偏差△P△Q电压偏差△V的平方分别为:\n',a);disp(FX);enddisp('收敛后的节点电压用极坐标表示为:');EF=E+j*F;fori=1:n-1fprintf('%d号节点电压的幅值为:',i)disp(abs(EF(i)));fprintf('%d号节点电压的相角度数为',i)disp(angle(EF(i))*180/pi);endPPH=0;fori=1:nPPH=PPH+EF(n)*conj(Y(n,i))*conj(EF(i));endfprintf('平衡节点的功率');disp(PPH);运行结果：运行结果复制如下：第1次迭代后的节点电压分别为:1.0034-0.1306i1.0285-0.1019i1.0339-0.0973i1.0000-0.0017i1.0000第1次迭代后功率偏差△P△Q电压偏差△V的平方分别为:-0.8055-0.180000.5000-0.37200.24740.38750第2次迭代后的节点电压分别为:0.9836-0.1290i1.0126-0.1038i1.0183-0.0995i1.0000-0.0035i1.0000第2次迭代后功率偏差△P△Q电压偏差△V的平方分别为:0.0512-0.0222-0.04030.0002-0.1012-0.0219-0.0099-0.0000第3次迭代后的节点电压分别为:0.9831-0.1290i1.0122-0.1038i1.0180-0.0996i1.0000-0.0035i1.0000第3次迭代后功率偏差△P△Q电压偏差△V的平方分别为:0.0008-0.0003-0.0005-0.0001-0.0021-0.0004-0.0003-0.0000第4次迭代后的节点电压分别为:0.9831-0.1290i1.0122-0.1038i1.0180-0.0996i1.0000-0.0035i1.0000第4次迭代后功率偏差△P△Q电压偏差△V的平方分别为:1.0e-005*0.0280-0.0083-0.0164-0.0121-0.1085-0.0199-0.0135-0.0005收敛后的节点电压用极坐标表示为:1号节点电压的幅值为:0.99161号节点电压的相角度数为-7.47482号节点电压的幅值为:1.01752号节点电压的相角度数为-5.85483号节点电压的幅值为:1.02293号节点电压的相角度数为-5.58644号节点电压的幅值为:1.00004号节点电压的相角度数为-0.2021平衡节点的功率0.4968-10.3280i