电力系统数字仿真作业1
电力系统分析仿真实验报告模板
电力系统分析仿真实验报告模板.doc
电力系统分析仿真实验报告模板电力系统综合分析实验一数据录入和单线图的绘制
一、实验目的通过输入9点系统的基础数据并绘制其单线图熟悉PSASP软件的使用，并未今后的分析提供基础数据。
二、实验内容1.进入PSASP主画面点击开始\程序\电力软件
\PSASP6.2\电力系统综合分析程序，弹出PSASP封面后，按任意键，即进入PSASP主画面，在该画面中完成如下工作1指定数据目录第一次可以通过“创建数据目录”按钮，建立新目录，以后可通过“选择数据目录”按钮，选择该目录。
2给定系统基准容量在系统基准容量项中，键入该系统基准容量100MV A。
建立了数据之后，该数不要轻易改动。
2.文本方式下的数据输入在PSASP主画面中点击“文本支持环境”按钮，进入文本支持环境。
点击“数据”，下拉菜单中选择“基础数据”，下拉第二级菜单1建立母线数据点击“母线”，弹出母线数据录入窗口，在窗口中依次录入该系统的母线数据。
其中母线名和基准电压必填。
该系统的母线数据如下母线名基准电压区域号电压。
PSC的电力系统仿真大作业
电力系统分析课程报告姓名学院自动化与电气工程学院专业控制科学与工程班级指导老师二〇一六年五月十三一、同步发电机三相短路仿真1、仿真模型的建立选取三相同步发电机模型,以三相视图表示;励磁电压和原动机输入转矩Ef 与Tm均为定常值,且发电机空载;当运行至时,发电机发生三相短路故障;同步发电机三相短路实验仿真模型如图1所示;图1 同步发电机三相短路实验仿真模型2、发电机参数对仿真结果的影响及分析衰减时间常数Ta对于直流分量的影响三相短路电流的直流分量大小不等,但衰减规律相同,均按指数规律衰减,衰减时间常数为Ta,由定子回路的电阻和等值电感决定大约;pscad同步发电机模型衰减时间常数Ta对应位置如图3所示当前Ta=;图3 同步发电机模型参数Ta对应位置1Ta=时,直流分量的衰减过程以励磁电流作为分析如图4所示;图4 Ta=发生短路If波形2Ta=时,直流分量的衰减过程以励磁电流作为分析如图5所示;图5 Ta=发生短路If波形短路时刻的不同对短路电流的影响由于短路电流的直流分量起始值越大,短路电流瞬时值就越大,而直流分量的起始值于短路时刻的电流相位有关,即直流分量是由于短路后电流不能突变而产生的;Pscad模型中对短路时刻的设置如图6所示图6 Pscad对于短路时刻的设置1当在t=时发生三相短路,三相短路电流波形如图7所示;图7 t=时三相短路电流波形2）当在t=时发生三相短路,三相短路电流波形如图8所示;图8 t=6时三相短路电流波形Xd、Xd`、Xd``对短路电流的影响1 Xd的影响Pscad中对于Xd的设置如图9所示：图9 Pscad对于D轴同步电抗Xd的设置下面验证不同Xd时A相短路电流的稳定值;i.Xd=标幺制,下同时,仿真波形如图10所示图10 Xd=时A相短路电流波形ii.Xd=10时,仿真波形如图11所示图11 Xd=时A相短路电流波形2Xd`的影响在Pscad中暂态电抗Xd`的设置如图13所示：图13 Pscad对于暂态电抗Xd的设置下面验证不同Xd`时A相短路电流的暂态过程;i.Xd`=时A相短路电流的波形如图14所示：图14 Xd`=时A相短路电流波形ii.Xd`=1时A相短路电流的波形如图15所示：图15 Xd``=1时A相短路电流波形3Xd``的影响这里次暂态电抗Xd``与暂态电抗Xd`相似,Xd``影响的是短路后的次暂态过程;在Pscad中次暂态电抗Xd``的设置如图16所示：图16 Pscad对于次暂态电抗Xd``的设置下面验证不同Xd``时A相短路电流的暂态过程;i.Xd``=时A相短路电流的波形如图17所示：图17 Xd``=时A相短路电流波形ii.Xd``=时A相短路电流的波形如图18所示：图18 Xd``=时A相短路电流波形Td`、Td``对短路电流的影响在Pscad中衰减时间常数的设置如图19所示：图19 Pscad对于衰减时间常数的设置1）下面验证不同Td`时A相短路电流暂态交流分量衰减速度;i.Td`=时短路励磁电流的波形如图20所示图20 Td`=时短路励磁电流的波形ii.Td`=时短路励磁电流的波形如图21所示图21 Td`=时短路励磁电流的波形2）下面验证不同Td``时A相短路电流暂态交流分量衰减速度;i.Td``=时短路励磁电流的波形如图22所示图22 Td``=时短路励磁电流的波形ii.Td``=时短路励磁电流的波形如图23所示图23 Td``=时短路励磁电流的波形二、简单电力网络的线路故障仿真1、仿真模型的建立仿真模型预览图根据题目要求,建立如图24所示的仿真模型;图24 简单电力网络仿真模型i.三相电源参数、输电线路π型等值电路参数对应在仿真模型中的设置位置如图25ab所示;图25a 电源参数设置位置注：题目要求中两侧电源电势夹角为15°,即两侧Phase值相差15°;图25b π型等值电路参数设置位置ii.变压器参数设置如图26所示,变压器采用星三角连接且不接地,零序电流不流通;图26 变压器模型参数设置位置2、输出通道的设计由于题目要求三相电压电流的瞬时波形和向量图,所以需要设计不同的输出通道;瞬时波形用覆盖图形和分离图形显示如图27所示；相量图需要将各相电气量经过FFT模块处理如图28所示,采用其幅值相位模式,用处理完成后的基波表示如图29所示,用向量仪显示效果如图30所示;图27 瞬时波形显示方式图28 FFT处理模块图29 FFT模式选择图20 相量仪输出幅相特征3、不同故障类型下的仿真结果线路发生三相故障由于是对称故障且系统两侧基本相同,这里只看M侧A相;故障前幅相特性故障后幅相特性线路发生A相接地故障主要看特殊相A相的电气量变化情况M侧A相母线电压波形M侧母线ABC相电流波形短路点A相电压波形短路点ABC三相电流波形故障前重合;故障前幅相特性故障后幅相特性线路发生BC两相故障仍然主要看特殊相A的电气量变化母线M侧A相电压波形短路点ABC三相电压波形故障前重合,故障后BC两相电压减小短路点ABC三相电流波形短路前的幅相特性短路后的幅相特性4课程学习心得通过对于电力系统分析这门课一个学期的学习,我更加深入的了解了电力系统分析中主要考虑的问题;本学期主要是对电力系统的暂态的学习,通过刘益青老师的耐心讲解和课下用PSCAD仿真软件对课程的巩固,使我自己在故障分析的能力上又上了一个台阶;对于课程的建议：可能由于课时的原因,课程进行的略微紧张,基础稍欠缺的地方有些跟不上;因为了解到老师有做过故障保护相关工作的经验,所以希望老师以后能为学生拓展一些电力系统其他方面的知识;在此再次感谢老师在学习上的支持和帮助;。
电力系统建模仿真作业1
无穷大功率电源供电系统仿真假设无穷大功率电源供电系统，在0.02s时刻变压器低压母线发生三相短路故障，仿真其短路电流周期分量幅值和冲击电流的大小。
线路参数L=50km，x1=0.4Ω/km，r1=0.17Ω/km；变压器Sn=20MV·A，短路电压Us%=10.5，短路损耗ΔPs=135kw，空载损耗ΔP0=22kw，空载电流I0%=0.8，变比kT=110/11，高低压绕组均为Y行联接；并设供电点电压为110KV。
其对应的Simulink仿真模型如图1-1所示。
图1-1 无穷大功率电源供电系统的Simulink仿真图表1-1 图1-1仿真电路中各模块名称及提取路径模块名提取路径无穷大功率电源Three-Phase Source SimPowerSystems/Eletrical Sources三相并联RLC负荷模块5MW SimPowerSystems/Elements串联RLC支路Three-phaseParallelRLCBranch SimPowerSystems/Elements三相故障模块Three-phase-Fault SimPowerSystems/Elements三相电压电流测量模块V-I-M SimPowerSystems/Measurements示波器模块Scope Simulink/Sinks电力系统图形用户界面Poweigui SimPowerSystems双绕组变压器模块Three-PhaseTransformer SimPowerSystems/Elements图1-2 电源模块的参数设置变压器T 采用“Three-PhaseTransformer （Two Windings ）”模型。
根据给定的数据，计算折算到110kv 侧的参数如下：变压器的电阻为2233221351101010 4.0820000s N T N PU R S ∆⨯=⨯=⨯Ω=Ω 变压器的电抗为22332%10.5110101063.5310010020000s N T N U U X S ⨯=⨯=⨯Ω=Ω⨯ 变压器的漏感：63.53/(2)0.2022 3.1450T T L X f H H π===⨯⨯变压器的励磁电阻为2233301101010 5.51022N m U R P =⨯=⨯Ω=⨯Ω∆ 变压器的励磁电抗为22330100100110101075625%0.820000N m N U X I S ⨯=⨯=⨯Ω=Ω⨯ 变压器的励磁电感为75625/(2)240.82 3.1450m m L X f H H π===⨯⨯变压器模块中的参数采用有名值则设置如图1-3所示图1-3采用有名值时变压器模块的参数设置如果要采用标幺值，则在Similink 的三相变压器模型中，一次、二次绕组漏感和电阻的标幺值以额定功率和一次、二次侧各自的额定线电压为基准值，励磁电阻和励磁电感以额定功率和一次额定线电压为基准值。
(完整版)电力系统仿真
电力系统分析与设计例子：电力系统可视化仿真介绍EXAMPLE1-1：题目：双总线电力系统初始条件：总线1电压为16kV，总线2为15.75KV，负载功率为5MW，发电机功率为5.1MW。
总线1与总线2之间由一条传输线连接。
实验步骤：保持其他参数不变，依次调节负载功率参数，观察其他参数的变化。
实验现象：①当负载功率为5MW时，发电机的输出功率为5.1MW。
②当负载功率调整为6MW时，发电机的输出功率为6.1MW。
③当负载功率调整为4MW 时，发电机的输出功率为4.0MW。
实验结论：在双总线电力系统中，当其他线路装置参数不变时，负载功率增大时，发电机的输出功率相应增大，负载功率减小时，发电机的输出功率相应减小。
EXAMPLE1-2：题目：植入新的总线初始条件：在上图中保持其他条件不变，植入新的总线”Bus3”。
实验步骤：在powerworld选择edit mode，在Draw中选择Network---bus，将”Bus”放置图中，双击”Bus”，将对话框中的名称改为”Bus3”，电压改为16kV。
实验结果：如下图所示EXAMPLE1-3：题目：三总线电力系统初始条件：在EXAMPLE2的基础上，通过传输线路将Bus1和Bus2与Bus3连接在一起。
实验步骤：在edit mode下，选择draw选项，选择Network中的transmission line，单击Bus1，然后将线路连接到Bus3，双击完成连接。
并调节字体大小和线路的颜色。
在Network中选择load选项，选择load的大小。
最后把系统名字改为Three Bus Powr system。
实验结果：如下图所示②对新系统进行调节参数实验：实验步骤：⑴调节新总线Bus3下负载参数，观察对其它参数的影响：①当负载功率为11MW时，如图②当负载功率为9MW时，如图实验结论：当Bus3下负载功率增大时，Bus2和Bus3上的电压降低，发电机的输出功率增大；当Bus3下负载功率减小时，Bus2和Bus3上的电压增大，发电机的输出功率变小。
电力系统建模仿真作业
风电并网后静态电压稳定性分析的建模与仿真电力系统经常采用P-V曲线分析法来分析有关静态电压稳定性的问题,P代表穿越传输断面传送的功率或者一个区域的总负荷,V代表代表性节点或关键节点的电压。
P-V曲线分析法即是建立一个区域负荷或者传输界面潮流和节点电压之间的关系曲线,从电力系统当前的稳定运行点开始,通过不断增加P,使用潮流计算,描出代表节点的电压变化曲线,用P-V曲线的拐点来表示区域负荷或者传输界面功率的增加导致整个系统临界电压崩溃的程度,即系统静态电压稳定极点。
在把P-V曲线法用于研究风电的接入对电压静态稳定性的影响时,P代表的是风电场输出的有功功率,V为机端电压、风电接入点电压（PCC电压）等其他需要监测的母线电压。
实际上，P-V曲线法是在静态情况下,研究风速变化导致的风电场输出有功功率的变化对电网电压的影响。
用风电输出的有功功率引起的电压水平的变化及当前运行点到电压崩溃点的“距离”,反映风电接入的电网的电压稳定裕度。
在求取风电接入系统的P-V曲线时 ,除了系统平衡节点外,一般不考虑网内其他常规机组的有功功率的变化以及网内负荷的变化情况。
综上,电网基于静态电压稳定性的风电接纳能力,即是以电网的静态电压稳定性作为约束条件,在保证电网静态电压稳定的基础上尽可能多接入风电。
通常系统静态电压越限临界点所接入的风电容量即为系统可接纳的最大风电并网容量。
1算例本文通过IEEE14节点标准测试系统作为算例，风电场通过变压器和110 kV 线路接入IEEEl4节点标准测试系统的14号节点，使用以上算法对基于静态电压稳定性下的一风电场的并网功率极限进行计算。
风电场IEEE14节点系统110kv线路图2.2 风电场接入IEEE14系统图图中变压器标幺变比取1（在实际运行中，可以通过改变变压器的分接头来调控特定节点的电压），风电场接入系统的线路参数为12.6+j24.96Ω。
本文基于双馈感应风机的风电场进行电压静态稳定约束下接纳能力计算。
电力系统潮流计算数字仿真实验实验报告（PASAP软件的使用）
电力系统潮流计算数字仿真实验实验报告（PASAP软件的使用）电力系统潮流计算数字仿真实验实验报告班级姓名学号五节点潮流计算（原图） (3)不同负荷水平时刻的潮流计算——负荷加倍 (6)各节点负荷减半作为谷时 (8)两种不同运行方式的潮流计算 (10)两种不同调度结果的潮流计算 (12)2一、五节点潮流计算（原图）结果导出：①物理母线单位：kA\kV\MW\Mvar区域分区厂站全网全网全网母线名称电压幅值电压相角-------- -------- --------* b 103.48250 -8.0217g 217.28740 -5.1551x 106.31940 -10.2418Ⅰ 242.00000 0Ⅱ 116.39210 -2.35793②发电机单位：kA\kV\MW\Mvar区域分区厂站全网全网全网发电机名称母线名类型有功发电无功发电功率因数---------- ------ ---- -------- -------- -------- Gen_1 Ⅰ Vθ 202.4800 155.9800 0.79220 Gen_2 Ⅱ PQ 40.0000 30.0000 0.80000 Gen_3 g PQ0.0000 10.0000 0.00000③负荷结果报表单位：kA\kV\MW\Mvar区域分区厂站全网全网全网负荷名称母线名类型有功负荷无功负荷功率因数-------- ------ ---- -------- -------- -------- Load_1 b PQ 50.0000 30.0000 0.85749 Load_3 x PQ 180.0000 100.00000.87416④交流线结果报表单位：kA\kV\MW\Mvar区域分区全网全网交流线名称 I侧母线 J侧母线 I侧电压 I侧有功 I侧无功---------- ------- -------- ------- ------- ------- AC_1 Ⅰ g 242.00000 180.4800 162.8100 AC_3 Ⅱ b 116.39210 61.9700 22.1800 AC_9 bx 103.48250 6.7700 -15.8100 J侧电压 J侧有功 J侧无功------- ------- -------217.28740 174.5300 131.0300103.48250 56.7700 14.1900106.31940 6.3200 -16.50004⑤两绕组变压器结果报表单位：kA\kV\MW\Mvar区域分区厂站全网全网全网两绕组变压器名称 I侧母线 J侧母线 I侧有功---------------- ------- ------- -------T2w_1 ⅠⅡ 21.9900 T2w_2 g x 174.5300I侧无功 J侧有功 J侧无功------- ------- --------6.8200 21.9600 -7.8200141.0300 173.6800 116.50005二、不同负荷水平时刻的潮流计算——负荷加倍问题：如果不增加无功功率，负荷在正常运行基础上加倍时潮流计算不成功，因为根据电压损耗公式，负荷加倍后电压损耗变大，部分母线电压过低，所以潮流计算不成功。
电力系统仿真软件介绍(作业)
电力系统分析软件介绍1 EMTDC/PSCADEMTDC是一种世界各国广泛使用的电力系统仿真软件，PSCAD是其用户界面，一般直接将其称为PSCAD。
使得用户能更方便地使用EMTDC进行电力系统分析，使电力系统复杂部分可视化成为可能。
PSCAD/EMTDC基于dommel电磁暂态计算理论，适用于电力系统电磁暂态仿真。
EMTDC（Electro Magnetic Transient in DC System）即可以研究交直流电力系统问题，又能完成电力电子仿真及其非线性控制的多功能工具。
PSCAD由Manitoba HVDC research center开发。
2 PSAPACPSAPAC由美国EPRI开发，是一个全面分析电力系统静态和动态性能的软件工具。
其包含多个模块，其中部分模块可以单独使用。
模块和功能如下：DYNRED（Dynamic Reduction Program）：网络化简与系统的动态等值，保留需要的节点。
LOADSYN（Load Synthesis Program）：模拟静态负荷模型和动态负荷模型。
IPFLOW（Interactive Power Flow Program）：采用快速分解法和牛顿－拉夫逊法相结合的潮流分析方法，由电压稳态分析工具和不同负荷、事故及发电调度的潮流条件构成。
TLIM（Transfer Limit Program）：快速计算电力潮流和各种负荷、事故及发电调度的输电线的传输极限。
DIRECT：直接法稳定分析软件弥补了传统时域仿真工作量大、费时的缺陷，并且提供了计算稳定裕度的方法，增强了时域仿真的能力。
LTSP（Long Term Stability Program）：LTSP是时域仿真程序，用来模拟大型电力系统受到扰动后的长期动态过程。
为了保证仿真的精确性，提供了详细的模型和方法。
VSTAB（V oltage Stability Program）：该程序用来评价大型复杂电力系统的电压稳定性，给出接近于电压不稳定的信息和不稳定机理。
电力系统实验十二大电网数字仿真
实验十二 大电网数字仿真1、实验目的：本实验针对我国某省级实际大电网的全数字动态模型来进行。
通过实验研究对比不同故障对大电网暂态稳定性的影响情况；对特定故障，研究继电保护动作时延对大电网暂态稳定性的影响，并用试探法求取故障极限切除时间。
2、实验内容：1)在指定的同一条220kV输电线路上，分别设置三相短路和两相短路，故障延时相同，观察系统各发电机的摇摆曲线图有何区别。
2)针对指定的5条不同的220kV输电线路，分别设置三相短路，故障延时相同，观察系统各发电机的摇摆曲线图有何区别。
3)用试探法求取指定220kV输电线路发生三相短路故障时的故障极限切除时间。
3、实验步骤：1.不同故障类型对大电网暂态稳定性的影响a) 启动我国某省级实际大电网的全数字动态仿真系统，建立初始运行方式和潮流。
然后，在不自动模拟继电保护的情况下，针对指定的220kV输电线路，故障设置为三相短路，故障延时设置为100ms，故障位置设置为距线路首端50%，启动动态仿真，查看各发电机的摇摆曲线，并保存之。
b) 仿真系统返回到初始状态。
针对同一条输电线路，故障设置为AB两相短路，故障延时设置为100ms，故障位置设置为距线路首端50%，启动动态仿真，查看各发电机的摇摆曲线，并保存之。
2.不同故障地点对大电网暂态稳定性的影响在不自动模拟继电保护的情况下，分别针对指定的5条不同的220kV输电线路，故障设置为三相短路，故障延时设置为100ms，故障位置设置为距线路首端50%，在仿真系统返回到初始状态后启动动态仿真，分别查看各发电机的摇摆曲线，并保存之。
3.测定故障极限切除时间a) 仿真系统返回到初始状态，不自动模拟继电保护。
针对指定的220kV输电线路，故障设置为三相短路，故障延时设置为50ms，故障位置设置为距线路首端50%，启动动态仿真，查看各发电机的摇摆曲线，并保存之。
b) 仿真系统返回到初始状态，针对指定的220kV输电线路，设置同样的故障类型和位置，只是故障延时增加一个Δt，比如Δt=50ms，察看各发电机的摇摆曲线，并保存之。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
电力系统数字仿真作业电力系统及其自动化研硕10-13 韩暘1 . Matlab潮流计算作业本程序的功能是用牛顿——拉夫逊法进行潮流计算。
其中：B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳；5、支路的变比；6、支路首端处于K 侧为1，1侧为0B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值；4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量；6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；3为PV节点；程序如下：n=input('请输入节点数:n=');nl=input('请输入支路数:nl=');isb=input('请输入平衡母线节电号:isb=');pr=input('请输入误差精度:pr=');B1=input('请输入由支路参数形成的矩阵:B1=');%变压器侧为1，否则为0B2=input('请输入各节点参数形成的矩阵:B2=');Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);% % %---------------------------------------------------for i=1:nl %支路数if B1(i,6)==0 %左节点处于1侧p=B1(i,1);q=B1(i,2);else %左节点处于K侧p=B1(i,2);q=B1(i,1);endY(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元Y(q,p)=Y(p,q); %非对角元Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2; %对角元K侧Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; %对角元1侧end%求导纳矩阵disp('导纳矩阵Y=');disp(Y)e(1)=1.00;e(2)=1.00;e(3)=1.10;e(4)=1.05;f(1)=0;f(2)=0;f(3)=0;f(4)=0;G=real(Y);B=imag(Y);%设S=P+Qj；S(1)=-0.30-0.18j;S(2)=-0.55-0.13j;S(3)=0.50;P=real(S);Q=imag(S);V(3)=1.10;V(4)=1.05+0j;k=0;precision=1;N1=3;while precision > 0.00001for m=1:N1for n=1:4Pt(n)=(e(m)*(G(m,n)*e(n)-B(m,n)*f(n))+f(m)*(G(m,n)*f(n)+B(m,n)*e(n))); enddP(m)=P(m)-sum(Pt);endfor m=1:3for n=1:4Qt(n)=(f(m)*(G(m,n)*e(n)-B(m,n)*f(n))-e(m)*(G(m,n)*f(n)+B(m,n)*e(n))); enddQ(m)=Q(m)-sum(Qt);endVt=e(3)*e(3)+f(3)*f(3);dV=(V(3)^2)-Vt;for m=1:3for n=1:4Bi(n)=G(m,n)*f(n)+B(m,n)*e(n);Ai(n)=G(m,n)*e(n)-B(m,n)*f(n);endendfor m=1:3for n=1:3if m==nH(m,n)=sum(Ai)-G(m,n)*e(m)-B(m,n)*f(m);N(m,n)=-sum(Bi)+B(m,n)*e(m)-G(m,n)*f(m);J(m,n)=sum(Bi)+B(m,n)*e(m)-G(m,n)*f(m);L(m,n)=-sum(Ai)+G(m,n)*e(m)+B(m,n)*f(m);J(3,3)=-2*e(3);L(3,3)=-2*f(3);elseH(m,n)=-G(m,n)*e(m)-B(m,n)*f(m);N(m,n)=B(m,n)*e(m)-G(m,n)*f(m);J(m,n)=B(m,n)*e(m)-G(m,n)*f(m);L(m,n)=G(m,n)*e(m)+B(m,n)*f(m);J(3,1)=0;L(3,1)=0;J(3,2)=0;L(3,2)=0;endendend%生成雅克比矩阵。
%对角线上的元素for m=1:3Ja(2*m-1,2*m-1)=H(m,m);Ja(2*m-1,2*m)=N(m,m);Ja(2*m,2*m-1)=J(m,m);Ja(2*m,2*m)=L(m,m);end%非对角线上的元素for m=1:3for n=1:3if m==nelseJa(2*m-1,2*n-1)=H(m,n);Ja(2*m-1,2*n)=N(m,n);Ja(2*m,2*n-1)=J(m,n);Ja(2*m,2*n)=L(m,n);endendendPQ=[dP(1);dQ(1);dP(2);dQ(2);dP(3);dV];dU=-inv(Ja)*PQ;precision=max(abs(dU));for n=1:3e(n)=e(n)+dU(2*n-1);f(n)=f(n)+dU(2*n);endk=k+1;for n=1:4Un(n)=e(n)+f(n)*j;enddisp('第')k-1disp('次迭代')disp('各节点电压')Undisp('迭代过程中节点不平衡量的变化情况') PQenddisp('计算平衡节点的功率')for m=1:4In(m)=Y(m,4)*Un(m);endS=Un(4)*sum(conj(In))输入数据为：B1=[1 2 0.1+0.4i 0.03056i 1 0;2 4 0.08+0.4i 0.02826i 1 0;3 1 0.3i 0 0.909 1 0;4 1 0.12+0.5i 0.0384i 1 0]B2=[1 2 0.1+0.4i 0.03056i 1 0;2 4 0.08+0.4i 0.02826i 1 0;3 1 0.3i 0 1. 1 1;4 1 0.12+0.5i 0.0384i 1 0]仿真结果为：请输入节点数:n=4请输入支路数:nl=>> 4请输入平衡母线节电号:isb=4请输入误差精度:pr=0.00001请输入由支路参数形成的矩阵:B1=[1 2 0.1+0.4i 0.03056i 1 0;2 4 0.08+0.4i 0.02826i 1 0;3 1 0.3i 0 0.9091 0;4 1 0.12+0.5i 0.0384i 1 0]请输入各节点参数形成的矩阵:B2=[1 2 0.1+0.4i 0.03056i 1 0;2 4 0.08+0.4i 0.02826i 1 0;3 1 0.3i 0 1.1 1;4 1 0.12+0.5i 0.0384i 1 0]导纳矩阵Y=1.0421 - 8.2428i -0.5882 +2.3529i 0 +3.6666i -0.4539 + 1.8911i-0.5882 + 2.3529i 1.0690 - 4.7274i 0 -0.4808 + 2.4038i0 + 3.6666i 0 0 - 3.3333i 0-0.4539 + 1.8911i -0.4808 + 2.4038i 0 0.9346 - 4.2616i第ans =次迭代各节点电压Un =0.9938 - 0.0118i 0.9771 - 0.1127i 1.1000 + 0.1234i 1.0500迭代过程中节点不平衡量的变化情况PQ =-0.2773-0.0509-0.52600.01960.5000第ans =1次迭代各节点电压Un =0.9848 - 0.0106i 0.9590 - 0.1088i 1.0929 + 0.1249i 1.0500 迭代过程中节点不平衡量的变化情况PQ =0.0003-0.00460.0020-0.06500.0028-0.0152第ans =2次迭代各节点电压Un =0.9844 - 0.0085i 0.9585 - 0.1083i 1.0924 + 0.1289i 1.0500 迭代过程中节点不平衡量的变化情况PQ =0.0011-0.0010-0.0026-0.00010.0065-0.0001第ans =3次迭代各节点电压Un =0.9846 - 0.0086i 0.9587 - 0.1084i 1.0924 + 0.1289i 1.0500 迭代过程中节点不平衡量的变化情况PQ =-0.00040.0015-0.00020.00040.0004-0.0000第ans =4次迭代各节点电压Un =0.9846 - 0.0086i 0.9587 - 0.1084i 1.0924 + 0.1290i 1.0500迭代过程中节点不平衡量的变化情况PQ =1.0e-004 *-0.2215-0.02110.3567-0.31260.0816-0.0000第ans =5次迭代各节点电压Un =0.9846 - 0.0086i 0.9587 - 0.1084i 1.0924 + 0.1290i 1.0500迭代过程中节点不平衡量的变化情况PQ =1.0e-005 *0.00140.0954-0.29090.88590.0524-0.0000计算平衡节点的功率S =0.3679 + 0.2647i2 . Simulink仿真作业TLf供电系统如图所示，各元件参数如下：线路L：长50km，x=0.4Ω/km；变压器T：假定供电点电压为106.5kV，保持恒定，空载运行时变压器低压母线发生三相短路。
搭建仿真电路图如下：仿真结果如下：图示为短路电流的波形图。