电力系统课程设计scut

1、电力系统运行方式分析和计算设计报告课程名称：题目：电力系统运行方式分析和计算专业：XXXXXXX班级：XXXXXXX姓 名: XXXXXX扌旨导教师： XXXXXX目录课程设计任务书 3第一章 绪论 6第二章 运行方式与 powerworld 建模 7第三章 “N1 ”潮流校核11第四章 短路计算 20第五章 课程设计总结与心得 22第六章 附录 23附录一 潮流计算程序 23附录二 短路计算程序 2929课程设计任务书课程设计题目A2:电力系统运行方式分析和计算姓名： 乞 指导教师：_x、一个220kV分网结构和参数如下：500kV站（#1 ）的220kV母线视为无穷大母线，电压恒定在 23

2、0kV图中，各变电站参数如下表：编号类型220kV最大负荷，MVA#1500kV 站平衡节点#2220kV 站230+j20#3220kV 站350+j90#4220kV 站250+j85#5220kV 站410+j90各变电站负荷曲线基本一致。日负荷曲线主要参数为: 日负荷率：0.88，日最小负荷系数：0.64各线路长度如图所示。所有线路型号均为LGJ-2*300,基本电气参数为:正序参数：r = 0.054Q/km, x = 0.308Q/km, C = 0.0116 卩 F/km零序参数：r0 = 0.204Q /km, x0 = 0.968Q 7/km, Co = 0.0078卩 Fj

3、/km40oC长期运行允许的最大电流：1190A。燃煤发电厂G有三台机组，均采用单元接线。电厂220kV侧采用双母接线。发电机组主 要参数如下表(在 PowerWorld中选择GENTRA模型)：机 组 台数单台容量(MW)额定电 压(EV)功 率 因 数升压变容量MVAXdXdXqTd0TJ=2Ha,2t/(MW 2 h)ai,1t/(MWh)ai,0t/hPmax(MW)Pmin(MW)130010.50.853501.80.181.2870.000040.29810.22300120130010.50.853501.80.181.2870.000030.30510.32300120125

4、010.50.853002.10.21.5760.000020.3219.38250100升压变参数均为Vs%=10.5%变比10.5kV/242kV不计内阻和空载损耗。稳定仿真中发电机采用无阻尼绕组的凸极机模型。不考虑调速器和原动机模型。不考虑 电力系统稳定器模型。励磁系统模型为：RegulatorStabilizer该模型在PowerWorld中为BPA_EG莫型，主要参数如下：KA=40 TA=0.1 TA仁 0.1KF=0.05TF=0.7 VRmax=3.7VRmi n=0.0发电厂按PV方式运行，高压母线电压定值为1.05Vn。考虑两种有功出力安排方式: 满发方式： 开机三台，所有

5、发电机保留10%的功率裕度； 轻载方式： 仅开250MW机组，且保留10%的功率裕度；发电厂厂用电均按出力的7%考虑。二、设计的主要内容：1、根据负荷变化和机组出力变化，拟定至少两种典型运行方式；2、进行参数计算和标幺化，形成两种典型运行方式的潮流计算参数;3、用Matlab编制潮流计算程序，可任选一种潮流计算方法；4、用所编制的潮流程序完成典型运行方式的潮流计算，进行电压和网损分析；5、用PowerWorld软件进行潮流计算并与自己编制的软件计算结果进行校核和分析；6、用所编制的潮流程序完成大方式的“ N-1 ”潮流校核，进行线路载流能力和电压水平 分析；7、用Matlab编制三相短路的短路

6、容量计算程序；8对主要220kV母线进行三相短路容量测算，并与 PowerWorld的计算结果进行校核；9、自行选择2-3种故障方案，用PowerWorld进行稳定计算，给出摇摆曲线，并计算故 障的极限切除时间。10、假定电网公司下发给燃煤发电厂 G的日发电计划曲线如下图，按照等微增率准则或 对三台机组进行经济负荷分配，同时采用 matlab中的quadprog函数对三台机组进行负 荷优化分配，并对两种分配结果进行分析比较。要求给出三台机组的日发电计划曲线。11、编制课程设计报告P/MW920 780 |II；540III0IIIIIIJ5t/h21 24600 ：Ihlik400 丨:IIi

