电力系统下课程设计-短路电流计算
电力系统的短路电流计算方法
电力系统的短路电流计算方法在电力系统的运行过程中，短路事故是一种常见的故障形式。
短路电流的计算是电力系统设计和运行中重要的一部分，对于确保电力系统的稳定和安全运行至关重要。
本文将介绍电力系统的短路电流计算方法。
一、短路电流的概念和意义短路电流是指在系统中发生短路故障时产生的电流。
短路故障是指两个或多个系统元件之间的短接，导致电流异常增加。
短路电流的大小直接关系到系统设备的安全运行和保护装置的选择。
因此，准确计算短路电流对于系统的设计和运行至关重要。
二、对称短路电流的计算方法对称短路电流是指发生对称型短路故障（如三相短路故障）时的电流。
对称短路电流的计算方法主要有两种：解析法和数值法。
1. 解析法解析法是通过应用基本的电路理论和计算公式来计算短路电流。
首先需要确定短路电流的路线，然后根据系统参数和电路拓扑关系计算短路电流。
这种方法的优点是计算结果准确，但对于复杂的系统结构和参数较多的情况下，计算过程繁琐。
2. 数值法数值法是通过建立系统的模型，根据短路电流计算方程和计算程序进行计算。
数值法的优点是计算过程简单，适用于复杂系统结构和参数较多的情况。
常用的数值法有潮流法、有限差分法和外推法等。
这些方法在复杂系统中具有较大的优势，得到了广泛应用。
三、非对称短路电流的计算方法非对称短路电流是指发生非对称型短路故障时的电流。
由于非对称故障导致的电流不对称，计算方法相对复杂。
1. 正序、负序和零序分量法正序、负序和零序分量法是计算非对称短路电流的常用方法之一。
该方法将非对称电流分解为三个分量，即正序、负序和零序分量。
通过计算各个分量的电流值，再结合系统的参数和拓扑关系进行计算。
这种方法在非对称分析和保护装置选择中应用广泛。
2. 矩阵法矩阵法是一种基于复数计算的方法，通过建立节点矩阵和支路矩阵，求解节点电压和支路电流的未知量。
这种方法具有较强的适应性，能够计算各种复杂情况下的非对称短路电流。
四、短路电流计算中的注意事项在进行短路电流计算时，还需注意以下几个方面：1. 系统参数的准确性系统参数对于计算结果的准确性具有重要影响。
短路电流的计算方法
短路电流的计算方法短路电流是电力系统中的一种重要电气特性，在电路中出现故障时会产生短路电流，对设备、线路和人员等产生威胁。
因此，计算短路电流是电气系统设计和运行中必不可少的一项任务。
本文将介绍短路电流的计算方法。
1.短路电流的定义。
短路电流，也称为故障电流，是指在电路中发生短路时，电源输出电流超过额定电流的情况。
在电气系统设计时，短路电流是评估系统安全性的重要参数之一、计算短路电流的目的是为了保证系统能承受故障时的电流，从而达到系统安全运行的目的。
2.短路电流的计算方法。
（1）简单短路电流的计算方法。
简单短路电流指的是在发生短路时，电路中只有一个源和一个负载的情况。
在这种情况下，短路电流的计算方法如下：Isc = E / Z。
其中，Isc表示短路电流；E表示电源的电动势；Z表示短路阻抗。
在实际应用中，Z是根据电路的图纸和电气参数计算得出的。
因此，短路电流的计算只需知道电源电动势即可。
（2）对称分量法。
对称分量法是计算三相电路短路电流的常用方法。
它将三相电路分解为正序、负序和零序三部分，分别计算其短路电流，再根据三者合成得到总短路电流。
在正常情况下，三相电路中的电流包含正、负、零三种分量。
而在短路情况下，正、负分量的相位角发生变化，但其大小仍然相等，而零序分量则减小为0。
这些特点是对称分量法计算短路电流的基础。
对于一个三相电路，它的短路电流按对称分量法计算的步骤如下：1）分解正、负、零序。
三相电路中，正、负、零序分量的计算方法分别如下：正序分量：Ia0 = Ia, Ib0 = Ibe某p(-2π/3i), Ic0 = Ibe某p(2π/3i)。
负序分量：Ia2 = Ia, Ib2= Ibe某p(2π/3i), Ic2 = Ibe某p(-2π/3i)。
零序分量：I0=(Ia+Ib+Ic)/3。
其中，i为虚数单位。
2）计算短路电流。
在计算正、负、零序分量短路电流前，需先确定短路点的相序。
短路点的相序为零序时：I0sc = 3E / Z。
电力系统短路电流计算
电力系统短路电流计算电力系统短路电流计算是电力系统设计和运行中非常重要的一项工作。
短路电流是指在系统发生故障时电流的最大值，通常由短路电流计算来确定。
短路电流的计算对于保护设备的选择、电路设计和系统运行状态的分析都具有重要意义。
短路电流计算主要分为对称分量法和非对称分量法两种方法。
下面将对这两种方法进行详细介绍。
1.对称分量法：对称分量法是一种传统的短路电流计算方法，它将三相电流分解为正序、负序和零序三个对称分量，然后再计算每个分量的短路电流。
对称分量法的计算步骤如下：a.首先需要确定系统的短路电流初始值。
可以通过测量系统的各个节点电压和电流来获得。
一般来说，短路电流初始值取系统额定电流的2-3倍。
b.将系统的正常运行条件下的三相电流表示为复数形式：iA,iB和iC。
c.计算三相电流的正序分量：I1=(iA+α^2*iB+α*iC)/3，其中α=e^(j2π/3)，j为虚数单位。
d.计算三相电流的负序分量：I2=(iA+α*iB+α^2*iC)/3e.计算三相电流的零序分量：I0=(iA+iB+iC)/3f.计算每个分量的短路电流。
可以使用短路电流公式和阻抗矩阵来计算。
例如，正序分量的短路电流I1'=Z1*I1，其中Z1为正序阻抗。
g.将三个分量的短路电流叠加得到总的短路电流。
2.非对称分量法：非对称分量法是一种更加准确的短路电流计算方法，它考虑了系统故障时的非对称特性，可以更好地反映系统的短路电流分布。
非对称分量法的计算步骤如下：a.获取系统正常运行条件下的三相电流。
b. 将三相电流转换为abc坐标系下的矢量形式。
c.计算叠加故障电流矢量。
d. 将叠加故障电流矢量转换为dq0坐标系的正序、负序和零序分量。
e.根据正、负、零序分量计算短路电流。
非对称分量法相比于对称分量法更加准确，但在计算过程中需要考虑更多的参数和细节，计算复杂度较高。
需要注意的是，短路电流计算是在假设系统中所有设备均采用理想的电气参数的情况下进行的。
电力系统短路电流计算方法
电力系统短路电流计算方法短路电流计算呢，是电力系统分析里挺重要的一事儿。
有一种比较简单的方法叫欧姆法。
这个就有点像咱们以前学的简单电路计算，根据电路里的电阻、电抗这些参数来算。
比如说，知道了电源的电动势，再把电路里的电阻、电抗啥的都找出来，按照欧姆定律的那些变形公式就能算出短路电流啦。
不过这个方法呢，在一些复杂的电力系统里就有点吃力喽，因为电力系统一复杂，电路结构就变得超级复杂，找那些参数就像在一团乱麻里找线头一样。
还有一种方法是标幺制法。
这名字听起来是不是有点酷？这个方法就很巧妙啦。
它是把各个物理量都转化成相对值来计算。
就好比把不同单位的东西都换算成一种特殊的相对单位，这样在计算的时候就方便很多。
在标幺制法里呢，我们要先选一个基准值，这个基准值就像是一个标准的尺子，然后把所有的参数都按照这个尺子来换算成标幺值。
这样做的好处就是不管电力系统多复杂，只要按照这个套路来，计算就会变得比较清晰。
在计算短路电流的时候啊，我们还得考虑很多因素呢。
像发电机的电抗，线路的电抗这些都是影响短路电流大小的关键因素。
发电机要是电抗小，短路的时候电流可能就会很大，就像一个大力水手突然发力一样。
线路的电抗也类似，如果线路短，电抗小，那短路电流也容易变大。
另外呀，计算的时候还要注意短路的类型。
有三相短路、两相短路这些不同的情况呢。
三相短路是最严重的一种短路情况，就像三个小伙伴同时出了问题一样。
它产生的短路电流往往是最大的。
不同类型的短路，计算方法也会有一些小区别。
宝子，电力系统短路电流计算虽然有点复杂，但是只要掌握了这些方法，多做几道题，也就慢慢熟悉啦。
就像学骑自行车，刚开始有点难，摔几跤，慢慢就骑得顺溜啦。
短路电流计算方法
短路电流计算方法短路电流是指在电路中出现短路时所产生的电流。
短路电流的计算对于电路的设计和保护具有重要意义。
正确计算短路电流可以帮助我们选择合适的电器设备和保护装置，从而确保电路的安全运行。
下面将介绍一些常见的短路电流计算方法。
首先，我们需要了解短路电流的定义。
短路电流是指在电路中出现短路时，电流突然增大的现象。
短路电流的大小取决于电路的阻抗、电压和负载等因素。
在进行短路电流计算时，我们需要考虑这些因素，并采用适当的方法进行计算。
一种常见的短路电流计算方法是采用对称分量法。
对称分量法是一种基于对称分量理论的电路分析方法，通过将三相电路中的不对称系统转化为对称系统，简化了电路的分析和计算过程。
在使用对称分量法进行短路电流计算时，我们需要先将电路转化为正序、负序和零序对称分量，然后分别计算它们的短路电流，最后将它们合成为总的短路电流。
另一种常用的短路电流计算方法是采用复功率法。
复功率法是一种基于复功率理论的电路分析方法，通过将电路中的各个元件的功率表示为复数形式，简化了电路的分析和计算过程。
在使用复功率法进行短路电流计算时，我们需要先将电路中各个元件的复功率表示出来，然后利用复功率的运算规则进行计算，最终得到短路电流的大小和相位。
除了对称分量法和复功率法，还有一些其他的短路电流计算方法，如有限元法、潮流法等。
这些方法各有特点，适用于不同类型的电路和不同的计算要求。
在实际工程中，我们可以根据具体的情况选择合适的方法进行短路电流计算。
总的来说，短路电流的计算对于电路的设计和保护具有重要意义。
正确计算短路电流可以帮助我们选择合适的电器设备和保护装置，从而确保电路的安全运行。
在进行短路电流计算时，我们可以采用对称分量法、复功率法等不同的方法，根据具体的情况选择合适的方法进行计算。
希望本文介绍的短路电流计算方法对大家有所帮助。
电力系统的短路电流计算
电力系统的短路电流计算电力系统的短路电流计算是电力工程中一个非常重要的环节，它可以帮助工程师确保电力系统的运行安全和稳定。
短路电流计算通常涉及到电力系统的拓扑结构、电气设备的参数以及电源的特性等多个方面，本文将详细介绍短路电流计算的方法和步骤。
一、短路电流计算的目的短路电流计算的主要目的是确定电力系统中的各个节点、支路以及设备上出现短路时所产生的电流大小，从而判断设备和电气系统是否能够承受这些电流并确保系统的正常运行。
通过短路电流计算，我们可以评估电力系统的稳定性、选择合适的保护设备以及确定设备参数和系统结构等重要工作。
二、短路电流计算的方法1. 传统短路电流计算法传统的短路电流计算法主要通过手工计算实现，通常包括以下几个步骤：首先，需要确定电力系统的拓扑结构，包括各个节点的连线关系和支路连接情况；其次，需要收集系统中各个设备的参数，如电流互感器、变压器、发电机等的额定值以及阻抗等参数；然后，根据短路电流计算公式，对各个节点进行计算，并确定电流的大小和方向；最后，通过对计算结果的分析，判断系统的稳定性和是否需要采取相应的措施进行改进。
2. 计算软件辅助短路电流计算法随着计算机技术的不断发展，短路电流计算方法也得到了很大的改进。
现在，我们可以利用专业的电力系统计算软件来辅助进行短路电流的计算。
这些软件可以根据用户输入的电力系统拓扑结构和设备参数，自动进行计算并输出结果。
相比传统的手工计算方法，计算软件的优势在于可以大大提高计算效率和准确性，并且可以处理更加复杂的电力系统结构和参数。
三、短路电流计算的步骤无论是传统的手工计算方法还是计算软件辅助计算方法，短路电流计算的步骤大体上是相似的，下面是一个典型的短路电流计算的步骤：1. 收集系统参数：包括电力系统的拓扑结构、设备参数以及电源特性等信息。
2. 建立短路电流模型：根据系统参数，建立电力系统的等值电路模型，主要包括发电机、线路、变压器、负荷等元件。
短路计算课程设计
课程设计(论文)题目名称短路计算课程设计课程名称电力系统暂态分析学生姓名学号系、专业电气工程系、电气工程及其自动化指导教师年月日邵阳学院课程设计（论文）任务书年级专业级电气工程及其自动化学生姓名学号题目名称短路计算课程设计设计时间课程名称电力系统暂态分析课程编号121202307设计地点校内一、课程设计（论文）目的1.通过课程设计, 使学生掌握电力系统三相短路计算的基本原理与方法；2.掌握并能熟练运用PSCAD/MATLAB仿真软件；3.采用PSCAD/MATLAB软件，做出系统三相接线图。
二、已知技术参数和条件测试系统由三个控制区域组成，区域 1 是一种典型供电系统，总装机容量为5700MV A、最大负荷为5000MW，大部分装机容量距离负荷相对较近，即接于母线3容量为4400MV A的发电厂。
另外，1300MV A是离负荷较远的核电机组，通过长距离500kV 线路向负荷送电。
区域2 代表附近地区的总装机容量和总负荷，此系统总装机容量为60000MV A，最大负荷为40000MW。
区域2通过2条500kV线路与区域1相连，即图1中的线路A和线路B。
区域3是一个大规模相邻系统，装机容量为70000MV A ，最大负荷为50000MW，此区域也通过 2 条500kV线路与区域1 相连，即图1.1中线路F和G。
注：1．此表由指导教师填写，经系、教研室审批，指导教师、学生签字后生效；2．此表1式3份，学生、指导教师、教研室各1份。
指导教师（签字）：学生（签字）：邵阳学院课程设计（论文）评阅表学生姓名学号系电气工程系专业班级电气工程及其自动化题目名称短路计算课程设计课程名称电力系统暂态分析一、学生自我总结二、指导教师评定注：1、本表是学生课程设计（论文）成绩评定的依据，装订在设计说明书（或论文）的“任务书”页后面；2、表中的“评分项目”及“权重”根据各系的考核细则和评分标准确定。
摘要随着经济的快速增长和社会的不断进步,供电系统的安全性和可靠性问题日益成为人们关注的焦点。
电力系统短路计算课程设计
电力系统短路计算课程设计1. 引言电力系统短路计算是电力系统工程中的重要内容之一。
它用于确定电力系统中各个组件（如发电机、变压器、线路、开关等）的短路电流以及电力系统的短路容量。
本课程设计旨在帮助学生深入理解电力系统短路计算的基本原理和方法，培养学生的问题分析和解决能力。
2. 实验目的•掌握电力系统短路计算的基本原理和方法；•学习使用电力系统短路计算软件进行短路计算；•培养学生的实际操作能力和数据处理能力。
3. 实验内容本次课程设计包括以下实验内容：1.了解电力系统短路计算的基本原理和方法；2.学习使用PSS/E软件进行短路计算；3.对示例电力系统进行短路计算，并绘制短路电流分布图；4.分析短路电流对电力系统设备的影响。
4. 实验步骤4.1 实验准备安装PSS/E软件并了解其基本操作。
4.2 系统建模•根据实验要求，选择合适的电力系统进行建模；•绘制电力系统的单线图。
4.3 数据采集•从电力系统实际运行数据中采集所需的电气参数；•对采集到的数据进行整理和校验。
4.4 短路计算•使用PSS/E软件对电力系统进行短路计算；•分析计算结果，得到各个节点的短路电流。
4.5 短路电流分布•根据计算结果，绘制电力系统的短路电流分布图；•分析电力系统中短路电流的分布规律。
4.6 设备影响分析•根据短路电流分布图，分析短路电流对电力系统设备的影响；•提出相应的设备保护措施。
5. 结果与分析根据实际操作和数据处理的结果，对电力系统短路计算进行结果分析。
可以对不同节点的短路电流进行比较，并针对计算结果进行讨论和总结。
6. 实验总结本次课程设计通过实际操作和数据处理，加深了对电力系统短路计算基本原理和方法的理解。
同时，培养了学生的实际操作能力和问题分析能力。
通过分析电力系统的短路电流分布，提出了针对电力系统设备的保护措施。
本次课程设计对于提高学生的专业素养和解决实际工程问题具有一定的指导意义。
7. 参考文献1.电力系统短路计算教程2.PSS/E软件使用手册。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
电力系统下课程设计-短路电流计算《电力系统分析》课程设计报告题目： 3G9bus短路电流计算系别电气工程学院专业班级10级电气四班学生姓名学号指导教师提交日期 2012年12月10日目录一、设计目的 (4)二、短路电流计算的基本原理和方法 (4)2.1电力系统节点方程的建立 (4)2.2利用节点阻抗矩阵计算短路电流 (5)三、3G9bus短路电流在计算机的编程 (7)3.1、三机九节点系统 (7)3.3输出并计算结果 (14)四．总结 (16)一、设计目的1.掌握电力系统短路计算的基本原理；2.掌握并能熟练运用一门计算机语言（MATLAB 语言或FORTRAN 或C 语言或C++语言）；3.采用计算机语言对短路计算进行计算机编程计算。
二、短路电流计算的基本原理和方法2.1电力系统节点方程的建立利用节点方程作故障计算，需要形成系统的节点导纳（或阻抗）矩阵。
一般短路电流计算以前要作电力系统的潮流计算，假定潮流计算的节点导纳矩阵已经形成，在此基础上通过追加支路的方式形成电力短路电流计算的节点导纳矩阵YN 。
1）对发电机节点在每一发电机节点增加接地有源支路 i E 与i i i Z R jX =+串联 求短路稳态解:i QiE E =i i qiZ R jX =+求短路起始次暂态电流解:i i E E ''= i i i Z R jX ''=+一般情况下发电机定子绕组电阻忽略掉，并将i E 与i i i Z R jX =+的有源支路转化成电流源i i i I E Z =与导纳1i i i i i Y G B R jX =+=+并联的形式2）负荷节点的处理负荷节点在短路计一算中一般作为节点的接地支路，并用恒定阻抗表示，其数值由短路前瞬间的负荷功率和节点实际电压算出，即首先根据给定的电力系统运行方式制订系统的等值电路，并进行各元件标么值参数的计算，然后利用变压器和线路的参数形成不含发电机和负荷的节点导纳矩阵 YN 。
2ˆk LDk LDk LDkLDkV Z R jX S=+=2ˆLDkLDk LDk LDkk SY G jB V =+=3）如果短路故障发生在输电线路内，在短路点要增加新节点，将一条输电线分成两段，并形成短路电流计算的节点导纳矩阵Y，最后形成包括所有发电机支路和负荷支路的节点方程如下：YV I=一般，节点导纳矩阵 Y 阵与 YN 阶次相同。
其差别只在于 YN 阵不含发电机和负荷的阻抗；节点注入电流向量 I中只有发电机端节点的电流不为零。
有非零电流源注人的节点称为有源节点。
系统中的同步调相机可按发电机处理。
任进行起始次暂态电流计算时，大型同步电动机、感应电动机以及以电动机为主要成分的综合负荷，特别是在短路点近处的这些负荷，必要时也可以用有源支路表示，并仿照发电机进行处理。
必须指出，在计算机已普遍应用的情况下，如果有必要的话，只要能提供短路计算所需的准确的原始数据，对短路进行更精确的计算并不存在什么障碍。
2.2利用节点阻抗矩阵计算短路电流假定系统中的节点f经过渡阻抗zf发生短路。
过渡阻抗渡阻抗zf，不参与形成网络的节点导纳（或阻抗）矩阵。
图6-3中方框内的有源网络代表系统正常状态的单相等值网络。
现在我们保持故障处的边界条件不变，把网络的原有部分同故障支路分开（见图6-3）。
容易看出，对于正常状态的网络而言，发生短路相当于在故障节点f增加了一个注人电流一 If （短路电流以流出故障点为正，节点电流则以注入为正）。
因此，网络中任一节点 i的电压可表示为：(1,2,..,,.,)i ij j if fj GV Z I Z I i f n∈=-=∑（6-3）由式（6-3）可见，任一节点 i 的电压郁由两项叠加而成。
第一项表示当 If=0时由网络内所有电源在节点i产生的电压，也就是短路前瞬间正常运行状态一 F 的节点电压，这是节点电压的正常分量，记为(0)i V是电网的潮流解。
第二项是当网络中所有电流源都断开，仅仅由短路电流If 在节点i 产生的电压，这就是节点电压的故障分童。
上述两个分量的叠加，就等于发生短路后节点，的实际电压，即(0)(1,2,..,,.,)i i if fV V Z I i f n =-=(0)f f if fV V Z I =- （6-4）式中，(0)i V 是短路前故障点的正常电压；当i=f 时，ff Z 是故障节点 f 的自阻抗，也称输入阻抗。
f f fV z I =（边界条件方程）带入(0)f f if fV V Z I =-可得：(0)f f ff fV I Z z =+(0)(0)(1,2,..,,.,)ifi i f ff fZ V V V i f n Z Z =-=+p qpq pqkV V I Z -=对于非变压器支路，令 k=1。
从计算公式（6-7）和（6-8）可以看到，式中所用到的阻抗矩阵元素都带有列标f 。
这就是说，如果网络在正常状态下的节点电压为已知，为了进行短路计算，只须利用节点阻抗矩阵中与故障点f 对应的一列元素。
因此，尽管是采用了阻抗型的节点方程，但是并不需要作出全部阻抗矩阵。
在短路的实际计一算中，一般只需形成网络的节点导纳矩阵．并根据具体要求，用第四章所讲的方法求出阻抗矩阵的某一列或某几列元素即可。
在应用节点阻抗矩阵进行短路计算时，我们都将采用这种算法。
1） 解潮流计算，(0)(1,2,..,,.,)i V i f n =2） 修正解潮流的YN 形成Y 3） 指定短路点f4） 计算节点阻抗矩阵第k 列111211212222121200100n k n k n n nn nk f k ,,,n Y Y Y Z YY Y Z Y Y Y Z =⎡⎤⎢⎥⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=←⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦第元素,()5）计算短路电流：(0)ffff fVIZ z=+6）计算节点电压：(0)(1,2,..,,.,)i i if fV V Z I i f n =-=7）计算支路电流：p q pqpqkV V IZ-=8）输出计算结果三、3G9bus短路电流在计算机的编程3.1、三机九节点系统29836415(3)f7图1 三机九节点系统表1 九节点系统支路参数支路R（p.u.）X（p.u.）B/2（p.u.）1~4 0 0.0576 1.02~7 0 0.0625 1.03~9 0 0.0586 1.04~5 0.01 0.085 0.0884~6 0.017 0.092 0.0795~7 0.032 0.161 0.1536~9 0.039 0.17 0.1797~8 0.0085 0.072 0.07458~9 0.0119 0.1008 0.1045表2 九节点系统发电机参数发电机编号节点类型P G（p.u.）V G（p.u.）dX''（p.u.）E''（p.u.）1 Vδ 1.04 0.3 1.1372 PV 1.63 1.025 0.3 1.2113 PV 0.85 1.025 0.3 1.043表3 九节点系统负荷参数节点编号节点类型Pi（p.u.）Qi（p.u.）4 PQ 0 05 PQ 1.25 0.56 PQ 0.9 0.37 PQ 0 08 PQ 1 0.359 PQ 0 03.2 程序设计主函数Sbase_MV A=100.fid=fopen('Nodedata.txt');N=textscan(fid, '%s %u %d %f %f %f %f %f %f')fclose(fid);busnumber=size(N{1},1)for i=1:busnumberBus(i).name=N{1}(i);Bus(i).type=N{2}(i);Bus(i).no=i;Bus(i).Base_KV=N{3}(i);Bus(i).PG=N{4}(i);Bus(i).QG=N{5}(i);Bus(i).PL=N{6}(i);Bus(i).QL=N{7}(i);Bus(i).pb=N{8}(i);Bus(i).V=1.0;Bus(i).angle=0;endfid=fopen('Aclinedata.txt');A=textscan(fid, '%s %s %f %f %f %f')fclose(fid);aclinenumber=size(A{1},1)for i=1:aclinenumberAcline(i).fbname=A{1}(i);Acline(i).tbname=A{2}(i);Acline(i).Base_KV=A{3}(i);Acline(i).R=A{4}(i);Acline(i).X=A{5}(i);Acline(i).hB=A{6}(i);for k=1:busnumberif strcmp(Acline(i).fbname, Bus(k).name)Acline(i).fbno=Bus(k).no;endif strcmp(Acline(i).tbname, Bus(k).name)Acline(i).tbno=Bus(k).no;endendendfid=fopen('Transdata.txt');T=textscan(fid, '%s %f %f %s %f %f %f %f') fclose(fid);tansnumber=size(T{1},1)for i=1:tansnumberTrans(i).fbname=T{1}(i);Trans(i).fbBase_KV=T{2}(i);Trans(i).fbrated_KV=T{3}(i);Trans(i).tbname=T{4}(i);Trans(i).tbBase_KV=T{5}(i);Trans(i).tbrated_KV=T{6}(i);Trans(i).R=T{7}(i);Trans(i).X=T{8}(i);for k=1:busnumberif strcmp(Trans(i).fbname, Bus(k).name)Trans(i).fbno=Bus(k).no;endif strcmp(Trans(i).tbname, Bus(k).name)Trans(i).tbno=Bus(k).no;endendTrans(i).k=Trans(i).tbrated_KV*Trans(i).fbBase_KV/Trans(i).fbrated_KV/Trans(i).tbBase_ KV;tempx=Trans(i).fbrated_KV^2/Trans(i).fbBase_KV^2;Trans(i).X=tempx*Trans(i).X;Trans(i).R=tempx*Trans(i).R;end%N=0%Trans(1)%Trans(2)% for Y=G+jB matrix[G,B,B2]=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber);%B:=B';B2:=B"dlmwrite('Gmatrix.txt', G, 'delimiter', '\t','precision', 6);dlmwrite('Bmatrix.txt', B, 'delimiter', '\t','precision', 6);GBB2pause[JP,JQ]=FormJPQmatrix(Bus,B,B2,busnumber);JPiJP=-inv(JP)JQiJQ=-inv(JQ)pause%maxiteration=0for i=1:busnumberNodeV(i)=Bus(i).V;Nodea(i)=Bus(i).angle;VX(i)=Bus(i).V*cos(Bus(i).angle);VY(i)=Bus(i).V*sin(Bus(i).angle);dQGQL(i)=Bus(i).QG-Bus(i).QL;dPGPL(i)=Bus(i).PG-Bus(i).PL;endNodeV=NodeV'Nodea=Nodea'%VX=VX'%VY=VY'dQGQL=dQGQL'dPGPL=dPGPL'pause%for nointer=1:10maxdP=1.;maxdQ=1.;epsilon=0.000001;noiteration=0;while (maxdP>epsilon)&(maxdP>epsilon)[deltaP,deltaQ,maxdP,maxdQ]=FormdPQvector(Bus,NodeV,Nodea,dQGQL,dPGPL,B, G,busnumber);deltaP;deltaQ;maxdP;maxdQ;da=iJP*deltaP;dV=iJQ*deltaQ;Nodea=Nodea+da;NodeV=NodeV+dV;noiteration=noiteration+1;if noiteration>20breakendendfor i=1:busnumberBus(i).V=NodeV(i);NodeV(i)=NodeV(i)*Bus(i).Base_KV;Bus(i).angle=Nodea(i);Nodea(i)=Nodea(i)*180/pi;endnoiterationNodea=Nodea'NodeV=NodeV'Clear子函数%生成G、B矩阵function[G,B,X]=FormYmatrix(Bus,busnumber,Acline,aclinenumber,Trans,tansnumber) Y=zeros(busnumber);X=zeros(busnumber);for i=1:busnumberY(i,i)=Y(i,i)+Bus(i).pb*j;endfor i=1:aclinenumberf=Acline(i).fbno;t=Acline(i).tbno;Y(f,f)=Y(f,f)+Acline(i).hB*j+1/(Acline(i).R+Acline(i).X*j);Y(t,t)=Y(t,t)+Acline(i).hB*j+1/(Acline(i).R+Acline(i).X*j);Y(f,t)=Y(f,t)-1/(Acline(i).R+Acline(i).X*j);Y(t,f)=Y(t,f)-1/(Acline(i).R+Acline(i).X*j);X(f,f)=X(f,f)-1/Acline(i).X;X(t,t)=X(t,t)-1/Acline(i).X;X(f,t)=1/Acline(i).X;X(t,f)=1/Acline(i).X;endfor i=1:tansnumberf=Trans(i).fbno;t=Trans(i).tbno;Y(f,f)=Y(f,f)+1/(Trans(i).R+Trans(i).X*j);Y(t,t)=Y(t,t)+1/(Trans(i).R+Trans(i).X*j)/Trans(i).k^2;Y(f,t)=Y(f,t)-1/(Trans(i).R+Trans(i).X*j)/Trans(i).k;Y(t,f)=Y(t,f)-1/(Trans(i).R+Trans(i).X*j)/Trans(i).k;X(f,f)=X(f,f)-1/Trans(i).X;X(t,t)=X(t,t)-1/Trans(i).X;X(f,t)=1/Trans(i).X;X(t,f)=1/Trans(i).X;endG=real(Y);B=imag(Y);end%生成JP、JQ矩阵function [JP,JQ]=FormJPQmatrix(Bus,B,B2,busnumber)JP=B;JQ=B2;for i=1:busnumberif Bus(i).type==1for k=1:busnumberJQ(i,k)=0.;JQ(k,i)=0.;JP(i,k)=0.;JP(k,i)=0.;endJQ(i,i)=1.;JP(i,i)=1.;endif Bus(i).type==3for k=1:busnumberJQ(i,k)=0.;JQ(k,i)=0.;endJQ(i,i)=1.;endendend%计算偏节点PQ差量function[deltaP,deltaQ,maxdP,maxdQ]=FormdPQvector(Bus,NodeV,Nodea,dQGQL,dPGPL,B,G,bu snumber)deltaQ=dQGQL;deltaP=dPGPL;maxdP=0.;maxdQ=0.;for i=1:busnumberif Bus(i).type==1deltaQ(i)=0.;deltaP(i)=0.;endif Bus(i).type==3deltaQ(i)=0.;%y1=0;%y2=0;y3=0;for k=1:busnumberif (B(i,k)~=0|G(i,k)~=0)%y1=y1+(G(i,k)*VX(k)-B(i,k)*VY(k));%y2=y2+(G(i,k)*VY(k)+B(i,k)*VX(k));y3=y3+NodeV(k)*(G(i,k)*cos(Nodea(i)-Nodea(k))+B(i,k)*sin(Nodea(i)-Nodea(k)));endenddeltaP(i)=deltaP(i)-y3*NodeV(i);%deltaP2(i)=(deltaP2(i)-(y1*VX(i)+y2*VY(i)))/Bus(i).V;endif Bus(i).type==2%y1=0;%y2=0;y3=0;y4=0;for k=1:busnumberif (B(i,k)~=0|G(i,k)~=0)%y1=y1+(G(i,k)*VX(k)-B(i,k)*VY(k));%y2=y2+(G(i,k)*VY(k)+B(i,k)*VX(k));y3=y3+NodeV(k)*(G(i,k)*cos(Nodea(i)-Nodea(k))+B(i,k)*sin(Nodea(i)-Nodea(k))); y4=y4+NodeV(k)*(G(i,k)*sin(Nodea(i)-Nodea(k))-B(i,k)*cos(Nodea(i)-Nodea(k)));endenddeltaP(i)=deltaP(i)-y3*NodeV(i);%deltaP2(i)=(deltaP2(i)-(y1*VX(i)+y2*VY(i)))/Bus(i).V;deltaQ(i)=deltaQ(i)-y4*NodeV(i);%deltaQ2(i)=(deltaQ2(i)-(y1*VY(i)-y2*VX(i)))/Bus(i).V;endif maxdP<abs(deltaP(i));maxdP=abs(deltaP(i));endif maxdQ<abs(deltaQ(i));maxdQ=abs(deltaQ(i));enddeltaP(i)=deltaP(i)/NodeV(i);deltaQ(i)=deltaQ(i)/NodeV(i);endend3.3输出并计算结果由以上程序可以得出：V：（1）进行系统正常运行状态的潮流计算，求得(0)iV＝1-（Z if/(Z ff+Z f)）因为金属性短路时Z f＝0公式为己知公式(0)i(0)V＝1-（Z if/Z ff）再根据公式可得，iV1&(0) =0.2774 ；V2&(0) ＝0.6770；V3&(0)＝0.6770；V4&(0) ＝0；V5&(0)＝0.1870；V6&(0)＝0.1934；V7&(0) ＝0.54197；V8&(0)＝0.5456；V9&(0)＝0.5511（2）不含发电机和负荷的节点导纳矩阵Y N：(3)形成包括所有发电机支路和负荷支路的节点导纳矩阵Y ，即在Y N 中的发电机节点和负荷节点的自导纳上分别增加发电机导纳i y 和负荷导纳,LD i y ：（4）1Z Y -=，计算节点阻抗矩阵，从而得到阻抗矩阵中的第f 列；（5）计算短路电流因为Z f ＝0，所以短路电流公式为：I f &(0)=1/ Z ff母线对地短路电流标幺值：I f = - 9.1239i母线对地短路电流有名值：I fymzh＝2.2903（6）计算各支路的短路电流己知公式I pq&(0)=(k*V p&(0)- V q&(0))/Z pq所以I14＝4.85i；I27＝2.16i；I39＝2.15i；I45＝-0.7399+6.28907i；I46＝-0.37562+2.03277i；I57＝-0.421564+2.12099i；I69＝-0.458575+1.998i；I78＝-0.00587+0.04972i；I89＝-0.00635+0.05381i；四．总结MATLAB是一个数学建模和计算软件，通过编程能够将复杂的矩阵计算转移到计算机上进行演算，极大的降低了人手计算的强度，给计算工作带来了极大的便利。