电力系统暂态上机计算课程设计报告

1、 课程设计报告( 20122013年度第二学期)名 称：电力系统暂态上机 院 系：电气与电子工程学院 成 绩： 日期： 2013年 7月1日一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。二、设计正文（详细内容见附录,用a4纸，页数不限）1 对称短路计算过程流程图和计算结果2 不对称短路计算过程流程图和计算结果3 静态稳定性计算过程流程图和计算结果4 暂态稳定性计算过程流程图和计算结果5、 思考题6、 实验程序三、课程设计总结或结论本次电力系统暂

2、态上机主要讨论了计算各种类型短路故障下系统网络中的电压电流分布以及电力系统的静态稳定、暂态稳定问题。在短路电流的计算过程中，先通过计算系统各个节点的导纳矩阵，再通过判断系统的故障节点号和故障类型代入相应的公式进行计算。其中，将故障后的电气量分成故障分量和正常分量，通过故障类型判断故障分量和正常分量的值，相加后通过对称分量法的矩阵转换即可得出所求的电气量的数值。在静态稳定性的计算中，计算出发电机中的各个参数值和放大倍数，根据相应的对应关系通过特征值判据判断稳定性，进而计算pe的值，从而判断其最大功率。在暂态稳定性的计算中，设置了相应的切除时间和分段时间，每一个分段时间进行一次计算，在切除故障前和

3、后分别进行计算，进而得出相应曲线，并通过多次尝试的除最大切除时间大约为2s。通过此次学习，了解到电力系统的故障时，大部分电磁量将随时间变化，描述其特性的是微分方程，这给分析计算带来一定困难。在分析过程中通常尽量避免对微分方程直接求解，而是采用一定的工具和假设使问题得以简化，即把“微分方程代数化，暂态分析稳态化”。在分析不对称故障时，各相之间电磁量的耦合使问题的分析更为复杂，此时常用的分析方法是采用对称分量法将不对称问题转化为对称问题来求解。另外，通过此次课程设计，我学习并掌握了一定的matlab知识，对简单的选择语句、循环语句、数组的使用等有了一些了解，为今后学习、工作中再次学习使用此软件打下

4、了基础，获益匪浅。四、参考文献1. 电力系统暂态分析，李光琦，中国电力出版社，2006年，第三版；2. 电力系统分析（上、下），何仰赞，华中科技大学出版社，1998年，第二版。3. 电力系统故障的计算机辅助分析 重庆大学出版社 米麟书等 4. 电力系统潮流计算 天津大学出版社 宋文南等 5. 电力系统故障分析 清华大学出版社 周荣光6. 短路电流实用计算方法 电力工业出版社 西安交通大学等7. 精通matlab6.5北京航空航天大学出版社，张志涌附录（设计流程图、计算结果、思考题答案）1、对称短路计算流程图及计算结果l 过程流程图输入数据根据fault值选择故障类型用节点导纳阵求逆得到节点阻抗

5、矩阵选择短路点f计算各序电流的故障分量计算故障点各序电流计算故障点各序电压的故障分量计算故障后各节点电压值形成节点导纳矩阵确定fault值开始计算支路电流计算计及发电机端点的各节点电压结束l 计算结果 节点导纳矩阵 图1 y1、y2的计算结果 图2 y0的计算结果 图3 yy1、yy2的结算结果(2) 节点3发生三相短路故障a、故障点三序电流、三相电流：i1= 0 - 9.8592i i2=0 i0=0 ; ia=0 - 9.8592i ib=-8. 5383+ 4.9296i ic=8.5383 + 4.9296ib.各节点三序电压、三相电压：节点u1u2u0uaubuc10.4507000

6、.4507-0.2254- 0.3903i-0.2254+ 0.3903i20.5352000.5352-0.2676- 0.4635i-0.2676 + 0.4635i300000040.6338000.5489 + 0.3169i0.0000 - 0.6338i-0.5489+ 0.3169i50.6091000.5977+ 0.3451i0.0000 - 0.6901i-0.5977 + 0.3451ic.各支路三序电流、三相电流：支路i1i2i0iaibic1-20 + 0.8451i000 + 0.8451i0.7319 - 0.4225i-0.7319 - 0.4225i1-30