7、IIIHPIII4iIIIII*IH0 6 8 18设计要求和设计成果：1、每两位同学为一组，自行分工，但任务不能重复；2、每位同学对自己的设计任务编写课程设计说明书一份；3、一组同学共同完成一份完整的设计报告；4、设计说明和报告应包含：以上设计任务每一部分的计算过程和结果分析；所编制的潮流和短路源程序（主要语句应加注释）；潮流计算结果（潮流图）稳定计算的功角曲线等；第一章 绪论1.1 课程设计目的 通过区域电力网络电力系统运行方式分析和计算， 巩固和加深对课程所学理 论知识的理解， 培养独立分析和解决问题的能力， 掌握电力工程运行计算设计的 基本方法和基本步骤。设计任务是根据给出的数据及要求

8、，利用 Matlab 软件及 Powerword 软件进 行程序编写及潮流计算、“ N-1”潮流校核、线路载流能力和电压水平分析、短路 容量测算等任务。并编写各自参与项目说明书，及总的报告一份。1.2 设计主要内容1 、根据负荷变化和机组出力变化，拟定至少两种典型运行方式；2、进行参数计算和标幺化，形成两种典型运行方式的潮流计算参数；3、用 Matlab 编制潮流计算程序，可任选一种潮流计算方法。4、用所编制的潮流程序完成典型运行方式的潮流计算，进行电压和网损分析；5、用 PowerWorld 软件进行潮流计算并与自己编制的软件计算结果进行校核和分 析；6用所编制的潮流程序完成大方式的“ N-

9、1 ”潮流校核，进行线路载流能力和电 压水平分析；7、用 Matlab 编制三相短路的短路容量计算程序；8、对所有220kV母线进行三相短路容量测算，并与PowerWorld的计算结果进行 校核；9、自行选择 2-3 种故障方案，用 PowerWorld 进行稳定计算，给出摇摆曲线，并 计算故障的极限切除时间。本次课程设计我负责的任务主要是使用powerworld进行该220KV系统的模型建立，以及在借鉴前人的基础上使用 matlab 编制三相短路程序，以及与 powerworld 的计算结果进行校核等。第二章 运行方式与 powerworld 建模根据负荷变化和机组出力，我们拟定了两种典型的

10、运行方式，分别是发电机满发时负荷最小和发电机轻载时负荷最大这两种运行方式。（已知最大负荷,可由日最小负荷系数求出最小负荷）2.1、节点根据题目要求建立并编号6个节点：发电厂汇总母线1个（母线6基准值220KV）,变电站母线4个（母线2/3/4/5基准值220KV），平衡节点母线1个（母线1基准值220KV ）。在母线1的信息菜单（Bus In formation Dialog）中勾选系统平衡母线选项（systemslack bus）以设置平衡节点。另外，需要为母线设置电压表（标幺值），选择节点后点击相应母线。馆圉册里器指示*H列霧网貉集合背昱帯F插入对象2.2、发电机建立3个发电机组，其中：1

11、.两个容量为300MW的发电机组联接在 & 9母线上，并如下图设置相应数据雄电机进项 养tfiE晋挾1我态书点容称9爭営傾匸1豐3D1畔*恥1VFKI labS机划粪型LJHOimwQ-显丁信目 功丰融圧控制I咸本1足路壷数I序有却IEH盼电站I agM I鮭 I功菊鉀1有功出力251.100叵启用出GC瑩与因子10.00最小有功附0.J.限最大肖办INlflM-OMl无龙出力130.跻节自fi号9最小无功斷-9900.000宕用A压电压设宦值LOWC使用君蚩抽夭囲方制节100.0MAT最小Mvar1时11 1X取消1?斜助其中，353=300085, 249=300*0.87节戌案号gu*按

