电力系统不对称故障分析与计算

1、复杂网络N-R法潮流分析与计算的设计？摘 要随着电力事业的快速发展，电力电子新技术得到了广泛应用；出于技术、经济等方面的考虑，500kV及以上的超高压输电线路普遍不换位，再加上大量非线性元件的应用，电力系统的不对称问题日益严重。因此电力系统不对称故障分析与计算显得尤为重要。 基于对称分量法的基本理论，对称分量法采取的具体方法之一是解析法，即把该网络分解为正，负，零序三个对称序网，这三组对称序分量可分别按对称的三相电路分解。计算机程序法。通过计算机形成三个序网的节点导纳矩阵，然后利用高斯消去法通过相应公式对他们进行数据运算，即可求得故障端点的等值阻抗。最后根据故障类型选取相关公式计算故障处各序电

2、流，电压，进而合成三相电流电压。进行了参数不对称电网故障计算方法的研究。通过引计算机算法，系统介绍电网参数不对称的计算机算法方法。根据断相故障和短路故障的特点，通过在故障点引入计算机算法，给出了各种断相故障和短路故障的仿真计算。此方法以将故障电网分为对称网络和不网络两部分，在程序法则下建立起不对称电网故障计算统一模型，根据线性电路的基本理论，并借助于相序参数变换技术完成故障计算。 关键词: 参数不对称; 电网; 故障计算目 录前 言11.电力系统短路故障的基本知识21.1 短路故障的概述21.2 标幺制42 对称分量法在不对称短路计算中的应用22.1 不对称三相量的分解32.2 对称分量法在不

3、对称短路计算中的应用43 简单不对称短路的分析与计算43.1 单相（a相）接地短路3.2 两相（b，c相）短路 73.3 两相（b相和c相）短路 接地74 简单不对称短路的分析与计算计算机计算程序法84.1 简单故障的计算程序原理94.2 网络节点方程的形成105 电力系统不对称短路计算实例115.1 单相接地短路和两相短路不对称故障分析与计算115.2 两种计算方法的对比18结语19参考资料19附录：不对称短路电流计算程序20前 言电力系统的安全、稳定、经济运行无疑是历代电力工作者所致力追求的，但是从电力系统建立之初至今电力系统就一直伴随着故障的发生而且电力系统的故障类型多样。所以电力系统故

4、障分析计算方法年来一直是学术研究的热点，其为预防及消除电力系统的故障准备，必要的理论知识。引起电网不对称的原因，超高压架空输电线路不换位.变压器结构不对称.交直流变换器的存在，系统负荷不平衡，系统中存在非线性元件。建立用于故障分析的网络的数学模型是计算机编程计算的基础选取合适的数学模型直接关系到故障分析计算的优劣。当前一般采用节点导纳矩阵或节点阻抗矩阵作为原网数学模型节点导纳矩阵描述的是网络的短路参数，只包含网络的局部信息，因而具有较好的稀疏性。因此，节点导纳矩阵主要运用于解题规模为主要矛盾的情况,节点阻抗矩阵要运用于进行大量故障计算的情况。随着计算机技术的快速发展，特别是计算机内存容量的扩大

5、和计算速度的提高，节点阻抗矩阵的解题规模和解题速度都有了很大的提高，一般可以满足大多数电网故障分析计算的需要。电网中的故障可分为两大类：简单故障和复杂故障。复杂故障一般可以由两种或两种以上的简单故障组合而成，简单故障又分为对称型故障和不对称型故障；而不对称故障又可以分为短路故障(横向故障)和断路故障(纵向故障)。故障计算方法总体上可以划分为两种：一种解析法，另一种计算机程序法。1.电力系统短路故障的基本知识1.1短路故障的概述 在电力系统运行过程中，时常发生故障，其中大多数是短路故障。所谓短路：是指电力系统正常运行情况以外的相与相之间或相与地（或中性线）之间的连接。除中性点外，相与相或相与地之

6、间都是绝缘的。电力系统短路可分为三相短路，单相接地短路。两相短路和两相接地短路等。三相短路的三相回路依旧是对称的，故称为不对称短路。其他的几种短路的三相回路均不对称，故称为不对称短路。电力系统运行经念表明，单相短路占大多数，上述短路均是指在同一地点短路，实际上也可能在不同地点同时发生短路，例如两相在不同地点接地短路。依照短路发生的地点和持续时间不同，它的后果可能使用户的供电情况部分地或全部地发生故障。当在有由多发电厂组成的电力系统发生端来了时，其后果更为严重，由于短路造成电网电压的大幅度下降，可能导致并行运行的发电机失去同步，或者导致电网枢纽点电压崩溃，所有这些可能引起电力系统瓦解而造成大面积

