“三相短路故障分析与计算的算法设计”电力系统专业本科课程设计论文学生姓名：学号：专业：电力系统及其自动化日期：摘要电力系统发生三
相短路故障造成的危害性是最大的。作为电力系统三大计算之一，分析与计算三相短路故障的参数更为重要。本设计是通过两种不同的方法进行分析
与计算三相短路故障的各参数，进一步提高短路故障分析与计算的精度和速度，为电力系统的规划设计、安全运行、设备选择、继电保护等提供重
要依据。基于Matlab最重要的组件之一Simulink中的电力元件库（SimPowerSystems）构建电力系统仿真模型，在M
atlab的平台下仿真电力系统为工程设计和维修提供依据，同时也为电力研究带来大大的便利，利用Simulink中的画图工具搭建电力系
统模型也是进行电力系统故障分析的常用方法，它让电力研究者从大量繁琐的理论分析及复杂的矩阵计算中解脱出来，让庞大的电力系统很直观的呈
现在研究者的面前，从而将庞大的电力网搬进了计算机，为研究带来了巨大的便利。关键词：三相短路计算，Matlab，SimPowerS
ystems，仿真目录1.绪论1.1电力系统三大计算1.2电力系统短路故障概述2.短路电流计算（解析法）2.1参数数据2.2电
抗标幺值定义2.3短路次暂态电流（功率）标幺值计算2.4各元件电抗标幺值2.5三相短路电流及短路功率3.导纳矩阵形成与计算（Y矩
阵法）3.1导纳矩阵等值电路3.2导纳计算公式3.3变压器变比的定义3.4导纳矩阵的形成3.5三相短路电流及短路功率4.计算机
算法设计4.1计算机编程语言4.2程序流程4.3程序清单及说明4.4程序结果及分析5.基于Matlab电力系统模型搭建与仿真5.
1Matlab简介5.2总体设计5.2.1设计内容概述5.2.2设计方案简介5.3详细设计5.3.1电力元件设计5.3.2
电力系统模型的搭建5.3.3电力系统模型三相短路仿真6．结论7．致谢8．参考文献9．附录9.1系统等值电路图9.2计算机算法设
计流程图9.3计算机算法设计程序清单9.4Matlab电力系统模型仿真图绪论1.1电力系统三大计算1.潮流计算研究电力系统稳态
运行情况的一种基本电气计算，常规潮流计算的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）
、网络中的功率分布以及功率损耗等。潮流计算的结果是电力系统稳定计算和故障分析的基础。2.短路故障计算主要研究电力系统中发生故障（包
括短路、断线和非正常操作）时，故障电流、电压及其在电力网中的分布。短路电流计算是故障分析的的主要内容。短路电流计算的目的，是确定短
路故障的严重程度，选择电气设备参数。整定继电保护，分析系统中负序及零序电流的分布，从而确定其对电气设备和系统的影响。3.稳定性计算
给定运行条件下的电力系统，在受到扰动后，重新回复到运行平衡状态的能力。系统中的多数变量可维持在一定的范围，使整个系统能稳定运行。根
据性质的不同，电力系统稳定性可分为功角稳定、电压稳定和频率稳定三类。1.2电力系统短路故障概述电力系统的短路故障是严重的，而又是发
生几率最多的故障。一般说来，最严重的短路是三相短路。当发生短路时，其短路电流可达数万安以至十几万安，它们所产生的热效应和电动力效
应将使电气设备遭受严重破环。为此，当发生短路时，继电保护装置必须迅速切除故障线路，以避免故障部分继续遭受危害，并使非故障部分从不正
常运行情况下解脱出来，这要求电气设备必须有足够的机械强度和热稳定度，开关电气设备必须具备足够的开断能力，即必须经得起“可能最大短路
的侵扰”而不致损坏。因此，电力系统短路电流计算是电力系统运行分析，设计计算的重要环节，许多电业设计单位和个人倾注极大精力从事这一工
作的研究。由于电力系统结构复杂，随着生产发展，技术进步系统日趋扩大和复杂化，短路电流计算工作量也随之增大，采用计算机辅助计算势在必
行。短路故障（横向故障）：电力系统正常运行情况以外的相与相之间或相与地（或中性线）之间的连接。短路类型：4种。单相接地短路（发
生概率最大的短路类型）两相短路两相短路接地三相短路（对称短路、最严重的短路类型）短路原因及危害原因：客观（绝缘破坏：架空线绝缘子
表面放电，大风、冰雹、台风）、主观（误操作）。危害：短路电流大（热效应、电动效应）、故障点附件电压下降、功率不平衡失去稳定、不对
称故障产生不平衡磁通影响通信线路。解决措施：继电保护快速隔离、自动重合闸、串联电抗器等。短路电流计算（解析法）2.1参数数据设基
准容量SB=100MV·A；基准电压UB=UavkV。（1）SB的选取是为了计算元件参数标幺值计算方便，取SB-100MV·A，可
任意设值但必须唯一值进行分析与计算。（2）UB的选取是根据所设计的题目可知系统电压有110kV、6kV、10kV，而平均额定电压分
别为115、6.3、10.5kV。平均电压Uav与线路额定电压相差5%的原则，故取UB=Uav。（3）为次暂态短路电流有效值，短路
电流周期分量的时间t等于初值（零）时的有效值。满足产生最大短路电流的三个条件下的最大次暂态短路电流作为计算依据。（4）为冲击电流，
即为短路电流的最大瞬时值（满足产生最大短路电流的三个条件及时间=0.01s）。一般取冲击电流=××=2.55。（5）为短路电流冲击
系数，主要取决于电路衰减时间常数和短路故障的时刻。其范围为1≤≤2，高压网络一般冲击系数=1.8。2.2电抗标幺值定义（1）发电机
电抗标幺值公式①式中——发电机电抗百分数，由发电机铭牌参数的；——已设定的基准容量（基值功率），；——发电机的
额定有功功率，MW——发电机额定有功功率因数。（2）负载电抗标幺值公式②式中U——元件所在网络的电压标幺值；——负载容
量标幺值；——负载无功功率标幺值。（3）变压器电抗标幺值公式③变压器中主要指电抗，因其电抗，即可忽略，由变压器电抗有名值
推出变压器电抗标幺值为公式④式中%——变压器阻抗电压百分数；——基准容量，MVA、——变压器铭牌参数给定额定容量，
MVA、额定电压，kV；——基准电压取平均电压，kV。（4）线路电抗标幺值公式⑤式中——线路单位长度电抗；——线路
长度，km；——基准容量，MVA；——输电线路额定平均电压，基准电压，kV。输电线路的等值电路中有四个参数，一般电抗，故0
。由于不做特殊说明，故电导、电纳一般不计，故而只求电抗标幺值。（5）电动机电抗标幺值（近似值）cos公式⑥式中——设定的
基准容量，MVA；——电动机额定的有功功率，MW；cos——电动机额定有功功率因数。2.3短路次暂态电流（功率）标幺值计算
（1）短路次暂态电流标幺值（）（取）（kA）公式⑦基准容量；基准电压(kV)。（2）冲击电流（）的计算（kA）公式⑧（3
）短路容量的计算（）公式⑨2.4各元件电抗标幺值1.电力系统等值电路如图2各元件电抗标幺值的计算。设基准容量；基准电压。（1
）发电机电抗标幺值由公式①得;（2）变压器电抗值标幺值由公式③得；；（3）线路电抗标幺值由公式④得；；（4）负载电抗标幺值由公式
②得（5）电动机电抗标幺值由公式⑥得3.等值简化电路图。（1）等值电路简化过程如图2和图3所示。（2）考虑电动机的影响后，短路点
的等值电抗为2.5三相短路电流及短路功率短路次暂态电流标幺值：短路次暂态电流有名值：冲击电流：短路功率：导纳矩阵形成与计算（
Y矩阵法）3.1导纳矩阵等值电路导纳矩阵等值电路如图4所示，节点数为⑥，电抗标幺值参考图2。3.2导纳计算公式公式⑩式中
3.3变压器变比的定义式中变压器变比标幺值3.4导纳矩阵的形成对地电纳Y=3.5三相短路电流及短路功率短路点
的电抗标幺值为短路点次暂态短路电流为短路点次暂态短路电流有名值为(kA)短路点冲击电流为(kA)短路点短路功率为(MVA)两种算法
的次暂态短路电流比较误差为ΔI=10.08-9.22=0.86(kA)计算机算法设计4.1计算机编程语言C是一门通用语言，应用的面
比较广，而MatLab是一门用于特殊用途的语言，提供了专业水平的符号计算，文字处理，http://baike.baidu.com/
view/69230.htm可视化建模http://baike.baidu.com/view/239632.htm仿真和http:
//baike.baidu.com/view/635721.htm实时控制等功能。?Matlab的基本数据单位是http://ba
ike.baidu.com/view/10337.htm矩阵，它的指令表达式与数学,工程中常用的形式十分相似,故用MATLAB来解
算问题要比用C,FORTRAN等语言完相同的事情简捷得多.在新的版本中也加入了对http://baike.baidu.com/vi
ew/10075.htmC,http://baike.baidu.com/view/36402.htmFORTRAN,http:/
/baike.baidu.com/view/824.htmC++,http://baike.baidu.com/view/29.
htmJAVA的支持.可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用。Matlab是一个高
度集成的系统，集科学计算、图象处理、声音处理于一体，具有极高的编程效率。近年来，Matlab已经从最初的“矩阵实验室”，渗透到科学
与工程计算的多个领域，在自动控制、信号处理、神经网络、模糊逻辑、小波分析等多个方向，都有着广泛的应用。4.2程序流程见附录。4.3
程序清单及说明见附录。4.4程序结果及分析见附录。基于Matlab电力系统模型搭建与仿真5.1Matlab简介MATLAB是由
美国Mathworks公司开发的大型软件。在MATLAB软件中，包括了两大部分：数学计算和工程仿真。其数学计算部分提供了强大的
矩阵处理和绘图功能。在工程仿真方面，MATLAB提供的软件支持几乎遍布各个工程领域，并且不断加以完善。在电力系统反面有专门的工具
箱——SimPowerSystems，给从事电力研究的工作人员带来了巨大的便利。Simulink是http://baike.ba
idu.com/view/10598.htmMATLAB最重要的组件之一，它提供一个动态系统建模、仿真和综合分析的集成环境。在该环
境中，无需大量书写程序，而只需要通过简单直观的鼠标操作，就可构造出复杂的系统。Simulink具有适应面广、结构和流程清晰及仿真精
细、贴近实际、效率高、灵活等优点，并基于以上优点Simulink已被广泛应用于控制理论和数字信号处理的复杂仿真和设计。同时有大量的
第三方软件和硬件可应用于或被要求应用于Simulink。其中电力行业的专用工具箱SimPowerSystems也在电力系统分析与计
算中发挥了重大作用。电力系统元件库SimPowerSystems包括10类库元件分别是电源元件(ElectricalSource
s)、线路元件(Elements)、电力电子元件(PowerElectronics)、电机元件(Machines)、连接器元件(
Connectors)、电路测量仪器(Measurements)、附加元件(Extras)、演示教程(Demos)、电力图形用户
接口(powergui)和电力系统元件库模型(powerlib_models)。5.2总体设计5.2.1设计内容概述用Matla
b中Simulink组件的SimPowerSystems工具箱构建设计要求所给的电力系统模型，然后在电动机母线上设置三相短路，运行
该系统，通过测量元件获得该母线三相短路的短路电流。设计方案简介步骤一：打开Matlab中Simulink组件的SimPowerSy
stems工具箱；步骤二：在工具箱中找到设计要求所需的电气元件复制到Simulink界面中;步骤三：双击所复制元件打开参数设置选项
框，按设计要求设置相关参数；步骤四：将相关元件放置在适合的位置并连线构成一个完整的系统。步骤五：加入测量元件并设置节点的三相短路，
然后运行电路，获得短路电流。5.3详细设计在Matlab的SimPowerSystems工具箱中利用相关的电力元件搭建仿真模型并
运行系统，设置短路点得出访真结果。5.3.1电力元件设计三相电源（1）元件模型三相电源(3-PhaseSource)元件的模型如
图5-1所示，见附录。图5-1三相电源（2）功能说明三相电源是电路设计中常见的电路元件。三相电源是电路中最重要的元件，它的运行
特性对电力系统的运行状态起决定性的影响。三相电源元件提供了带有串联RLC支路的三相电源。（3）参数设置双击三相电源元件，则弹出三相
电源参数设置对话框，如图5-2所示，见附录。图5-2三相电源参数设置对话框变压器元件在电力系统电路中，变压器元件用来实现各种电
路中的变压器。在变压器元件中包括6种元件，在本次仿真计算中采用三相变压器（双绕组）元件。(1)元件模型，如图5-3所示，见附
录。图5-3三相变压器（双绕组）元件(2)参数设置双击三相变压器（双绕组）元件，则弹出参数设置对话框，如图5-4所示，见附录。图
5-4三相变压器（双绕组）元件输电线路输电线路有两种包括双回线和单回线。(1)元件模型，如图5-5所示，见附录。图5-5输电线路元
件(2)参数设置双击三相输电线路元件，在三相输电线路元件对话框中进行设置，如图5-6所示，见附录。图5-6三相输电线路元件参数设
置选项框5.3.2电力系统模型的搭建从电源元件库选择三相电源(3-PhaseSource)元件，复制后粘贴在电路图中。步骤一：
将交流电流源元件依次编号G1、G2、G3。步骤二：双击交流电流源元件，在交流电流源元件参数对话框中如下设置：峰值振幅(Peak
Amplitude)：10KV初始相位(Phase)：0频率(Frequency)：50内部连接方式（internalco
nnection）:Y测量(Measurements)选项：选择不测量单击OK按钮完成对三相电源(3-PhaseSource
)元件的设置。从线路元件(Elements)库中选择变压器（双绕组）元件，复制后粘贴在电路图中。步骤一：将交流变压器元件依次编号
T1、T2、T3。步骤二：双击变压器元件，在变压器元件参数对话框中如下设置：容量（Nominal）:根据系统选择频率(Frequ
ency)：50各绕组参数设置：根据系统选择双绕组内部连接方式（internalconnection）:Yn(ABCn)D1
1(abc)测量(Measurements)选项：选择不测量单击OK按钮完成对变压器元件的设置。从线路元件(Elements
)库中选择输电线路（双回线型及单回线型）元件，复制后粘贴在电路图中。步骤一：将输电线路元件元件依次编号L120、L150、L185
。步骤二：双击输电线路元件，在输电线路元件参数对话框中如下设置：频率(Frequency)：50单位长度阻抗（Resistan
ceperunitlength）：[0.012730.3864]单位长度电感(Inductanceperunitl
ength)：[0.9337e-34.1264e-3]单位长度电容(Capacitanceperunitlength
)：[12.74e-97.751e-9]线路长度：100Km。在本设计中，三相π型输电线路元件的参数如下设置：频率：50正序阻
抗和零序阻抗：[0.012730.3684]正序电感和零序电感：[0.9337e-34.1264e-3]正序电
容和零序电容：[12.74e-97.751e-9]线路长度(LineLength)：100Km。测量(Measurem
ents)选项：选择不测量。下面对测量元件进行设计。从电路测量仪器(Measurements)库中选择电流计(CurrentMe
asurement)元件，复制后粘贴在电路图中。将电流计元件名称改为I，表示对电源电流进行测量。双击电流计元件，在弹出的电流计
参数设置对话框中看到，输出的信号为电流幅值。选择接地(Ground)元件、节点等，进行合理放置，对上述已经放置好的电气元件进行接线
即可完成电路图的绘制。注意在接线时，接线端点的提示，如果接线错误，提示颜色显示为红色。电路连线完成后所得电路图及各子系统接线图见附
录所示（阴影部分表示含有子系统）。5.3.3电力系统模型三相短路仿真当电路图设计完成后，对其进行仿真，以达到观察交流电路中暂态变化
情况。1.在电路图的菜单选项中，选择仿真(Simulation)菜单，激活仿真参数(Simulationarameters)命
令，弹出仿真参数对话框，如图5-7所示，见附录。图5-7仿真参数对话框根据对暂态过程时间的估算，对仿真参数进行如下设置。开始时
间(Starttime)：0s停止时间(Stoptime)：0.0001s求解程序类型(Type)选项：可变步长(Vari
able－step)，ode23(Bogacki－Shampine)最大步长(Maxstepsize)选项：自动(auto)
最小步长(Minstepsize)选项：自动(auto)初始步长(Intialstepsize)选项：自动(auto)
相对容差(Relativetolerance)选项：1e-32.仿真结果加入测量元件并设置三相短路点运行系统，仿真结果如图
5-8所示，见附录。图5-8节点三相短路仿真结果6．结论1．解析法短路点的电抗标幺值为短路点的次暂态短路电流为2.Y矩阵法短
路点的电抗标幺值为短路点的导纳标幺值为短路点的次暂态短路电流为3.优缺点（1）解析法误差大，每一短路处需要逐一分析与计算。（
2）Y矩阵计算时考虑对地电容，变压器实际变比，则误差小；Y矩阵对角元素将各节点的等值短路电抗（阻抗）均求出；使分析其他点的短路故障
提供了更容易更直观的参数值；Y矩阵程序通用性强等特点。（3）两种分析与计算三相短路故障的各参数结果，如图6所示。7．致谢为期半个月
的电力系统分析课程设计，我以本设计题目入手，一开始没有任何思路和方向。但通过老师的指点，我明白了设计的关键点在于短路电流的计算，难
点在于程序的设计和仿真。依据老师推荐的方法，从图书馆借阅了有关于电力系统计算、Matlab电力系统建模以及仿真的书籍，认真研究并通
过Matlab做出了电力系统初步模型，但对于系统仿真还存在着许多问题。通过老师进一步的讲解以及对于Matlab电力系统元件库的进一
步学习，更加熟悉了电力系统模型的原理及仿真的过程。通过本次设计使我掌握了电力系统综合知识的运用，了解了现代电力系统分析的软件运行及
计算机建模仿真。在此，对老师在本次课程设计中对我的指导与鼓励表示衷心的感谢！8．参考文献著作：[1]作者、译者、书名、版本、出版
地、出版社、出版时间。期刊：[序号]作者、译者、文章题目、期刊名、年份、卷号（期数）、引用部分起止页。9．附录9.1系统等值
电路图9.2计算机算法设计流程图↓输入数据↓形成节点导纳矩阵/支路阻抗矩阵选择故障点f↓用公式计算短路电流
↓用公式计算个节点电压↓用公式计算指定支路的电流↓↓输出结果9.3计算机算法设计程序清单程序1清单：（C语言）本
程序的功能是计算各元件电抗标幺值。#include#include#defineSB1
00#defineUav115voidmain(){doubleXdG,SnG,CosG,X1,X2,UkT,SnT,X
3,X4,X5,P,Q,Uf,X9,X10,Pm,CosM,x,b,l,X6,X7,X8,Y1,Y2,Y3,Sl,Ufb,Slb;
//定义各支路电抗X1-X10,发电机次暂态电抗XdG,额定容量SnG,变压器百分数UkT,额定容量SnT,负荷功率P,Q,等指电
压Uf,电动机功率Pm,功角因数CosM,线路电阻x,电容b,长度l,对地电纳Y1-Y3.//发电机G2电抗标幺值：printf(
"\nPleaseinputtheparameterofgenerator\n请输入发电机G2参数:XdG,SnG
,CosG\n\n1:");//电抗标幺值定义，发电机电抗。scanf("%lf,%lf,%lf",&XdG,&SnG,&Co
sG);X1=XdGSB/(SnG/CosG);//发电机G3电抗标幺值：printf("\npleaseinputthe
parameterofgenerator\n请输入发电机G3参数:XdG,SnG,CosG\n2:");//电抗标幺
值定义，发电机电抗。scanf("\n%lf,%lf,%lf",&XdG,&SnG,&CosG);X2=XdGSB/(SnG/C
osG);//变压器T1-T3电抗标幺值：printf("\nPleaseinputtheparameteroftra
nsformer\n请输入变压器T1参数:UkT,SnT\n\n1:");scanf("%lf,%lf",&UkT,&SnT)
;X3=(UkT/100)SB/SnT;printf("\n2:");scanf("%lf,%lf",&UkT,&SnT);
X4=(UkT/100)SB/SnT;printf("\n3:");scanf("%lf,%lf",&UkT,&SnT);X5
=(UkT/100)SB/SnT;//负载电抗标幺值：printf("\nPleaseinputtheparamete
rofload\n请输入负载参数:Uf,PQ\n\n");scanf("%lf,%lf,%lf",&Uf,&P,&Q)
;Sl=sqrt(PP+QQ);Ufb=Uf/Uav;Slb=Sl/SB;X9=(UfbUfb)/(SlbSlb)(Q/
SB);//负载电动机电抗标幺值：printf("\nPleaseinputtheMofload\n请输入电动机负
载参数:Pm,CosM\n\n");scanf("%lf,%lf",&Pm,&CosM);X10=1/6.5100/(Pm/C
osM);//线路L1-L3电抗标幺值及对地电纳Y1-Y3：printf("\nPleaseinputtheparamete
rofline\n请输入线路参数:x,b,l\n");printf("\nForexample:x=0.4,b=2.7
4e-6,l=100");printf("\nL1:");scanf("%lf,%lf,%lf",&x,&b,&l);X6=xl
SB/(UavUav);Y1=(double)1/2bl(UavUav)/SB;printf("\n2:");scan
f("%lf,%lf,%lf",&x,&b,&l);X7=xlSB/(UavUav);Y2=(double)1/2bl
(UavUav)/SB;printf("\n3:");scanf("%lf,%lf,%lf",&x,&b,&l);X8=xl
SB/(UavUav);Y3=(double)1/2bl(UavUav)/SB;printf("\n各元件电抗标幺
值计算结果如下：");//输出发电机电抗标幺值X1,X2：printf("\nX1=%lf.X2=%lf\n",X1,X2);/
/输出变压器电抗标幺值X3,X4,X5printf("\nX3=%lf,X4=%lf,X5=%lf\n",X3,X4,X5)
;//输出线路电抗标幺值x6,x7,x8：printf("\nX6=%lf,X7=%lf,X8=%lf\n",X6,X7,X8)
;//输出负载电抗标幺值X9,X10:printf("\nX9=%lf,X10=%lf\n",X9,X10);//输出对地电纳Y
1,Y2,Y3printf("\nY1=%lf,Y2=%lf,Y3=%lf\n",Y1,Y2,Y3);printf("\nTh
eEnd!输入任意键退出\n\n");getchar();}程序1结果分析：程序2清单：（Matlab语言）>>%本程序的
功能是计算三相短路的短路电流，各节点电压，各支路电流NF=input(''请输入短路点的数目：NF='');n=input(''请输
入独立节点数：n='');nl=input(''请输入支路数：nl='');B=input(''请输入由线路参数形成的矩阵：B='');
V0=input(''请输入由各节点的初电压标幺值形成的列矩阵：V0='');D=input(''请输入由短路号，短路点阻抗组成的矩
阵:D='');m=0;Z=zeros(n);V=zeros(n);I=zeros(nl);fork1=1:nlp=B(k
1,1);q=B(k1,2);ifB(k1,6)==0k=1./B(k1,5);elsek=B(k1,5);
endifp==0ifq>mZ(q,q)=B(k1,3);m=m+1;elsefori=1:mZ(i
,m+1)=-Z(i,q);Z(m+1,i)=-Z(q,i);endZ(m+1,m+1)=Z(q,q)+B(k1,3);
fori=1:mforj=1:mZ(i,j)=Z(i,j)-Z(i,m+1)Z(m+1,j)./Z(m+1,m+1
);endZ(i,m+1)=0;endfori=1:m+1Z(m+1,i)=0;endendel
seifq>mfori=1:mZ(i,q)=Z(i,p)k;Z(q,i)=Z(p,i)k;endZ(q,
q)=k^2Z(p,p)+k^2B(k1,3);m=m+1;elsefori=1:mZ(i,m+1)=kZ
(i,p)-Z(i,q);Z(m+1,i)=kZ(p,i)-Z(q,i);endZ(m+1,m+1)=k^2Z(p
,p)+Z(q,q)-2kZ(p,q)+k^2B(k1,3);fori=1:mforj=1:mZ(i,j)=
Z(i,j)-Z(i,m+1)Z(m+1,j)./Z(m+1,m+1);endZ(i,m+1)=0;endfor
i=1:m+1Z(m+1,i)=0;endendendendfork=1:NFI(D(k,1),D(k
,1))=V0(D(k,1),1)./(Z(D(k,1),D(k,1))+D(k,2));ft=num2str(D(k,1))
;ts1=(''点短路时'');ts2=(''电流的标幺值If='');dn=strcat(ft,ts1,ts2);dis
p(dn);disp(I(D(k,1),D(k,1)));fori=1:nV(i,i)=V0(i,1)-I(D(k,1),D(k,1))Z(i,D(k,1));endfori=1:nlifB(i,6)==0k=B(i,5);elsek=1./B(i,5);endp=B(i,1);q=B(i,2);ifp==0e=0;b=B(i,3);I(i,i)=(e-V(q,q)./k)./b;elseI(i,i)=(V(p,p)-V(q,q)./k)./B(i,3);endenddisp(''各节点的电压标幺值U（节点号从小到大）依次为：'');fori=1:ndisp(V(i,i));enddisp(''各支路电流标幺值I(顺序同B)依次为：'');fori=1:nldisp(I(i,i));endend程序2结果分析：请输入短路点的数目：NF=1请输入独立节点数：n=6请输入支路数：nl=9请输入由线路参数形成的矩阵：B=[016.62010;030.3030.018310;040.208010;060.204010;121.05010;230.2980.009310;340.165010;350.3090.016710;560.656010;]请输入由各节点的初电压标幺值形成的列矩阵：V0=[1;1;1;1;1;1]请输入由短路号，短路点阻抗组成的矩阵:D=[1,0]1点短路时电流的标幺值If=1.1034各节点的电压标幺值U（节点号从小到大）依次为：011111各支路电流标幺值I(顺序同B)依次为：0-3.3003-4.8077-4.9020-0.95240000>>9.4Matlab电力系统模型仿真图2 |
|
