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

1、 课程设计报告( 20142015年度第二学期)名 称：电力系统暂态上机计算院 系： 电气与电子工程学院 班 级： 电气1211 学 号： 学生姓名： 郝阳 指导教师： 陈艳波 设计周数： 两周 成 绩： 日期： 2015年 7月4日一、课程设计的目的与要求巩固电力系统暂态分析理论知识，使学生掌握采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础。二、设计正文（详细内容见附录,用A4纸，页数不限）1 对称短路计算过程流程图和计算结果2 不对称短路计算过程流程图和计算结果3 静态稳定性计算过程流程图和计算结果4 暂

2、态稳定性计算过程流程图和计算结果5 思考题三、课程设计总结或结论本次电力系统暂态上机主要讨论了计算各种类型短路故障下系统网络中的电压电流分布以及电力系统的静态稳定、暂态稳定问题。通过本次课程设计，本人对电力系统故障分析有了更深刻的理解，电力系统的故障时，大部分电磁量将随时间变化，描述其特性的是微分方程，这给分析计算带来一定困难。在分析过程中通常尽量避免对微分方程直接求解，而是采用一定的工具和假设使问题得以简化，即把“微分方程代数化，暂态分析稳态化”。在分析不对称故障时，各相之间电磁量的耦合使问题的分析更为复杂，此时常用的分析方法是采用对称分量法将不对称问题转化为对称问题来求解。同时我对用来分析

3、电力系统静态稳定的试探法，用来分析电力系统暂态稳定的改进欧拉法有了一些使用心得。这与手算系统短路电流时使用的网络化简方法大大不同。在学习中，参照潮流程序，我加深了对节点导纳矩阵建立方法的理解与学习，巩固了不同类型短路的短路电流计算方法，和序电压、序电流，相电压、相电流基于matlab软件的计算。利用小干扰分析法判别静态稳定可以使用劳斯判据，也可以使用特征根判别的方法。二者都能够判断，在书中的例题中使用的是劳斯判据，但是，特征根判别放法更适合在matlab软件的环境下使用，因为该软件提供了计算矩阵特征根的函数，使用起来简便易懂。利用改进欧拉法计算最大切除角或切除时间虽然比较方便，但是由于计算机有

4、效位数的限制而引起的舍入误差随步长的减小以及运算次数的增多而增大。在上机编程中体会比较深的是自己犯的几个错误，比如说在算短路电流时，发电机之路的电流和电压是否要经过相位的变换，在静态稳定计算中的角度是弧度制，在计算中使用的都是标幺值等等，这些看似细微的地方，往往考验的就是自己对电力系统分析中的知识的基础，虽然小但是也非常的重要，往往就是这些错误会令程序得不到正确的结果，值得我注意和反思。总的来说，结合我们上学期潮流上机编程的经验，虽然这次用的是matlab，但是明显感觉到对计算机处理电力系统问题的简便，程序的编写也显得更加的轻松和得心应手，体会到matlab软件的强大。经过两个星期的回顾与学习

5、，我巩固了电力系统暂态分析理论知识，在编程的过程中很深刻的感受到要想得出正确的结果，就必须认真地理解课本上对应的例题，知道这些原理，不能仅仅照着书抄公式；同时我也初步掌握了采用计算机对电力系统电磁暂态过程和机电暂态过程进行计算的基本方法，并进一步巩固计算机编程能力，为将来从事相关的技术工作打下必要的基础，更好地建设我国的电力系统。四、参考文献1. 电力系统暂态过程，常鲜戎、赵书强编，机械工业出版社，2010年1月，第一版；2. 电力系统分析基础，李庚银，机械工业出版社，2011年，第一版；附录（设计流程图、计算结果、思考题答案）1. 对称短路计算过程流程图和计算结果流程图：输入数据选择故障类型

6、用节点导纳阵求逆得到节点阻抗矩阵计算各节点的序电压、相电压根据故障类型计算故障点序电流、相电流计算各支路的序电流、相电流计算发电机节点的序电压、相电压形成节点导纳矩阵选择短路点开始计算发电机支路的序电流、相电流结束计算结果：（1） 导纳矩阵：正序阻抗：负序阻抗：零序阻抗：（2） 节点3发生三相短路故障a. 故障点三序电流、三相电流有效值：I1 = 0.0000 - 9.8592iIa = 9.8592I2 = 0Ib = 9.8592I0 = 0Ic = 9.8592电压b. 各节点三序电压、三相电压有效值：节点U1U2U0UaUbUc10.4507000.45070.45070.450720