7、的停电事故，这是最危险的后果。产生短路的原因很多主要有如下几个方面（1）原件损坏，例如绝缘材料的自然老化，设计，安装及维护不良所带来的设备缺陷发展成短路，(2)气象条件恶化，例如雷电造成的闪络放电或避雷针动作，架空线路由于大风或导线覆冰引起电杆倒塌等。(3)违规操作，例如运行人员带负荷拉刀闸(4)其他，例如挖沟损伤电缆。在电力系统和电气设备的设计和运行中，短路计算是解决一系列技术问题所不可缺少的基本计算，比如在选择发电厂和电力系统的主接线时为啦比较不同方案接图，进行电力系统暂态稳定计算，研究短路对用户的影响。合理配置各种继电保护和自动装置并正确整定其参数都必须进行短路的计算和分析.1.2 标幺

8、制 在短路计算中，各电气量如电流.电压.电阻.电抗.功率等数值可以用名值表示，也可以用标幺值表示，为了计算方便，通常在1KV以下的低压系统用有名值，而高压等级中有多个电压等级，存在电抗换算问题，宜采用标幺值。1.21 标幺制概念所谓标幺制，就是把各个物理量均用标幺值来表示的一种相对单位制。某一物理量的标幺值,等于它的实际值A与所选定 的基准值的比值，即： （1-1） 在进行标幺值计算时，首先选定基准值。基准值原则上可以任意选定，但因物理量之间有内在的必然联系，所以并非所有的基准值都可以任意选取。在短路计算中经常用到的四个物理量容量S，电压U，电流I，和电抗X。通常选定基准容量和基准电压则基准电

9、流和基准电抗分别为： （1-2） （1-3）为了计算方便，常取基准容量Sd=100MV.A，基准电压用各级线路的平均额定电压，即。所谓线路平均额定电压，是指线路始端最大额定电压与线路末端最小额定电压的平均值，如表1-1所示。表1-1 线路的额定电压与平均额定电压额定电压0.220.3836103560110220330平均额定电压0.230.43.156.310.53763115230345在产品样本中，电力系统中各电气设备如发电机，变压器，电抗器等所给出的标幺值，都是以其本身额定值为基准的标幺值或百分值。由于各电气设备的额定值往往不尽相同，基准值不同的标幺值是不能直接进行运算的，因此，必须把

10、不同基准值的标幺值换算成统一基准值的标幺值。 电力系统各原件电抗标幺值的计算（1）发电机电抗标幺值： （1-4）式中-发电机电抗百分数，有发电机名牌参数的 -已设定的基准容量（基准功率），MV.A ； -发电机额定容量。（2）负载电抗标幺值： （1-5） -原件所在网络的电压标幺值；-负载容量标幺值；-负载无功率标幺值 （3） 变压器电抗标幺值： （1-6）在变压器中电抗，即忽略，因此在变压器中阻抗主要是指电抗，由变压器电抗有名值推出变压器电抗标幺值为： （1-7）式中 -变压器阻抗百分数；-基准容量，MV.A; -变压器铭牌参数给定额定容量,MV.A; -基准电压取平均电压 。4 线路电抗标

11、幺值： （1-8）式中： -线路单位长度电抗；L-线路单位长度，KM。-基准容量.MV.A -线路额定平均电压。基准电压输电线路的等值电路有四个参数，一般电抗故。由于不做特殊说明，故电导.电纳一般不计，故而只求电抗标幺值。2. 对称分量法在不对称短路计算中的应用2.1不对称三相量的分解 人们在长期的实践中发现，在三相电路中，任意一组不对称的三相相量（电压或电流），可以分解为三组三相对称的相量分量，式（2-1）。在线性电路中，可以用叠加原理对这三组对称分量按照三相电路去解，然后将其结果叠加起来。就是不对称三相电路的解答，这个方法就叫做对称分量法。设，为三相系统中任意一组不对称的三相量，可以分解为

