新版电力系统暂态上机计算课程设计报告完整无误版附程序.doc

课程设计汇报
( —第二学期)

名 称：电力系统暂态上机计算
院 系： 电气与电子工程学院
班 级： 电气1211
学 号： 05
学生姓名： 郝阳
指导教师： 陈艳波
设计周数： 两周

成 绩：

日期： 7月4日
一、课程设计目与规定
巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算基本措施，并深入巩固计算机编程能力，为未来从事有关技术工作打下必要基础。

二、设计正文（详细内容见附录,用A4纸，页数不限）
1． 对称短路计算过程流程图和计算成果
2． 不对称短路计算过程流程图和计算成果
3． 静态稳定性计算过程流程图和计算成果
4． 暂态稳定性计算过程流程图和计算成果
5． 思索题

三、课程设计总结或结论
本次电力系统暂态上机重要讨论了计算多种类型短路故障下系统网络中电压电流分布以及电力系统静态稳定、暂态稳定问题。通过本次课程设计，本人对电力系统故障分析有了更深刻理解，电力系统故障时，大部分电磁量将随时间变化，描述其特性是微分方程，这给分析计算带来一定困难。在分析过程中一般尽量防止对微分方程直接求解，而是采用一定工具和假设使问题得以简化，即把“微分方程代数化，暂态分析稳态化”。在分析不对称故障时，各相之间电磁量耦合使问题分析更为复杂，此时常用分析措施是采用对称分量法将不对称问题转化为对称问题来求解。同步我对用来分析电力系统静态稳定试探法，用来分析电力系统暂态稳定改善欧拉法有了某些使专心得。这与手算系统短路电流时使用网络化简措施大大不一样。
在学习中，参照时尚程序，我加深了对节点导纳矩阵建立措施理解与学习，巩固了不一样类型短路短路电流计算措施，和序电压、序电流，相电压、相电流基于matlab软件计算。运用小干扰分析法鉴别静态稳定可以使用劳斯判据，也可以使用特性根鉴别措施。两者都可以判断，在书中例题中使用是劳斯判据，不过，特性根鉴别放法更适合在matlab软件环境下使用，由于该软件提供了计算矩阵特性根函数，使用起来简便易懂。运用改善欧拉法计算最大切除角或切除时间虽然比较以便，不过由于计算机有效位数限制而引起舍入误差随步长减小以及运算次数增多而增大。在上机编程中体会比较深是自己犯几种错误，例如说在算短路电流时，发电机之路电流和电压与否要通过相位变换，在静态稳定计算中角度是弧度制，在计算中使用都是标幺值等等，这些看似细微地方，往往考验就是自己对电力系统分析中知识基础，虽然小不过也非常重要，往往就是这些错误会令程序得不到对成果，值得我注意和反思。
总来说，结合我们上学期时尚上机编程经验，虽然这次用是matlab，不过明显感觉到对计算机处理电力系统问题简便，程序编写也显得愈加轻松和得心应手，体会到matlab软件强大。通过两个星期回忆与学习，我巩固了电力系统暂态分析理论知识，在编程过程中很深刻感受到要想得出对成果，就必须认真地理解书本上对应例题，懂得这些原理，不能仅仅照着书抄公式；同步我也初步掌握了采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算基本措施，并深入巩固计算机编程能力，为未来从事有关技术工作打下必要基础，更好地建设我国电力系统。

四、参照文献
1. 《电力系统暂态过程》，常鲜戎、赵书强编，机械工业出版社，1月，第一版；
2. 《电力系统分析基础》，李庚银，机械工业出版社，，第一版；

附录（设计流程图、计算成果、思索题答案）
1. 对称短路计算过程流程图和计算成果
流程图：
输入数据
选择故障类型
用节点导纳阵求逆得到节点阻抗矩阵
计算各节点序电压、相电压

根据故障类型计算故障点序电流、相电流
计算各支路序电流、相电流

计算发电机节点序电压、相电压
形成节点导纳矩阵

选择短路点
开始
计算发电机支路序电流、相电流
结束

计算成果：
（1） 导纳矩阵：
正序阻抗：

负序阻抗：
零序阻抗：

（2） 节点3发生三相短路故障
a. 故障点三序电流、三相电流有效值：
I1 = 0.0000 - 9.8592i Ia = 9.8592
I2 = 0 Ib = 9.8592
I0 = 0 Ic = 9.8592
电压
b. 各节点三序电压、三相电压有效值：
节点