7、.5352000.53520.53520.535230.0000000.00000.00000.000040.633800.63380.63380.633850.690100.69010.69010.6901c. 各支路三序电流、三相电流有效值：支路I1I2I0IaIbIc1-20.0000 + 0.8451i000.84510.84510.84511-30.0000 - 4.5070i004.50704.50704.50702-30.0000 - 5.3521i005.35215.35215.35211-40.0000 + 3.6620i003.66203.66203.66202-50.00

8、00 + 6.1972i006.19726.19726.1972G10.0000 - 3.6620i03.66203.66203.6620G20.0000 - 6.1972i06.19726.19726.1972电流2. 不对称短路计算过程流程图和计算结果（1） 节点3发生A相短路接地故障：a. 故障点三序电流、三相电流有效值：I1 = 0.0000 - 3.1167iIa = 9.3500I2 = 0.0000 - 3.1167iIb = 0.0000I0 = 0.0000 - 3.1167iIc = 0.0000节点电压b. 各节点三序电压、三相电压有效值：U1U2U0UaUbUc10.8

9、264-0.1736-0.06860.58410.95180.951820.8531-0.1469-0.04360.66250.95260.952630.6839-0.3161-0.36780.00001.02681.026840.8842-0.11580.83241.00000.832450.9020-0.09800.85731.00000.8573c. 各支路三序电流、三相电流有效值：支路I1I2I0IaIbIc1-20.0000 + 0.2671i0.0000 + 0.2671i0.0000 + 0.1247i0.65890.14250.14251-30.0000 - 1.4248i0.

10、0000 - 1.4248i0.0000 - 1.4960i4.34550.07120.07122-30.0000 - 1.6919i0.0000 - 1.6919i0.0000 - 1.6207i5.00450.07120.07121-40.0000 + 1.1576i0.0000 + 1.1576i0.0000 + 1.3713i3.68660.21370.21372-50.0000 + 1.9590i0.0000 + 1.9590i0.0000 + 1.7453i5.66340.21370.2137G10.0000 - 1.1576i0.0000 - 1.1576i2.00500.000

11、02.0050G20.0000 - 1.9590i0.0000 - 1.9590i3.39320.00003.3932电流（2） 节点3发生A相经10电阻接地故障a. 故障点三序电流、三相电流有效值：I1 = 0.9793 - 2.7705iIa = 8.8154I2 = 0.9793 - 2.7705iIb = 0.0000I0 = 0.9793 - 2.7705iIc = 0.0000节点电压b. 各节点三序电压、三相电压有效值：U1U2U0UaUbUc10.8456 - 0.0546i-0.1544 - 0.0546i-0.0610 - 0.0215i0.64370.92690.9867

12、20.8694 - 0.0462i-0.1306 - 0.0462i-0.0388 - 0.0137i0.70800.92810.986930.7190 - 0.0993i-0.2810 - 0.0993i-0.3269 - 0.1156i0.33331.03751.010040.8971 - 0.0364i-0.1029 - 0.0364i0.81491.00000.888850.9129 - 0.0308i-0.0871 - 0.0308i0.84321.00000.9043电流c. 各支路三序电流、三相电流有效值：支路I1I2I0IaIbIc1-2-0.0839 + 0.2375i-0.

13、0839 + 0.2375i-0.0392 + 0.1108i0.62130.13430.13431-30.4477 - 1.2665i0.4477 - 1.2665i0.4701 - 1.3298i4.09710.06720.06722-30.5316 - 1.5040i0.5316 - 1.5040i0.5093 - 1.4406i4.71830.06720.06721-4-0.3638 + 1.0290i-0.3638 + 1.0290i-0.4309 + 1.2190i3.47580.20150.20152-5-0.6156 + 1.7414i-0.6156 + 1.7414i-0.5

14、484 + 1.5515i5.33960.20150.2015G10.3638 - 1.0290i0.3638 - 1.0290i1.89040.00001.8904G20.6156 - 1.7414i0.6156 - 1.7414i3.19920.00003.1992（3） 节点3发生b、c两相相间短路故障a. 故障点三序电流、三相电流有效值：I1 = 0.0000 - 4.9296iIa = 0I2 = 0.0000 + 4.9296iIb = 8.5383I0 = 0Ic = 8.5383节点电压b. 各节点三序电压、三相电压有效值：U1U2U0UaUbUc10.72540.274601

