三相短路接地实验.docVIP

三相短路的计算的计算机算法 NF input 请输入短路点的数目：NF ; n input 请输入独立节点数：n ; nl input 请输入支路数：nl ; B input 请输入由线路参数形成的矩阵：B ; V0 input 请输入由各节点的初电压标么值形成的列矩阵:V0 ; D input 请输入由短路号，短路点阻抗组成的矩阵D ; m 0;Z zeros n ;V zeros n ;I zeros nl ; for k1 1:nl p B k1,1 ;q B k1,2 ; if B k1,6 0 k 1./B k1,5 ; else k B k1,5 ; end if p 0 if q m Z q,q B k1,3 ;m m+1; else for i 1:m, Z i,m+1 -Z i,q ;Z m+1,i -Z q,i ; end Z m+1,m+1 Z q,q +B k1,3 ; for i 1:m for j 1:m Z i,j Z i,j -Z i,m+1 *Z m+1,j ./Z m+1,m+1 ; end Z i,m+1 0; end for i 1:m+1 Z m+1,i 0; end end else if q m for i 1:m Z i,q Z i,p *k;Z q,i Z p,i *k; end Z q,q k^2*Z p,p +k^2*B k1,3 ; m m+1; else for i 1:m Z i,m+1 k*Z i,p -Z i,q ; Z m+1,i k*Z p,i -Z q,i ; end Z m+1,m+1 k^2*Z p,p +Z q,q -2*k*Z p,q +k^2*B k1,3 ; for i 1:m for j 1:m Z i,j Z i,j -Z i,m+1 *Z m+1,j ./Z m+1,m+1 ; end Z i,m+1 0; end for i 1:m+1 Z m+1,i 0; end end end end for k 1:NF 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 ; disp I D k,1 ,D k,1 ; for i 1:n V i,i V0 i,1 -I D k,1 ,D k,1 *Z i,D k,1 ; end for i 1:nl if B i,6 0 k B i,5 ; else k 1./B i,5 ; end p B i,1 ;q B i,2 ; if p 0 e 0;b B i,3 ; I i,i e-V q,q ./k ./b; else I i,i V p,p -V q,q ./k ./B i,3 ; end end disp 各节点的电压标么值U为（节点号从小到大排）: ; for i 1:n disp V i,i ; end disp 各支路短路电流的标么值I为（顺序同您输入B时一样）: ; for i 1:nl disp I i,i ; end end 请输入短路点的数目：NF 1 请输入独立节点数：n 4 请输入支路数：nl 5 请输入由线路参数形成的矩阵：B [0 1 2i 0 1 0;0 2 2.1i 0 1 0;1 3 1.51i 0 1 0;2 3 0.5i 0 1 0;3 4 1.4i 0 1 0] 请输入由各节点的初电压标么值形成的列矩阵:V0 [1;1;1;1] 请输入由短路号，短路点阻抗组成的矩阵D [4,0] 4点短路时电流的标么值If 0 - 0.3456i 各节点的电压标么值U为（节点号从小到大排）: 0.7059 0.5831 0.4838 0 各支路短路电流的标么值I为（顺序同您输入B时一样）: 0 + 0.3529i 0 + 0.2777i 0 - 0.1471i 0 - 0.1985i 0 - 0.3456i 电气F1203 徐晖 座机电话号码0803