电力系统简单不对称故障的计算机算法(共27页).doc

精选优质文档-倾情为你奉上长 沙 学 院课程设计说明书题目电力系统简单不对称故障的计算机算法系(部)电子与通信工程系专业(班级)电气工程及其自动化02姓名扶泗洲学号指导教师张丹，冯婉起止日期2012.12.22-2012.12.26电力系统分析课程设计任务书系（部）：电子与通信工程系 专业：电气工程及其自动化 指导教师：冯婉课题名称电力系统简单不对称故障的计算机算法设计内容及要求 当电力系统发生简单不对称故障时，用计算机编写程序的方式完成电力网络中短路参数的计算。要求：1、根据所给的电力系统，手算方式完成计算；2、编制计算程序，通过计算机调试，机算方式完成计算；3、提交一份完整的课程设计资料，包括原理、手算过程、程序框图、源程序、符号说明表、算例及计算结果、参考文献、设计总结等。 电力网络具体见电力系统分析P340 16.4图设计工作量1、 形成所给网络的节点导纳矩阵，准备原始数据，并手工计算短路电流；2、 复习系统三相短路的基本原理，建立数学模型；3、 确定合适的数值计算方法；4、 上机编程调试，分析；5、 仿真系统的短路参数，并与手工计算比较。进度安排起止日期（或时间量）设计内容（或预期目标）备注第1天课题介绍，答疑，收集材料第2天学习MATLAB第3天手算短路参数，编写计算机程序第4天编写计算机程序并调试，比较手算与机算结果第5天编写设计说明书教研室意见年 月 日系（部）主管领导意见年 月 日长沙学院课程设计鉴定表姓名扶泗洲学号专业电气工程及其自动化班级2设计题目指导教师指导教师意见：评定成绩： 教师签名： 日期： 答辩小组意见：评定成绩：答辩小组长签名：日期：教研室意见：最终评定等级：教研室主任签名：日期：说明课程设计成绩分“优秀”、“良好”、“中等”、“及格”、“不及格”五等。目录 摘 要在电力系统运行过程中，时常会发生故障，且大多是短路故障。短路通常分为三相短路、单相接地短路、两相短路和两相接地短路。其中三相短路为对称短路，后三者为不对称短路。电力运行经验指出单相接地短路占大多数，因此分析与计算不对称短路具有非常重要意义。分析计算不对称短路方法很多，目前实际最常用的方法是对称分量法。而以对称分量法为核心的计算方法又可有解析法和计算机程序算法等，本论文的主要工作即介绍这两种计算方法。解析法，是将微分方程代数化、暂态分析稳态化、不对称转化为对称并叠加完成不对称故障的分析与计算。计算机程序算法是在形成三个序网的节点导纳矩阵后，对其应用高斯消去法求得故障端点等值阻抗，根据故障类型选用相应公式计算各序电流、电压，进而合成三相电流、电压。 MATLAB的基本数据单位是矩阵，它的指令表达式与数学,工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完相同的事情简捷得多.在新的版本中也加入了对C,FORTRAN,c+ ,JAVA的支持.可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用。 Matlab是一个高度集成的系统，集科学计算、图象处理、声音处理于一体，具有极高的编程效率。近年来，Matlab已经从最初的“矩阵实验室”，渗透到科学与工程计算的多个领域，在自动控制、信号处理、神经网络、模糊逻辑、小波分析等多个方向，都有着广泛的应用。本文根据电力系统路的特点，建立了合理的短路的数学模型，在此基础上，形成电力系统短路电流实用汁算方法编制了对任意一个电力系统在任意点发生短路时短路电流及其分布的通用计算程序从一个侧面展示了计算机应用于电力系统的广阔前景  根据所给的电力系统，编制短路电流计算程序，通过计算机进行调试，最后完成一个切实可行的电力系统计算应用程序。 电力系统短路电流计算是电力系统运行分析，设计计算的重要环节，许多电业设计单位和个人倾注极大精力从事这一工作的研究。由于电力系统结构复杂，随着生产发展，技术进步系统日趋扩大和复杂化，短路电流计算工作量也随之增大，采用计算机辅助计算势在并行。关键词: 单相接地短路，两相短路，两相接地短路，对称分量法，节点导纳矩阵一前言电力系统分析是一门介绍电力系统稳态运行分析、故障分析和暂态过程分析的课程。电力系统分析的基础为电力系统潮流计算、短路故障计算和稳定计算。在电力系统运行过程中，时常会发生故障，其中大多数是短路故障（简称短路）。所谓短路，是指电力系统正常运行情况以外的相与相之间或相与地（或中性线）之间的连接。产生短路的主要原因是电气设备载流部分的相间绝缘或相对地绝缘损坏。此外运行人员在短路检修后未拆除地线就加电压等误操作也会引起短路故障。短路问题是电力技术方面的基本问题之一。在发电厂、变电站以及整个电力系统的设计和运行工作中，都必须事先进行短路计算，以此作为合理选择电气接线、确定限制短路电流措施等的重要依据。为此计算短路时各种运行参量（电流、电压等）是非常必要的。 短路通常分为三相短路、单相接地短路、两相短路和两相接地短路。其中三相短路为对称短路，后三者为不对称短路。电力运行经验指出，单相接地短路占大多数。因此分析与计算不对称短路具有非常重要意义。二.电力系统短路故障的基本知识2.1短路故障的概述在电力系统运行过程中，时常发生故障，其中大多数是短路故障。所谓短路：是指电力系统正常运行情况以外的相与相之间或相与地（或中性线）之间的连接。除中性点外，相与相或相与地之间都是绝缘的。电力系统短路可分为三相短路，单相接地短路。两相短路和两相接地短路等。三相短路的三相回路依旧是对称的，故称为不对称短路。其他的几种短路的三相回路均不对称，故称为不对称短路。电力系统运行经念表明，单相短路占大多数，上述短路均是指在同一地点短路，实际上也可能在不同地点同时发生短路，例如两相在不同地点接地短路。依照短路发生的地点和持续时间不同，它的后果可能使用户的供电情况部分地或全部地发生故障。当在有由多发电厂组成的电力系统发生端来了时，其后果更为严重，由于短路造成电网电压的大幅度下降，可能导致并行运行的发电机失去同步，或者导致电网枢纽点电压崩溃，所有这些可能引起电力系统瓦解而造成大面积的停电事故，这是最危险的后果。产生短路的原因很多主要有如下几个方面（1）原件损坏，例如绝缘材料的自然老化，设计，安装及维护不良所带来的设备缺陷发展成短路，(2)气象条件恶化，例如雷电造成的闪络放电或避雷针动作，架空线路由于大风或导线覆冰引起电杆倒塌等。(3)违规操作，例如运行人员带负荷拉刀闸(4)其他，例如挖沟损伤电缆。在电力系统和电气设备的设计和运行中，短路计算是解决一系列技术问题所不可缺少的基本计算，比如在选择发电厂和电力系统的主接线时为啦比较不同方案接图，进行电力系统暂态稳定计算，研究短路对用户的影响。合理配置各种继电保护和自动装置并正确整定其参数都必须进行短路的计算和分析.三.对称分量法及其在计算中的应用3.1不对称三相量的分解人们在长期的实践中发现，在三相电路中，任意一组不对称的三相相量（电压或电流），可以分解为三组三相对称的相量分量，式（2-1）。在线性电路中，可以用叠加原理对这三组对称分量按照三相电路去解，然后将其结果叠加起来。就是不对称三相电路的解答，这个方法就叫做对称分量法。设，为三相系统中任意一组不对称的三相量，可以分解为三组对称的三序分量如下： （2.1）三相序分量如图2-1所示： 正序分量 负序分量 零序分量 图2-1 三序分量正序分量：三相的正序分量大小相等，彼此相位相差，与系统正常对称运行对称运行方式下的相序相同，达到最大值，在电机内部产生正转磁场，这就是正序分量。此正序分量为一平衡的三相系统，因此有： 负序分量：三相的负序分量大小相等，彼此相位相差，与系统正常对称运行对称运行方式下的相序相反，达到最大值，在电机内部产生反转磁场，这就是正序分量。此正序分量为一平衡的三相系统，因此有： 零序分量：三相的零序分量大小相等，相位相同，三相的零序分量同时达到最大值，在电机内部产生漏磁，其合成磁场为零，这就是零序分量。 如果以A相为基准相，各序分量有如下关系：正序分量： 负序分量： 零序分量： 其中： ； ；于是有： （2-2） （2-3）则有： 其逆关系式为： （2-4） 这样根据式（2-3）可以把三组三相对称向量合成三个不对称向量，而根据式（2-4）可以把三个不对称向量分解成三组对称量。四简单不对称短路的分析与计算电力系统发生不对称故障时，无论是单相接地短路，还是两相短路，两相短路接地，只是短路点的电压，电流出现不对称，利用对称分量法将不对称的电流电压分解为三组对称的序分量，由于每一序系统中三相对称，则在选好一相为基准后，每一序只需要计算一相即可，用对称分量法计算电力系统的不对称故障。其大概步骤如下：（1）计算电力系统各个 原件的序阻抗；（2）制定电力系统的各序网络；（3）由各序网络和故障列出对应方程；（4）从联立方程组解出故障点电流和电压的各序分量，将相对应的各序分量相加，以求得故障点的各相电流和电压；（5）计算各序电流和各序电压在网络中的分布，进而求得各指定支路的各相电流和指定节点的各相电压。4.1单相（a相）接地短路单相接地短路时的系统接线图如图3-1所示，故障处的三相的边界条件： 图3-1 单相接地短路时的等值接线图 用对称分量表示为： （3-1）即有： （3-2）根据单相接地短路时的边界条件（式3-1，式3-2）连接复合网，如图3-2所示。由复合网图可以写出各序分量：图3-2单相接地短路时的复合序网 （3-3）于是可以用对称分量法得到短路点的各相电流电压： （3-4）短路点电流，电压的相量图如图3-3所示。这里按纯电感性电路画的，电流滞后电压，若不是纯电感电路，则电流与电压角度由的阻抗角决定，一般小于。在相量图中，将每相的序分量相加，得各相电流电压的大小和相位。 图3-3 单相接地时短路处的电流，电压的相量图4.2 两相（b,c相）短路 两相短路的系统接线如图3-4所示，在k点发生b，c两相短路。短路点的边界条件：图3-4两相短路等值接线图 用对称分量表示为： （3-5）于是有： (3-6)由式（3-6）可知，故障点不与大地相连，零序电流无通路，因此无零序网络。复合网络是正负序网并联后的网络。如图3-5所示：图3-5两相短路的复合序网图从复合序网中可以直接求出电流，电压的各序分量： （3-7）由对称分量法可求得短路点各相电流和电压为： （3-8）短路点电压和电流的相量图。如图3-6所示，这里任然是纯电感电路。电流滞后电压 图3-6两相接地短路点电流和电压相量图4.3两相（b相和c相）短路 接地 两相短路接地时系统接线图如图3-7所示，短路点的边界条件为：图3-7两相短路接地的等值接线图 用对称分量表示为： （3-9） 由（3-9）可以画出两相短路接地的复合序网图是三个序网并联，如图3-8所示，根据复合序网可求出电流，电压各序分量： 图3-8两相短路接地的复合序网图 （3-10）用对称分量法合成各相电流电压为： （ 3-11）短路点流入地中的电流为： （3-12）短路点电压，电流的相量图，如图3-9所示。这里任然是纯电感电路。电流滞后电压 图3-9两相接地短路处的电压电流相量图五简单不对称短路的计算机算法5.1简单故障的计算程序原理 1如果要求准确计算故障前的运行情况，则需要进行潮流计算，在近似使用计算中，对短路故障可假设节点UI0I.均为1IGI2这里采用形成节点导纳矩阵的方法。发电机的正序电抗用Xdn,可计算故障后瞬时量。发电机的负序的电抗近似等于Xdn.当计算中不计负何影响时，在正，负序网络中不接负何阻值。如果记及负何影响时，负何的正序阻抗可通过其额定功率和电压计算。负何阻抗很难却确定，一般取X2=0.35以负何额定功率为基准。负何的中性点一般不接地，零序无通路3形成三个序网的节点导纳矩阵后，利用公式3-19和高斯消去法可求得 ,即为Z1f Znf .对于短路故障，只需令 =1 （其余节电电流均为零），分别应用高斯消去法求解一次所得电压，即为三序网和f点有关节点阻抗（4）根据不同短路故障，可分别应用 表3.1中 相应公式计算故障处各序电流，电压，进而合成得三相电流，电压。（5）计算网络中任一点电压时，负序和零电压只须计算由故障点电流引起的电压。对于正序则还需要加上正常运行时的电压表4-1不对出故障处的各序电流电压计算公式。故障类型故障端口各序电流公式故障端口各序电压公式单相接地短路两相短路两相接地短路六设计数据及记录6.1设计要求 当电力系统发生简单不对称故障时，用计算机编写程序的方式完成电力网络中短路参数的计算。要求：1、根据所给的电力系统，手算方式完成计算；2、编制计算程序，通过计算机调试，机算方式完成计算；6.2.手算数据一、解：取=100MVA，（1）各元件序电抗标幺值的计算正序电抗：发电机 变压器T1和T2 线路L1 线路L2 正序总电抗 正序网络化简过程如图a所示。图a 正序网络及简化过程负序电抗：因发电机和变压器的正负序电抗相等，所以正序网络化简过程如图b所示。图b 负序网络及简化过程零序电抗：变压器零序电抗与正序电抗相等，即T1与T2的正序和零序电抗相等；线路L1 线路L2 变压器T3 零序总电抗 零序网络化简过程如图c所示。图c 零序网络及化简过程（2）单相短路电流计算附加电抗：计算电抗：查运算曲线得时，次暂态短路电流 所以KA 因为， 所以KA（3）两相短路电流计算附加电抗：计算电抗：查运算曲线得时，次暂态短路电流 所以KA 因为， 所以KA（4）两相接地短路电流计算附加电抗：计算电抗：查运算曲线得时，次暂态短路电流 所以KA 因为， 所以KA七.总结电力系统发生不对称短路后，由于短路点对地故障支路的不对称，使得整个网络电流电压三相不对称.解决不对称短路的问题核心是对称分量法。根据对称分量法采取的具体方法之一是解析法，即把该网络分解为正，负，零序三个对称序网，这三组对称序分量可分别按对称的三相电路分解。然后将其结果叠加起来。求解不对称短路。首先应该计算各原件的序参数和画出等值电路。然后制定各序网络。根据不同的故障类型，确定出以相分量表示的边界条件，进而列出以序分量表示的边界条件，按边界条件将三个序网联合成复合网，由复合网求出故障处各序电流和电压，进而合成三相电流电压。二是计算机程序法。通过计算机形成三个序网的节点导纳矩阵，然后利用高斯消去法通过相应公式对他们进行数据运算，即可求得故障端点的等值阻抗。最后根据故障类型选取相关公式计算故障处各序电流电压，进而合成三相电流电压。通过上面两种方法放入计算进行对比可以得出，计算机程序法比较解析法具有计算过程及结果更精确，通用性更强的优点。作为一名大三的电气类学生，电力系统这门专业课对于我们学习的并不算好，主要是投入的时间不够，以至于拿到课题之后不知道应该怎么去做，我们这组通过分工协作开始做这次设计，开始的确很难，让我们无从下手。但是通过查阅资料书本，咨询老师同学，渐渐地有了方向和目标，然后进行了手算，并得出了最终结果。但是最后调试软件时候，虽然查阅了许多资料也询问了同学，但是还是没有得出最后结果，是本次课程设计得最大遗憾。但我相信经过这次电力系统分析课程设计，大家会对它进一步加深学习，这对我们以后找工作也会有一定的好处，总之，要学好它，必须下功夫，将来才能更好地应用它。任何事情在实施之前均要有详细的计划，无论是做查资料，还是编写软件，而且在实施之前都要对所做的事做到心中有数，在调试软件的过程中遇到了很多困难，但是通过查找资料和做实物有学到了很多书本上学不到的东西，既丰富了自己的知识，又加强了自己的动手能力。在做本次课程设计的过程中，Matlab程序设计语言制在整个设计过程中都用到了，我又感觉到了Matlab功能的强大，尤其是它的操作界面，给人一种强大的震撼力，希望以后用到的时候会进一步学习它，由于水平有限，错误在所难免，希望老师批评指正。再次感谢冯婉老师对我们这组不辞辛苦的辅导，顺便在此祝老师新年快乐！参考文献：1 祝淑萍 等.电力系统分析课程设计综合实验M.中国电力出版社;20072 陈珩.电力系统稳态分析M.中国电力出版社;2007:220-2503 苏金明，王永利.MATLAB实用指南M.电子工业出版社;20044 孟祥萍，高嬿.电力系统分析M.高等教育出版社; 2010:201-221附录短路电流计算程序：%本程序是对简单不对称故障的计算NF=input('请输入短路点的数目：NF=');for i1=1:NFclear;n1=input('请输入正序节点的数目：n1=');n2=input('请输入负序节点的数目：n2=');n0=input('请输入零序节点的数目：n0=');nl1=input('请输入正序支路数：nl1=');nl2=input('请输入负序支路数：nl2=');nl0=input('请输入零序支路数：nl0=');Lf=input('请输入短路的类型：Lf=');f=input('请输入短路节点号：f=');If=input('请输入纵横故障标志：If=');zf=input('请输入接地阻抗：zf=');zg=input('请输入接地阻抗：zg=');B1=input('请输入正序支路参数：B1=');B2=input('请输入负序支路参数：B2=');B0=input('请输入零序支路参数：B0=');A1=input('请输入不对称系数：A1=');Vcs=input('请输入由各节点的初始电压标幺值形成的列矩阵：Vcs=');Z1=zeros(n1);Z2=zeros(n2);Z0=zeros(n0);Y1=zeros(n1);Y2=zeros(n2);Y0=zeros(n0);V1=zeros(n1,1);V2=zeros(n2,1);V0=zeros(n0,1);I1=zeros(nl1,1);I2=zeros(nl2,1);I0=zeros(nl0,1);for m1=1:3 m=0; if m1=1 n1=nl1;B=B1;Z=Z1;n=n1; elseif m1=2 n1=nl2;B=B2;Z=Z2;n=n2; elseif m1=3 n1=nl0;B=B0;Z=Z0;n=n0; end m=0;Z=zeros(n);V=zeros(n);I=zeros(n1);for k1=1:n1 p=B(k1,1);q=B(k1,2); if B(k1,6)=0 k=1./B(k1,5); else k=B(k1,5); end if p=0 if q>m %追加接地树支 Z(q,q)=B(k1,3);m=m+1; else %追加接地连支 for i=1:m, Z(i,m+1)=-Z(i,q);Z(m+1,i)=-Z(q,i); end Z(m+1,m+1)=Z(q,q)+B(k1,3); for i=1:m for j=1:m Z(i,j)=Z(i,j)-Z(i,m+1)*Z(m+1,j)./Z(m+1,m+1); end Z(i,m+1)=0; end for i=1:m+1 Z(m+1,i)=0; end end else if q>m %追加不接地树支 for i=1:m Z(i,q)=Z(i,p)*k;Z(q,i)=Z(p,i)*k; end Z(q,q)=k2*Z(p,p)+k2*B(k1,3) m=m+1; else for i=1:m %追加不接地连支 Z(i,m+1)=k*Z(i,p)-Z(i,q); Z(m+1,i)=k*Z(p,i)-Z(q,i); end Z(m+1,m+1)=k2*Z(p,p)+Z(q,q)-2*k*Z(p,q)+k2*B(k1,3); for i=1:m for j=1:m Z(i,j)=Z(i,j)-Z(i,m+1)*Z(m+1,j)./Z(m+1,m+1); end Z(i,m+1)=0; end for i=1:m+1 Z(m+1,i)=0; end end endend if m1=1 Z1=Z; elseif m1=2 Z2=Z; elseif m1=3 Z0=Z; endend %以上程序求出正、负、零序对应的阻抗阵if If=0 %求出横向故障的阻抗阵 Z1(f,f)=Z1(f,f); Z2(f,f)=Z2(f,f); Z0(f,f)=Z0(f,f);else %求出纵向故障的阻抗阵 Z1(f,f)=Z1(f,f)+Z1(1,1)-2*Z1(f,1); Z2(f,f)=Z2(f,f)+Z2(1,1)-2*Z2(f,1); Z0(f,f)=Z0(f,f)+Z0(1,1)-2*Z0(f,1);endif Lf=1 %求出正负零序短路节点的电流标幺值 Z4=Z2(f,f)+Z0(f,f)+3*zf;k2=1;k0=1;I1(f)=Vcs(f)./(Z1(f,f)+Z4);I2(f)=k2*I1(f); I0(f)=k0*I1(f);a1=I1(f);a2=I2(f);a0=I0(f);elseif Lf=2 Z4=Z2(f,f)+2*zf;k2=-1;k0=0; I1(f)=Vcs(f)./(Z1(f,f)+Z4);I2(f)=k2*I1(f);I0(f)=k0*I1(f); a1=I1(f);a2=I2(f);a0=I0(f);elseif Lf=3 Z4=zf+(Z2(f,f)+zf)*(Z0(f,f)+zf+3*zg)./(Z2(f,f)+Z0(f,f)+2*zf+3*zg); k2=-(Z0(f,f)+zf+3*zg)./(Z2(f,f)+Z0(f,f)+2*zf+3*zg); k0=-(Z2(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf+3*zg); I1(f)=Vcs(f)./(Z1(f,f)+Z4); I2(f)=k2*I1(f); I0(f)=k0*I1(f); a1=I1(f); a2=I2(f); a0=I0(f);elseif Lf=4 Z4=zf+(Z2(f,f)+zf)*(Z0(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf); k2=-(Z0(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf); k0=-(Z2(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf); I1(f)=Vcs(f)./(Z1(f,f)+Z4); I2(f)=k2*I1(f); I0(f)=k0*I1(f); a1=I1(f); a2=I2(f); a0=I0(f);elseif Lf=5 Z4=Z2(f,f)+Z0(f,f)+3*zf; k2=1; k0=1; I1(f)=Vcs(f)./(Z1(f,f)+Z4); I2(f)=k2*I1(f); I0(f)=k0*I1(f); a1=I1(f); a2=I2(f); a0=I0(f);enddisp('故障点处的正序电流标幺值为：');disp(a1);disp('故障点处的负序电流标幺值为：');disp(a2);disp('故障点处的零序电流标幺值为：');disp(a0);for i=1:n1 %求各节点各序的电压标幺值 V1(i)=Vcs(i)-I1(f)*Z1(i,f);endfor i=1:n2 V2(i)=-I2(f)*Z2(i,f);endfor i=1:n0 V0(i)=-I0(f)*z0(i,f);enddisp('各节点的正序电压标幺值为（节点号从小到大排）：');disp(V1);disp('各节点的负序电压标幺值为（节点号从小到大排）：');disp(V2);disp('各节点的零序电压标幺值为（节点号从小到大排）：');disp(V0);for i=1:3 %求各支路各序的电流标幺值 if i=1 B=B1;V=V1;n1=nl1; elseif i=2 B=B2;V=V2;n1=nl2; elseif i=3 B=B0;V=V0;n1=nl0; end for j=1:n1 if B(j,6)=0 k=B(j,5); else k=1./B(j,5); end p=B(j,1);q=B(j,2); if p=0 e=0;b=B(j,3); I(j)=(e-V(q)./k)./b; else I(j)=(V(p)-V(q)./k)./B(j,3); end end if i=1 I1=I; elseif i=2 I2=I; elseif i=3 I0=I; endenddisp('各支路的正序电流标幺值为（顺序同您输入B时一样）：');for i=1:n1 disp(I1(i);enddisp('各支路的负序电流标幺值为（顺序同您输入B时一样）：');for i=1:n1 disp(I2(i);enddisp('各支路的零序电流标幺值为（顺序同您输入B时一样）：');for i=1:n1 disp(I0(i);endfor i=1:3 if i=1 s0=a0;s1=a1;s2=a2; elseif i=2 s0=V0;s1=V1;s2=V2; else s0=I0;s1=I1;s2=I2; end for j=1:3 CA=s0*A1(j,1)+s1*A1(j,2)+s2*A1(j,3); if j=1 A=CA; elseif j=2 B=CA; elseif j=3 C=CA; end end if i=1 Iadl=A;Ibdl=B;Icdl=C; elseif i=2 Vajy=A;Vbjy=B;Vcjy=C; else Iazl=A;Ibzl=B;Iczl=C; endenddisp('故障点处A相电流Iadl为：');disp(Iadl);disp('故障点处B相电流Ibdl为：');disp(Ibdl);disp('故障点处C相电流Icdl为：');disp(Icdl);disp('各节点处A相电压Vajy分别为（节点号从小到大排）：');disp(Vajy);disp('各节点处B相电压Vbjy分别为（节点号从小到大排）：');disp(Vbjy);disp('各节点处C相电压Vcjy分别为（节点号从小到大排）：');disp(Vcjy);disp('各支路A相电流Iazl分别为（顺序同您输入B时一样）：');for i=1:n1 disp(Iazl(i);enddisp('各支路B相电流Ibzl分别为（顺序同您输入B时一样）：');for i=1:n1 disp(Ibzl(i);enddisp('各支路C相电流Iczl分别为（顺序同您输入B时一样）：');for i=1:n1 disp(Iczl(i);endend6.5程序运算结果:单相短路电流为：    2.1941两相短路电流为：    2.0849两相短路接地电流为： 2.3157 专心-专注-专业