基于牛顿拉夫逊法潮流计算的matlab实验报告

一、实验目的

应用MATLAB语言编写具有一定通用性的牛 顿-拉夫逊法潮流计算程序。 要求：

（1）潮流计算方法为牛顿-拉夫逊法。 （2）编程语言为MATLAB。 （3）程序具有较强通用性。 二、程序流程图

所用公式

j?n???Pi?Pi??[ei(Gijej?Bijfj)?fi(Gijfj?Bijej)]j?1?j?n????Qi?Qi??[fi(Gijej?Bijfj)?ei(Gijfj?Bijej)]

j?1???U2?U2?(e2?f2)iii?i??

i?j

?Pi?H??ij?f??Bijei?Gijfii???PiN??ij?e?Gijei?Bijfii???QiJ???Gijei?Bijfi??Nij?ij?fi? ??Q?Lij?i??Bijei?Gijfi?Hij?ei??2?R??Ui?0?ij?fi??Ui2??Sij??e?0i?

i?j?Pi??Hii??f??Bijei?Gijfi?biii???PiN??ii?e?Gijei?Bijfi?aiii???QiJ???Gijei?Bijfi?aii?ii?fi???Lii??Qi??Bijei?Gijfi?bii

?ei??2?R??Ui?2fi?ii?fi??Ui2??Sii??e?2eii?

其中

j?n??aii?Giiei?Biifi??(Gijej?Bijfj)j?1?j?i??j?n?b?Gf?Be?(Gf?Be)?iiiiiiijjijj?iij?1?j?i?

开始 输入原始数据 形成节点导纳矩阵 设节点电压ei(0)fi(0)，i=1,2?,n,i?s 置迭代次数k?0 按式（3-2）计算PQ节点的?P(k)，?Q(k)，PV节点的?P(k)，?U(k)2 iiii 置节点号i=1 雅克比矩阵是否已全部形成？ 按式（3-3），（3-4）计算雅克比矩阵元素 迭代次数 k=k+1 i=i+1 计算各节点电压 的新值： 求解修正方程式，得?ei(k)，?fi(k) ei(k?1)?e(k)??e(k) fi(k?1)?f(k)??f(k) 求|?e(k)|max，|?f(k)|max |?e(k)|max,|?f(k)|max??？ 计算平衡节点及PV节点功率 潮流计算完成 三、求解问题及其结果 IEEE-美国新英格兰10机39节点测试系统

一、 系统单线图

G2130G3725262738318171615212829G39GG4578131210192034332235142436239611G

31G32GGG 二、系统参数

1）系统容量基准值为100MVA。 2) 负荷数据见表D-1

表D-1 负荷 数据

节点号 BUS- 3 BUS- 4 BUS- 7 BUS- 8 BUS- 12 BUS- 15 BUS- 16 BUS- 18 BUS- 20 有功（MW） 无功（Mvar） 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 有功（MW） 无功（Mvar） 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 发电机数据 发电机 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 BUS-31 有功（MW） 250.00 650.00 632.00 508.00 650.00 560.00 540.00 830.00 1000.00 无功（Mvar） 175.90 103.35 96.88 0.0(电压幅角) 电压（p.u.） 1.04750 1.01230 1.04930 1.02780 1.02650 1.03000 1.0（幅值）

4）线路参数见表D-3

表D-3 线路参数 序号 LN1 LN2 LN3 LN4 LN5 LN6 LN7 LN8 LN9 节点I BUS-2 BUS-39 BUS-3 BUS-25 BUS-4 BUS-18 BUS-5 BUS-14 BUS-6 节点J BUS-1 BUS-1 BUS-2 BUS-2 BUS-3 BUS-3 BUS-4 BUS-4 BUS-5 R(p.u.) 0.00350 0.00100 0.00130 0.00700 0.00130 0.00110 0.00080 0.00080 0.00020 X(p.u.) 0.04110 0.02500 0.01510 0.00860 0.02130 0.01330 0.01280 0.01290 0.00260 B/2(p.u.) 0.34935 0.3750 0.12860 0.07300 0.11070 0.10690 0.06710 0.06910 0.02170

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库基于牛顿拉夫逊法潮流计算的matlab实验报告在线全文阅读。