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);endY2=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);endYY2=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）a2 a 1a a2 1;%-第一步： 计算短路点的序电流，相 电流-%计根据故障 类 型选择不同的计算公式，计算故障点各序电流if Fault_Type=0I_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=2I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node)+30*60/1152);I_Fault2=I_Fault1;I_Fault0=I_Fault1; elseif Fault_Type=3I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node);I_Fault2=-I_Fault1;I_Fault0=0; elseif Fault_Type=4I_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_NodeIfault_node1(m,1)=-1*I_Fault1;Ifault_node2(m,1)=-1*I_Fault2;Ifault_node0(m,1)=-1*I_Fault0; elseIfault_node1(m,1)=0;Ifault_node2(m,1)=0;Ifault_node0(m,1)=0;endendUfault_node1=zeros(bus_Num1,1);Ufault_node2=zeros(bus_Num1,1);Ufault_node0=zeros(bus_Num1,1);Ufault_node1=Y1Ifault_node1;Ufault_node2=Y2Ifault_node2;Ufault_node0=Y0Ifault_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) %有效值%-第三步： 计算各个支路的序电流，相 电流-%for k=1:branch_Num1m=branch1_No1(k);n=branch1_No2(k); Ibranch1(k)=-(Uall_node1(m)-Uall_node1(n)*Y1(m,n);Ibranch2(k)=-(Uall_node2(m)-Uall_node2(n)*Y1(m,n);Ibranch0(k)=-(Uall_node0(m)-Uall_node0(n)*Y0(m,n); %或者不加符号，而是除以z0_branch(k),像第五步那样子mn(1)=m;mn(2)=n;mn %显示支路号str= 该支路的序电流Ibranch1(k)Ibranch2(k)Ibranch0(k) str= 该支路的相电流Ibranchabc=T*Ibranch1(k) Ibranch2(k) Ibranch0(k). %相量Ibranchabc_effective=abs(Ib