7、- 4.5070i000 - 4.5070i-3.9032 + 2.2535i3.9032 + 2.2535i2-30 - 5.3521i000 - 5.3521i-4.6351 + 2.6761i4.6351 + 2.6761i1-4 0 + 3.6620i000+3.6620i3.1714-1.8310i-3.1714-1.8310i2-5 0 + 6.1972i000+6.1972i5.3669-3.0968i-5.3669-3.0968i发电机1-1.8310 + 3.1714i00-1.8310 + 3.1714i3.6620-1.8310 -3.1714i发电机23-3.0986

8、+ 5.3669i00-3.0986 + 5.3669i6.1972 - 0.0000i-3.0986 -5.3669i2、不对称短路计算的流程图及计算结果l 流程图输入数据根据fault值选择故障类型用节点导纳阵求逆得到节点阻抗矩阵选择短路点f计算各序电流的故障分量计算故障点各序电流计算故障点各序电压的故障分量计算故障后各节点电压值形成节点导纳矩阵确定fault值开始计算支路电流计算计及发电机端点的各节点电压结束l 计算结果(1) 节点3发生a相短路接地故障a. 故障点三序电流、三相电流：i1= 0 - 3.1080i i2= 0 - 3.1080i i0= 0 - 3.1080i ; ia

9、= 0 - 9.3241i ib= 0 - 0.0000i ic= 0 - 0.0000ib. 各节点三序电压、三相电压：节点u1u2u0uaubuc10.8268-0.1732-0.07150.5822-0.3983 - 0.8660i -0.3983 + 0.8660i20.8535-0.1465-0.04590.6610-0.3994 - 0.8660i-0.3994 + 0.8660i30.6848-0.3152-0.36950.0000 -0.5543 - 0.8660i -0.5543 + 0.8660i40.8846-0.115400.6661 + 0.5000i0.0000 -

10、 1.0000i-0.6661 + 0.5000i50.9023 -0.097700.6968 + 0.5000i0.0000 - 1.0000i-0.6968 + 0.5000ic. 各支路三序电流、三相电流：支路i1i2i0iaibic1-20 + 0.2664i0 + 0.2664i0 + 0.1280i0 + 0.6608i0 - 0.1384i0 - 0.1384i1-30 - 1.4208i0 - 1.4208i0 - 1.4900i0 - 4.3517i0.0000 - 0.0692i-0.0000 - 0.0692i2-30 - 1.6872i0 - 1.6872i0 - 1.

11、6180i0 - 4.9924i0.0000 + 0.0692i-0.0000 + 0.0692i1-40 + 1.1544i0 + 1.1544i 0 + 1.3621i0 + 3.6709i0.0000 + 0.2077i-0.0000 + 0.2077i2-50 + 1.9536i0 + 1.9536i 0 + 1.7459i0 + 5.6532i 0.0000 - 0.2077i- 0.0000 - 0.2077i发电机1 -0.5772 + 0.9997i 0.5772 + 0.9997i0-0.0000 + 1.9995i0.0000 - 0.0000i-0.0000 - 1.99

12、95i发电机23-0.9768 + 1.6919i0.9768 + 1.6919i0 -0.0000 + 3.3838i0.0000 - 0.0000i -0.0000 -3.3838i(2) 节点3发生a相经10电阻接地故障a. 故障点三序电流、三相电流：i1=1.1152 - 2.6363i i2= 1.1152 - 2.6363i i0= 1.1152 - 2.6363i ; ia= 3.3456-7.9088i ib= 0 ic= 0b. 各节点三序电压、三相电压：节点u1u2u0uaubuc10.8531 - 0.0621i-0.1469 - 0.0621i-0.0607 - 0.0

13、257i0.6456 - 0.1499i -0.4138 - 0.8296i-0.4138 + 0.9025i2 0.8757 - 0.0526i-0.1243 - 0.0526i-0.0389 - 0.0165i0.7125 - 0.1216i-0.4147 - 0.8299i -0.4147 + 0.9021i 30.7326 - 0.1131i-0.2674 - 0.1131i-0.3134 - 0.1326i 0.1518 - 0.3588i -0.5460 - 0.8855i-0.5460 + 0.8466i40.9021 - 0.0414i-0.0979 - 0.0414i00.6

