基于牛顿拉夫逊法潮流计算的matlab实验报告(含源程序和结果).doc
《基于牛顿拉夫逊法潮流计算的matlab实验报告(含源程序和结果).doc》由会员分享，可在线阅读，更多相关《基于牛顿拉夫逊法潮流计算的matlab实验报告(含源程序和结果).doc（24页珍藏版）》请在冰豆网上搜索。
基于牛顿拉夫逊法潮流计算的matlab实验报告
一、实验目的和要求
1.学习掌握matlab的基本用法
2.应用MATLAB语言编写具有一定通用性的牛顿-拉夫逊法潮流计算程序。
要求：
（1）潮流计算方法为牛顿-拉夫逊法。
（2）编程语言为MATLAB。
（3）程序具有较强通用性。
二、程序流程图
1.程序流程图
开始
形成节点导纳矩阵
输入原始数据，节点重新编号
设节点电压初值，i=1,2…,n,is
置迭代次数P=0
置节点号i=1
计算雅克比矩阵元素
按公式计算节点的，，节点的，
求解修正方程式，得，
雅克比矩阵是否已全部形成？
节点编号逆回，算平衡节点及PV节点功率
求，
迭代次数P=P+1
i=i+1
？
潮流计算完成
计算各节点电压的新值：
三、求解问题及其结果
1.求解问题：
IEEE-美国新英格兰10机39节点测试系统
1）系统单线图
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
PV
PV
PQ
PV
PV
PV
BUS-30
BUS-32
BUS-33
BUS-34
BUS-35
BUS-36
BUS-37
BUS-38
BUS-39
250.00
650.00
632.00
508.00
650.00
560.00
540.00
830.00
1000.00
175.90
103.35
96.88
1.04750
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
LN33
LN34
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-29
BUS-29
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
BUS-26
BUS-28
0.00350
0.00100
0.00130
0.00700
0.00130
0.00110
0.00080
0.00080
0.00020
0.00080
0.00060
0.00070
0.00040
0.00230
0.00100
0.00040
0.00040
0.00090
0.00180
0.00090
0.00070
0.00160
0.00080
0.00030
0.00070
0.00130
0.00080
0.00060
0.00220
0.00320
0.00140
0.00430
0.00570
0.00140
0.04110
0.02500
0.01510
0.00860
0.02130
0.01330
0.01280
0.01290
0.00260
0.01120
0.00920
0.00820
0.00460
0.03630
0.02500
0.00430
0.00430
0.01010
0.02170
0.00940
0.00890
0.01950
0.01350
0.00590
0.00820
0.01730
0.01400
0.00960
0.03500
0.03230
0.01470
0.04740
0.06250
0.01510
0.34935
0.3750
0.12860
0.07300
0.11070
0.10690
0.06710
0.06910
0.02170
0.07380
0.05650
0.06945
0.03900
0.19020
0.60000
0.03645
0.03645
0.08615
0.18300
0.08550
0.06710
0.15200
0.12740
0.03400
0.06595
0.16080
0.12825
0.09230
0.18050
0.25650
0.11980
0.39010
0.51450
0.12450
LN35:
BUS-4接有并联电容器，B4＝1.0000
LN36:
BUS-5接有并联电容器，B4＝2.0000
5）变压器参数见表D-4
表D-4变压器参数
序号
节点I
节点J
R（p.u.）
X（p.u.）
变比（%）
TR37
TR38
TR39
TR40
TR41
TR42
TR43
TR44
TR45
TR46
TR47
TR48
BUS-11
BUS-13
BUS-30
BUS-31
BUS-32
BUS-34
BUS-33
BUS-35
BUS-36
BUS-37
BUS-38
BUS-20
BUS-12
BUS-12
BUS-2
BUS-6
BUS-10
BUS-20
BUS-19
BUS-22
BUS-23
BUS-25
BUS-29
BUS-19
0.00160
0.00160
0.00000
0.00000
0.00000
0.00090
0.00070
0.00000
0.00050
0.00060
0.00080
0.00070
0.04350
0.04350
0.01810
0.02500
0.02000
0.01800
0.01420
0.01430
0.02720
0.02320
0.01560
0.01380
100.6
100.6
102.5
107.0
107.0
100.9
107.0
102.5
100.0
102.5
102.5
106.0
%IEEE-美国新英格兰10机39节点测试系统
%123456
%busvoltanglepqtype
bus=[11.00000.000.000.001
21.00000.000.000.001
31.00000.00-3.22-0.0241
41.00000.00-5.00-1.841
51.00000.000.000.001
61.00000.000.000.001
71.00000.00-2.338-0.841
81.00000.00-5.22-1.761
91.00000.000.000.001
101.00000.000.000.001
111.00000.000.000.001
121.00000.00-0.085-0.881
131.00000.000.000.001
141.00000.000.000.001
151.00000.00-3.20-1.531
161.00000.00-3.29-0.3231
171.00000.000.000.001
181.00000.00-1.58-0.301
191.00000.000.000.001
201.00000.00-6.80-1.031
211.00000.00-2.74-1.151
221.00000.000.000.001
231.00000.00-2.475-1.151
241.00000.00-3.08-0.9221
251.00000.00-2.24-0.4721
261.00
升级会员 