《电力系统暂态上机节点导纳矩阵.doc》由会员分享，可在线阅读，更多相关《电力系统暂态上机节点导纳矩阵.doc（4页珍藏版）》请在金锄头文库上搜索。
1、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
2、;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
3、; 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;%矩阵赋初值为0for p=1:N1 for q=1:N1 if(p=q) %求自导纳 for a=1:N1 i
4、f 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
5、)*0; %矩阵赋初值为0ZZ1(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
6、%计算各序电流的故障分量%求节点导纳矩阵，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;Y1Y2Y0YY1YY2