U1
U2
U0
Ua
Ub
Uc
1
0.4507
0
0
0.4507
0.4507
0.4507
2
0.5352
0
0
0.5352
0.5352
0.5352
3
0.0000
0
0
0.0000
0.0000
0.0000
4
0.6338
0
－
0.6338
0.6338
0.6338
5
0.6901
0
－
0.6901
0.6901
0.6901
c. 各支路三序电流、三相电流有效值：

支路

I1
I2
I0
Ia
Ib
Ic
1-2
0.0000 + 0.8451i
0
0
0.8451
0.8451
0.8451
1-3
0.0000 - 4.5070i
0
0
4.5070
4.5070
4.5070
2-3
0.0000 - 5.3521i
0
0
5.3521
5.3521
5.3521
1-4
0.0000 + 3.6620i
0
0
3.6620
3.6620
3.6620
2-5
0.0000 + 6.1972i
0
0
6.1972
6.1972
6.1972
G1
0.0000 - 3.6620i
0
—
3.6620
3.6620
3.6620
G2
0.0000 - 6.1972i
0
－
6.1972
6.1972
6.1972
电流
2. 不对称短路计算过程流程图和计算成果
（1） 节点3发生A相短路接地故障：
a. 故障点三序电流、三相电流有效值：
I1 = 0.0000 - 3.1167i Ia = 9.3500
I2 = 0.0000 - 3.1167i Ib = 0.0000
I0 = 0.0000 - 3.1167i Ic = 0.0000
节点
电压
b. 各节点三序电压、三相电压有效值：

U1
U2
U0
Ua
Ub
Uc
1
0.8264
-0.1736
-0.0686
0.5841
0.9518
0.9518
2
0.8531
-0.1469
-0.0436
0.6625
0.9526
0.9526
3
0.6839
-0.3161
-0.3678
0.0000
1.0268
1.0268
4
0.8842
-0.1158
－
0.8324
1.0000
0.8324
5
0.9020
-0.0980
－
0.8573
1.0000
0.8573
c. 各支路三序电流、三相电流有效值：
支路

I1
I2
I0
Ia
Ib
Ic
1-2
0.0000 + 0.2671i
0.0000 + 0.2671i
0.0000 + 0.1247i
0.6589
0.1425
0.1425
1-3
0.0000 - 1.4248i
0.0000 - 1.4248i
0.0000 - 1.4960i
4.3455
0.0712
0.0712
2-3
0.0000 - 1.6919i
0.0000 - 1.6919i
0.0000 - 1.6207i
5.0045
0.0712
0.0712
1-4
0.0000 + 1.1576i
0.0000 + 1.1576i
0.0000 + 1.3713i
3.6866
0.2137
0.2137
2-5
0.0000 + 1.9590i
0.0000 + 1.9590i
0.0000 + 1.7453i
5.6634
0.2137
0.2137
G1
0.0000 - 1.1576i
0.0000 - 1.1576i
—
2.0050
0.0000
2.0050
G2
0.0000 - 1.9590i
0.0000 - 1.9590i
－
3.3932
0.0000
3.3932
电流
（2） 节点3发生A相经10Ω电阻接地故障
a. 故障点三序电流、三相电流有效值：
I1 = 0.9793 - 2.7705i Ia = 8.8154
I2 = 0.9793 - 2.7705i Ib = 0.0000
I0 = 0.9793 - 2.7705i Ic = 0.0000
节点
电压
b. 各节点三序电压、三相电压有效值：

U1
U2
U0
Ua
Ub
Uc
1
0.8456 - 0.0546i
-0.1544 - 0.0546i
-0.0610 - 0.0215i
0.6437
0.9269
0.9867
2
0.8694 - 0.0462i
-0.1306 - 0.0462i
-0.0388 - 0.0137i
0.7080
0.9281
0.9869
3
0.7190 - 0.0993i
-0.2810 - 0.0993i
-0.3269 - 0.1156i
0.3333
1.0375
1.0100
4
0.8971 - 0.0364i
-0.1029 - 0.0364i
－
0.8149
1.0000
0.8888
5
0.9129 - 0.0308i
-0.0871 - 0.0308i
－
0.8432
1.0000
0.9043
电流
c. 各支路三序电流、三相电流有效值：
支路

