基于牛顿拉夫逊法潮流计算的matlab实验报告(含源程序和结果).docxVIP

形成节点导纳矩阵输入原始数据，节点重新编号ef设 形成节点导纳矩阵 输入原始数据，节点重新编号 e f 设节点电压初值 (0) (0) i i ，i=1,2…,n,i ?s 置迭代次数 P=0 置节点号 i=1 计算雅克比矩阵元素 ( k ) 按公式计算 PQ 节点的 ?Pi ( k ) ， ?Qi ， PV 节点的 ?Pi ( k ) ， ?U ( k )2 i 迭代次数 P=P+1 i=i+1 一、 基于牛顿拉夫逊法潮流计算的 matlab 实验报告 实验目的和要求 学习掌握 matlab 的基本用法 应用 MATLAB 语言编写具有一定通用性的牛顿-拉夫逊法潮流计算程序。 要求： 潮流计算方法为牛顿-拉夫逊法。 编程语言为 MATLAB。 程序具有较强通用性。 二、程序流程图 1.程序流程图 开始 雅克比矩阵是否已全部形成？ 3、求解问题及其结果 1. 求解问题：IEEE-美国新英格兰 10 机 39 节点测试系统 1）系统单线图 32541451239 3 25 4 14 5 12 39 G  30 G  37 26 28 29 2 27 1 G  18  17 38  G 16 21 15 G 24 36 9  6  13  19  23 7 11 20 22 8  31 G G  32 10  34 G G  33  G  35 2）系统参数 1）系统容量基准值为 100MVA。 2) 负荷数据见表 D-1 表 D-1 负荷 数据 节点号 有功 （MW） 无功 （Mvar） 节点号 有功 （MW） 无功 （Mvar） BUS- 3 BUS- 4 BUS- 7 BUS- 8 BUS- 12 BUS- 15 BUS- 16 BUS- 18 BUS- 20 322.0 500.0 233.8 522.0 8.5 320.0 329.0 158.0 680.0 2.4 184.0 84.0 176.0 88.0 153.0 32.3 30.0 103.0 BUS- 21 BUS- 23 BUS- 24 BUS- 25 BUS- 26 BUS- 27 BUS- 28 BUS- 29 BUS- 39 274.0 247.5 308.0 224.0 139.0 281.0 206.0 283.5 1104.0 115.0 84.6 -92.2 47.2 17.0 75.5 27.6 26.9 3）发电机数据见表 D-2 表 D-2 发电机数据 发电机 节点号 有功（MW） 无功（Mvar） 电压（p.u.） PV PQ PQ BUS-30 BUS-32 BUS-33 250.00 650.00 632.00 175.90 103.35 1.04750 PV PV PQ PV PV PV BUS-34 BUS-35 BUS-36 BUS-37 BUS-38 BUS-39 508.00 650.00 560.00 540.00 830.00 1000.00 96.88 1.01230 1.04930 1.02780 1.02650 1.03000 平衡节点 BUS-31 0.0(电压幅角) 1.0（幅值） 4）线路参数见表 D-3 表 D-3 线路参数 序号 节点 I 节点 J R(p.u.) X(p.u.) B/2(p.u.) LN1 LN2 LN3 LN4 LN5 LN6 LN7 LN8 LN9 LN10 LN11 LN12 LN13 LN14 LN15 LN16 LN17 LN18 LN19 LN20 LN21 LN22 LN23 LN24 LN25 LN26 LN27 LN28 LN29 LN30 LN31 LN32 BUS-2 BUS-39 BUS-3 BUS-25 BUS-4 BUS-18 BUS-5 BUS-14 BUS-6 BUS-8 BUS-7 BUS-11 BUS-8 BUS-9 BUS-39 BUS-11 BUS-13 BUS-14 BUS-15 BUS-16 BUS-17 BUS-19 BUS-21 BUS-24 BUS-18 BUS-27 BUS-22 BUS-23 BUS-24 BUS-26 BUS-27 BUS-28 BUS-1 BUS-1 BUS-2 BUS-2 BUS-3 BUS-3 BUS-4 BUS-4 BUS-5 BUS-5 BUS-6 BUS-6 BUS-7 BUS-8 BUS-9 BUS-10 BUS-10 BUS-13 BUS-14 BUS-15 BUS-16 BUS-16 BUS-16 BUS-16 BUS-17 BUS-17 BUS-21 BUS-22 BUS-23 BUS-25 BUS-26 BUS-26 0.00350