15、.00000.63430.634320.76760.232401.00000.68180.681830.50000.500001.00000.50000.500040.81690.18310.92220.63380.922250.84510.15490.93220.69010.9322电流c. 各支路三序电流、三相电流有效值：支路I1I2I0IaIbIc1-20.0000 + 0.4225i0.0000 + 0.4225i00.00000.73190.73191-30.0000 - 2.2535i0.0000 - 2.2535i00.00003.90323.90322-30.0000 - 2.

16、6761i0.0000 - 2.6761i00.00004.63514.63511-40.0000 + 1.8310i0.0000 + 1.8310i00.00003.17143.17142-50.0000 + 3.0986i0.0000 + 3.0986i00.00005.36695.3669G10.0000 - 1.8310i0.0000 - 1.8310i1.83103.66201.8310G20.0000 - 3.0986i0.0000 - 3.0986i3.09866.19723.0986（4） 节点3发生b、c两相接地短路故障a. 故障点三序电流、三相电流有效值：I1 = 0.00

17、00 - 6.4114iIa = 0.0000I2 = 0.0000 + 3.4478iIb = 9.6262I0 = 0.0000 + 2.9636iIc = 9.6262节点电压b. 各节点三序电压、三相电压有效值：U1U2U0UaUbUc10.64280.19210.06520.90010.52580.525820.69770.16250.04150.90180.60490.604930.34970.34970.34971.04910.00000.000040.76190.12810.83330.63380.833350.79850.10840.85780.69010.8578电流c.

18、各支路三序电流、三相电流有效值：支路I1I2I0IaIbIc1-20.0000 + 0.5495i0.0000 - 0.2955i0.0000 - 0.1185i0.13550.77190.77191-30.0000 - 2.9309i0.0000 + 1.5761i0.0000 + 1.4225i0.06774.43224.43222-30.0000 - 3.4805i0.0000 + 1.8717i0.0000 + 1.5411i0.06775.19475.19471-40.0000 + 2.3814i0.0000 - 1.2806i0.0000 - 1.3040i0.20323.6737

19、3.67372-50.0000 + 4.0300i0.0000 - 2.1672i0.0000 - 1.6596i0.20325.95965.9596G10.0000 - 2.3814i0.0000 + 1.2806i2.06433.66202.0643G20.0000 - 4.0300i0.0000 + 2.1672i3.49346.19723.49343. 静态稳定性计算过程流程图和计算结果流程图：否否否是是是开始清屏设初值：Xd_all,Xdd_all,Xq,Xq_all,Eq0,UG0,Xe,Xd,XddKe=0.1Ke_result=Ke,delta_result=delta,P_r

20、esult=Eq*U*sin(delta*pi/180)/Xd_all, delta=49求A阵特征值求Eq,Eqq求UGd,UGq,UG求K1K6Ke5.7delta=110(V10)&(V20)&(V30)结束打印Ke_result，delta_result, P_resultKe=Ke+1delta=delta+1计算结果：Ke=0., delta=92, P=1.Ke=0., delta=93, P=1.Ke=0., delta=94, P=1.Ke=0., delta=96, P=1.Ke=0., delta=97, P=1.Ke=0., delta=98, P=1.Ke=0., d

21、elta=99, P=1.Ke=0., delta=100, P=1.Ke=0., delta=101, P=1.Ke=1., delta=102, P=1.Ke=1., delta=102, P=1.Ke=1., delta=103, P=1.Ke=1., delta=104, P=1.Ke=1., delta=104, P=1.Ke=1., delta=105, P=1.Ke=1., delta=106, P=1.Ke=1., delta=106, P=1.Ke=1., delta=106, P=1.Ke=1., delta=103, P=1.Ke=2., delta=100, P=1.K

22、e=2., delta=98, P=1.62837Ke=2., delta=95, P=1.Ke=2., delta=93, P=1.Ke=2., delta=90, P=1.Ke=2., delta=88, P=1.Ke=2., delta=86, P=1.Ke=2., delta=84, P=1.Ke=2., delta=82, P=1.Ke=2., delta=81, P=1.Ke=3., delta=79, P=1.Ke=3., delta=77, P=1.Ke=3., delta=76, P=1.Ke=3., delta=74, P=1.Ke=3., delta=73, P=1.Ke