12、謂号登拱狀态节点窘称弓1 丨i就-二囲絢51媾料类型UdmwiT【帚，“ 1nnlate 呂机组夷型UH (UfiknoMdri)矗示信罷I M車书砲压觀1丨騙本 理蹈變故丨斯有皆4费讨区悭电站丨自走次曙定掠存*蜩|?笨肋卄压芟压盘1轨端和磁茎救 Lr赵岌电机 100.0世刑AjHhab此PUvflbesed usng device base qF 1烫,业 MvA千 |Rao.aooo :咖 0.1600斗Tdope.aooo :5(1.0)0.0000 土H3.0000 :5(1.2)0.0000 *DD.dooo :Af 0.2500 二Xdl.aoao ；X1.2OT0 :上图数据按照

13、给出参数输入KaW.OOOOt0.1QQ0TaiC. 1ODDEfdmna DODOE fid max3.7000Kf0.0500Tf 0.7000上图数据按照300MW发电机组给出的参数输入2. 一个250MW的发电机组联接在母线8上，参数设置如下图世电现远琐节曲n号匚列册号直直3杭粧节点电称 6I扶屯解迈撫一IID II 査找1Unlonown&s.rw3 labelsUN CUnknawin) 圭保有X取沽?帮肋其中，294=250085 207=250*0.87按照250MW的发电机组给出的参数输入2.3线路线路的参数需与线路的长度对应，信息对话框中有相应的换算工具其中，电阻和电抗值题

14、中已给出，电容值需计算后输入（B）零序相应数据编号区域拒称基准电压ii百号向r点11(1)220.0标釜 r labels1220.0I 按编号蛋找 按名毬直找 查找显示锻SH踏參数断有者/区分区糜电站自宦丈垂序阻抗R:C.00752SX:0.029510C :0.047600零序井联导纳苜端G;0.000000ttsHa：0.000000末韓G:o.oooaoo末端日：0.0000000注意：谨甌置仅:夬定克压器绩组的接地方式 必须输入移相器用度作为我压蓋腔制埶拥0?帮助X取消2.4负荷根据所选定的运行方式选择不同的负荷，最终形成下图:V10S0 MW-23 Mvar1.0355 pu1.0

15、317 pu1,0500 pu350 MW &0 Mvar230 MW20 Mvar410 MW90 Mvgr1.0455 pu250 MW85 Mv日1.0340 pu249 MW9D Mvar249 MW90 Mvar208 MW374 MVar第三章 “ N-1”潮流校核3.1 matlab 计算N-1的处理方法：本课程设计“ N-1”指是在发电机满发方式，负荷重载方 式进行的断开一回线路。运用 matlab 编织的潮流程序（如附录 1 所示），再采用如下支路数据，得出 最终的结果。1-2：B1=12(0.054+0.308i)*27 /4840.00000364*27*4841 023

16、(0.054+0.308i)*17 /484/20.00000364*17*484*210;24(0.054+0.308i)*6 /484/20.00000364*6*484*210;36(0.054+0.308i)*18 /484/20.00000364*18*484*210;45(0.054+0.308i)*5 /484/20.00000364*5*484*21056(0.054+0.308i)*24 /484/20.00000364*24*484*210; 2-3B1=12(0.054+0.308i)*27 /484/20.00000364*27*484*210;23(0.054+0.3

17、08i)*17/4840.00000364*17*4841 0;24(0.054+0.308i)*6 /484/20.00000364*6*484*210;36(0.054+0.308i)*18 /484/20.00000364*18*484*210;45(0.054+0.308i)*5 /484/20.00000364*5*484*210;56(0.054+0.308i)*24 /484/20.00000364*24*484*210; 2-4B1=12(0.054+0.308i)*27 /484/20.00000364*27*484*210;23(0.054+0.308i)*17/484/2

18、0.00000364*17*484*210;24(0.054+0.308i)*6 /484 0.00000364*6*4841 0;36(0.054+0.308i)*18 /484/20.00000364*18*484*210;45(0.054+0.308i)*5 /484/20.00000364*5*484*210;56(0.054+0.308i)*24 /484/20.00000364*24*484*210; 3-6B1=12(0.054+0.308i)*27 /484/20.00000364*27*484*210;23(0.054+0.308i)*17/484/20.00000364*1

