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

《电力系统暂态上机节点导纳矩阵》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

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)

2、=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)=Z

3、Z0(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;

4、Y2=eye(N1)*0;Y0=eye(N1)*0;%矩阵赋初值为0forp=1:N1forq=1:N1if(p==q)%求自导纳fora=1:N1ifZZ1(p,a)~=0Y1(p,q)=Y1(p,q)+1/ZZ1(p,a);endifZZ2(p,a)~=0Y2(p,q)=Y2(p,q)+1/ZZ2(p,a);endifZZ0(p,a)~=0Y0(p,q)=Y0(p,q)+1/ZZ0(p,a);endendelseifZZ1(p,q)~=0Y1(p,q)=-1/ZZ1(p,q);endifZZ2(p,q)~=0Y2(p,q)=-1/

5、ZZ2(p,q);endifZZ0(p,q)~=0Y0(p,q)=-1/ZZ0(p,q);endendendendYY1=eye(N2)*0;YY2=eye(N2)*0;%矩阵赋初值为0ZZ1(1,1)=0;ZZ2(1,1)=0;ZZ1(2,2)=0;ZZ2(2,2)=0;forp=1:N2forq=1:N2if(p==q)%求自导纳fora=1:N2ifZZ1(p,a)~=0YY1(p,q)=YY1(p,q)+1/ZZ1(p,a);endifZZ2(p,a)~=0YY2(p,q)=YY2(p,q)+1/ZZ2(p,a);endend

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

7、Y2=[];%正序节点阻抗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