23、=3., delta=71, P=1.Ke=3., delta=70, P=1.Ke=3., delta=69, P=1.Ke=3., delta=67, P=1.Ke=3., delta=66, P=1.Ke=4., delta=65, P=1.Ke=4., delta=64, P=1.Ke=4., delta=63, P=1.Ke=4., delta=62, P=1.Ke=4., delta=61, P=1.Ke=4., delta=60, P=1.Ke=4., delta=59, P=1.Ke=4., delta=58, P=1.Ke=4., delta=57, P=1.Ke=4., d

24、elta=56, P=1.Ke=5., delta=55, P=1.Ke=5., delta=54, P=1.Ke=5., delta=53, P=1.Ke=5., delta=52, P=1.Ke=5., delta=52, P=1.Ke=5., delta=51, P=1.Ke=5., delta=50, P=1.Ke=5., delta=49, P=1.分析结果：最终选择放大倍数Ke=2.1特性曲线：4. 暂态稳定性计算过程流程图和计算结果流程图：是否开始结束输入系统数据，赋初值（delta,omega,步长h,时段长度）Duration等）等输入故障切除时间CutTimei=0i=ro

25、und(CutTime/h)求时段初功角、角速度变化率求时段末功角、角速度估计值求时段末功角、角速度变化率求平均功角、角速度变化率求时段末功角、角速度计算值i=i+1显示此时功角同理求故障切除后的功角变化过程画出摇摆曲线摇摆曲线：(1) 0.15s时切除故障：(2) 0.25s时切除故障：分析结果：根据摇摆曲线判断：0.15s时切除故障系统暂态稳定，0.25s时切除故障系统失稳。5. 思考题(1) 计算短路电流，书中给出的手算方法与计算机编程方法有何区别？答：手算短路电流时，要将用各序等值电抗表示的复杂网络简化为一个电抗值，并根据边界条件计算出三序网络的关系做出复合序网图，然后计算出正负零各序

26、电流，之后转化为abc三相短路电流。而用计算机时可以不用简化网络,求出各序导纳矩阵后利用I=YU来计算各节点正负零各序电流后再转换成三相短路电流。(2) 用计算机方法，进行电力系统潮流计算和短路电流计算，有哪些区别？答：潮流计算时所求的导纳矩阵为正序的导纳矩阵，而短路计算时要求出正负零三序的导纳矩阵。潮流计算时需要迭代计算出结果，而短路计算一次就可以计算出短路电流。潮流计算需要选定PQ和PV 节点，短路计算不需要区分节点类型。(3) 如果交给你一个任务，请你用已有的短路电流计算软件计算某地区电网某点发生短路后的短路电流，你应该做哪些工作？（收集哪些数据，如何整理，计算结果如何分析）答：1.网络

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

28、用串联电容补偿三、 提高线路标称等级四、 改善系统的结构和采用中间补偿设备提高暂态稳定的措施：一、改变制动功率1、故障快速切除和应用自动重合闸2、对发电机进行强行励磁3、电气制动4、变压器中性点经小电阻接地5、输电线路设置开关站6、输电线路采用强行串联电容补偿二、改变原动功率1、快速的自动调速系统或快速关闭汽门2、联锁切除部分发电机3、合理选择远距离输电系统的运行接线三、系统失去稳定后的措施1、设置解列点2、短期异步运行和再同步3、做好系统“黑启动”方案6. 程序源代码（1） formAdmitanceMatrix.mfunction Z1,Z0=formAdmitanceMatrix(nbu

29、s_1,nbus_0,nbranch_1,nbranch_0,branch_1,branch_0)if( nargin3 ) 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

30、_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; endenddisp(正序阻抗);Y1Y2=Y1; %负序等于正序disp(负序阻抗);Y2Y0=zeros(nbus_0);for m=1:nbranch_0 from

31、_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;

32、 endenddisp(零序阻抗);Y0%-计算节点阻抗矩阵-%Z1=inv(Y1);Z2=inv(Y2);Z0=inv(Y0);return;（2） shortcircuit.mfunction 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

33、(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

34、,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*

35、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)

36、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） a2 a 1 a a2 1;%-第一步：计算短路点的序电流，相电流-%计根据故障类型选择不同的计算公式，计算故障点各序电流if Fault_Type=0 %三相短路 I_Fault1=1/Z1(Fault_Node,Fa

37、ult_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 her

38、e Zf=10*50/(1152); 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); %正序阻抗和负序阻抗并联