12、三组对称的三序分量如下： （2-1） 三相序分量如图2-1所示： 正序分量 负序分量 零序分量 图2-1 三序分量正序分量：三相的正序分量大小相等，彼此相位相差，与系统正常对称运行对称运行方式下的相序相同，达到最大值，在电机内部产生正转磁场，这就是正序分量。此正序分量为一平衡的三相系统，因此有： 负序分量：三相的负序分量大小相等，彼此相位相差，与系统正常对称运行对称运行方式下的相序相反，达到最大值，在电机内部产生反转磁场，这就是正序分量。此正序分量为一平衡的三相系统，因此有： 零序分量：三相的零序分量大小相等，相位相同，三相的零序分量同时达到最大值，在电机内部产生漏磁，其合成磁场为零，这就是零

13、序分量。 如果以A相为基准相，各序分量有如下关系：正序分量： 负序分量： 零序分量： 其中： ； ；于是有： （2-2） （2-3）则有： 其逆关系式为： （2-4） 这样根据式（2-3）可以把三组三相对称向量合成三个不对称向量，而根据式（2-4）可以把三个不对称向量分解成三组对称量。 2.2 对称分量法在不对称短路计算中的应用电力系统的正常运行一般是对称的，它的三相电路的参数相同，各相的电流，电压对称，这就是说只有正序分量存在。当电力系统的某一点发生不对称故障时，三相电路的对称条件受到破坏，三相对称电路就成为不对称的了。此时，可用对称分量法，将实际的故障系统变成三个互相独立的序分量系统，而每

14、个序分量系统本身又是三相对称的，从而就可以用进行电路计算了。 图2-2 简单系统单相接地故障图 如图2-2所示的简单系统发生单相接地短路故障。应用对称分量法，可绘出三序网图（三序等值电路图），如图2-3所示为最简化的三序网图，三序网的参数可分为正序，负序，零序。图中分为正序阻抗，负序阻抗，零序阻抗。图2-3简化三序网图列出电压方程： （2-5） 由此可见，应用对称分量法进行不对称故障计算时，其关键问题是先求出各序网络的等效电抗（即要求出系统中各主要原件（发电机，变压器，线路等）的各序电抗值），然后根据短路的类型，边界条件，把正，负，零序网连接成串，并联的形式，从而可求出电流，电压的各序分量，再

15、应用对称分量法进而可求出各相电流和电压等.3 简单不对称短路的分析与计算（解析法）电力系统发生不对称故障时，无论是单相接地短路，还是两相短路，两相短路接地，只是短路点的电压，电流出现不对称，利用对称分量法将不对称的电流电压分解为三组对称的序分量，由于每一序系统中三相对称，则在选好一相为基准后，每一序只需要计算一相即可，用对称分量法计算电力系统的不对称故障。其大概步骤如下：（1）计算电力系统各个 原件的序阻抗；（2）制定电力系统的各序网络；（3）由各序网络和故障列出对应方程；（4）从联立方程组解出故障点电流和电压的各序分量，将相对应的各序分量相加，以求得故障点的各相电流和电压；（5）计算各序电流

16、和各序电压在网络中的分布，进而求得各指定支路的各相电流和指定节点的各相电压。【4】3.1单相（a相）接地短路单相接地短路时的系统接线图如图3-1所示，故障处的三相的边界条件： 图3-1 单相接地短路时的等值接线图 用对称分量表示为： （3-1）即有： （3-2）根据单相接地短路时的边界条件（式3-1，式3-2）连接复合网，如图3-2所示。由复合网图可以写出各序分量：图3-2单相接地短路时的复合序网 （3-3）于是可以用对称分量法得到短路点的各相电流电压： （3-4）短路点电流，电压的相量图如图3-3所示。这里按纯电感性电路画的，电流滞后电压，若不是纯电感电路，则电流与电压角度由的阻抗角决定，一

17、般小于。在相量图中，将每相的序分量相加，得各相电流电压的大小和相位。图3-3 单相接地时短路处的电流，电压的相量图3.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所示

18、，这里任然是纯电感电路。电流滞后电压图3-6两相接地短路点电流和电压相量图3.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两相接地短路处的电压电流相量图4

19、 简单不对称短路的分析与计算计算机计算程序法4.1简单故障的计算程序原理框图 1如果要求准确计算故障前的运行情况，则需要进行潮流计算，在近似使用计算中，对短路故障可假设节点UI0I.均为1IGI2这里采用形成节点导纳矩阵的方法。发电机的正序电抗用Xdn,可计算故障后瞬时量。发电机的负序的电抗近似等于Xdn.当计算中不计负何影响时，在正，负序网络中不接负何阻值。如果记及负何影响时，负何的正序阻抗可通过其额定功率和电压计算。负何阻抗很难却确定，一般取X2=0.35以负何额定功率为基准。负何的中性点一般不接地，零序无通路3形成三个序网的节点导纳矩阵后，利用公式3-19和高斯消去法可求得 ,即为Z1f

