其中计算模块包含一个子流程图：
主要数据变量说明： 程序中：
n 节点数
f 短路节点序号
z 任意节点间阻抗矩阵 g_z 发电机电抗矩阵 g s 发电厂额定容量矩阵 byqy 双绕组变压器参数矩阵 byqy3 三绕组变压器参数矩阵 sz 系统参数矩阵 Z 网络节点阻抗矩阵 Y 网络节点导纳矩阵 If 短路点电流
V 网络节点电压矩阵 I 网络支路电流矩阵
6
四、关键程序及程序注释
for i=1:n %以下生成节点导纳矩阵 for j=1:n if z(i,j)~=0
Y(i,i)=Y(i,i)+1/z(i,j); end end end
%考虑发电机参数影响 for pp=1:p i=gz(pp,1);
Y(i,i)=Y(i,i)+1/sqrt(-1)/gz(pp,2); end
%考虑系统参数影响 for ss=1:s i=sz(ss,1);
Y(i,i)=Y(i,i)+1/sqrt(-1)/sz(ss,2); end
for i=1:n %计算非对角线元素 for j=i+1:n if z(i,j)==0 Y(i,j)=0; else Y(i,j)=-1/z(i,j); end Y(j,i)=Y(i,j); end end
%考虑双绕组变压器参数影响 for qq=1:q2 i=byqy(qq,1); j=byqy(qq,2);
7
Y(i,i)=Y(i,i)+1/sqrt(-1)/byqy(qq,4); Y(j,j)=Y(j,j)+1/sqrt(-1)/byqy(qq,4); Y(i,j)=Y(i,j)-1/sqrt(-1)/byqy(qq,4); Y(j,i)=Y(i,j); end
%考虑三绕组变压器参数影响 for qq3=1:q3 i=byqy3(qq3,1); j=byqy3(qq3,2); p=byqy3(qq3,3); q=byqy3(qq3,4);
Y(i,i)=Y(i,i)+1/sqrt(-1)/byqy3(qq3,8); Y(j,j)=Y(j,j)+1/sqrt(-1)/byqy3(qq3,9); Y(p,p)=Y(p,p)+1/sqrt(-1)/byqy3(qq3,10);
Y(q,q)=Y(q,q)+(1/byqy3(qq3,8)+1/byqy3(qq3,9)+1/byqy3(qq3,10))/1i; Y(i,q)=Y(i,q)-1/sqrt(-1)/byqy3(qq3,8); Y(j,q)=Y(j,q)-1/sqrt(-1)/byqy3(qq3,9); Y(p,q)=Y(p,q)-1/sqrt(-1)/byqy3(qq3,10); Y(q,i)=Y(i,q); Y(q,j)=Y(j,q); Y(q,p)=Y(p,q); end
Z=inv(Y);
%计算短路点电流 If=1/Z(f,f);
%计算各节点电压 for i=1:n
V(i)=1-If*Z(i,f);
8
end
for i=1:n %计算各支路电流 for j=1:n
I(i,j)=(V(i)-V(j))/z(i,j); end end
五、输入参数，以及输出计算结果与分析
算例接线图
1．输入
节点总数n=10 短路点序号f=4 发电机台数p=4
双绕组变压器数q2=2 系统数s=1
9
三绕组变压器数q3=1 线路数x=5
具体参数输入
2．输出数据及分析
短路电流大小：If= -6.412600
Y矩阵如下
32.619048 -33.333333 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
-33.333333 28.797508 2.666667 1.869159 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
0.000000 2.666667 -2.666667 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000
0.000000 1.869159 0.000000 -8.790706 6.921547 0.000000 0.000000 0.000000 0.000000 0.000000
0.000000 0.000000 0.000000 6.921547 -20.112996 6.736451 6.454998 0.000000 0.000000 0.000000
0.000000 0.000000 0.000000 0.000000 6.736451 -98.981476 82.245025 0.000000 0.000000 0.000000
0.000000 0.000000 0.000000 0.000000 6.454998 82.245025 -95.298337 6.598314 0.000000 0.000000
0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 6.598314 -10.234677 1.818182
10