I1
I2
I0
Ia
Ib
Ic
1-2
-0.0839 + 0.2375i
-0.0839 + 0.2375i
-0.0392 + 0.1108i
0.6213
0.1343
0.1343
1-3
0.4477 - 1.2665i
0.4477 - 1.2665i
0.4701 - 1.3298i
4.0971
0.0672
0.0672
2-3
0.5316 - 1.5040i
0.5316 - 1.5040i
0.5093 - 1.4406i
4.7183
0.0672
0.0672
1-4
-0.3638 + 1.0290i
-0.3638 + 1.0290i
-0.4309 + 1.2190i
3.4758
0.
0.
2-5
-0.6156 + 1.7414i
-0.6156 + 1.7414i
-0.5484 + 1.5515i
5.3396
0.
0.
G1
0.3638 - 1.0290i
0.3638 - 1.0290i
—
1.8904
0.0000
1.8904
G2
0.6156 - 1.7414i
0.6156 - 1.7414i
－
3.1992
0.0000
3.1992
（3） 节点3发生b、c两相相间短路故障
a. 故障点三序电流、三相电流有效值：
I1 = 0.0000 - 4.9296i Ia = 0
I2 = 0.0000 + 4.9296i Ib = 8.5383
I0 = 0 Ic = 8.5383
节点
电压
b. 各节点三序电压、三相电压有效值：

U1
U2
U0
Ua
Ub
Uc
1
0.7254
0.2746
0
1.0000
0.6343
0.6343
2
0.7676
0.2324
0
1.0000
0.6818
0.6818
3
0.5000
0.5000
0
1.0000
0.5000
0.5000
4
0.8169
0.1831
－
0.9222
0.6338
0.9222
5
0.8451
0.1549
－
0.9322
0.6901
0.9322
电流
c. 各支路三序电流、三相电流有效值：
支路

I1
I2
I0
Ia
Ib
Ic
1-2
0.0000 + 0.4225i
0.0000 + 0.4225i
0
0.0000
0.7319
0.7319
1-3
0.0000 - 2.2535i
0.0000 - 2.2535i
0
0.0000
3.9032
3.9032
2-3
0.0000 - 2.6761i
0.0000 - 2.6761i
0
0.0000
4.6351
4.6351
1-4
0.0000 + 1.8310i
0.0000 + 1.8310i
0
0.0000
3.1714
3.1714
2-5
0.0000 + 3.0986i
0.0000 + 3.0986i
0
0.0000
5.3669
5.3669
G1
0.0000 - 1.8310i
0.0000 - 1.8310i
—
1.8310
3.6620
1.8310
G2
0.0000 - 3.0986i
0.0000 - 3.0986i
－
3.0986
6.1972
3.0986
（4） 节点3发生b、c两相接地短路故障
a. 故障点三序电流、三相电流有效值：
I1 = 0.0000 - 6.4114i Ia = 0.0000
I2 = 0.0000 + 3.4478i Ib = 9.6262
I0 = 0.0000 + 2.9636i Ic = 9.6262

节点
电压
b. 各节点三序电压、三相电压有效值：

U1
U2
U0
Ua
Ub
Uc
1
0.6428
0.1921
0.0652
0.9001
0.5258
0.5258
2
0.6977
0.1625
0.0415
0.9018
0.6049
0.6049
3
0.3497
0.3497
0.3497
1.0491
0.0000
0.0000
4
0.7619
0.1281
－
0.8333
0.6338
0.8333
5
0.7985
0.1084
－
0.8578
0.6901
0.8578

电流
c. 各支路三序电流、三相电流有效值：
支路

I1
I2
I0
Ia
Ib
Ic
1-2
0.0000 + 0.5495i
0.0000 - 0.2955i
0.0000 - 0.1185i
0.1355
0.7719
0.7719
1-3
0.0000 - 2.9309i
0.0000 + 1.5761i
0.0000 + 1.4225i
0.0677
4.4322
4.4322
2-3
0.0000 - 3.4805i
0.0000 + 1.8717i
0.0000 + 1.5411i
0.0677
5.1947
5.1947
1-4
0.0000 + 2.3814i
0.0000 - 1.2806i
0.0000 - 1.3040i
0.2032
3.6737
3.6737
2-5
0.0000 + 4.0300i
0.0000 - 2.1672i
0.0000 - 1.6596i
0.2032
5.9596
5.9596
G1
0.0000 - 2.3814i
0.0000 + 1.2806i
—
2.0643
3.6620
2.0643
G2
0.0000 - 4.0300i
0.0000 + 2.1672i
－
3.4934
6.1972
3.4934

