5节点电力系统潮流计算程序

input('系统数据=');
Y=zeros(n);
%把发电机电抗加入自阻抗.....................................................
for i=1:n
if X(i,2)~=0;
p=X(i,1);
Y(p,p)=1./X(i,2);
end
end
%算节点阻抗矩阵............................................................
[l m]=size(A);
for i=1:l;
p=A(i,1);
q=A(i,2);
Y(p,q)=-1./A(i,3);
Y(q,p)=Y(p,q);
Y(p,p)=Y(p,p)-Y(p,q);
Y(q,q)=Y(q,q)-Y(p,q);
end
%节点导纳矩阵...............................................................
disp('短路分量网络的导纳矩阵Y=');
disp(Y);
%节点阻抗矩阵...............................................................
Z=Y^(-1);
disp('短路分量网络的阻抗矩阵Z=');
disp(Z);
%故障点f的选定，及其对应故障分量网络的节点注入电流矩阵I.......................
f=3;
%故障分量网络节点电压方程，阻抗形式..........................................
Zf=zeros(n,1);
for i=1:n
Zf(i)=Z(i,f);
end %把短路支路取一列阻抗...................................................
If3=Uf0(3)/Zf(f);
disp('短路点电流If3=');
disp(If3);
dU=Zf*If3;
disp('短路时电压故障分量dU=');
disp(dU);
%运用叠加原理求个节点残余电压...................................................
U=Uf0+dU;
disp('各节点残余电压U=');
disp(U);
%发电机电流................................................................
IG1=(-dU(1))/X(1,2);
disp('发电机G1的短路电流IG1=');
disp(IG1);
IG2=0.5*(-dU(2))/X(2,2);
disp('发电机G1的短路电流IG2=');
disp(IG2);
%线路电流...................................................................
Iij=zeros(n);
for i=1:n
for j=1:n;
Iij(i,j)=(dU(i)-dU(j))*Y(i,j);
end
end
disp('各支路电流矩阵Iij=');
disp(Iij);













%不计负荷，短路前所有电动势，电压均为1，电流为零
% n 为节点数
n=5;
% A为支路参数形成的矩阵;
% (1)为该支路首端节点号;
% (2)为该支路首端节点号;
% (3)为支路阻抗;
%(4)为支路变压器变比;
A=[2,3,0.145+0.581i,0;
1,5,0.104+0.581i,0;
3,5,0.031+0.248i,0;
3,4,0.082+0.427i,0.95;
1,2,0.163+0.754i,0];
% X 为发电机节点号阻抗参数矩阵
% X(i,1)为节点号;
% X(i,2)为发电机电抗;
X=[1,0.15i;
2,0.075i;
3,0;
4,0;
5,0];
%短路前各节点电压Uf0,不考虑负荷
Uf0=[1;1;1;1;1];