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

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);%请同学们求正序节点导纳矩阵for m=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(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);%请同学们求零序节点导纳矩阵for m=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))=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);%请同学们求正序节点导纳矩阵for m=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_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_branch(m);end YY2=YY1;%负序等于正序 %---------------第三步：计算节点阻抗矩阵----------%Z1=inv(Y1);Z2=inv(Y2);Z0=inv(Y0);clc;clear; %数据来源于教材《电力系统暂态分析》P77例（3-4），P143例（5-7） Matrix %计算节点导纳矩阵、节点阻抗矩阵，形成全局变量 Fault_Node=input('输入短路点编号 ;\n Fault_Node='); Fault_Type=input('输入短路类型 ;\n(1)Fault_Type=0为三相短路;\n(2)Fault_Type=1为a相接地短路;\n(3)Fault_Type=2为a相经10欧姆电阻接地短路\n(4)Fault_Type=3为bc两相相间短路\n(5)Fault_Type=4为bc两相短路接地\nFault_Type=') a=-0.5+j*sqrt(3)/2;T=[1 1 1 % T为对称分量法的变换矩阵，见P87公式（4-4） a^2 a 1 a a^2 1]; %---------------第一步：计算短路点的序电流，相电流---------------%%计根据故障类型选择不同的计算公式，计算故障点各序电流if Fault_Type==0 I_Fault1=1/Z1(Fault_Node,Fault_Node); I_Fault2=0; I_Fault0=0; elseif Fault_Type==1 I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node)); I_Fault2=I_Fault1; I_Fault0=I_Fault1; elseif Fault_Type==2 I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node)+30*60/115^2); I_Fault2=I_Fault1; I_Fault0=I_Fault1; elseif Fault_Type==3 I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)); I_Fault2=-I_Fault1; I_Fault0=0; elseif Fault_Type==4 I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)*Z0(Fault_Node,Fault_Node)/(Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node))); I_Fault2=-I_Fault1*Z0(Fault_Node,Fault_Node)/(Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node)); I_Fault0=-I_Fault1*Z2(Fault_Node,Fault_Node)/(Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node)); end; str='短路电流';strI_Fault1I_Fault2I_Fault0 %故障电流序分量Iabc=T*[I_Fault1 I_Fault2 I_Fault0].' %相分量相量Iabc_effective=abs(Iabc) %有效值 %---------------第二步：计算各个节点的序电压，相电压---------------%%故障分量 Ifault_node1=zeros(bus_Num1,1);Ifault_node2=zeros(bus_Num1,1);Ifault_node0=zeros(bus_Num1,1);for m=1:bus_Num1 if m==Fault_Node Ifault_node1(m,1)=-1*I_Fault1; Ifault_node2(m,1)=-1*I_Fault2; Ifault_node0(m,1)=-1*I_Fault0; else Ifault_node1(m,1)=0; Ifault_node2(m,1)=0; Ifault_node0(m,1)=0; endend Ufault_node1=zeros(bus_Num1,1);Ufault_node2=zeros(bus_Num1,1);Ufault_node0=zeros(bus_Num1,1);Ufault_node1=Y1\Ifault_node1;Ufault_node2=Y2\Ifault_node2;Ufault_node0=Y0\Ifault_node0; %正常分量 Unormal_node1=ones(bus_Num1,1);Unormal_node2=zeros(bus_Num1,1);Unormal_node0=zeros(bus_Num1,1); %故障后电压值=故障分量+正常分量Uall_node1=Ufault_node1+Unormal_node1Uall_node2=Ufault_node2+Unormal_node2Uall_node0=Ufault_node0+Unormal_node0 %序分量str='各节点相电压';strUabc=T*[Uall_node1 Uall_node2 Uall_node0].' %相量Uabc_effective=abs(Uabc) %有效值 %---------------第三步：计算各个。