19、7*484*210;24(0.054+0.308i)*6 /484/20.00000364*6*484*210;36(0.054+0.308i)*18 /4840.00000364*18*4841 0;45(0.054+0.308i)*5 /484/20.00000364*5*484*210;56(0.054+0.308i)*24 /484/20.00000364*24*484*210; 4-5B1=12(0.054+0.308i)*27 /484/20.00000364*27*484*210;23(0.054+0.308i)*17/484/20.00000364*17*484*210;24(

20、0.054+0.308i)*6 /484/20.00000364*6*484*210;36(0.054+0.308i)*18 /484/20.00000364*18*484*210;45(0.054+0.308i)*5 /484 0.00000364*5*48410;56(0.054+0.308i)*24 /484/20.00000364*24*484*210; 5-6B1=12(0.054+0.308i)*27 /484/20.00000364*27*484*210;23(0.054+0.308i)*17/484/20.00000364*17*484*210;24(0.054+0.308i)

21、*6 /484/20.00000364*6*484*210;36(0.054+0.308i)*18 /484/20.00000364*18*484*210;45(0.054+0.308i)*5 /484/20.00000364*5*484*210;56(0.054+0.308i)*24 /4840.00000364*24*48410; B2=0 01.04551.04550 1;0 (230+20i)/100100 2;0 (350+90i)/100100 2;0 (250+85i)/100100 2;0 (410+90i)/100100 2;( 1+0.6197i) *711.45/1000

22、1 1.050 3;运行时输入参数： chaoliujisuan1是否在进行 n-1 校验，是为 1，否为 0: 1 请输入 n-1 校验断线序号，从 1-6: 1 chaoliujisuan1是否在进行 n-1 校验，是为 1，否为 0: 1 请输入 n-1 校验断线序号，从 1-6: 2 chaoliujisuan1是否在进行 n-1 校验，是为 1，否为 0: 1 请输入 n-1 校验断线序号，从 1-6: 3 chaoliujisuan1是否在进行 n-1 校验，是为 1，否为 0: 1 请输入 n-1 校验断线序号，从 1-6: 4 chaoliujisuan1是否在进行 n-1 校

23、验，是为 1，否为 0: 1 请输入 n-1 校验断线序号，从 1-6: 5 chaoliujisuan1是否在进行 n-1 校验，是为 1，否为 0: 1 请输入 n-1 校验断线序号，从 1-6: 6最终结果如下：n-1 1-2 各节点的电压大小 V 为 （节点号从小到大排列 ）：1.0455 1.0332 1.0373 1.0320 1.0330 1.0500 各节点的电压相角 sida 为 （节点号从小到大排列 ）：0 -4.9250 -4.9544 -5.2359 -5.2940 -3.9780n-12-3各节点的电压大小 V为（节点号从小到大排列）：1.04551.03531.03

24、931.03371.03441.0500各节点的电压相角sida为（节点号从小到大排列）：0-2.5405-2.6057-2.8603-2.9248-1.6195n-12-4各节点的电压大小 V为（节点号从小到大排列）：1.04551.03631.03881.03331.03411.0500各节点的电压相角sida为（节点号从小到大排列）：0-2.5592-2.6782-3.1444-3.1869-1.7769n-13-6各节点的电压大小 V为（节点号从小到大排列）：1.04551.03331.03311.03201.03301.0500各节点的电压相角sida为（节点号从小到大排列）：0-2

25、.5290-2.7567-2.7757-2.7785-1.1796n-14-5各节点的电压大小 V为（节点号从小到大排列）：1.04551.03581.03851.03401.03511.0500各节点的电压相角sida为（节点号从小到大排列）：0-2.5545-2.6142-2.8661-2.9827-1.6508n-15-6各节点的电压大小 V为（节点号从小到大排列）：1.04551.03271.03691.02971.02931.0500各节点的电压相角sida为（节点号从小到大排列）：0-2.5202-2.3105-2.9206-3.0514-1.06463.2使用powerworld

26、进行仿真断开线路：12表10结点电压节点123456电压幅值1.04551.03391.0381.03261.03361.0500电压相角/ 0-4.980-5.016-5.293-5.353-4.044由powerworld的实例信息可导出下面的最大相电流，如下表所示:表11 一回线的最大相电流线路122324364556电流/kA1.37040.0980.39150.48780.11800.4917250 MW85 Mv打 f41.0455 pu1B1.0339 pu1*038 puL0500 pu6230 MW20 Mvar50 IMWMvar546 MW-4 Mvar200 MW249