14、964 + 0.4283i0.0000 - 1.0000i-0.6964 + 0.5717i50.9171 - 0.0350i-0.0829 - 0.0350i0 0.7225 + 0.4393i0.0000 - 1.0000i-0.7225 + 0.5607ic. 各支路三序电流、三相电流：支路i1i2i0iaibic1-2-0.0956 + 0.2260i -0.0956 + 0.2260i-0.0459 + 0.1085i-0.2371 + 0.5605i0.0497 - 0.1174i 0.0497 - 0.1174ii1-3 0.5098 - 1.2052i0.5098 - 1.20

15、52i0.5346 - 1.2639i1.5542 - 3.6742i0.0248 - 0.0587i0.0248 - 0.0587i2-30.6054 - 1.4311i0.6054 - 1.4311i0.5806 - 1.3724i1.7913 - 4.2346i-0.0248 + 0.0587i -0.0248 + 0.0587i1-4-0.4142 + 0.9792i-0.4142 + 0.9792i-0.4887 + 1.1553i-1.3172 + 3.1137i-0.0745 + 0.1761i-0.0745 + 0.1761i2-5-0.7010 + 1.6571i-0.701

16、0 + 1.6571i -0.6265 + 1.4809i -2.0284 + 4.7951i0.0745 - 0.1761i0.0745 - 0.1761i发电机1-0.8483 + 0.6409i0.1309 + 1.0551i0-0.7174 + 1.6960i0.0000 - 0.0000i0.7174 - 1.6960i发电机23 -1.4356 + 1.0846i 0.2215 + 1.7856i0-1.2141 + 2.8702i-0.0000 - 0.0000i1.2141 - 2.8702i(3) 节点3发生b、c两相短路故障a. 故障点三序电流、三相电流：i1= 0 - 4

17、.9296i i2= 0 + 4.9296i i0= 0 ; ia= 0 ib= -8.6078 - 0.0000i ic= 8.6078 + 0.0000ib.各节点三序电压、三相电压：节点u1u2u0uaubuc1 0.72540.274601.0000-0.5000 - 0.3903i-0.5000 + 0.3903i20.76760.232401.0000-0.5000 - 0.4635i-0.5000 + 0.4635i30.50000.500001.0000-0.5000-0.500040.81690.183100.8660 + 0.3169i0.0000 - 0.6338i-0.

18、8660 + 0.3169i50.84510.15490 0.8660 + 0.3451i0.0000 - 0.6901i-0.8660 + 0.3451ic.各支路三序电流、三相电流：支路i1i2i0iaibic1-20 + 0.4225i0 - 0.4225i000.7319 -0.7319 1-30 - 2.2535i0 + 2.2535i00-3.9032 3.9032 2-30 - 2.6761i0 + 2.6761i00-4.63514.63511-40 + 1.8310i0 +-1.8310i003.1714 3.17142-5 0 + 3.0986i 0 -3.0986i005

19、.3669-5.3669发电机1-0.9155 + 1.5857i-0.9155 - 1.5857i0-1.83103.6620-1.8310发电机2、3-1.5493 + 2.6835i-1.5493 - 2.6835i0-3.09866.1972-3.0986(4) 节点3发生b、c两相短路接地故障a. 故障点三序电流、三相电流：i1= 0 - 6.4036i i2= 0 + 3.4556i i0= 0 + 2.9480i ; ia= 0 ib= -8.5383 + 4.4220i ic= 8.5383+ 4.4220ib. 各节点三序电压、三相电压：节点u1u2u0uaubuc10.64

20、320.19250.06780.9036 -0.3500 - 0.3903i-0.3500 + 0.3903i20.69810.16290.04360.9046-0.3870 - 0.4635i-0.3870 + 0.4635i30.35050.35050.35051.05150040.76220.128300.7712 + 0.3169i0.0000 - 0.6338i-0.7712 + 0.3169i50.79870.108600.7858 + 0.3451i0.0000 - 0.6901i-0.7858 + 0.3451ic. 各支路三序电流、三相电流：支路i1i2i0iaibic1-2