3. 静态稳定性计算过程流程图和计算成果
流程图：
否
否
否
是
是
是
开始
清屏
设初值：Xd_all,Xdd_all,Xq,Xq_all,Eq0,UG0,Xe,Xd,Xdd
Ke=0.1
Ke_result=Ke,delta_result=delta,P_result=Eq*U*sin(delta*pi/180)/Xd_all,
delta=49
求A阵特性值
求Eq,Eqq
求UGd,UGq,UG
求K1~K6
Ke<5.7
delta<=110
(V1<0)&(V2<0)&(V3<0)

结束
打印Ke_result，delta_result, P_result

Ke=Ke+1
delta=delta+1

计算成果：
Ke=0.100000, delta=92, P=1.345052 Ke=0.00, delta=93, P=1.363051
Ke=0.300000, delta=94, P=1.380840 Ke=0.400000, delta=96, P=1.398218
Ke=0.500000, delta=97, P=1.415607 Ke=0.600000, delta=98, P=1.432731
Ke=0.700000, delta=99, P=1.449575 Ke=0.800000, delta=100, P=1.466124
Ke=0.900000, delta=101, P=1.482360 Ke=1.000000, delta=102, P=1.498265
Ke=1.100000, delta=102, P=1.514127 Ke=1.00, delta=103, P=1.529419
Ke=1.300000, delta=104, P=1.544342 Ke=1.400000, delta=104, P=1.559111
Ke=1.500000, delta=105, P=1.573400 Ke=1.600000, delta=106, P=1.587273
Ke=1.700000, delta=106, P=1.601053 Ke=1.800000, delta=106, P=1.614344
Ke=1.900000, delta=103, P=1.623638 Ke=2.000000, delta=100, P=1.626568
Ke=2.100000, delta=98, P=1.62837 Ke=2.00, delta=95, P=1.620900
Ke=2.300000, delta=93, P=1.615423 Ke=2.400000, delta=90, P=1.598429
Ke=2.500000, delta=88, P=1.586369 Ke=2.600000, delta=86, P=1.571759
Ke=2.700000, delta=84, P=1.554738 Ke=2.800000, delta=82, P=1.535446
Ke=2.900000, delta=81, P=1.527090 Ke=3.000000, delta=79, P=1.504406
Ke=3.100000, delta=77, P=1.479812 Ke=3.00, delta=76, P=1.468345
Ke=3.300000, delta=74, P=1.440938 Ke=3.400000, delta=73, P=1.427839
Ke=3.500000, delta=71, P=1.397958 Ke=3.600000, delta=70, P=1.383436
Ke=3.700000, delta=69, P=1.368427 Ke=3.800000, delta=67, P=1.335646
Ke=3.900000, delta=66, P=1.319478 Ke=4.000000, delta=65, P=1.302918
Ke=4.100000, delta=64, P=1.285990 Ke=4.00, delta=63, P=1.268718
Ke=4.300000, delta=62, P=1.251124 Ke=4.400000, delta=61, P=1.233229
Ke=4.500000, delta=60, P=1.215054 Ke=4.600000, delta=59, P=1.196618
Ke=4.700000, delta=58, P=1.177939 Ke=4.800000, delta=57, P=1.159036
Ke=4.900000, delta=56, P=1.139925 Ke=5.000000, delta=55, P=1.120623
Ke=5.100000, delta=54, P=1.101144 Ke=5.00, delta=53, P=1.081504
Ke=5.300000, delta=52, P=1.061715 Ke=5.400000, delta=52, P=1.061829
Ke=5.500000, delta=51, P=1.041864 Ke=5.600000, delta=50, P=1.021780
Ke=5.700000, delta=49, P=1.001590

分析成果：
最终选择放大倍数Ke=2.1

特性曲线：

4. 暂态稳定性计算过程流程图和计算成果
流程图：
是
否
开始
结束
输入系统数据，赋初值（delta,omega,步长h,时段长度）Duration等）等

输入故障切除时间CutTime
i=0
i<=round(CutTime/h)

求时段初功角、角速度变化率
求时段末功角、角速度估计值
求时段末功角、角速度变化率
求平均功角、角速度变化率
求时段末功角、角速度计算值
i=i+1
显示此时功角
同理求故障切除后功角变化过程
画出摇摆曲线

