电力系统暂态上机节点导纳矩阵.doc

clc;%清空屏幕显示
clear;%清空所有变量

N1=3; %输入网络的节点数
N2=5; %输入包括所有发电机节点的网络的节点数
k=3; %输入短路点的节点号

ZZ1=eye(N2)*0;
ZZ2=eye(N2)*0;
ZZ0=eye(N1)*0;

%输入支路各序阻抗，ZZ1表示支路正序阻抗，ZZ2表示支路负序阻抗，ZZ0表示支路零序阻抗
ZZ1(1,2)=j*0.1;
ZZ1(1,3)=j*0.1;
ZZ1(2,3)=j*0.1;
ZZ1(2,1)=ZZ1(1,2);
ZZ1(3,1)=ZZ1(1,3);
ZZ1(3,2)=ZZ1(2,3);
ZZ1(1,1)=j*0.15;
ZZ1(2,2)=j*0.075;
ZZ1(3,3)=0;

ZZ2(1,2)=j*0.1;
ZZ2(2,3)=j*0.1;
ZZ2(1,3)=j*0.1;
ZZ2(2,1)=ZZ2(1,2);
ZZ2(3,1)=ZZ2(1,3);
ZZ2(3,2)=ZZ2(2,3);
ZZ2(1,1)=j*0.15;
ZZ2(2,2)=j*0.075;
ZZ2(3,3)=0;

ZZ0(1,2)=j*0.2;
ZZ0(2,3)=j*0.2;
ZZ0(1,3)=j*0.2;
ZZ0(2,1)=ZZ0(1,2);
ZZ0(3,1)=ZZ0(1,3);
ZZ0(3,2)=ZZ0(2,3);
ZZ0(1,1)=j*0.0525;
ZZ0(2,2)=j*0.0263;
ZZ0(3,3)=0;

ZZ1(4,1)=j*0.05;
ZZ1(1,4)=ZZ1(4,1);
ZZ1(2,5)=j*0.025;
ZZ1(5,2)=ZZ1(2,5);
ZZ1(4,4)=j*0.10;
ZZ1(5,5)=j*0.05;

ZZ2(4,1)=j*0.05;
ZZ2(1,4)=ZZ2(4,1);
ZZ2(2,5)=j*0.025;
ZZ2(5,2)=ZZ2(2,5);
ZZ2(4,4)=j*0.10;
ZZ2(5,5)=j*0.05;

Y1=eye(N1)*0;
Y2=eye(N1)*0;
Y0=eye(N1)*0;
%矩阵赋初值为0
for p=1:N1
for q=1:N1
if(p==q) %求自导纳
for a=1:N1
if ZZ1(p,a)~=0
Y1(p,q)=Y1(p,q)+1/ZZ1(p,a);
end
if ZZ2(p,a)~=0
Y2(p,q)=Y2(p,q)+1/ZZ2(p,a);
end
if ZZ0(p,a)~=0
Y0(p,q)=Y0(p,q)+1/ZZ0(p,a);
end
end
else
if ZZ1(p,q)~=0
Y1(p,q)=-1/ZZ1(p,q);
end
if ZZ2(p,q)~=0
Y2(p,q)=-1/ZZ2(p,q);
end
if ZZ0(p,q)~=0
Y0(p,q)=-1/ZZ0(p,q);
end
end
end
end

YY1=eye(N2)*0;
YY2=eye(N2)*0;
%矩阵赋初值为0
ZZ1(1,1)=0;
ZZ2(1,1)=0;
ZZ1(2,2)=0;
ZZ2(2,2)=0;

for p=1:N2
for q=1:N2
if(p==q) %求自导纳
for a=1:N2
if ZZ1(p,a)~=0
YY1(p,q)=YY1(p,q)+1/ZZ1(p,a);
end
if ZZ2(p,a)~=0
YY2(p,q)=YY2(p,q)+1/ZZ2(p,a);
end
end
else
if ZZ1(p,q)~=0
YY1(p,q)=-1/ZZ1(p,q);
end
if ZZ2(p,q)~=0
YY2(p,q)=-1/ZZ2(p,q);
end
end
end
end
%计算各序电流的故障分量
%求节点导纳矩阵，Y1表示不计发电机节点的正序网络节点导纳阵，Y2表示不计发电机节点的负序网络节点导纳阵，Y0表示不计发电机节点的零序网络节点导纳阵，
% Y1=[];
% Y2=[];
% Y0=[];
%节点导纳矩阵参考结果，YY1表示计及发电机节点的正序网络节点导纳阵，YY2表示计及发电机节点的负序网络节点导纳阵
%矩阵赋初值为0

% YY1=[];
% YY2=[];
%正序节点阻抗
Z1=eye(N1)/Y1;%eye(N1)是3阶的单位矩阵,求节点导纳矩阵Y1的逆阵,即节点阻抗矩阵Z%Z1=inv(Y1)与上式等效%负序节点阻抗
Z2=eye(N1)/Y2;%零序节点阻抗
Z0=eye(N1)/Y0;

Z11=eye(N2)/YY1;
Z22=eye(N2)/YY2;
Y1
Y2
Y0
YY1
YY2