20、 Znf .对于短路故障，只需令 =1 （其余节电电流均为零），分别应用高斯消去法求解一次所得电压，即为三序网和f点有关节点阻抗（4）根据不同短路故障，可分别应用 表3.1中 相应公式计算故障处各序电流，电压，进而合成得三相电流，电压。（5）计算网络中任一点电压时，负序和零电压只须计算由故障点电流引起的电压。对于正序则还需要加上正常运行时的电压表4-1不对出故障处的各序电流电压计算公式。故障类型故障端口各序电流公式故障端口各序电压公式单相接地短路两相短路两相接地短路下图4-1即为不对称短路计算原理图： 输入数据形成网络节点导纳矩阵选择短路点短路点注入单位电流时应用高斯消去法解各点电压即 应用式

21、表（4-1）中各种故障对应电流公式计算应用节点电压公式计算任一节点电压应用节点电压公式计算任一节点电压结束图4-1 不对称短路计算程序框图4.2网络节点方程的形成计算机编程计算中，考虑了对地电容和变压器的实际变比7导纳计算公式如下：对角线元素： （ jI, Ui=1,Uj=0,ij） （4-1）非对角线元素 ： （Ui=1,Uj=0,ij） （4-2）式中i， j-第i , j节点 Ii， Ij ，Ui ，Uj-第i,j 节点电流电压变压器支路的导纳变化 （高压侧） （4-3） （4-4） 4） 网络节点方程的形成矩阵表达式 网络节点方程： （4-5） （4-6）5 电力系统不对称短路计算实例

22、5.1 单相接地短路和两相短路不对称故障分析与计算 已知量如图5-1中所示数据，现求取节点K1单相接地短路时和节点K2两相短路时的短路电流，冲击电流以及短路容量。图5-1系统电路图 （一）：解析法分析与计算单相接地短路和两相短路时情况下：短路电流，冲击电流以及短路容量：取基准容量： 100MV.A; 基准电压 ： （kv）等值电路图如图5-2：图5-2等值电路图各原件电抗标幺值计算 ： 取： 100MV.A, （1）发电机。发电机G1电机电抗标幺值由式（1-4）得： 正序电抗标幺值等于负序电抗标幺值等于零序电抗即： =发电机G2的正序电抗标幺值由式（1-4）得到 正序电抗标幺值等于负序电抗标幺

23、值等于负序电抗即：=0.416 .对下标的说明：各原件下标括号外的数字表示为第几个原件，括号内的数字1，2，0分别为正序，负序，零序。例如，X第一个原件（负序） （2）负载。负载电抗标幺值由式（1-5）得;负载（负何）正序电抗标幺值;负序电抗标幺值;零序电抗标幺值. （3）变压器。由式（1-6）得：T1变压器正序电抗标幺值X2: ; 正序电抗标幺值等于负序电抗标幺值等于零序电抗标幺值。即 。变压器正序电抗标幺值： ;即（每台） 。变压器正序电抗标幺值: ; .（4）线路。由式（1-8）得：L1线路的正序电抗标幺值X3: ; ; ；L2线路的正序电抗标幺值X6: ; ; L3 线路的正序电抗标幺

24、值X8: ; ; . 2不对称短路电流的计算：各序等值电路图，正序等值电路图如图5-3所示图5-3 正序等值电路图短路点K1正序标幺值为：）/）=0.221短路点K2正序标幺值为： ）/ ）=0.339（1） 负序等值电路图：负载的负序电抗标幺值；其他原件的负序电抗标幺值等于正序标幺值，如图5-4所示：图 5-4负序等值电路图短路点K1负序标幺值为：短路点K2负序标幺值为：=0.313零序等值电路图：负载的零序电阻抗标幺值；线路的零序电抗标幺值等于分别为：; ; .其他原件的零序阻抗标幺值等于正序标幺值，如图5-5所示图5-5零序等值电路图短路点K1零序标幺值为:短路点K1零序标幺值为:K1点