摇摆曲线：
(1) 0.15s时切除故障：

(2) 0.25s时切除故障：

分析成果：
根据摇摆曲线判断：0.15s时切除故障系统暂态稳定，0.25s时切除故障系统失稳。
5. 思索题
(1) 计算短路电流，书中给出手算措施与计算机编程措施有何区别？
答：手算短路电流时，要将用各序等值电抗表达复杂网络简化为一种电抗值，并根据边界条件计算出三序网络关系做出复合序网图，然后计算出正负零各序电流，之后转化为abc三相短路电流。而用计算机时可以不用简化网络,求出各序导纳矩阵后运用I=YU来计算各节点正负零各序电流后再转换成三相短路电流。

(2) 用计算机措施，进行电力系统时尚计算和短路电流计算，有哪些区别？
答：时尚计算时所求导纳矩阵为正序导纳矩阵，而短路计算时规定出正负零三序导纳矩阵。
时尚计算时需要迭代计算出成果，而短路计算一次就可以计算出短路电流。
时尚计算需要选定PQ和PV 节点，短路计算不需要辨别节点类型。

(3) 假如交给你一种任务，请你用已经有短路电流计算软件计算某地区电网某点发生短路后短路电流，你应当做哪些工作？（搜集哪些数据，怎样整顿，计算成果怎样分析）
答： 1.网络中各个发电机容量、额定电压、功率因数和正负序阻抗；变压器容量、额定电压、额定变比和短路电压百分数；各个线路电抗。
2.计算各元件等值电抗，求出系统等值电路；
3.短路前运行状况分析计算；
4.修改程序中参数，计算故障分量。

(4) 静态稳定概念，暂态稳定概念？
答：静态稳定：指电力系统受到小干扰后，不发生非周期性失步或自发振荡，自动恢复到初始运行状态能力。
暂态稳定：指电力系统受到大干扰后，各发电机组保持同步运行并过渡到新或恢复到本来稳定运行状态能力。

(5) 提高静态稳定、暂态稳定措施？
答：提高静态稳定措施：
一、 采用自动调整励磁装置
二、 减小元件电抗
1、 采用分裂导线
2、 采用串联电容赔偿
三、 提高线路标称等级
四、 改善系统构造和采用中间赔偿设备
提高暂态稳定措施：
一、变化制动功率
1、故障迅速切除和应用自动重叠闸
2、对发电机进行强行励磁
3、电气制动
4、变压器中性点经小电阻接地
5、输电线路设置开关站
6、输电线路采用强行串联电容赔偿
二、变化原动功率
1、迅速自动调速系统或迅速关闭汽门
2、联锁切除部分发电机
3、合理选择远距离输电系统运行接线
三、系统失去稳定后措施
1、设置解列点
2、短期异步运行和再同步
3、做好系统“黑启动”方案

6. 程序源代码
（1） formAdmitanceMatrix.m
function [Z1,Z0]=formAdmitanceMatrix(nbus_1,nbus_0,nbranch_1,nbranch_0,branch_1,branch_0)

if( nargin<3 )
disp('The input parameters are not enough for calculation! please check function formAdmitanceMatrix!');
end

%---------------计算节点导纳矩阵----------%
Y1=zeros(nbus_1);
for m=1:nbranch_1
from_bus=branch_1(m).from;
to_bus=branch_1(m).to;
%write your programme here
No1=from_bus;
No2=to_bus;
if No2==0
Y1(No1,No1)=Y1(No1,No1)+1/branch_1(m).x;
else
Y1(No1,No1)=Y1(No1,No1)+1/branch_1(m).x;
Y1(No2,No2)=Y1(No2,No2)+1/branch_1(m).x;
Y1(No1,No2)=-1/branch_1(m).x;
Y1(No2,No1)=-1/branch_1(m).x;
end
end
disp('正序阻抗');
Y1
Y2=Y1; %负序等于正序
disp('负序阻抗');
Y2
Y0=zeros(nbus_0);
for m=1:nbranch_0
from_bus=branch_0(m).from;
to_bus=branch_0(m).to;
%write your programme here
No1=from_bus;
No2=to_bus;
if No2==0
Y0(No1,No1)=Y0(No1,No1)+1/branch_0(m).x0;
else
Y0(No1,No1)=Y0(No1,No1)+1/branch_0(m).x0;
Y0(No2,No2)=Y0(No2,No2)+1/branch_0(m).x0;
Y0(No1,No2)=-1/branch_0(m).x0;
Y0(No2,No1)=-1/branch_0(m).x0;
end
end
disp('零序阻抗');
Y0
%---------------计算节点阻抗矩阵----------%
Z1=inv(Y1);
Z2=inv(Y2);
Z0=inv(Y0);

