电力系统各种短路故障的MATLAB算法.docVIP

function [sxdldl]=sxdldl() zoko=input(总控矩阵[功率基准值SB 发电机数n1 双绕组变压器数n2 a类三绕组变压器数n3a b类三绕组变压器数n3b 线路数n4 节点总数n]=) %不经变压器直接连在公共母线上的取一台，且电抗百分值去等值电抗 、三绕组变不在短路网络的节点不计入总数（即忽略的带负荷支路） %a类为考虑三绕组变在短路网络中只有两个绕组接入，即有一个绕组由于带负荷忽略掉；b类为三个电压等级都在短路网络中 sb=zoko(1,1); n1=zoko(1,2); n2=zoko(1,3); n3a=zoko(1,4); n3b=zoko(1,5); n4=zoko(1,6); n=zoko(1,7); %以上为取控制量 zg=input(发电机支路[所在节点 额定功率 额定功率因数 电抗]=); %发电机支路原始参数 zt2=input(双绕组变压器支路控制阵[节点一 节点二 额定容量 电抗百分值]=); %双绕组变压器支路原始参数 if n3a~=0 zt3a=input(三绕组变压器支路控制阵a[节点一 节点二 额定容量 短路网络包含绕组（12 13 23 ） 1-2 1-3 2-3]=); %三绕组变压器支路原始参数,小号在前！ %短路网络包含绕组：12 13 23 分别指包含1绕组2绕组、1绕组3绕组、2绕组3绕组 end if n3b~=0 zt3b=input(三绕组变压器支路控制阵b[节点一 节点二 节点三 额定容量 1-2 1-3 2-3]=); %三绕组变压器支路原始参数,小号在前！ end zl=input(线路支路控制阵（节点号小号在前！）[节点一 节点二 长度 阻抗/每公里 电压等级（平均额定电压）]=); %线路支路原始参数，小号在前！ %其他支路 n5=input(电抗器数n5=); %电抗器支路原始数据 if n5~=0 dkq=input(电抗器参数矩阵[节点一 节点二 额定电压 额定电流 电抗百分数 所在支路电压等级（平均额定电压）]=); end %等值系统支路 n6=input(等值系统数n6=); if n6~=0 dzxt=input(等值系统参数[所在节点 短路电流（有名值） 短路电流（标幺值） 电压等级(平均额定电压)]=); %以短路电流形式给出阻抗参数 end %%%%%%%%%%%%%%%%%%%形成系统阻抗%%%%%%%%%%%%%%%% z=zeros(n,n); %存放系统阻抗矩阵 %%%%%发电机支路形成系统阻抗%%%%%% for i=1:n1 z(zg(i,1),zg(i,1))=zg(i,4)*sb*zg(i,3)/zg(i,2) %发电机电抗标幺值 end %%%%双绕组变压器支路形成系统矩阵%%%%% for i=1:n2 z(zt2(i,1),zt2(i,2))=zt2(i,4)*sb/zt2(i,3) %双绕组变压器电抗标幺值 end %%%%计算三绕组变压器的阻抗(并形成系统阻抗)%%%% if n3a~=0 %若有a类三绕组或者自耦变压器 ZT3=zeros(n3a,3); for i=1:n3a ZT3(i,1)=(zt3a(i,5)+zt3a(i,6)-zt3a(i,7))*0.5*sb/zt3a(i,3) %1绕组电抗标幺值 ZT3(i,2)=(zt3a(i,5)+zt3a(i,7)-zt3a(i,6))*0.5*sb/zt3a(i,3) %2绕组电抗标幺值 ZT3(i,3)=(zt3a(i,6)+zt3a(i,7)-zt3a(i,5))*0.5*sb/zt3a(i,3) %3绕组电抗标幺值 if zt3a(i,4)==12 z(zt3a(i,1),zt3a(i,2))=ZT3(i,1)+ZT