短路故障计算.docxVIP

一、三相短路故障计算MATLAB程序n=4; %独立节点数nl=3; %支路数B1=[1 3 0.51i 0 1 0;2 3 0.59i 0 1 0;3 4 1.43i 0 1 0]; %线路参数形成的矩阵X=[1 0.2i;2 4i;3 0;4 0]; %对地阻抗形成的矩阵V0=[1;1;1;1]; %各节点的初电压标幺值形成的列矩阵D=[4 0]; %短路点阻抗组成的矩阵NF=1; %短路点数目B=[0;0;0;1]; %常数项矩阵(短路点对应为1,其余节点对应为0)Y=zeros(n); %初始化节点导纳矩阵for i=1:n %初始化节点是否非直接接地 if X(i,2)~=0; %判断i节点是否非直接接地 p=X(i,1); Y(p,p)=1./X(i,2); %算自导中接地阻抗的导纳 endend%形成节点导纳矩阵for i=1:nl if B1(i,6)==0 %判断变压器的变比是否在低压侧(it=0说明i侧为低压侧；it=1说明i侧为高压侧) p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); endY(p,q)=Y(p,q)-1./B1(i,3)*B1(i,5); %节点互导纳 Y(q,p)=Y(p,q); %节点互导纳Y(p,p)=Y(p,p)+1./B1(i,3)*B1(i,5)^2 +B1(i,4)./2; %节点导纳=自导纳+X修正量Y(q,q)=Y(q,q)+1./B1(i,3)+B1(i,4)./2;enddisp(导纳矩阵);disp(Y); %输出导纳矩阵A=Y; [n,m]=size(A);%解线性方程组，形成矩阵因子表Afor i=1:n A(i,i)=1./A(i,i); for j=i+1:n A(i,j)=A(i,j)*A(i,i); %矩阵规格化 end for k=i+1:n for j=i+1:n A(k,j)=A(k,j)-A(k,i)*A(i,j); end endenddisp(矩阵A的因子表为);disp(A);B=input(请输入常数项矩阵(短路点对应为1,其余节点对应为0),:B=);%利用因子表对常数项进行的前代过程（修正矩阵B）for i=1:n B(i)=B(i)*A(i,i); for j=i+1:n B(j)=B(j)-A(j,i)*B(i); endenddisp(利用因子表对常数项进行会带的结果为: B=) ;disp(B)%利用因子表的上三角回带过程for i=n-1:-1:1 for j=i+1:-1:2 B(j-1)=B(j-1)-A(j-1,i+1)*B(i+1); endenddisp(在因子表的基础上求解线性方程组的解为：X=);disp(B); %解出各节点电压即得到Zik，此时B为列向量V0=input(请输入由各节点的初始电压标幺值形成的列矩阵：V0=);D=input(请输入由短路号、短路点阻抗组成的矩阵：D=);NF=input(请输入短路点的数目：NF=);Z=zeros(n);V=zeros(n);I=zeros(nl);for k=1:NF %求各短路点的电流标幺值 disp(短路点与其他各节点的互阻抗); for i=1:n Z(i,D(k,1))=B(i,1); disp(Z(i,D(k,1))); end I(D(k,1),D(k,1))=V0(D(k,1),1)./(Z(D(k,1),D(k,1))+D(k,2)); ft=num2str(D(k,1)); ts1=(点短路时); ts2=(电流的标幺值If=); dn=strcat(ft,ts1,ts2); disp(dn)