牛顿拉夫逊法潮流计算matlab程序

牛顿拉夫逊潮流计算程序

电力系统潮流计算的计算机算法(algorithm )提示 1、概 述

2、潮流计算(power-flow calculation )的基本方程3、牛顿-拉夫逊法(Newton-Laphson method )潮流计算 4、P Q分解法潮流计算

小

结

电力系统分析

牛顿拉夫逊潮流计算程序

本章提示 节点分类的概念； 潮流计算的基本方程式； 牛顿—拉夫逊法潮流计算的计算机算法；

P—Q分解法潮流计算的计算机算法。

电力系统分析

牛顿拉夫逊潮流计算程序

1概 述

类型：导纳法

阻抗法牛顿-拉夫逊法(N—R法) 快速分解法( PQ分解法)

电力系统分析

牛顿拉夫逊潮流计算程序

2 潮流计算的基本方程

2.1节点的分类 2.2基本方程式

电力系统分析

牛顿拉夫逊潮流计算程序

2.1节点(bus)的分类根据电力系统中各节点性质的不同，可把节点分成三种类型。 1.PQ节点 事先给定的是节点功率(P、Q),待求的是节点电压向量 (U、θ)。通常变电所母线都是PQ节点，当某些发电机的 出力P、Q给定时，也可作为PQ节点。PQ节点上的发电机称 之为PQ机(或PQ给定型发电机)。在潮流计算中，系统大 部分节点属于PQ节点。

电力系统分析

牛顿拉夫逊潮流计算程序

2.PU节

%本程序的功能是用牛顿——拉夫逊法进行潮流计算

% B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 % 5、支路的变比；6、支路首端处于K侧为1，1侧为0

% B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值 % 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 % 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点； % 3为PV节点； clear;

n=10;%input('请输入节点数：n='); nl=10;%input('请输入支路数：nl=');

isb=1;%input('请输入平衡母线节点号：isb='); pr=0.00001;%input('请输入误差精度：pr=');

