clc清空屏幕显示clear清空所有变量N13输入网络的节点数N25输入包括所有发电机节点的网络的节点数k3输入短路点的节点号ZZ1eyeN20ZZ2eyeN20ZZ0eyeN10输入支路各序阻抗，ZZ1表示支路正序阻抗，ZZ2表示支路负序阻抗，ZZ0表示支路零序阻抗ZZ112j01ZZ113j01ZZ123j01ZZ121ZZ112ZZ131ZZ113ZZ132ZZ123ZZ111j015ZZ122j0075ZZ1330
ZZ212j01ZZ223j01ZZ213j01ZZ221ZZ212ZZ231ZZ213ZZ232ZZ223ZZ211j015ZZ222j0075ZZ2330
ZZ012j02ZZ023j02ZZ013j02ZZ021ZZ012ZZ031ZZ013ZZ032ZZ023ZZ011j00525ZZ022j00263ZZ0330
fZZ141j005ZZ114ZZ141ZZ125j0025ZZ152ZZ125ZZ144j010ZZ155j005
ZZ241j005ZZ214ZZ241ZZ225j0025ZZ252ZZ225ZZ244j010ZZ255j005
Y1eyeN10Y2eyeN10Y0eyeN10矩阵赋初值为0forp1N1forq1N1ifpq求自导纳fora1N1ifZZ1pa0Y1pqY1pq1ZZ1pae
difZZ2pa0Y2pqY2pq1ZZ2pae
difZZ0pa0Y0pqY0pq1ZZ0pae
de
delseifZZ1pq0Y1pq1ZZ1pqe
difZZ2pq0Y2pq1ZZ2pqe
difZZ0pq0Y0pq1ZZ0pq
fe
de
de
de
dYY1eyeN20YY2eyeN20矩阵赋初值为0ZZ1110ZZ2110ZZ1220ZZ2220forp1N2forq1N2ifpq求自导纳fora1N2ifZZ1pa0YY1pqYY1pq1ZZ1pae
difZZ2pa0YY2pqYY2pq1ZZ2pae
de
delseifZZ1pq0YY1pq1ZZ1pqe
difZZ2pq0YY2pq1ZZ2pqe
de
de
de
d计算各序电流的故障分量求节点导纳矩阵，Y1表示不计发电机节点的正序网络节点导纳阵，Y2表示不计发电机节点的负序网络节点导纳阵，Y0表示不计发电机节点的零序网络节点导纳阵，Y1Y2Y0节点导纳矩阵参考结果，YY1表示计及发电机节点的正序网络节点导纳阵，YY2表示计及发电机节点的负序网络节点导纳阵矩阵赋初值为0
fYY1YY2正序节点阻抗Z1eyeN1Y1eyeN1是3阶的单位矩阵求节点导纳矩阵Y1的逆阵即节点阻抗矩阵ZZ1i
vY1与上式等效负序节点阻抗Z2eyeN1Y2零序节点阻抗Z0eyeN1Y0
Z11eyeN2YY1Z22eyeN2YY2Y1Y2Y0YY1YY2
fr