基于matlab的牛顿拉夫逊法潮流计算

姓名

学号

指导教师

专业与班级

所在学院

1

目录

一、 实验目的和要求 ....................................................................................................................................... 4 二、 实验原理与流程图.................................................................................................................................... 4

1.实验原理 ........................................................................................................................................................ 4

1.1 输入初始数据，节点编号，形成节点导纳矩阵 ............................................................................. 4

(0)(0)(0)(0)?QU??P1.2 对于初始值和，求出功率偏差和 ................................................................ 4

(0)J1.3 依据公式算出雅可比矩阵 ......................................................................................................... 4

1.4 应用高斯消元法解线性方程组，得到 、 ......................................................................... 5

1.5 依据如下公式计算修正量，并形成迭代公式 ................................................................................. 5

1.6 计算线路功率分布，PV节点无功功率和平衡节点注入功率 ....................................................... 6 2.程序流程图 .................................................................................................................................................... 6 三、 算例及结果 ............................................................................................................................................... 8

1.算例：IEEE-美国新英格兰10机39节点测试系统 ................................................................................... 8

1.1 系统单线图 ........................................................................................................................................ 8 1.2 系统参数 ............................................................................................................................................ 8 1.3发电机数据 ......................................................................................................................................... 9 1.4线路参数 ............................................................................................................................................. 9 2.计算结果 ...................................................................................................................................................... 13

2.1 节点计算结果................................................................................................................................... 13 2.2线路计算结果.................................................................................................................................... 14 3. 实验结果分析 ............................................................................................................................................ 15 附录： ..................................................................................................................................................................... 17

1.实验源程序 .................................................................................................................................................. 17

1.1主程序——nr.m ................................................................................................................................ 17 1.2 读取节点与线路数据——fileopen.m ............................................................................................. 20 1.3 对节点重新编号——num.m ........................................................................................................... 21 1.4 求节点导纳阵——y.m ..................................................................................................................... 22 1.5 计算功率偏差——deltaPQ.m .......................................................................................................... 23 1.6求雅克比矩阵——form_jac.m ......................................................................................................... 24 1.7 高斯消元法解线性方程组——gs.m ............................................................................................... 26 1.8 节点编号逆回——renum.m ............................................................................................................ 27 1.9 计算线路潮流及损耗——flow.m ................................................................................................... 28

作业一 ..................................................................................................................................................................... 30

1.写一个 MATLAB 函数pifun.m 来计算下列级数 ............................................................................. 30 2.一个平面上的椭圆可以表示成下列方程式： ..................................................................................... 31 2.1 程序——y.m ..................................................................................................................................... 31

作业二 ..................................................................................................................................................................... 32

1.试用牛顿－拉夫逊法求解如下的非线性方程组，要求精度达到10－10以上 ..................................... 32

1.1 计算结果 .......................................................................................................................................... 32 1.2 程序 .................................................................................................................................................. 34

2

2. 解方程组 .................................................................................................................................................... 35

2.1 计算结果 .......................................................................................................................................... 36 2.2 程序 .................................................................................................................................................. 37

3