资源描述：
《短路故障计算》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、一、三相短路故障计算MATLAB程序n=4;%独立节点数nl=3;%支路数B1=[130.51i010;230.59i010;341.43i010];%线路参数形成的矩阵X=[10.2i;24i;30;40];%对地阻抗形成的矩阵V0=[1;1;1;1];%各节点的初电压标幺值形成的列矩阵D=[40];%短路点阻抗组成的矩阵NF=1;%短路点数目B=[0;0;0;1];%常数项矩阵(短路点对应为1,其余节点对应为0)Y=zeros(n);%初始化节点导纳矩阵fori=1:n%初始化节点是否非直接接地ifX(i,2)~=0;%判断i节点是否非直接接地p=
2、X(i,1);Y(p,p)=1./X(i,2);%算自导中接地阻抗的导纳endend%形成节点导纳矩阵fori=1:nlifB1(i,6)==0%判断变压器的变比是否在低压侧(it=0说明i侧为低压侧；it=1说明i侧为高压侧)p=B1(i,1);q=B1(i,2);elsep=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修正
3、量Y(q,q)=Y(q,q)+1./B1(i,3)+B1(i,4)./2;enddisp('导纳矩阵');disp(Y);%输出导纳矩阵A=Y;[n,m]=size(A);%解线性方程组，形成矩阵因子表Afori=1:nA(i,i)=1./A(i,i);forj=i+1:nA(i,j)=A(i,j)*A(i,i);%矩阵规格化endfork=i+1:nforj=i+1:nA(k,j)=A(k,j)-A(k,i)*A(i,j);endendenddisp('矩阵A的因子表为');disp(A);B=input('请输入常数项矩阵(短路点对应为1,其余节点
4、对应为0),:B=');%利用因子表对常数项进行的前代过程（修正矩阵B）fori=1:nB(i)=B(i)*A(i,i);forj=i+1:nB(j)=B(j)-A(j,i)*B(i);endenddisp('利用因子表对常数项进行会带的结果为:B=');disp(B)%利用因子表的上三角回带过程fori=n-1:-1:1forj=i+1:-1:2B(j-1)=B(j-1)-A(j-1,i+1)*B(i+1);endenddisp('在因子表的基础上求解线性方程组的解为：X=');disp(B);%解出各节点电压即得到Zik，此时B为列向量V0=inp
5、ut('请输入由各节点的初始电压标幺值形成的列矩阵：V0=');D=input('请输入由短路号、短路点阻抗组成的矩阵：D=');NF=input('请输入短路点的数目：NF=');Z=zeros(n);V=zeros(n);I=zeros(nl);fork=1:NF%求各短路点的电流标幺值disp('短路点与其他各节点的互阻抗');fori=1:nZ(i,D(k,1))=B(i,1);disp(Z(i,D(k,1)));endI(D(k,1),D(k,1))=V0(D(k,1),1)./(Z(D(k,1),D(k,1))+D(k,2));ft=num
6、2str(D(k,1));ts1=('点短路时');ts2=('电流的标幺值If=');dn=strcat(ft,ts1,ts2);disp(dn);disp(I(D(k,1),D(k,1)));fori=1:nV(i,i)=V0(i,1)-I(D(k,1),D(k,1))*Z(i,D(k,1));%求各节点的电压标幺值电压end%求非接地支路的短路电流标幺值fori=1:nlifB1(i,6)==0%判断该支路是否为含有变压器支路且变比在低压侧k=B1(i,5);elsek=1./B1(i,5);endp=B1(i,1);q=B1(i,2);I(i,
7、i)=(V(p,p)-V(q,q)./k)./B1(i,3);%各支路短路电流=支路残压/支路阻抗(标幺值)enddisp('各节点的电压标幺值U为(节点号从小到大排)：');fori=1:ndisp(V(i,i));enddisp('各接地支路短路电流的标幺值I为(顺序同您输入B时一样');fori=1:nifX(i,2)~=0;%判断i节点是否为直接接地e=0;b=X(i,2);Ii0=(e-V(i,i))./b;disp(Ii0);endenddisp('各非接地支路短路电流的标幺值I为(顺序同您输入B时一样)：');fori=1:nldisp(
8、I(i,i));endend二、程序输出结果导纳矩阵0-6.9608i00+1.9608i00
