一、实验目的应用 MATLAB 语言编写具有一定通用性的牛顿-拉夫逊法潮流计算程序。要求：（1）潮流计算方法为牛顿-拉夫逊法。（2）编程语言为 MATLAB。（3）程序具有较强通用性。二、程序流程图所用公式 其中 开始形成节点导纳矩阵输入原始数据设节点电压， i=1,2…,n,is置迭代次数置节点号 i=1按式（ 3-3 ），（ 3-4 ）计算雅克比矩阵元素按式（ 3-2 ）计算节点的，，节点的，求解修正方程式，得，雅 克 比 矩 阵 是 否 已 全 部 形成？计算平衡节点及 PV 节点功率求，迭代次数 k=k+1i=i+1？潮流计算完成计算各节点电压的新值：三、求解问题及其结果 IEEE-美国新英格兰 10 机 39 节点测试系统一、系统单线图GGGGGGGGGG303912253729172693381654182728362435222120342319331011131415831126327 二、系统参数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.0224.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.03000平衡节点BUS-310.0(电压幅角)1.0（幅值）4）线路参数见表 D-3表 D-3 线路参数序号节点 I节点 JR(p.u.)X(p.u.)B/2(p.u.)LN1LN2LN3LN4LN5BUS-2BUS-39BUS-3BUS-25BUS-4BUS-1BUS-1BUS-2BUS-2BUS-30.003500.001000.001300.007000.001300.041100.025000.015100.008600.021300.349350.37500.128600.073000.11070LN6LN7LN8LN9LN10LN11LN12LN13LN14LN15LN16LN17LN18LN19LN20LN21LN22LN23LN24LN25LN26LN27LN28LN29LN30LN31LN32LN33LN34BUS-18BUS-5BUS-14BUS-6BUS-8BUS-7BUS-11BUS-8BUS-9BUS-39BUS-11BUS-13BUS-14BUS-15BUS-16BUS-17BUS-19BUS-21BUS-24BUS-18BUS-27BUS-22BUS-23BUS-24BUS-26BUS-27BUS-28BUS-29BUS-29BUS-3BUS-4BUS-4BUS-5BUS-5BUS-6BUS-6BUS-7B...
