电力系统三相短路计算的MATLAB代码

(k,2);Zt=line(k,3)+li*line(k,4);Yt=l/Zt;Ym=line(k,5)+li*line(k,6);ifJ〜=0Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+Yt+Ym;Y(I,J)=Y(I,J)-Yt;Y(J,I)=Y(I,J);endifJ==0Y(IZI)=Y(I,I)+Ym;endendendYmatrix3.mfunctionY=Ymatrix3(bus,line)%考虑变压器变比为1,采用支路参数近似1[nb,mb]=size(bus);[nl,ml]=size(line);Y=zeros(nb,nb);fork=l:nlI=line(k,1);J=line(k,2);Zt=line(k,3)+li*line(k,4);Yt=imag(1/Zt);Ym=imag(line(k,5)+li*line(k,6));ifJ~=0Y(IzI)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+Yt+Ym;Y(I,J)=Y(I,J)-Yt;Y(JZI)=Y(IZJ);endifJ==0Y(I,I)=Y(I,I)+Ym;endendendYmatrix4.mfunctionY=Ymatrix4(bus,line)%变压器变比近似为1,采用支路等效参数2[nb,mb]=size(bus);[nl,ml]=size(line);Y=zeros(nbznb);fork=l:nlI=line(k,1);J=line(k,2);Zt=li*line(k,4);Yt=l/Zt;Ym=li*line(k,6);ifJ〜=0Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+Yt+Ym;Y(IZJ)=Y(IZJ)-Yt;Y(J,I)=Y(IZJ);endifJ==0Y(I,I)=Y(I,I)+Ym;endendend