25、单相接地短路；K2点两相短路： ；（二）：计算机程序法析与计算单相接地短路和两相短路时情况下：短路电流，冲击电流以及短路容量：1. 由式, 得到：；2. 正序导纳标幺值Y矩阵Y矩阵的形成：正序电抗标幺值和对地电容标幺值的等值电路图如图5-6所示：图5-6正序电抗标幺值和对地电容标幺值的等值电路图正序导纳矩阵如下：对地导纳式如下（5-1）： 其中，是图5-6节点13得对地导纳之半， ，分别为节点35，36的导纳之半【9】.Y矩阵的元素为：; ； ；; （5-1）3 负序阻抗Y矩阵图5-7负序阻抗标幺值和对地电纳标幺值的等值电路图 同理的得：如图5-7负序阻抗标幺值和对地电纳标幺值的等值电路图，及

26、负序导纳矩阵如下式如下（5-2）： （5-2）3. 零序列Y矩阵：同理零序列电抗标幺值的等值电路图如图5-8所示：图5-8所示零序列电抗标幺值的等值电路图及 负序导纳矩阵如下如下式（5-3）： （5-3）4 求不对称短路电流及相关参数：由式3-19及3-21节相关知识，运用高斯消去程序【10】解得：K1点单相接地短路时的转移阻抗为如下（5-4）： （5-4）K2点两相短路时的转移阻如下式（5-5）： （5-5）由式表（3-1）相关公式得：K1点单相接地短路；;K2点两相短路:；;； 程序流程图如下： 输入数据形成网络节点导纳矩阵选择短路点短路点注入单位电流时应用高斯消去法解各点电压即应用式表（

27、4.1）中各种故障对应电流公式计算IF. 应用公式 计算任意点节点电压应用公式，计算任一点节点电流应用公式1-9，1-13计算次暂态电流和容量结束5-10不对称短路计算程序框图5.2两种计算方法的对比下面将解析法与计算机程序法进行简单比较。将例一和例二结果如表5-1下：表5-1解析法与计算机程序法数据比较项目短路类型故障处短路电流短路冲击电流（KA）短路容量（MV.A） 1 2 1 2 1 2单相接地短路5.8035.8687.42827.5047508.3568.3两相短路2.6572.6483.4023.389265.7264.81-代表解析法结果；2-代表计算机程序法结果。由例一和例二计

28、算过程以及上面的表格数据及比较可得：解析法误差大，计算过程数据大多需要进行四舍五入。计算机程序法中，利用Y矩阵计算时考虑对地电容，变压器实际变比误差小，通过变换式（3-19）中的电流的取值可以将网络各个节点转移阻抗求出，使分析其他点的短路故障提供了更容易更直观的参数值，因此计算机程序法有通用性强的特点。 结束语电力系统发生不对称短路后，由于短路点对地故障支路的不对称，使得整个网络电流电压三相不对称：。本论文解决不对称短路的问题核心是对称分量法。根据对称分量法采取的具体方法之一是解析法，即把该网络分解为正，负，零序三个对称序网，这三组对称序分量可分别按对称的三相电路分解。然后将其结果叠加起来。求

29、解不对称短路。首先应该计算各原件的序参数和画出等值电路。然后制定各序网络。根据不同的故障类型，确定出以相分量表示的边界条件，进而列出以序分量表示的边界条件，按边界条件将三个序网联合成复合网，由复合网求出故障处各序电流和电压，进而合成三相电流电压。方法二是计算机程序法。通过计算机形成三个序网的节点导纳矩阵，然后利用高斯消去法通过相应公式对他们进行数据运算，即可求得故障端点的等值阻抗。最后根据故障类型选取相关公式计算故障处各序电流电压，进而合成三相电流电压。通过上面两种方法放入计算进行对比可以得出，计算机程序法比较解析法具有计算过程及结果更精确，通用性更强的优点。参考资料1 孙丽华电力工程基础北京

30、：机械工业出版社，20062 李光琦电力系统暂态分析(第二版)北京:中国电力出版社,19953 于永源，杨绮雯电力系统分析(第二版)北京:中国电力出版社,20044 刘万顺. 电力系统故障分析北京：中国电力出版社，19995 刘从爱电力工程基础济南：山东科学技术出版社，19976 电气工程师手册第二版编辑委员会电气工程师手册北京：机械工业出版社，20007 西安交通大学，电力工业部西北电力设计院，电力工业部西北勘测设计院短路电流实用计算方法.北京：电力工业出版社，19827 何仰赞 温增银，电力系统分析(第三版).武汉 华中科技大学出版社， 20019 刘万顺电力系统故障分析北京：中国电力出版