21、0 + 0.5489i0 - 0.2962i0 - 0.1214i0 + 0.1304i0.7319 - 0.2477i-0.7319 - 0.2477i1-30 - 2.9274i0 + 1.5797i0 + 1.4133i0 + 0.0657i-3.9032 + 2.0872i3.9032 + 2.0872i2-30 - 3.4762i0 + 1.8759i0 + 1.5347i0 - 0.0657i-4.6351 + 2.3349i4.6351 + 2.3349i1-4 0 + 2.3785i0 - 1.2835i0 - 1.2920i0 - 0.1970i3.1714 - 1.8394

22、i -3.1714 - 1.8394i2-50 + 4.0251i0 - 2.1721i0 - 1.6561i0 + 0.1970i5.3669 - 2.5826i-5.3669 - 2.5826i发电机1-1.1892 + 2.0598i-0.6417 - 1.1115i0-1.8310 + 0.9483i3.6620 -1.8310 - 0.9483i发电机2、3-2.0126 + 3.4858i-1.0860 - 1.8811i0-3.0986 + 1.6048i6.1972 + 0.0000i-3.0986 - 1.6048i3静态稳定性计算过程流程图和计算结果l 过程流程图否否否是是

23、是开始清屏设初值：xd_all,xdd_all,xq,xq_all,eq0,ug0,xe,xd,xddke=0.1ke_result=ke,delta_result=delta, k1_result= k1,kemax_result=kemaxp_result=eq*u*sin(delta*pi/180)/xd_all,delta=delta+0.01delta=49.00求kemax,kemin求eq,eqq求ugd,ugq,ug求k1k6ke5.7delta=110.00特征根实部是否小于0结束打印ke_result，delta_result, k1_result,kemax_result

24、l 计算结果ke=0.100000, delta=91, p=1.345117 ke=0.200000, delta=92, p=1.363000 ke=0.300000, delta=93, p=1.380693ke=0.400000, delta=95, p=1.398409ke=0.500000, delta=96, p=1.415747ke=0.600000, delta=97, p=1.432842ke=0.700000, delta=98, p=1.449679ke=0.800000, delta=99, p=1.466244ke=0.900000, delta=100, p=1.

25、482520ke=1.000000, delta=101, p=1.498487ke=1.100000, delta=101, p=1.513985ke=1.200000, delta=102, p=1.529380ke=1.300000, delta=103, p=1.544429ke=1.400000, delta=103, p=1.558878ke=1.500000, delta=104, p=1.573336ke=1.600000, delta=105, p=1.587402ke=1.700000, delta=105, p=1.600902ke=1.800000, delta=105

26、, p=1.613924ke=1.900000, delta=102, p=1.621480ke=2.000000, delta=99, p=1.622734ke=2.100000, delta=97, p=1.623379ke=2.200000, delta=94, p=1.614349ke=2.300000, delta=92, p=1.607801ke=2.400000, delta=89, p=1.589402ke=2.500000, delta=87, p=1.576380ke=2.600000, delta=85, p=1.560855ke=2.700000, delta=83,

27、p=1.542968ke=2.800000, delta=81, p=1.522859ke=2.900000, delta=80, p=1.514020ke=3.000000, delta=78, p=1.490598ke=3.100000, delta=76, p=1.465314ke=3.200000, delta=75, p=1.453440ke=3.300000, delta=73, p=1.425420ke=3.400000, delta=72, p=1.411962ke=3.500000, delta=70, p=1.381541ke=3.600000, delta=69, p=1

28、.366704ke=3.700000, delta=68, p=1.351397ke=3.800000, delta=66, p=1.318171ke=3.900000, delta=65, p=1.301745ke=4.000000, delta=64, p=1.284942ke=4.100000, delta=63, p=1.267785ke=4.200000, delta=62, p=1.250299ke=4.300000, delta=61, p=1.232503ke=4.400000, delta=60, p=1.214420ke=4.500000, delta=59, p=1.19

