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

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