计算机辅助分析课程设计.doc

课 程 设 计课程名称： 电力系统分析 设计题目：基于Matlab计算程序的电力系统运行分析学 院： 电力工程学院 专 业： 电气工程自动化 年 级： 学生姓名： 指导教师： 日 期： 教 务 处 制目录前 言1第1章 参数计算 2 第一节 原始数据 2 第二节 电网模型的建立 3第二章 潮流计算 6 第一节 系统参数的设置 6 第二节 程序的调试 7第三节 网损和潮流分析 14第三章 短路故障的分析计算16第一节 序网分析16 第二节 短路计算19心得体会35参考文献36前 言电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷.各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。在电力系统中可能发生的各种故障中，危害最大且发生概率较高的首推短路故障。产生短路故障的主要原因是电力设备绝缘损坏。短路故障分为三相短路、两相短路、单相接地短路及两相接地短路。其中三相短路时三相电流仍然对称，其余三类短路统成为不对称短路。短路故障大多数发生在架空输电线路。电力系统设计与运行时，要采取适当的措施降低短路故障的发生概率。短路计算可以为设备的选择提供原始数据。38第1章 参数计算第一节 原始数据一、 目标电网接线图图11二、系统参数表1. 线路参数表线路编号线路型号线路长度（km）线路电阻/km线路正序电抗/km线路容纳之半S/km4-5LGJ-240/301130.0470.41.784-6LGJ-120/701200.0741.475-7LGJ-120/251650.0791.606-9LGJ-95/551660.0921.807-8LGJ-240/30920.0471.788-9LGJ-240/301220.0471.78说明：线路零序电抗为正序电抗3倍。表2. 变压器参数表线路编号变压器型号变压器变比（kV）短路电压百分数(%)2-7SSPL-22000024232.5%/2010.433-9SSPL-12000024232.5%/155.811-4SSPL-24000024232.5%/17.511.42说明：变压器零序电抗与正序电抗相等，且均为/Y0接法。表3. 发电机参数表1发电机额定功率MW额定电压kV额定功率因数120016.50.852180180.85310013.80.85表4. 发电机参数表2发电机母线名（S）（S）（S）1147.2800.320.130.210.218.962212.8001.930.261.870.436.000.535336.0201.510.211.450.298.590.60表5. 负荷数据表节点号有功负荷（MW）无功负荷（MVA）51355061003088035第二节 电网模型的建立本设计中，采用精确计算算法，选取=100MVA，=220KV，将所有支路的参数都折算到220KV电压等级侧；一、计算公式 （1）线路参数的计算：R=rlSBUB2 X=xlSBUB2 B2=blUB2SB（2）变压器参数的计算：（3）发电机参数的计算：（暂态分析时，只用到发电机的暂态电抗来代替其次暂态电抗，故只求出暂态电抗）（4）负荷节点的计算二、使用编程的结果如下：线路电抗X= 0.0934 0.0992 0.1364 0.1372 0.0760 0.1008线路电阻R= 0.0109 0.0183 0.0269 0.0316 0.0089 0.0118线路正序、负序阻抗Z1=Z2= Columns 1 through 5 0.0109 + 0.0937i 0.0183 + 0.0992i 0.0269 + 0.1364i 0.0316 + 0.1372i 0.0089 + 0.0760i Column 6 0.0118 + 0.1008i 线路容纳之半B/2= Columns 1 through 5 0 + 0.0974i 0 + 0.0854i 0 + 0.1280i 0 + 0.1450i 0 + 0.0793i Column 6 0 + 0.1050i线路零序阻抗Z0: Columns 1 through 5 0.0109 + 0.2802i 0.0183 + 0.2975i 0.0269 + 0.4091i 0.0316 + 0.4112i 0.0089 + 0.2281i Column 6 0.0118 + 0.3024i 变压器正、负序电抗值Xt： 0 + 0.0576i 0 + 0.0574i 0 + 0.0586i变压器零序电抗值Xt0： 0 + 0.0576i 0 + 0.0574i 0 + 0.0586i发电机暂态电抗Xdg= 0 + 0.0476i 0 + 0.0803i 0 + 0.1103i负荷节点阻抗Zl= 0.6514 - 0.2413i 0.9174- 0.2752i 1.049 - 0.4590i以上结果均采用程序计算。三、等值电路模型该等值电路图全部数值以标幺值表示，具体图12所示：图12 第二章 潮流计算第一节 系统参数的设置设计中要求所有结点电压不得低于1.0p.u.，也不得高于1.05p.u.，若电压不符合该条件，可采取下面的方法进行调压：（1） 改变发电机的机端电压（2） 改变变压器的变比（即改变分接头）（3） 改变发电机的出力（4） 在电压不符合要求的结点处增加无功补偿调压方式应属于逆调压。结点的分类：根据电力系统中各结点性质的不同，将结点分为三类：PQ结点、PV结点和平衡结点，在潮流计算中，大部分结点属于PQ结点，小部分结点属于PV结点，一般只设一个平衡结点。对于平衡结点，给定其电压的幅值和相位，整个系统的功率平衡由这一点承担。本设计中，选1号节点为平衡节点；2、3号节点为P、U节点；4、5、6、7、8、9号结点为P、Q节点。设计中，节点数:n=9，支路数:nl=9，平衡母线节点号:isb=1，误差精度:pr=0.0001。由支路参数形成的矩阵:B1= 1 4 0.0576i 0 1.0000 0; 2 7 0.0574i 0 1.0000 0; 3 9 0.0586i 0 1.0000 0; 4 5 0.0109+0.0934i 2*0.0974i 1.0000 0; 4 6 0.0183+0.0992i 2*0.0854i 1.0000 0; 5 7 0.0269+0.1364i 2*0.1280i 1.0000 0; 6 9 0.0316+0.1372i 2*0.1450i 1.0000 0; 7 8 0.0089+0.0760i 2*0.0793i 1.0000 0; 8 9 0.0118+0.1008i 2*0.1050i 1.0000 0; ;%支路参数矩阵由各节点参数形成的矩阵:B2=0 0 1.0000 1.0000 0 1.0000;1.8000 0 1.0000 1.0000 0 3.0000;1.0000 0 1.0000 1.0000 0 3.0000;0 0 1.0000 0 0 2.0000;0 1.3500+0.5000i 1.0000 0 0 2.0000;0 0.0000+0.3000i 1.0000 0 0 2.0000;0 0 1.0000 0 0 2.0000;0 0.8000+0.3500i 1.0000 0 0 2.0000;0 0 1.0000 0 0 2.0000;%节点参数矩阵由节点号及其对地阻抗形成的矩阵:X=1,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0; 9,0;第二节 程序的调试一 未调试前，原始参数运行结果如下：输入“var=1”时将选用牛顿-拉夫逊法来进行潮流计算；输入“var=其他数字时将选用PQ分解法来进行潮流计算”var=1迭代次数 4没有达到精确要求的个数 14 16 8 0各节点的实际电压标么值E为(节点号从小到大排列): Columns 1 through 5 1.0000 0.9759 + 0.2183i 0.9906 + 0.1368i 0.9840 - 0.0240i 0.9498 - 0.0429i Columns 6 through 9 0.9689 - 0.0388i 0.9863 + 0.1148i 0.9814 + 0.0641i 0.9969 + 0.0785i各节点的电压大小V为(节点号从小到大排列): 1.0000 1.0000 1.0000 0.9843 0.9508 0.9697 0.9929 0.9835 1.0000各节点的电压相角O为(节点号从小到大排列): 0 12.6101 7.8603 -1.3973 -2.5839 -2.2909 6.6372 3.7377 4.5008各节点的功率S为(节点号从小到大排列): Columns 1 through 5 0.4167 + 0.2781i 1.8000 + 0.2175i 1.0000 + 0.0298i 0.0000 - 0.0000i -1.3500 - 0.5000i Columns 6 through 9 -1.0000 - 0.3000i -0.0000 + 0.0000i -0.8000 - 0.3500i 0.0000 - 0.0000i各条支路的首端功率Si为(顺序同您输入B1时一样): 0.4167 + 0.2781i 1.8000 + 0.2175i 1.0000 + 0.0298i -0.2437 - 0.3983i -0.1703 - 0.1903i -1.1063 - 0.1017i -0.8297 - 0.1097i -0.6533 - 0.1056i -0.1467 - 0.2444i各条支路的末端功率Sj为(顺序同您输入B1时一样): -0.4167 - 0.2636i -1.8000 - 0.0288i -1.0000 + 0.0289i 0.2456 + 0.2320i 0.1711 + 0.0316i 1.1427 + 0.0445i 0.8528 - 0.0711i 0.6573 - 0.0157i 0.1472 + 0.0422i各条支路的功率损耗DS为(顺序同您输入B1时一样): 0 + 0.0145i 0.0000 + 0.1887i 0.0000 + 0.0587i 0.0019 - 0.1663i 0.0008 - 0.1587i 0.0364 - 0.0572i 0.0232 - 0.1808i 0.0039 - 0.1213i 0.0005 - 0.2022i以下是每次迭代后各节点的电压值(如图21所示)图21 由运行结果可知，节点4、5、6、7、8电压均不满足要求。故需进行调试，以期各结点电压均满足要求。下面用N-L法进行潮流计算。二 采用NL法进行潮流的计算和分析1）第一次调试将2号发电机电压初值由1.0000升高至1.0400，则修改结果如下：B2=0 0 1.0000 1.0000 0 1.0000;1.8000 0 1.0400 1.0000 0 3.0000;运行结果如下：各节点的电压大小V为(节点号从小到大排列):1.0000 1.0000 1.0000 0.9843 0.9508 0.9697 0.9929 0.9835 1.0000 节点4、5、6、7、8都不满足要在1.00001.0500范围内的要求，改变发电机机端电压的影响不大。2）第二次调试 将1号变压器变比初值由1.0000改至1.0250，则修改结果如下：B1=1.0000 4.0000 0.0576i 0 1.0250 0; 2.0000 7.0000 0.0574i 0 1.0000 0; 运行结果如下：各节点的电压大小V为(节点号从小到大排列): 1.0000 1.0000 1.0000 1.0022 0.9638 0.9823 0.9968 0.9875 1.0038节点5、6、7、8的值，虽然较第二次调试仍然不满足要在1.00001.0500范围内的要求，但电压有所升高。3）第三次调试 将1号变压器变比由1.0250改至1.0500，则修改结果如下：B1=1.0000 4.0000 0.0576i 0 1.0500 0; 2.0000 7.0000 0.0574i 0 1.0000 0; 运行结果如下：各节点的电压大小V为(节点号从小到大排列): 1.0000 1.0000 1.0000 1.0197 0.9765 0.9945 1.0006 0.9914 1.0075节点5、6、8的值，虽然较第三次调试仍然不满足要在1.00001.0500范围内的要求，但电压有所升高。4）第四次调试将1号变压器变比由1.0500改至1.0750，则修改结果如下：B1=1.0000 4.0000 0.0576i 0 1.0750 0; 2.0000 7.0000 0.0574i 0 1.0000 0; 运行结果如下：各节点的电压大小V为(节点号从小到大排列):1.0000 1.0000 1.0000 1.0366 0.9887 1.0064 1.0043 0.9951 1.0111节点5、8的值，不满足要在1.00001.0500范围内的要求，且变压器T1的变比以调至最大。应调节与节点5、8临近的变压器T2的变比。5）第五次调试将2号变压器变比由1.0000改至1.0250，则修改结果如下：B1=1.0000 4.0000 0.0576i 0 1.0750 0;2.0000 7.0000 0.0574i 0 1.0250 0;运行结果如下：各节点的电压大小V为(节点号从小到大排列):1.0000 1.0000 1.0000 1.0414 0.9998 1.0115 1.0225 1.0080 1.0160节点5的值，虽然较第五次调试仍然不满足要在1.00001.0500范围内的要求，但电压有所升高。6)第六次调试将2号变压器变比由1.0250改至1.0500，则修改结果如下：B1=1.0000 4.0000 0.0576i 0 1.0750 0;2.0000 7.0000 0.0574i 0 1.0500 0;运行结果如下：各节点的电压大小V为(节点号从小到大排列):1.0000 1.0000 1.0000 1.0460 1.0105 1.0164 1.0402 1.0204 1.0206所有节点的值，都满足在1.00001.0500范围内的要求。NL法 程序运行结果如下：输入“var=1”时将选用牛顿-拉夫逊法来进行潮流计算；输入“var=其他数字时将选用PQ分解法来进行潮流计算”var=1迭代次数 4没有达到精确要求的个数 15 16 11 0各节点的实际电压标么值E为(节点号从小到大排列): Columns 1 through 8 1.0000 0.9777 + 0.2101i 0.9908 + 0.1353i 1.0457 - 0.0256i 1.0096 - 0.0414i 1.0157 - 0.0373i 1.0342 + 0.1113i 1.0184 + 0.0640i Column 9 1.0175 + 0.0798i各节点的电压大小V为(节点号从小到大排列): 1.0000 1.0000 1.0000 1.0460 1.0105 1.0164 1.0402 1.0204 1.0206各节点的电压相角O为(节点号从小到大排列): 0 12.1309 7.7770 -1.4010 -2.3491 -2.1006 6.1444 3.5958 4.4856各节点的功率S为(节点号从小到大排列): Columns 1 through 8 0.4130 + 0.4739i 1.8000 + 0.2570i 1.0000 - 0.3234i -0.0000 + 0.0000i -1.3500 - 0.5000i -1.0000 - 0.3000i -0.0000 + 0.0000i -0.8000 - 0.3500i Column 9 -0.0000 - 0.0000i各条支路的首端功率Si为(顺序同您输入B1时一样): 0.4130 + 0.4739i 1.8000 + 0.2570i 1.0000 - 0.3234i -0.2287 - 0.4550i -0.1805 - 0.3573i -1.1213 - 0.0450i -0.8195 + 0.0573i -0.6416 - 0.2586i -0.1584 - 0.0914i各条支路的末端功率Sj为(顺序同您输入B1时一样): -0.4130 - 0.4511i -1.8000 - 0.0672i -1.0000 + 0.3881i 0.2307 + 0.2653i 0.1824 + 0.1858i 1.1546 - 0.0553i 0.8414 - 0.2633i 0.6454 + 0.1225i 0.1586 - 0.1248i各条支路的功率损耗DS为(顺序同您输入B1时一样): 0 + 0.0228i 0.0000 + 0.1898i 0.0000 + 0.0647i 0.0019 - 0.1897i 0.0019 - 0.1716i 0.0333 - 0.1003i 0.0219 - 0.2059i 0.0038 - 0.1361i 0.0003 - 0.2163i以下是每次迭代后各节点的电压值(如图所示)图22PQ分解法 程序运行结果如下：输入“var=1”时将选用牛顿-拉夫逊法来进行潮流计算；输入“var=其他数字时将选用PQ分解法来进行潮流计算”var=2迭代次数 8每次没有到达精度要求的有功功率个数为 5 8 8 8 7 4 2 0每次没有达到精度要求的无功功率个数为 6 6 6 6 5 4 2 0各节点的电压标幺值E为（节点号从小到大排）: Columns 1 through 8 1.0000 0.9777 + 0.2101i 0.9908 + 0.1353i 1.0457 - 0.0256i 1.0096 - 0.0414i 1.0157 - 0.0373i 1.0342 + 0.1113i 1.0184 + 0.0640i Column 9 1.0175 + 0.0798i各节点的电压U大小(节点号从小到大排)为: 1.0000 1.0000 1.0000 1.0460 1.0105 1.0164 1.0402 1.0204 1.0206各节点的电压相角O（节点号从小到大排）: 0 12.1308 7.7769 -1.4010 -2.3492 -2.1008 6.1443 3.5957 4.4854各节点的功率S（节点号从小到大排）为: Columns 1 through 8 0.4130 + 0.4739i 1.8000 + 0.2570i 1.0000 - 0.3234i 0.0000 - 0.0002i -1.3500 - 0.5000i -1.0000 - 0.3000i 0.0000 - 0.0000i -0.8000 - 0.3500i Column 9 -0.0000 + 0.0000i各支路的首段功率Si(顺序同您输入B1时一样）为: 0.4130 + 0.4739i 1.8000 + 0.2570i 1.0000 - 0.3234i -0.2288 - 0.4550i -0.1805 - 0.3573i -1.1213 - 0.0450i -0.8195 + 0.0573i -0.6417 - 0.2586i -0.1583 - 0.0914i各条支路的末端功率Sj（顺序同您输入B1时一样）为: -0.4130 - 0.4512i -1.8000 - 0.0672i -1.0000 + 0.3881i 0.2307 + 0.2653i 0.1824 + 0.1857i 1.1546 - 0.0553i 0.8414 - 0.2633i 0.6454 + 0.1225i 0.1586 - 0.1248i各条支路的功率损耗DS（顺序同您输入B1时一样）为: 0 + 0.0228i 0 + 0.1898i 0 + 0.0647i 0.0019 - 0.1897i 0.0019 - 0.1716i 0.0333 - 0.1003i 0.0219 - 0.2059i 0.0038 - 0.1361i 0.0003 - 0.2163i以下是每次迭代后各节点的电压值(如图所示)：图23三 对运行结果的分析：1）、为什么在用计算机对某网络初次进行潮流计算时往往是要调潮流，而并非任何情况下只一次送入初始值算出结果就行呢？要考虑什么条件？各变量是如何划分的？哪些可调？哪些不可调？ 答：潮流计算时功率方程是非线性，多元的具有多解。初始条件给定后得到的结果不一定能满足约束条件要求，要进行调整初值后才能满足。其约束条件有：，。负荷的PQ量为扰动变量，发电机的PQ为控制变量，各节点的V为状态变量。扰动变量是不可控变量，因而也是不可调节的，状态变量是控制变量的函数，因而状态变量和控制变量是可以调节的。所以，计算机对某网络初次进行潮流计算时往往是要调潮流的。2）牛顿拉夫逊法与PQ分解法有哪些联系？有哪些区别？二者的计算性能如何？答：在潮流计算中，牛顿拉夫逊法的实质 是一种逐步线性化的方法，它有很好的收敛性，但要求有合适的初值。 PQ分解法来源于N-L法极坐标形式，又充分考虑在交流高压电网中，变压器和线路的电抗远远大于电阻，且系统中母线有功功率的变化主要受电压相位的影响，而无功功率的变化主要受母线电压幅值的影响，在此基础之上对H、L阵进行简化。由于这种简化只涉及修正方程的系数矩阵，并为改变结点功率平衡方程和收敛判据，从NL法和PQ分解法程序运行结果可以看出不会降低计算结果的精度。由图22和图23比较可知，PQ分解法与N-L法相比大大提高了计算速度，但迭代次数却增多；PQ分解法为恒速率收敛，比起N-L法，其收敛速度变慢。第三节 网损和潮流分析一 选取NL法的数据来分析网损表21调整前后网络损耗比较节点编号1-42-73-94-54-65-76-97-88-9调整前0.0145i0.1887i0.0587i0.0019 - 0.1663i0.0008 - 0.1587i0.0364- 0.0572i0.0232 - 0.1808i0.0039 - 0.1213i0.0005 - 0.2022i牛拉法调整后0.0228i0.1898i0.0647i0.0019 - 0.1897i0.0019 - 0.1716i0.0333 - 0.1003i0.0219 - 0.2059i0.0038 - 0.1361i0.0003 - 0.2163i调节前有功损耗：P=（0+0+0+0.0019+0.0008+0.0364+0.0232+0.0039+0.0005）* 100=6.67MW调节后有功损耗：P= ( 0+0+0+0.0019+0.0019+0.0333+0.0291+0.0038+0.0003)* 100=6.31MW调节前无功损耗：Q=（0.0145+0.1887+0.0587+0.1663+0.1587+0.0572+0.1808+0.1213+0.2022）*100i=114.84i MVar调节后无功损耗：Q=（0.0228+0.1898+0.0647+0.1897+0.1716+0.1003+0.2059+0.1361+0.2163）X100i=129.74 Mvar由以上可知：调压前有功损耗为6.67MW，调节后有功损耗为6.31MW； 调节前无功损耗：114.84iMVar，调节后无功损耗为129.74iMVar。进行调压有功网损减少了，而无功网损却增加了，之所以会出现这种情况是由于调压后，许多节点电压升高了，对地电容的充电功率为：Q=1/2CU2，充电功率与网络电压的平方成正比，所以最终导致无功网损增加。 二 选取NL法的数据来分析潮流1 潮流分布如图24所示图242 根据图24进行分析1）线路有功潮流最有可能的流向是？线路无功潮流最有可能的流向是？发电机节点的注入无功为负值说明了什么？由潮流分布图可知，线路有功功率总是从电压相位超前的结点流向滞后的结点，无功功率的总是从电压幅值高的结点流向电压幅值低的结点。在潮流分布图中，3号发电机发出的无功功率为负值，说明3号发电机从系统吸收无功。2）潮流控制的主要手段有哪些？由于电网运行时负荷等的变化，引起系统中各节点电压的变化。由于发电机容量有限，而电压的调整又必须服从对电压质量的要求，调整的幅度不能太大，另一方面从电网运行经济性和安全等方面考虑，网络中的潮流往往需要控制。在实际的网络潮流控制中主要采用：串联电容（作用以容抗抵偿线路感抗）；串联电抗（作用在于限流）；附加串联加压器（作用在于产生一环流或强制循环功率使强制循环功率与自然分布功率的叠加可达到理想值）等手段控制潮流。第三章 短路故障的分析计算第一节 序网分析一 各序网参数计算设计中，采用近似计算算法，选取 =100MVA，=Uav，将所有支路的参数都用平均电压，计算过程及结果如下：1计算公式 （1）线路参数的计算：R=rlSBUB2 X=xlSBUB2 B2=blUB2SB（2）变压器参数的计算：（3）发电机参数的计算：（暂态分析时，只用到发电机的暂态电抗来代替其次暂态电抗，故只求出暂态电抗）（4）负荷节点的计算2 使用编程的结果如下：线路电抗X= 0.0854 0.0907 0.1248 0.1255 0.0696 0.0922线路电阻R= 0.0100 0.0168 0.0246 0.0289 0.0082 0.0108线路正序、负序阻抗Z1=Z2= Columns 1 through 5 0.0100+0.0854i 0.0168 + 0.0907i 0.0246 + 0.1248i 0.0289 + 0.1255i 0.0082 + 0.0696i Column 6 0.0108 + 0.0923i线路容纳之半B/2= Columns 1 through 5 0 + 0.1060i 0 + 0.0933i 0 + 0.1400i 0 + 0.1580i 0 + 0.0866i Column 6 0 + 0.1150i线路零序阻抗Z0: Columns 1 through 5 0.0100+0.2563i 0.0168 + 0.2722i 0.0246 + 0.3743i 0.0289 + 0.3766i 0.0082 + 0.2087i Column 6 0.0108 + 0.2767i变压器正、负序电抗值Xt： 0 + 0.0476i 0 + 0.0474i 0 + 0.0484i变压器零序电抗值Xt0： 0 + 0.0476i 0 + 0.0474i 0 + 0.0484i发电机暂态电抗Xdg= 0 + 0.0433i 0 + 0.0728i 0 + 0.1000i负荷节点阻抗Zl= 0.6514 - 0.2413i 0.9174- 0.2752i 1.049 - 0.4590i二 各序网图1.正序图图312.负序图图323.零序图图33第二节 短路计算一、三相短路1、三相短路电流的计算原理:利用结点阻抗矩阵和导纳矩阵都可以计算短路电流，其算法有所不同。利用结点阻抗阵时，只要形成了阻抗阵，计算网络中任意一点的对称短路电流和网络中电流、电压的分布非常方便，计算工作量小，但是，形成阻抗阵的工作量大，网络变化时的修改也比较麻烦，而且结点阻抗矩阵是满阵，需要计算机存储量较大。利用结点导纳矩阵计算短路电流，实质是先用它计算出与指定的短路点有关的结点阻抗矩阵的相应列的元素，然后利用公式（为连接节点i和j的支路阻抗）进行短路电流的计算。然而，导纳阵是对称、稀疏阵，极易形成，且网络结构变化时也易于修改。程序运行步骤及对变量的解释： （1） 请输入短路点的数目NF（2） 请输入节点数n（3） 请输入支路数nl（4） 输入各支路参数矩阵B 矩阵B的每行是由下列参数构成的： 某支路的首端号P； 末端号Q，且PQ； 支路的阻抗（R+jX）; 支路的对地电纳； 支路的变比K； 折算到哪一侧的标志（如果支路的首端P处于高压侧请输入“1”，否则请输入“0”）。（5） 输入由短路点号，短路点阻抗形成的矩阵D（6） 请输入由各节点的初电压标幺值形成的列矩阵V0（7） 形成节点阻抗矩阵Z（8） 求短路点电流（9） 求网络中各节点的电压（10）求网络中各支路电流2、7节点、8节点发生三相短路的程序结果：（1） 7节点发生三相短路： 1）短路点的电流：标幺值： 0.8838 -16.9406i有名值（kA）： 0.2218- 402521i 2）各节点的电压：表31节点号电压表幺值有名值（KV）1116.520.33856.13113.841.046240.585 0.5790 + 0.0303i 133.18 +6.9693i6 0.7592 + 0.0302i 174.61 +6.9387i7-2.22E-1608 0.1557 + 0.0446i 35.8192 +10.2614i9 0.4134 + 0.0385i 95.0725 + 8.8533i有名值与标幺值的转化公式： 3）各支路电流：表32支路编号电流表幺值电流有名值（KA）01 0 +23.0947i 0 + 5.7968i02 0 + 4.6497i 0 + 1.1671i03 0 +10.0000i 0 + 2.5100i05 -0.7665 - 0.3304i -0.1924 - 0.0829i06 -0.7502 - 0.2579i -0.1883 - 0.0647i08 -0.1090 - 0.0902i -0.0274 - 0.0226i14 0 - 0.5667i 0 - 0.1423i27 0 - 7.1412i 0 - 1.7925i39 -0.7953 -12.1207i -0.1996 - 3.0423i45 -0.3548 - 5.4681i -0.0891 - 1.3725i46 -0.3326 - 3.1623i -0.0835 - 0.7937i57 0.2428 - 4.6396i 0.0609 - 1.1645i69 -0.0663 - 2.7555i -0.0166 - 0.6916i78 -0.6410 + 2.2376i -0.1609 + 0.5616i89 0.0663 + 2.7912i 0.0166 + 0.7006i3、8节点发生三相短路 1）短路点的电流：标幺值： 0.9693 -10.1289i有名值（kA）： 0.2433 - 205424i 2）各节点的电压：表33节点号电压表幺值有名值（KV）1116.52 0.6697 + 0.0027i12.053113.841.046240.585 0.7952 + 0.0169i 182.90+3.8891i6 0.7118 + 0.0199i 163.72 +4.5738i7 0.5209 + 0.0043i 119.80 +0.98763i8 0 +1.3878e-017i09 0.3036 + 0.0084i 69.8173 + 1.9383i 3 )各支路电流：表34支路编号电流表幺值电流有名值（KA）01 0 +23.0947i 0 + 5.7968i02 -0.0375 + 9.1998i -0.0094 + 2.3091i03 0 +10.0000i 0 + 2.5100i05 -1.0650 - 0.4205i -0.2673 - 0.1055i06 -0.7059 - 0.2334i -0.1772 - 0.0586i08 4.8570e-018 -1.1102e-017i 1.2191e-018 -2.7867e-018i14 0 - 0.5667i 0 - 0.1423i27 -0.0287 - 3.6639i -0.0072 - 0.9196i39 -0.1741 -14.3894i -0.0437 - 3.6117i45 -0.1980 - 2.9367i -0.0497 - 0.7371i46 -0.2193 - 3.6844i -0.0550 - 0.9248i57 0.1011 - 2.1981i 0.0254 - 0.5517i69 0.0913 - 3.2531i 0.0229 - 0.8165i78 0.0617 - 7.4839i 0.015