29、6069ke=4.600000, delta=58, p=1.177469ke=4.700000, delta=57, p=1.158638ke=4.800000, delta=56, p=1.139594ke=4.900000, delta=55, p=1.120353ke=5.000000, delta=54, p=1.100930ke=5.100000, delta=53, p=1.081340ke=5.200000, delta=52, p=1.061599ke=5.300000, delta=51, p=1.041718ke=5.400000, delta=51, p=1.04179

30、2ke=5.500000, delta=50, p=1.021746ke=5.600000, delta=49, p=1.001589ke=5.600000, delta=49, p=1.001589因此，ke=2.1时，p最大，静态稳定性最好4、暂态稳定性计算过程流程图和计算结果l 过程流程图 是否开始结束赋初值（delta,omega,步长h,时段长度duration等）等输入故障切除时间cuttimei=0i=round(cuttime/h)求时段初变化率求时段末估计值求时段末变化率求平均变化率求时段末计算值i=i+1显示此时功角同理求故障切除后的功角变化过程画出摇摆曲线l 计算结果 0

31、.15s切除故障时的摇摆曲线 0.25s切除故障时的摇摆曲线0.15s切除故障时，系统为暂态稳定，最大切 此时系统出现失稳现象除时间大约为0.2s。5、思考题a. 计算短路电流，书中给出的手算方法与计算机编程方法有何区别？答：手算短路电流时，要将用各序等值电抗表示的复杂网络简化为一个电抗值，并根据边界条件计算出三序网络的关系做出复合序网图，然后计算出正负零各序电流，之后转化为abc三相短路电流。而用计算机时可以不用简化网络,求出各序导纳矩阵后利用i=yu来计算各节点正负零各序电流后再转换成三相短路电流。b. 用计算机方法，进行电力系统潮流计算和短路电流计算，有哪些区别？答：1、计算机潮流计算是

32、根据发电机和线路参数算出系统中各线路的功率和电压分布，所计算的数值实际各节点电压和各个线路流过的功率。而短路计算所涉及的电气量为电压和电流，需计算短路点的电压电流，以及故障后其他节点的电压和各个线路流过的电流。 2、. 潮流计算时所求的导纳矩阵为正序的导纳矩阵，而短路计算时要求出正负零三序的导纳矩阵。 3、潮流计算时需要迭代计算出结果，而短路计算一次就可以计算出短路电流。4、潮流计算需要选定pq和pv 节点，短路计算不需要区分节点类型。c.如果交给你一个任务，请你用已有的短路电流计算软件计算某地区电网某点发生短路后的短路电流，你应该做哪些工作？（收集哪些数据，如何整理，计算结果如何分析）答：1

33、.网络中各个发电机的容量、额定电压、功率因数和正负序阻抗；变压器的容量、额定电压、额定变比和短路电压百分数；各个线路的电抗。2.计算各元件等值电抗，求出系统的等值电路；3.短路前运行状况的分析计算；4.修改程序中的参数，计算故障分量。d.静态稳定的概念，暂态稳定的概念？答：静态稳定是指电力系统受到小干扰之后，不发生非周期性失步或自发震荡，自动恢复到初始运行状态的能力。 暂态稳定是指电力系统受到大的干扰后，各发电机组保持同步运行并且过渡到新的或恢复到原来稳定运行状态的能力。e.提高静态稳定、暂态稳定的措施？答：提高静态稳定性的措施有 采用自动调节励磁装置 减小元件的电抗 提高线路标称电压等级 改

34、善系统的结构和采用中间补偿设备 提高暂态稳定性的措施有n 改变制动功率1 故障的快速切除和自动重合闸装置的应用2 对发电机进行强行励磁3 电气制动4 变压器中性点经小电阻接地5 输电线路设置开关站6 输电线路采用强行串联电容补偿 n 改变原动功率1 快速的自动调速系统或者快速关闭进气门2 连锁切除部分发电机3 合理合理选择远距离输电系统的运行接线6、原始程序l matrixclc;clear; %-输入已知条件-%bus_num1=3; %节点数bus_num2=5; %包括发电机节点的节点数branch_num1=3; %线路数branch_num2=5; %包括发电机支路的支路数branc

