《基于牛顿拉夫逊法潮流计算的matlab实验报告》由会员分享,可在线阅读,更多相关《基于牛顿拉夫逊法潮流计算的matlab实验报告(26页珍藏版)》请在金锄头文库上搜索。
1、一、实验目的应用MATLAB语言编写具有一定通用性的牛顿-拉夫逊法潮流计算程序。要求:(1)潮流计算方法为牛顿-拉夫逊法。(2)编程语言为MATLAB。(3)程序具有较强通用性。二、程序流程图所用公式 其中开始形成节点导纳矩阵输入原始数据设节点电压,i=1,2,n,is置迭代次数置节点号i=1按式(3-3),(3-4)计算雅克比矩阵元素按式(3-2)计算节点的,节点的,求解修正方程式,得,雅克比矩阵是否已全部形成?计算平衡节点及PV节点功率求,迭代次数 k=k+1i=i+1?潮流计算完成计算各节点电压的新值: 3、 求解问题及其结果 IEEE-美国新英格兰10机39节点测试系统一、 系统单线图
2、 二、系统参数1)系统容量基准值为100MVA。2) 负荷数据见表D-1表D-1 负荷 数据节点号有功(MW)无功(Mvar)节点号有功(MW)无功(Mvar)BUS- 3BUS- 4BUS- 7BUS- 8BUS- 12BUS- 15BUS- 16BUS- 18BUS- 20322.0500.0233.8522.08.5320.0329.0158.0680.02.4184.084.0176.088.0153.032.330.0103.0BUS- 21BUS- 23BUS- 24BUS- 25BUS- 26BUS- 27BUS- 28BUS- 29BUS- 39274.0247.5308.02
3、24.0139.0281.0206.0283.51104.0115.084.6-92.247.217.075.527.626.93)发电机数据见表D-2表D-2 发电机数据 发电机节点号有功(MW)无功(Mvar)电压(p.u.)PVPQPQPVPVPQPVPVPVBUS-30BUS-32BUS-33BUS-34BUS-35BUS-36BUS-37BUS-38BUS-39250.00650.00632.00508.00650.00560.00540.00830.001000.00175.90103.3596.881.047501.012301.049301.027801.026501.0300
4、0平衡节点BUS-310.0(电压幅角)1.0(幅值)4)线路参数见表D-3表D-3 线路参数序号节点I节点JR(p.u.)X(p.u.)B/2(p.u.)LN1LN2LN3LN4LN5LN6LN7LN8LN9LN10LN11LN12LN13LN14LN15LN16LN17LN18LN19LN20LN21LN22LN23LN24LN25LN26LN27LN28LN29LN30LN31LN32LN33LN34BUS-2BUS-39BUS-3BUS-25BUS-4BUS-18BUS-5BUS-14BUS-6BUS-8BUS-7BUS-11BUS-8BUS-9BUS-39BUS-11BUS-13BU
5、S-14BUS-15BUS-16BUS-17BUS-19BUS-21BUS-24BUS-18BUS-27BUS-22BUS-23BUS-24BUS-26BUS-27BUS-28BUS-29BUS-29BUS-1BUS-1BUS-2BUS-2BUS-3BUS-3BUS-4BUS-4BUS-5BUS-5BUS-6BUS-6BUS-7BUS-8BUS-9BUS-10BUS-10BUS-13BUS-14BUS-15BUS-16BUS-16BUS-16BUS-16BUS-17BUS-17BUS-21BUS-22BUS-23BUS-25BUS-26BUS-26BUS-26BUS-280.003500.00
6、1000.001300.007000.001300.001100.000800.000800.000200.000800.000600.000700.000400.002300.001000.000400.000400.000900.001800.000900.000700.001600.000800.000300.000700.001300.000800.000600.002200.003200.001400.004300.005700.001400.041100.025000.015100.008600.021300.013300.012800.012900.002600.011200.0
7、09200.008200.004600.036300.025000.004300.004300.010100.021700.009400.008900.019500.013500.005900.008200.017300.014000.009600.035000.032300.014700.047400.062500.015100.349350.37500.128600.073000.110700.106900.067100.069100.021700.073800.056500.069450.039000.190200.600000.036450.036450.086150.183000.0
8、85500.067100.152000.127400.034000.065950.160800.128250.092300.180500.256500.119800.390100.514500.12450LN35: BUS-4接有并联电容器,B41.0000LN36: BUS-5接有并联电容器,B42.00005)变压器参数见表D-4表D-4 变压器参数序号节点I节点JR(p.u.)X(p.u.)变比(%)TR37TR38TR39TR40TR41TR42TR43TR44TR45TR46TR47TR48BUS-11BUS-13BUS-30BUS-31BUS-32BUS-34BUS-33BUS-3
9、5BUS-36BUS-37BUS-38BUS-20BUS-12BUS-12BUS-2BUS-6BUS-10BUS-20BUS-19BUS-22BUS-23BUS-25BUS-29BUS-190.001600.001600.000000.000000.000000.000900.000700.000000.000500.000600.000800.000700.043500.043500.018100.025000.020000.018000.014200.014300.027200.023200.015600.01380100.6100.6102.5107.0107.0100.9107.010
10、2.5100.0102.5102.5106.0%IEEE-美国新英格兰10机39节点测试系统 % 1 2 3 4 5 6 % bus volt angle p q typebus= 1 1.0000 0.00 0.00 0.00 1 2 1.0000 0.00 0.00 0.00 1 3 1.0000 0.00 -3.22 -0.024 1 4 1.0000 0.00 -5.00 -1.84 1 5 1.0000 0.00 0.00 0.00 1 6 1.0000 0.00 0.00 0.00 1 7 1.0000 0.00 -2.338 -0.84 1 8 1.0000 0.00 -5.22
11、 -1.76 1 9 1.0000 0.00 0.00 0.00 1 10 1.0000 0.00 0.00 0.00 1 11 1.0000 0.00 0.00 0.00 1 12 1.0000 0.00 -0.085 -0.88 1 13 1.0000 0.00 0.00 0.00 1 14 1.0000 0.00 0.00 0.00 1 15 1.0000 0.00 -3.20 -1.53 1 16 1.0000 0.00 -3.29 -0.323 1 17 1.0000 0.00 0.00 0.00 1 18 1.0000 0.00 -1.58 -0.30 1 19 1.0000 0.00 0.00 0.00 1 20 1.0000 0.00 -6.80 -1.03 1 21 1.0000 0.00 -2.74 -1.15 1 22 1.0000 0.00 0.00 0.00 1 23 1.0000 0.00 -2.475 -1.15 1 24 1.0000 0.00 -3.08 -0.922 1 25 1.0000 0.00 -2.24 -0.472 1