return;
（2） shortcircuit.m
function shortcircuit()
%数据来源于教材《电力系统分析基础》P77例（3-4），P182例（7-3）

%---------------输入已知条件----------%
bus_Num_1=5; %正负序网络节点数
bus_Num_0=3; %零序网络节点数
bus(1) = struct('idx',1,'shift',0);
bus(2) = struct('idx',2,'shift',0);
bus(3) = struct('idx',3,'shift',0);
bus(4) = struct('idx',4,'shift',pi/6);
bus(5) = struct('idx',5,'shift',pi/6);

branch_Num_1=7; %正负序网络支路数
branch_Num_0=5; %零序网络支路数
branch_1(1) = struct('from',1,'to',2,'x',j*0.1,'shift',0);
branch_1(2) = struct('from',1,'to',3,'x',j*0.1,'shift',0);
branch_1(3) = struct('from',2,'to',3,'x',j*0.1,'shift',0);
branch_1(4) = struct('from',1,'to',4,'x',j*0.05,'shift',0);
branch_1(5) = struct('from',2,'to',5,'x',j*0.025,'shift',0);
branch_1(6) = struct('from',4,'to',0,'x',j*0.1,'shift',pi/6);
branch_1(7) = struct('from',5,'to',0,'x',j*0.05,'shift',pi/6);

branch_0(1) = struct('from',1,'to',2,'x0',j*0.2);
branch_0(2) = struct('from',1,'to',3,'x0',j*0.2);
branch_0(3) = struct('from',2,'to',3,'x0',j*0.2);
branch_0(4) = struct('from',1,'to',0,'x0',j*0.05);
branch_0(5) = struct('from',2,'to',0,'x0',j*0.025);

%---------------计算正序和负序导纳矩阵----------%
[Z1,Z0]=formAdmitanceMatrix(bus_Num_1,bus_Num_0,branch_Num_1,branch_Num_0,branch_1,branch_0); %计算节点导纳矩阵、节点阻抗矩阵，形成全局变量
Z2=Z1;

Fault_Node=input('输入短路点编号：;\n Fault_Node=');
Fault_Type=input('输入短路类型：;\n(1)Fault_Type=0为三相短路;\n(2)Fault_Type=1为a相接地短路;\n(3)Fault_Type=2为a相经10欧姆电阻接地短路\n(4)Fault_Type=3为bc两相相间短路\n(5)Fault_Type=4为bc两相短路接地\nFault_Type=');

a=-0.5+j*sqrt(3)/2;
T=[1 1 1 % T为对称分量法变换矩阵，见P87公式（4-4）
a^2 a 1
a a^2 1];

%---------------第一步：计算短路点序电流，相电流---------------%
%计根据故障类型选择不一样计算公式，计算故障点各序电流
if Fault_Type==0 %三相短路
I_Fault1=1/Z1(Fault_Node,Fault_Node); %三相短路只有正序电流
I_Fault2=0;
I_Fault0=0;

elseif Fault_Type==1 %a相接地短路
%write your programme here
I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node)); %3个序网阻抗串联
I_Fault2=I_Fault1;
I_Fault0=I_Fault1;

elseif Fault_Type==2 %a相经10欧姆电阻接地短路
%write your programme here
Zf=10*50/(115^2);
I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node)+3*Zf); %3个序阻抗和电阻串联
I_Fault2=I_Fault1;
I_Fault0=I_Fault1;

elseif Fault_Type==3 %bc两相相间短路
%write your programme here
I_Fault1=1/(Z1(Fault_Node,Fault_Node)+Z2(Fault_Node,Fault_Node)); %正序阻抗和负序阻抗并联
I_Fault2=-I_Fault1;
I_Fault0=0;

elseif Fault_Type==4 %bc两相短路接地
%write your programme here
%三个序阻抗并联
I_Fault1=1/(Z1(Fault_Node,Fault_Node)+(Z2(Fault_Node,Fault_Node)*Z0(Fault_Node,Fault_Node)/(Z2(Fault_Node,Fault_Node)+Z0(Fault_Node,Fault_Node))));