35、h1_no1=1,1,2;branch1_no2=2,3,3;branch2_no1=1,1,2,1,2;branch2_no2=2,3,3,4,5; %输入支路各序阻抗，z1_branch表示支路正序阻抗，z2_branch表示支路负序阻抗，z0_branch表示支路零序阻抗z1_branch(1)=j*0.1;z1_branch(2)=j*0.1;z1_branch(3)=j*0.1; z2_branch=z1_branch; z0_branch(1)=j*0.2;z0_branch(2)=j*0.2;z0_branch(3)=j*0.2; %-第一步：不考虑发电机节点计算节点导纳矩阵-%

36、节点导纳矩阵，y1表示不计发电机节点的正序网络节点导纳阵，y2表示不计发电机节点的负序网络节点导纳阵，y0表示不计发电机节点的零序网络节点导纳阵，y1=zeros(bus_num1);y1(1,1)=1/(j*0.15);y1(2,2)=1/(j*0.075);%正序节点导纳矩阵for i = 1:branch_num1 a=1/z1_branch(i); b=branch1_no1(i); c=branch1_no2(i); y1(b,b)=y1(b,b)+a; y1(b,c)=y1(b,c)-a; y1(c,b)=y1(c,b)-a; y1(c,c)=y1(c,c)+a;end y2=y1

37、;%负序等于正序 y0=zeros(bus_num1);y0(1,1)=1/(j*0.0525);y0(2,2)=1/(j*0.0263);%零序节点导纳矩阵for i = 1:branch_num1 a=1/z0_branch(i); b=branch1_no1(i); c=branch1_no2(i); y0(b,b)=y0(b,b)+a; y0(b,c)=-a; y0(c,b)=-a; y0(c,c)=y0(c,c)+a;end %-第二步：考虑发电机节点计算节点导纳矩阵-%节点导纳矩阵，yy1表示计及发电机节点的正序网络节点导纳阵，yy2表示计及发电机节点的负序网络节点导纳阵z1_br

38、anch(4)=j*0.05;z1_branch(5)=j*0.025; z2_branch=z1_branch; yy1=zeros(bus_num2);yy1(4,4)=1/(j*0.1);yy1(5,5)=1/(j*0.05);%正序节点导纳矩阵for i=1:branch_num2 a=1/z1_branch(i); b=branch2_no1(i); c=branch2_no2(i); yy1(b,b)=yy1(b,b)+a; yy1(b,c)=-a; yy1(c,b)=-a; yy1(c,c)=yy1(c,c)+a;end yy2=yy1;%负序等于正序 z0_branch(4)=

39、j*0.1;z0_branch(5)=j*0.05;yy0=zeros(bus_num2);yy0(4,4)=1/(j*0.2);yy0(5,5)=1/(j*0.1);%零序节点导纳矩阵for i = 1:branch_num2 a=1/z0_branch(i); b=branch2_no1(i); c=branch2_no2(i); yy0(b,b)=yy0(b,b)+a; yy0(b,c)=-a; yy0(c,b)=-a; yy0(c,c)=yy0(c,c)+a;end %-第三步：计算节点阻抗矩阵-%z1=inv(y1);z2=inv(y2);z0=inv(y0);l shortcirc

40、uitclc;clear; %数据来源于教材电力系统暂态分析p77例（3-4），p143例（5-7） matrix %计算节点导纳矩阵、节点阻抗矩阵，形成全局变量 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*sq

