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

电力系统分析潮流计算大作业（源程序及实验报告）

源程序如下：

采用直角坐标系的牛顿-拉夫逊迭代

function chaoliujisuan()

m=3; %m=PQ节点个数

v=1;%v=PV节点个数

P=[-0.8055 -0.18 0]; %P=PQ节点的P值

Q=[-0.5320 -0.12 0]; %Q=PQ节点的Q值

PP=[0.5];%PP=PV节点的P值

V=[1.0];%V=PV节点的U值

E=[1 1 1 1.0 1.0]'; %E=PQ,PV,Vθ节点e的初值

F=[0 0 0 0 0]'; %F=PQ,PV,Vθ节点f的初值

G=[

6.3110 -3.5587 -2.7523 0 0;

-3.5587 8.5587 -5 0 0;

-2.7523 -5 7.7523 0 0;

0 0 0 0 0;

0 0 0 0 0

];

B=[

-20.4022 11.3879 9.1743 0 0;

11.3879 -31.00937 15 4.9889 0;

9.1743 15 -28.7757 0 4.9889;

0 4.9889 0 5.2493 0;

0 0 4.9889 0 -5.2493

];

Y=G+j*B;

X=[]; %X=△X

n=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);%最后的节点电压矩阵

while max(FX)>=10^(-5)

for i=1:m %PQ节点

FX(i)=P(i);%△P

FX(n+i-1)=Q(i);%△Q

for w=1:n

FX(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); %△P

FX(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); %△Q

end

end

for i=m+1:n-1 %PV节点

FX(i)=PP(i-m);%△P

FX(n+i-1)=V(i-m)^2-E(i)^2-F(i)^2; %△Q

for w=1:n

FX(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);%△P end

end

for i=1:m %PQ节点

for w=1:n-1

if i~=w

F1(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);

else

F1(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);

for k=1:n

F1(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);

end

end

end

end

for i=m+1:n-1 %PV节点

for w=1:n-1

if i~=w

F1(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;

else

F1(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);

for k=1:n

F1(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);

end

end

end

end

X=inv(F1)*(-FX) ;

for i=1:n-1

E(i)=E(i)+X(i);

F(i)=F(i)+X(n+i-1);

end

a=a+1;

fprintf('第%d次迭代后的节点电压分别为:\n',a);

disp(E+j*F);

fprintf('第%d次迭代后功率偏差△P △Q电压偏差△V的平方分别为:\n',a); disp(FX);

end

disp('收敛后的节点电压用极坐标表示为:');

EF=E+j*F;

for i=1:n-1

fprintf('%d号节点电压的幅值为:',i)

disp(abs(EF(i)));

fprintf('%d号节点电压的相角度数为',i)

disp(angle(EF(i))*180/pi);

end

PPH=0;

for i=1:n

PPH=PPH+EF(n)*conj(Y(n,i))*conj(EF(i));

end

fprintf('平衡节点的功率');

disp(PPH);

运行结果：