39、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); I_Fault2=-I_Fault1*Z0(Fault_Node,Fault_Node)/(Z2(Fault_Nod

40、e,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;disp(短路点序电流);I1=I_Fault1I2=I_Fault2I0=I_Fault0disp(短路电流有效值);Iabc=T*I_Fault1;I_Fault2;I_Fault0; %相量Iabc_effective=abs(Iabc) %有效值%-第二步：计算各个节点的序电压，相电压-%write yo

41、ur programme here%故障分量 Ifault_node1=zeros(bus_Num_1,1);Ifault_node2=zeros(bus_Num_1,1);Ifault_node0=zeros(bus_Num_0,1);for m=1:bus_Num_1 if m=Fault_Node Ifault_node1(m,1)=-1*I_Fault1; Ifault_node2(m,1)=-1*I_Fault2; else Ifault_node1(m,1)=0; Ifault_node2(m,1)=0; endendfor m=1:bus_Num_0 if m=Fault_Nod

42、e Ifault_node0(m,1)=-1*I_Fault0; else Ifault_node0(m,1)=0; endendUfault_node1=zeros(bus_Num_1,1);Ufault_node2=zeros(bus_Num_1,1);Ufault_node0=zeros(bus_Num_0,1);Ufault_node1=Z1*Ifault_node1;Ufault_node2=Z2*Ifault_node2;Ufault_node0=Z0*Ifault_node0;%正常分量 Unormal_node1=ones(bus_Num_1,1);Unormal_node2=

43、zeros(bus_Num_1,1);Unormal_node0=zeros(bus_Num_0,1);U_node1=zeros(bus_Num_1,1);U_node2=zeros(bus_Num_1,1);U_node0=zeros(bus_Num_0,1);U_node1=Ufault_node1+Unormal_node1; %正常分量加故障分量U_node2=Ufault_node2+Unormal_node2;U_node0=Ufault_node0+Unormal_node0;disp(各节点序电压);U1=U_node1U2=U_node2U0=U_node0Uabc=;fo

44、r m=1:bus_Num_1 if m=bus_Num_0 tmp=T*U1(m)*exp(j*bus(m).shift);U2(m)*exp(-j*bus(m).shift);U0(m); Uabc=Uabc,tmp; else tmp=T*U1(m)*exp(j*bus(m).shift);U2(m)*exp(-j*bus(m).shift);0; Uabc=Uabc,tmp; endenddisp(各节点相电压有效值);Uabc_effective=abs(Uabc)% -第三步：计算各个支路的序电流，相电流-%write your programme hereI_branch1=ze

45、ros(bus_Num_1,1);I_branch2=zeros(bus_Num_1,1);I_branch0=zeros(bus_Num_0,1);Iabc1=;for m=1:branch_Num_0 from_bus = branch_1(m).from; to_bus=branch_1(m).to; if(from_bus=0) I_branch1(m,1)=(0-U1(to_bus)/branch_1(m).x;%支路正序电流 I_branch2(m,1)=(0-U2(to_bus)/branch_1(m).x;%支路负序电流 elseif(to_bus=0) I_branch1(m

46、,1)=(U1(from_bus)-0)/branch_1(m).x;%支路正序电流 I_branch2(m,1)=(U2(from_bus)-0)/branch_1(m).x;%支路负序电流 else I_branch1(m,1)=(U1(from_bus)-U1(to_bus)/branch_1(m).x;%支路正序电流 I_branch2(m,1)=(U2(from_bus)-U2(to_bus)/branch_1(m).x;%支路负序电流 end from_bus = branch_0(m).from; to_bus=branch_0(m).to; if(from_bus=0) I_b

47、ranch0(m,1)=(0-U0(to_bus)/branch_0(m).x0;%支路零序电流 elseif(to_bus=0) I_branch0(m,1)=(U0(from_bus)-0)/branch_0(m).x0;%支路零序电流 else I_branch0(m,1)=(U0(from_bus)-U0(to_bus)/branch_0(m).x0;%支路零序电流 end %-考虑变压器支路造成的相移-% tmp=T*I_branch1(m,1)*exp(j*branch_1(m).shift);I_branch2(m,1)*exp(-j*branch_1(m).shift);I_branch0(m,1);%由序分量合成相分量 Iabc1=Iabc1,tmp;end;for m=branch_Num_0+1:branch_Num_1 from_bus = branch_1(m).from; I_branch1(m,1)=(1-U1(from_bus)/branch_1(m).x;%支路正序电流 I_branch2(m,1)=(0-U2(from_bus)/branch_1(m).x;%支路负序电流 %-考虑变压器支路造成的相移-% tmp=T*I_branch1(m,1)*exp(j*branch_1(m)