27、 MW 108 Mvar108 Mvar410 MW A90 NTvar*pII51.0327 pu1.0336 pu249 MW108 Mvar图9断开线路12电流电压分布断开线路：23表12结点电压节点123456电压幅值1.04551.03601.04001.03441.03511.0500电压相角/ 0-2.435-2.469-2.715-2.744-1.479表13 一回线的最大相电流线路122324364556电流/kA0.68080.09500.39240.46950.099350.4811250 MW 快、85 Mvar 11 4 一L410 MW/N90 Mvar541 MW

28、27 Ntvar1.0455 pu1,0344 puL0351 pu230 MW20 MvarL0360 pu3ha. _ 1.040 pu350 MW90 Mvar1.0500 pu249 MW89 Mvar208 MW249 MW 的 Mvar89 Mvar图10断开线路23电流电压分布断开线路：24表14结点电压节点123456电压幅值1.04551.03711.03901.03401.03471.0500电压相角/ 0-2.445-2.538-3.012-3.048-1.628表15 一回线的最大相电流线路122324364556电流/kA0.68040.07290.72370.446

29、20.07940.51011.0455 pu230 MW20 Mvar410 go21.0340 pu31.0500 pu6350 MW90 Mvar250 MW85 Mvar 斗-A1.0347 pu249 MW93 Mvar249 MW93 Mvar1.0371 pu1.039 pu208 MW93 Mvar图11断开线路24电流电压分布 断开线路：36表16结点电压节点123456电压幅值1.04551.03411.03401.03291.03371.0500电压相角/ 0-2.426-2.638-2.664-2.661-1.053表17 一回线的最大相电流线路122324364556电

30、流/kA0.68460.091830.30280.74660.073850.5680542 MW49 Mvar83 MvarB3 Mvar图12断开线路36电流电压分布断开线路：45表18结点电压节点123456电压幅值1.04551.03651.03901.03471.03571.0500电压相角/ 0-2.439-2.474-2.742-2.846-1.502表19 一回线的相电流线路122324364556电流/kA0.68040.08330.38900.47100.17380.48111.0455 pu21.0365 pu350 MW1,0347 pu0 Mvar1.039 pu362

31、08 MW90 Mvar410 MW90 Mvan541 MW21 Mvar230 MW20 Mvar250 MIV 785 Mvar 11 41.0357 puo o co249 MW90 Mvar249 MW90 Mvar1.0500 pu图13断开线路45电流电压分布断开线路：56表20结点电压节点123456电压幅值1.04551.03351.03701.03061.03011.0500电压相角/ 0-2.423-2.193-2.817-2.943-0.944表21 一回线的相电流线路122324364556电流/kA0.68610.14030.51910.58750.18880.70

32、22410 MW A90 Mvar1.0455 pu11.0335 pu1.0500 pu6350 MW90 Mvar250 MW85 Mvar 斗543 MW57 Mvar230 MW20 NTvar208 MW249 MW 82 Mvar82 Mvarli1.0306 puL037 puK0301 pu249 MW&2 Mvar图14断开线路56电流电压分布经过与powerworld仿真对比，matlab的所用潮流程序所得结果与powerworld仿真基本一致，故可说明 matlab的潮流程序基本正确第四章短路计算4.1 powerworld的短路容量计算本文电力系统暂态稳定是在发电机满发

33、方式，负荷重载方式进行的 故障：三相短路Fault NameASkdSolvedFault Object file Fermat)FaUtLotafionType for FadtlType for FaultsFaultResistaraFault ReactanceFault 1 CurrenlMagFault 1. CirrenlAng百l1Sus1!NO桂5BusI1SLG0.0M0,00022.7060266,441352Bus?NO性5Bus?SLGO.WOQ.OGO25.219105,山6船SBusJfNOYESBus?51G0,0000,00025,216035.7032斗Bu