31、社，199810 张俊才电力系统故障的分析计算一对称分量法与相分量法江西工业大学11 邱关源 电路 北京 ：高等教育出版社， 199912 罗建军 杨琦. MATLAB 教程 北京： 电子工业出版社 200513 潘小辉 陈强 MATLAB 5.1 全攻略宝典. 北京： 中国水利水电出版社， 200014 Billinton Roy.ALLAN RONALD (RONALD NORMAN ).RELIABILITY EVALUATION OFPOWER SYSTEMS;NEW YORE ;PLENUM PRESS ;199615 IJ.NAGRATH.D.P.K.OTHARI.MODERN.P

32、OWER SYSTEM ANALYSIS.NEW DELHI附录：不对称短路电流计算程序#include “stdio.h”#define SB 100#define Uav 115Main()int K,j,n,d,S Float x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,Y77,U7,I7,UK,SN,P,Q,X,b,L,Y1,Y2,Y3,K1,K2,K3,IK,IM,SD,PN; K1=K2=1.05;K3=0.96; Clear();电抗.对地导纳 Printf(“please input the patament of generater :Xd,

33、PN,cosan”); Scanf(“%f,%f,%f”,&Xd,&PN,&X); X1=Xd*Sn/(PN/X); Scanf(“%f,%f,%f”,&Xd,&PN,&X); X10=Xd*Sn/(PN/X);Printf(“X1=%f,X10=%fn”,X0,X10);Printf(“please input the putamenter of tramsfer :UK and SNn1:”);Scanf(“%f,%f,”,&UK,&SN);X2=UK*Sn/SN;Printf(“n2:n”);Scanf(“%f,%f,”,&

34、amp;UK,&SN);X3=UK*Sn/SN;Printf(“n3:n”);Scanf(“%f,%f,”,&UK,&SN);X4=UK*Sn/SN;Printf(“X2=%f,X3=%f,X4=%fn”,X2,X3,X4);Printf(“please input the putamenter of load:P,Qn”);Scanf(“%f,%f,”,&P,&Q);X8=1/(Q/100)*Q/100+(P/100)*(P/100)*Q/100);X11=0.35;X12=1.2;Printf(“ please input the M of load

35、:P,cosan”); Scanf(“%f,%f,”,&P,&X);X9=1/6.5*100/(P/X);Printf(“X8=%f,X9=%fn”,X8,X9);Printf(“please input the putamenter of line: n”);Printf(“for example:X=0.39,b=02.92e-6,i=100”);Printf(“n1:”);Scanf(“%f,%f,”,&X,b,&&L);X5=X*L*Sn/(Uav*Uav);Y1=(double)12*b*L*(Uav*Uav)/100;Printf(“n2:”

36、); Scanf(“%f,%f,”,&X,b,&&L); X6=X*L*Sn/(Uav*Uav);Y2=(double)12*b*L*(Uav*Uav)/100; Printf(“n3:”); Scanf(“%f,%f,”,&X,b,&&L);X7= X*L*Sn/(Uav*Uav);Y3=(double)12*b*L*(Uav*Uav)/100;Printf(“X5=%f,X6=%f,X7=%fn”,X5,X6,X7);Printf(“Y1=%f,Y2=%f,Y3=%fn”,Y1,Y2,Y3);Printf(“Y1=Y1,Y2=Y2,Y3=Y3

37、n”,Y1,Y2,Y3);Printf(“please input the K(7) to see Y array: n”);Scanf(“&d”,&K); For(n=0;n<K;n+) 导纳矩阵形成For(j=0;j<K;j+)Ynj=0;Y00=1/X1+1/X2;Y01=-1/(X2*K1);Y10= Y01;Y11=1/X5+Y1+1/(X2*K1*K1);Y12=-1/X5;Y21= Y12;Y22=Y1+1/X5+Y2+1/X7+(1/2*X3*K2*K2);Y23=-1/(1/2*X3*K2);Y24=-1/X6;Y25=-1/X7;Y32= Y23;Y33=1/(1/2*X3)+(1/2*X10);Y42= Y24;Y44=Y2+1/X8 +1/X6;Y55=Y3+1/X7+1/(X4*K3*K3);Y56=-1/(X4*K3);Y65= Y56;Y66= 1/X4;For(n=0;n<K;n+)For(j=0;j<K;j+) Printf(“j%,2ft”,Ynj); If(j=6) Printf(“n”)For(n=0;n<K;n+) 短路电流和功率 Un=0.0;Printf(“please input the number of