资源描述：
《华北电力大学电力系统暂态上机程序设计完整版》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、clc;clear;%---------------输入已知条件----------%bus_Num1=3;%节点数bus_Num2=5;%包括发电机节点的节点数branch_Num1=3;%线路数branch_Num2=5;%包括发电机支路的支路数branch1_No1=[1,1,2];branch1_No2=[2,3,3];branch2_No1=[1,1,2,1,2];branch2_No2=[2,3,3,4,5];%每条支路始末端编号%输入支路各序阻抗，z1_branch表示支路正序阻抗，z2_branch表示支路负序阻抗，z0_branch表示支路零序阻抗z1_branch(1)=j
2、*0.1;z1_branch(2)=j*0.1;z1_branch(3)=j*0.1;z2_branch=z1_branch;z0_branch(1)=j*0.2;z0_branch(2)=j*0.2;z0_branch(3)=j*0.2;%---------------第一步：不考虑发电机节点计算节点导纳矩阵----------%%节点导纳矩阵，Y1表示不计发电机节点的正序网络节点导纳阵，Y2表示不计发电机节点的负序网络节点导纳阵，Y0表示不计发电机节点的零序网络节点导纳阵，Y1=zeros(bus_Num1);%全为0矩阵（3*3）Y1(1,1)=1/(j*0.15);Y1(2,2)=1/
3、(j*0.075);%请同学们求正序节点导纳矩阵form=1:branch_Num1Y1(branch1_No1(m),branch1_No1(m))=Y1(branch1_No1(m),branch1_No1(m))+1/z1_branch(m);Y1(branch1_No1(m),branch1_No2(m))=Y1(branch1_No1(m),branch1_No2(m))-1/z1_branch(m);Y1(branch1_No2(m),branch1_No1(m))=Y1(branch1_No2(m),branch1_No1(m))-1/z1_branch(m);Y1(branch
4、1_No2(m),branch1_No2(m))=Y1(branch1_No2(m),branch1_No2(m))+1/z1_branch(m);endY2=Y1;%负序等于正序Y0=zeros(bus_Num1);Y0(1,1)=1/(j*0.0525);Y0(2,2)=1/(j*0.0263);%请同学们求零序节点导纳矩阵form=1:branch_Num1Y0(branch1_No1(m),branch1_No1(m))=Y0(branch1_No1(m),branch1_No1(m))+1/z0_branch(m);Y0(branch1_No1(m),branch1_No2(m))=
5、Y0(branch1_No1(m),branch1_No2(m))-1/z0_branch(m);Y0(branch1_No2(m),branch1_No1(m))=Y0(branch1_No2(m),branch1_No1(m))-1/z0_branch(m);Y0(branch1_No2(m),branch1_No2(m))=Y0(branch1_No2(m),branch1_No2(m))+1/z0_branch(m);end%---------------第二步：考虑发电机节点计算节点导纳矩阵----------%%节点导纳矩阵，YY1表示计及发电机节点的正序网络节点导纳阵，YY2表示
6、计及发电机节点的负序网络节点导纳阵z1_branch(4)=j*0.05;z1_branch(5)=j*0.025;z2_branch=z1_branch;YY1=zeros(bus_Num1);YY1(4,4)=1/(j*0.1);YY1(5,5)=1/(j*0.05);%请同学们求正序节点导纳矩阵form=1:branch_Num2YY1(branch2_No1(m),branch2_No1(m))=YY1(branch2_No1(m),branch2_No1(m))+1/z1_branch(m);YY1(branch2_No1(m),branch2_No2(m))=YY1(branch2
7、_No1(m),branch2_No2(m))-1/z1_branch(m);YY1(branch2_No2(m),branch2_No1(m))=YY1(branch2_No2(m),branch2_No1(m))-1/z1_branch(m);YY1(branch2_No2(m),branch2_No2(m))=YY1(branch2_No2(m),branch2_No2(m))+1/z1_b