41、rt(3)/2;t=1 1 1 % t为对称分量法的变换矩阵，见p87公式（4-4） a2 a 1 a a2 1; %-第一步：计算短路点的序电流，相电流-%计根据故障类型选择不同的计算公式，计算故障点各序电流if fault_type=0 %三相短路各序电流计算 i_fault1=1/z1(fault_node,fault_node); i_fault2=0; i_fault0=0; elseif fault_type=1 %a相接地短路各序电流计算 i_fault1=1/(z1(fault_node,fault_node)+z2(fault_node,fault_node)+z0(faul

42、t_node,fault_node); i_fault2=i_fault1; i_fault0=i_fault1; elseif fault_type=2 %a相经10欧姆电阻接地短路各序电流计算 i_fault1=1/(z1(fault_node,fault_node)+z2(fault_node,fault_node)+z0(fault_node,fault_node)+3*10*60/(1152); i_fault2=i_fault1; i_fault0=i_fault1;elseif fault_type=3 %bc两相相间短路各序电流计算 i_fault1=1/(z1(fault_n

43、ode,fault_node)+z2(fault_node,fault_node); i_fault2=-i_fault1; i_fault0=0;elseif fault_type=4 %bc两相短路接地各序电流计算 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); i_fault2=-i_fault1*z0(fault_node,fault_node)/

44、(z2(fault_node,fault_node)+z0(fault_node,fault_node); i_fault0=-i_fault1*z2(fault_node,fault_node)/(z2(fault_node,fault_node)+z0(fault_node,fault_node); end;str=短路序电流;stri1=i_fault1i2=i_fault2i0=i_fault0str=短路相电流;striabc=t*i_fault1 i_fault2 i_fault0. %相量iabc_effective=abs(iabc) %有效值 %-第二步：计算各个节点的序电压

45、，相电压-%故障分量 i1fault_node1=zeros(bus_num1,1);i1fault_node2=zeros(bus_num1,1);i1fault_node0=zeros(bus_num1,1);for m=1:bus_num1 if m=fault_node %故障点短路电流各序分量 i1fault_node1(m,1)=-1*i_fault1; i1fault_node2(m,1)=-1*i_fault2; i1fault_node0(m,1)=-1*i_fault0; else %非故障点各序电流各序分量 i1fault_node1(m,1)=0; i1fault_no

46、de2(m,1)=0; i1fault_node0(m,1)=0; endend%故障分量u1fault_node1=zeros(bus_num1,1);u1fault_node2=zeros(bus_num1,1);u1fault_node0=zeros(bus_num1,1);u1fault_node1=y1i1fault_node1;u1fault_node2=y2i1fault_node2;u1fault_node0=y0i1fault_node0; %正常分量 u1normal_node1=ones(bus_num1,1);u1normal_node2=zeros(bus_num1,

47、1);u1normal_node0=zeros(bus_num1,1); u1_node1=zeros(bus_num1,1);u1_node2=zeros(bus_num1,1);u1_node0=zeros(bus_num1,1);u1_node1=u1fault_node1+u1normal_node1;u1_node2=u1fault_node2+u1normal_node2;u1_node0=u1fault_node0+u1normal_node0;str=1.2.3节点电压stru1=u1_node1 u2=u1_node2u0=u1_node0uabc123=t*u1_node1

48、 u1_node2 u1_node0. %-第三步：计算各个支路的序电流，相电流-%i_branch1=zeros(bus_num1,1);i_branch2=zeros(bus_num1,1);i_branch0=zeros(bus_num1,1);for i=1:bus_num1 m=branch1_no1(i); n=branch1_no2(i); i_branch1(i,1)=(u1fault_node1(m,1)+u1normal_node1(m,1)-u1fault_node1(n,1)-u1normal_node1(n,1)/z1_branch(i); i_branch2(i,1

49、)=(u1fault_node2(m,1)+u1normal_node2(m,1)-u1fault_node2(n,1)-u1normal_node2(n,1)/z2_branch(i); i_branch0(i,1)=(u1fault_node0(m,1)+u1normal_node0(m,1)-u1fault_node0(n,1)-u1normal_node0(n,1)/z0_branch(i);endstr=各支路序电流;stri1= i_branch1i2= i_branch2i0=i_branch0str=各支路相电流;iabc1=t*i_branch1 i_branch2 i_branch0.iabc1_effective=abs(iabc1) %-第四步：计算发电机节点的序电压，相电压-%i2fault_node1=zeros(bus_num2,1);i2fault_node2=zeros(bus_num2,1);i2fault_node0=zeros(bus_num2,1);for m=1:bus_num2 if m=fault_node %故障点短路电流各序分量 i2fault_node1(m,1)=-1*i_fault1; i2fault_node2(m,1)=-1*i_fault2; i2fault_node0(m,1)