B1=[1 2 0.03512+0.08306i 0.13455i 1 0; 2 3 0.0068+0.18375i 0 1.02381 1; 1 4 0.05620+0.13289i 0.05382i 1

%本程序的功能是用牛顿——拉夫逊法进行潮流计算

% B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 % 5、支路的变比；6、支路首端处于K侧为1，1侧为0

% B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值 % 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 % 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点； % 3为PV节点； clear;

n=10;%input('请输入节点数：n='); nl=10;%input('请输入支路数：nl=');

isb=1;%input('请输入平衡母线节点号：isb='); pr=0.00001;%input('请输入误差精度：pr=');

B1=[1 2 0.03512+0.08306i 0.13455i 1 0; 2 3 0.0068+0.18375i 0 1.02381 1; 1 4 0.05620+0.13289i 0.05382i 1

安徽工程大学 本科生课程设计说明书

1

安徽工程大学课程设计说明书

目录

安徽工程大学课程设计任务书······················3 摘要··············································5 Abstract··········································5 第一章电力系统潮流计算概述························6 1.1电力系统概述···································6 1.2 电力系统潮流概述·····························7 1.3 潮流计算的目的································8 1.4电力系统的发展和分析计算·····················9 1.5、MATLAB软件的应用·····························10 第二章 牛顿—拉夫逊法潮流计算基本原理···········11 2.1牛顿—拉夫逊法潮流计算简介····················11 2.2牛顿—拉夫逊法潮流计算计算公式

安徽工程大学 本科生课程设计说明书

1

安徽工程大学课程设计说明书

目录

安徽工程大学课程设计任务书······················3 摘要··············································5 Abstract··········································5 第一章电力系统潮流计算概述························6 1.1电力系统概述···································6 1.2 电力系统潮流概述·····························7 1.3 潮流计算的目的································8 1.4电力系统的发展和分析计算·····················9 1.5、MATLAB软件的应用·····························10 第二章 牛顿—拉夫逊法潮流计算基本原理···········11 2.1牛顿—拉夫逊法潮流计算简介····················11 2.2牛顿—拉夫逊法潮流计算计算公式

电工专业必看

%简单潮流计算的牛顿拉夫逊程序，相关的原始数据数据数据输入格式如下： %B1是支路参数矩阵，第一列和第二列是节点编号。节点编号由小到大编写 %对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点 %编号，将变压器的串联阻抗置于低压侧处理。

%第三列为支路的串列阻抗参数。

%第四列为支路的对地导纳参数。

%第五列为含变压器支路的变压器的变比

%第六列为变压器是否含有变压器的参数，其中“1”为含有变压器，

%“0”为不含有变压器。

%B2为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列为节点负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数。

%X为节点号和对地参数矩阵。其中第一列为节点编号，第二列为节点对地 %参数。

n=input('请输入节点数:n=');

n1=input('请输入支路数:n1=');

isb=input('请输入平衡节点号:isb=');

pr=input('请输入误差精度:pr=');

B1=input('请输入支路参数:B1=');

B2=input('请输入节点参数:B2=');

X=input('节点号和对地参数:X=');

Y=zeros(n

电力系统潮流计算的MATLAB辅助程序设计

潮流计算，通常指负荷潮流，是电力系统分析和设计的主要组成部分，对系统规划、安全运行、经济调度和电力公司的功率交换非常重要。此外，潮流计算还是其它电力系统分析的基础，比如暂态稳定，突发事件处理等。现代电力系统潮流计算的方法主要:高斯法、牛顿法、快速解耦法和MATLAB的M语言编写的MATPOWER4.1，这里主要介绍高斯法、牛顿法和快速解耦法。高斯法的程序是lfgauss，其与lfybus、busout和lineflow程序联合使用求解潮流功率。lfybus、busout和lineflow程序也可与牛顿法的lfnewton程序和快速解耦法的decouple程序联合使用。(读者可以到MATPOWER主页下载MATPOWER4.1,然后将其解压到MATLAB目录下，即可使用该软件进行潮流计算)

一、高斯-赛德尔法潮流计算使用的程序：

高斯-赛德法的具体使用方法读者可参考后面的实例，这里仅介绍各程序的编写格式:

lfgauss：该程序是用高斯法对实际电力系统进行潮流计算，需要用到busdata和linedata两个文件。程序设计为输入负荷和发电机的有功MW和无功Mvar，以及节点电压标幺值和相角的角度值

电力系统计算潮流 14节点经典例子

ZL=20;

WC=0.01;

ZP=[1,2,0.01335,0.04211,0.0,0.0;

1,3,0.0,0.20912,0.0,0.0;

1,4,0.0,0.55618,0.0,0.0;

1,10,0.05811,0.17632,0.0,0.0340;

1,11,0.06701,0.17103,0.0,0.0128;

2,10,0.05695,0.17388,0.0,0.0346;

2,12,0.0,0.25202,0.0,0.0;

2,14,0.05403,0.22304,0.0,0.0492;

3,4,0.0,0.11001,0.0,0.0;

3,13,0.0,0.17615,0.0,0.0;

4,5,0.03181,0.08450,0.0,0.0;

4,9,0.12711,0.27038,0.0,0.0;

5,6,0.08205,0.19207,0.0,0.0;

6,12,0.09498,0.19890,0.0,0.0;

7,8,0.22092,0.19988,0.0,0.0;

7,12,0.12291,0.25581,0.0,0.0;

8,9,0.17093,0.34802,0.0,0.0;

8,12,0.06615,0.130

南华大学电气工程学院毕业设计

前言

电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性。可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和很基础的计算。

随着科学技术的发展，电力系统变得越来越复杂，电气工程师掌握一种好的能对电力系统进行仿真的软件是学习和研究的需要。文章简要介绍了MATLAB发展历史、组成和强大的功能，并用简单例子分别就编程和仿真两方面分析了MATIAB软件在电力系统研究中的具体应用。采取等效电路法，能对特殊、复杂地电力系统进行高效仿真研究，因此，掌握编程和仿真是学好MATLAB的基础。与众多专门的电力系统仿真软件相比，MATLAB软件具有易学、功能强大和开放性好，是电力系统仿真研究的有力工具。

第 1 页 共 49 页

南华大学电气工程学院毕业设计

1电力系统的基本概念

电力系统：发电机把机械能转化为电能，电能经变压器和电

%线路数据 Line=[1-节点编号 2-线路首端节点号 3-线路末端节点号 4-支路电阻 5-支路电抗 6-支路电纳(注意：此处取的是B/2)]

Line=[1 1 2 0.04 0.25 0.25 2 1 3 0.1 0.35 0 3 2 3 0.08 0.30 0.25 ];

%变压器数据 transform=[1-支路编号 2-支路首节点编号 3-支路末节点编号 4-支路电阻(p.u.) 5-支路电抗(p.u.) 6-变压器变比(p.u.)] transform=[ 1 2 4 0 0.015 1.05

2 3 5 0 0.03 1.05]; %数据预处理

Nbus=5; % 节点数 nline=size(Line,1); % 线路个数 ntrans=size(transform,1); % 变压器个数 slack=5; % 平衡节点号 Npq=4;