资源描述：
《短路计算程序.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、n=4;%节点数nl=3;%支路数B1=[130.51i010;230.59i010;341.43i010];X=[10.2i;24i;30;40];%线路参数形成的矩阵V0=[1;1;1;1];%电压矩阵D=[4,0];%短路号矩阵NF=1;%短路点个数B=[0;0;0;1];%常数项矩阵Y=zeros(n);%初始化节点导纳矩阵fori=1:nifX(i,2)~=0;%判断i节点是否非直接接地，若是算自导中接地阻抗的导纳p=X(i,1);Y(p,p)=1./X(i,2);%算出自导中接地支路导纳endend%形成节点导纳矩阵fori=1:nlifB1(i,6)==0%判断变压器的
2、变比是否在低压侧p=B1(i,1);q=B1(i,2);elsep=B1(i,2);q=B1(i,1);%从p→ｑ，变化为1:kendY(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;Y(q,q)=Y(q,q)+1./B1(i,3)+B1(i,4)./2;enddisp('导纳矩阵Y');%输出需要的字符串disp(Y);A=Y;%赋值Y等于A[n,m]=size(A);%用解线性方程组方法，形成矩阵A因子表fori=1:nA(i,i)=1./A(i
3、,i);%第一列元素规格化，即取倒数运算forj=i+1:n%逐列进行消去A(i,j)=A(i,j)*A(i,i);endfork=i+1:n%形成因子表forj=i+1:nA(k,j)=A(k,j)-A(k,i)*A(i,j);endendenddisp('矩阵A的因子表为');disp(A);%B=input('pleaseinput');%利用因子表对常数项进行的前代过程（按列消去）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)
4、;%利用因子表的回代过程fori=n-1:-1:1%做n-1次循环，即n-1次消去运算forj=i+1:-1:2B(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);l=zeros(nl);fork=1:NF%求各短路点的电流标幺值fori=1:nZ(i,D(k,1))=B(i,1);%将求得的Zik赋值
5、到阻抗相应位置disp(Z(i,D(k,1)));endl(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));%数转换成字符串并赋值给ftts1=('点短路时');ts2=('电流的标幺值lf=');dn=(strcat(ft,ts1,ts2));%将两字符串穿起来并赋值给dndisp(dn);disp(l(D(k,1),D(k,1)));fori=1:n%求各节点的电压标幺值V(i,i)=V0(i,1)-l(D(k,1),D(k,1))*Z(i,D(k,1));end%求非接地支路的短
6、路电流标幺值fori=1:nlifB1(i,6)==0%判断该支路是否为含有变压器支路且变比在低压侧k=B1(i,5);elsek=1./B1(i,5);endp=B1(i,1);q=B1(i,2);l(i,i)=(V(p,p)-V(q,q)./k)./B1(i,3);enddisp('各节点的电压标幺值U为（节点号从小到大）');fori=1:n%i从1到n每次循环for函数disp(V(i,i));enddisp('各非接地支路短路电流的标幺值I为（顺序与输入B一样）');fori=1:nldisp(l(i,i));endend%求接地支路短路电流标幺值disp('各接地支路短路
7、电流的标幺值I为（顺序与输入B一样）');fori=1:nifX(i,2)~=0;%判断i节点是否为直接接地e=0;b=X(i,2);l0=(e-V(i,i))./b;disp(l0);endend导纳矩阵Y0.0000-6.9608i0.0000+0.0000i0.0000+1.9608i0.0000+0.0000i0.0000+0.0000i0.0000-1.9449i0.0000+1.6949i0.0000+0.0000i0.0000+1.9608i