34、s4NOYESBus4SLG0.0000,00025,0994565画5i,IfESBus弓SIGQ.OOO0.00025.10384畐 289860U56MOYESBus 5SIG0.000O.OOO26,85433表三相短路短路电流及短路容量220kV母线#1#2#3#4#5#6短路电流（标幺值）22.706025.219125.2168:25.099425.103826.85433短路电流/kA5.69916.32776.32946.29996.30226.7404短路容量/ X04MVA2.27062.52192.52162.50992.51032.6854334.2 matlab的短

35、路容量计算Matlab所编制的参数如附录二所示。输入相应的参数矩阵，请输入各支路阻抗矩阵0.1672 + 0.0534i0.1657 + 0.0456i0.1651 + 0.0433i0.1656 + 0.0453i0.1654 + 0.0451i0.1646 + 0.0440i0.1657 + 0.0456i0.1657 + 0.0456i0.1650 + 0.0433i0.1655 + 0.0453i0.1653 + 0.0451i0.1646+ 0.0440i0.1651 + 0.0433i0.1650 + 0.0433i0.1650 + 0.0448i0.1649 + 0.0435i0

36、.1648 + 0.0436i0.1643+ 0.0442i0.1656 + 0.0453i0.1655 + 0.0453i0.1649 + 0.0435i0.1656 + 0.0466i0.1654 + 0.0463i0.1646+ 0.0446i0.1654 + 0.0451i0.1653 + 0.0451i0.1648 + 0.0436i0.1654 + 0.0463i0.1655 + 0.0473i0.1645+ 0.0451i0.1646 + 0.0440i0.1646 + 0.0440i0.1643 + 0.0442i0.1646 + 0.0446i0.1645 + 0.0451i

37、0.1644+ 0.0475i请输入各支路电纳矩阵00.10404879800.10404879800.0655122060.0231219550.0655122060.0693658660.0231219550.01926829600.0192682960.0924878210.06936586600.092487821请输入结点负荷0 2.3-02 3.50+09 2.50+0.85i 4.10+09 0请输入各节点电压1.00001.0111 - 0.0368i1.0399 - 0.0475i1.0144 - 0.0417i1.0191 -0.0424i1.0496 - 0.0285i可

38、得到各节点短路时的短路电流幅值22.938725.093926.902625.458126.849624.3111以上各个幅值都是标幺值，与上面的powerworld仿真基本一致，故而 matlab程序基本上是正确的。所以短路容量等如下表：220kV母线#1#2#3#4#5#6短路电流（标幺值）22.706025.219125.2168:25.099425.103826.85433短路电流/kA5.69916.32776.32946.29996.30226.7404短路容量/ X04MVA2.27062.52192.52162.50992.51032.685433第五章 课程设计总结与心得在本

39、次课程设计的过程中， 本人通过区域电力网络电力系统运行方式分析和 计算，巩固和加深对课程所学理论知识的理解， 培养独立分析和解决问题的能力， 掌握电力工程运行计算设计的基本方法和基本步骤。在潮流计算过程中，由于时间限制与水平差距， 许多地方尚且存在不足之处。 例如由于对电力系统分析的理解不是很透彻， 很多参数计算不是十分精确， 许多 矩阵的形成需要大量查阅资料。另外，由于 matlab 编程也不太熟悉，故而导致 我们只能参考前人的程序， 理解不是十分透彻， 导致程序运行中会偶尔出现问题， 在找出问题所在的时候花费了大量的时间。另外，在使用 powerworld 的过程中，由于刚开始接触，很多地

40、方不懂得如 何操作，许多参数也不知道如何输入， 以至于学习软件使用方法也花了很多的时 间，最 终通 过询问同学得到解 决。 我还发现 powerworld13 的模型 到了 powerworld17 ，发电机的模型会丢失， 导致刚开始暂态计算也无法进行， 幸好后 来发现问题并及时解决了。本小组的课程设计，在我们两个同学的不懈努力下终于完成， 由于水平有限， 依旧存在很多缺陷，希望老师批评指正end33第六章 附录附录一 潮流计算程序%本程序的功能是用牛顿拉夫逊法进行潮流计算n=6 ;nl=6 ; isb=1 ; pr=0.0001 B1=1J2(0.054+0.308i)*27/484/20.

