(精品word)华北电力大学电力系统暂态上机程序设计完整版.doc

(精品word)华北电力大学电力系统暂态上机程序设计完整版
(精品word)华北电力大学电力系统暂态上机程序设计完整版
(精品word)华北电力大学电力系统暂态上机程序设计完整版
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＊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/（j＊0。075）；
％请同学们求正序节点导纳矩阵
form=1:branch_Num1
Y1(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（branch1_No2（m)，branch1_No2（m）)=Y1（branch1_No2(m）,branch1_No2(m））+1/z1_branch（m)；
end

Y2=Y1;%负序等于正序

Y0=zeros(bus_Num1）；
Y0（1，1）=1/(j＊0.0525）;
Y0(2，2）=1/(j*0。0263)；
％请同学们求零序节点导纳矩阵
form=1:branch_Num1
Y0（branch1_No1（m）,branch1_No1(m))=Y0（branch1_No1（m),branch1_No1（m)）+1/z0_branch(m)；
Y0（branch1_No1（m)，branch1_No2（m)）=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表示计及发电机节点的负序网络节点导纳阵
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_Num2
YY1(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_No1（m），branch2_No2(m））-1/z1_branch（m）;
YY1（br