牛顿拉夫逊法潮流计算实验报告.docx

南昌大学实验报告

学生姓名：学号：专业班级：

实验类型：□验证□综合□设计□创新实验日期：实验成绩：

实验一牛顿-拉夫逊法潮流计算实验

一、实验目的

本实验针对IEEE37节点配电网结构特性，要求根据功率平衡方程的泰勒展开式近似求取系统状态变量，使系统潮流计算能够由计算机自行完成，即根据已知电力网络的数学模型（节点导纳矩阵）及各节点参数，由计算程序运行完成该电力系统的潮流计算。通过实验教学加深学生对复杂电力系统潮流计算方法的理解，学会运用电力系统的数学模型，掌握牛顿-拉夫逊法潮流计算的过程及其特点，熟悉应用软件MATLAB，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。

二、实验内容

编制调试牛顿-拉夫逊法潮流计算的计算机程序。程序要求根据已知的电网拓扑、线路参数及节点参数，完成电力系统的潮流计算，要求计算出节点电压、线路潮流等参数。

1.根据已知电网拓扑和系统参数生成系统导纳矩阵；

2.采用牛顿-拉夫逊法编制潮流程序；

3.对程序进行组织调试；

4.应用计算例题验证程序的计算效果；

5.对调试正确的计算程序进行存储、打印；

6.完成本次实验的实验报告。

三、程序清单

1.实验数据

节点

节点

长度(英尺)

欧姆/英里

欧姆/英里

有功负荷

无功负荷

额定电压4.8kV

1

2

1850

0.2926

0.1973

85kW

40kVar

选定S_b=100MVA

2

3

960

0.4751

0.2973

85kW

40kVar

Z_b=V_b*V_b/S_b

3

4

1320

0.4751

0.2973

85kW

40kVar

1英里=5280英尺

3

5

400

2.0952

0.7758

85kW

40kVar

节点1为平衡节点

3

6

360

1.2936

0.6713

85kW

40kVar

4

7

240

2.0952

0.7758

85kW

40kVar

4

8

600

1.2936

0.6713

85kW

40kVar

5

9

240

2.0952

0.7758

85kW

40kVar

5

10

320

2.0952

0.7758

85kW

40kVar

6

11

520

1.2936

0.6713

85kW

40kVar

7

12

280

1.2936

0.6713

85kW

40kVar

8

13

200

1.2936

0.6713

85kW

40kVar

11

14

80

2.0952

0.7758

85kW

40kVar

11

15

800

1.2936

0.6713

85kW

40kVar

12

16

200

2.0952

0.7758

85kW

40kVar

12

17

280

2.0952

0.7758

85kW

40kVar

13

18

600

1.2936

0.6713

85kW

40kVar

13

19

960

0.4751

0.2973

85kW

40kVar

13

20

320

1.2936

0.6713

85kW

40kVar

14

21

520

2.0952

0.7758

85kW

40kVar

15

22

600

1.2936

0.6713

85kW

40kVar

15

23

920

2.0952

0.7758

85kW

40kVar

20

24

320

2.0952

0.7758

85kW

40kVar

20

25

320

1.2936

0.6713

85kW

40kVar

22

26

280

2.0952

0.7758

85kW

40kVar

23

27

120

2.0952

0.7758

85kW

40kVar

23

28

760

2.0952

0.7758

85kW

40kVar

25

29

560

1.2936

0.6713

85kW

40kVar

29

30

520

2.0952

0.7758

85kW

40kVar

29

31

640

1.2936

0.6713

85kW

40kVar

30

32

200

2.0952

0.7758

85kW

40kVar

30

33

1280

2.0952

0.7758

85kW

40kVar

31

34

400

1.2936

0.6713

85kW

40kVar

34

35

400

1.2936

0.6713

85kW

40kVar

35

36

200

2.0952

0.7758

85kW

40kVar

35

37

400

1.2936

0.6713

85kW

40kVar

2.程序

clear;

clc;

n=37;%节