41、00000364*27*484*223(0.054+0.308i)*17/484/20.00000364*17*484*224(0.054+0.308i)*6/484/20.00000364*6*484*236(0.054+0.308i)*18/484/20.00000364*18*484*245(0.054+0.308i)*5/484/20.00000364*5*484*256(0.054+0.308i)*24/484/20.00000364*24*484*2clear;1111110;0;0;0;0;0; ;B2=001.04551.050 1;0(230+20i)/10011.0500(

42、350+90i)/10011.050 2;0(250+85i)/10011.050 2;0(410+90i)/10011.050 2(1 +0.6 1 97i)*209.2 5 /1 0001.051.050 3;% B1 矩阵：1 、支路首端号；2 、末端号； 3 、支路阻抗；4、支路对地电纳%5、支路的变比；6、支路首端处于K 侧为 1 ，1 侧为 0% B2 矩阵：1 、该节点发电机功率； 2、该节点负荷功率；3、节点电压初始值%4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量%6、节点分类标号：1为平衡节点（应为1号节点）;2为PQ节点;%3 为 PV 节点；isn_1=i

43、nput（ 是否在进行 n-1 校验，是为 1 ，否为 0:）;if isn_1=0isge=input（ 请输入发电机出力状态，满载为 1 ，轻载为 0:）;if isge=0 B2(6,1)=(1+0.6197i)*209.25/100;else if isge=1 B2(6,1)=(1+0.6197i)*711.45/100;end end elseif isn_1=1, B2(6,1)=(1+0.6197i)*711.45/100; nn_1=input( 请输入 n-1 校验断线序号，从 1-6:);B1(nn_1,3)=B1(nn_1,3)*2; B1(nn_1,4)=B1(nn_

44、1,4)/2;endY=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);% % %支路数%左节点处于 1 侧%左节点处于 K 侧for i=1:nl if B1(i,6)=0 p=B1(i,1);q=B1(i,2);elsep=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)A2)+B1(i,4)./ 2;%对角元

45、 K侧Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2;%对角元 1 侧end %求导纳矩阵disp(导纳矩阵Y=);disp(Y)%G=real(Y);B=imag(Y);%分解出导纳阵的实部和虚部for i=1:n%给定各节点初始电压的实部和虚部e(i)=real(B2(i,3);f(i)=imag(B2(i,3);V(i)=B2(i,4);%PV节点电压给定模值endfor i=1:nS(i)=B2(i,1)-B2(i,2);B(i,i)=B(i,i)+B2(i,5);end%给定各节点注入功率%i 节点注入功率 SG-SL%i 节点无功补偿量%=P=real(S)

46、;Q=imag(S);%分解出各节点注入的有功和无功功率ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; %迭代次数ICT1、a；不满足收敛要求的节点数 IT2 while IT2=0% N0=2*n雅可比矩阵的阶数；N=N0+1扩展列IT2=0;a=a+1;for i=1:nif i=isb%非平衡节点C(i)=0;D(i)=0;for j1=1:nC(i)=C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1);% 工(Gij*ej-Bij*fj)D(i)=D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1);% 工(Gij*fj+Bij*ej)endP

47、仁C(i)*e(i)+f(i)*D(i);% 节点功率 P 计算 eiS (Gij*ej-Bij*fj)+fi 工(Gij*fj+Bij*ej) Q仁C(i)*f(i)-e(i)*D(i);% 节点功率 Q 计算 fiS (Gij*ej-Bij*fj)-ei S (Gij*fj+Bij*ej)%求i节点有功和无功功率P,Q的计算值V2=e(i)A2+f(y2;% 电压模平方%=以下针对非 PV节点来求取功率差及 Jacobi矩阵元素 =if B2(i,6)=3%非 PV节点DP=P(i)-P1;%节点有功功率差DQ=Q(i)-Q1;%节点无功功率差%= 以上为除平衡节点外其它节点的功率计算%=求取 Jacobi 矩阵 =for j1=1:nif j1=isb&j1=i%非平衡节点 &非对角元X1=-G(i,j1)*e(i)-B(i,j1)*f(i);% dP/de=-dQ/dfX2=B(i,j1)*e(i)-G(i,j1)*f(i