现电力系统短路计算课程设计.docVIP

目 录 1手工计算………………………………………………………………1 2两相接地短路计算框图………………………………………………………5 3两相接地短路源程序…………………………………………………………6 4结果输出………………………………………………………………………11 设计总结………………………………………………………………………19 2两相接地短路计算框图 两相接地短路计算流程框图如图2.1所示。 图2.1 两相接地短路计算流程框图 3两相接地短路源程序 clear all; format long e Y1=[-j*40.000000 j*10.000000 j*10 j*20 0 j*10 -j*40 j*10 0 j*20.000000 j*10 j*10 -j*20 0 0 j*20.000000 0 0 -j*25 0 0 j*20.000000 0 0 -j*25.000000]; %输入正序负序导纳矩阵 Y0=[-j*86.666667 j*33.333333 j*33.333333 j*20.000000 j*33.333333 -j*86.666667 j*33.333333 0 j*33.333333 j*33.333333 -j*66.666667 0 j*20.000000 0 0 -j*27.142857 ]; %输入零序导纳矩阵、 N1=input(please input 网络节点数 ：); N2=input(please input 零序网络节点数 ：); N3=input(please input 短路节点号 ：); disp(fault表示输入短路类型 1表示单相a短路 2表示两相bc短路 3表示两相bc短路接地) fault=input(please input fault的值 ：); %求正序和负序因子表 YY1=zeros(N1,N1); for m=1:N1 for n=m :N1 if n==m YY1(m,m)=Y1(m,m); for k=1:m-1 YY1(m,m)=YY1(m,m)-YY1(k,m)*YY1(k,m)*YY1(k,k); end else YY1(m,n)=Y1(m,n); for k=1:m-1 YY1(m,n)=YY1(m,n)-YY1(k,m)*YY1(k,n)*YY1(k,k); %原来的程序在这里出错将n写成m了 end YY1(m,n)=YY1(m,n)/YY1(m,m); end end end for m=1:N1 YY1(m,m)=1/YY1(m,m);%原来的程序在这里出错将yy写出y了 end YY2=YY1; disp(正序因子表); YY1 disp(负序因子表); YY2 %求零序因子表 if fault~=2 YY0=zeros(N2,N2); for m=1:N2 for n=m :N2 if n==m YY0(m,m)=Y0(m,m); for k=1:m-1 YY0(m,m)=YY0(m,m)-YY0(k,m)*YY0(k,m)*YY0(k,k); end else YY0(m,n)=Y0(m,n); for k=1:m-1 YY0(m,n)=YY0(m,n)-YY0(k,m)*YY0(k,n)*YY0(k,k); %原来的程序在这里出错将n写成m了 end YY0(m,n)=YY0(m,n)/YY0(m,m);%原