山东交通学院

电力系统分析课程设计

报告书

院(部)别 班 级 学 号 姓 名 指导教师 时 间

课 程 设 计 任 务 书

题 目 电力系统分析课程设计 学 院 专 业

班 级 学生姓名 学 号

月 日至 月 日 共 周

指导教师(签字)

院 长(签字)

年 月 日

成绩评定表

摘 要

潮流计算，电力学名词，指在给定电力系统网络拓扑、元件参数和发电、负荷参量条件下，计算有功功率、无功功率及电压在电力网中的分布。潮流计算是根据给定的电网结构、参数和发电机、负荷等元件的运行条件，确定电力系统各部分稳态运行状态参数的计算。通常给定的运行条件有系统中各电源和负荷点的功率、枢纽点电压、平衡点的电压和相位角。待求的运行状态参量包括电网各母线节点的电压幅值和相角，以及各支路的功率分布、网络的功率损耗等。传统的潮流计算程序缺乏图形用户界面，结果显示不直观，难于与其他分析功能集成，网络原始数据输入工作量大且易于出错。结合电力系统的特点，对于复杂电力系统，根据定条件，应用牛顿-拉夫逊法进行计算，在手工计算中，由于涉及大量变量、微分方程、矩阵计算，求解很繁琐，计算不同系统时需要重新计算。运用 MATLAB软件进行仿真潮流计算，图形界面直观，运行稳定，计算准确，提高了计算速度，各个类的有效封装又使程序具有很好的模块性.可维护性和可重用性。

关键字： 潮流计算 牛拉法 Matlab

目录

摘 要 ............................................................................................................................................................ 6

1潮流计算...................................................................................................................................................... 1

1.1潮流计算概述 ..................................................................................................................... 1

1.2潮流计算的要求 ................................................................................................................. 2

1.3潮流计算的优势 ................................................................................................................. 2

1.4潮流计算的用途 ................................................................................................................. 3

2 MATLAB简介 ............................................................................................................................................ 4

2.1MATLAB概述 .................................................................................................................... 4

4.2牛顿拉夫逊法程序流程图 ............................................................................................... 13

4.3设计程序 ........................................................................................................................... 14

4.4运行结果 ........................................................................................................................... 16

心得体会....................................................................................................................................................... 22

参考文献....................................................................................................................................................... 23

1.1潮流计算概述 1潮流计算

潮流计算是研究电力系统稳态运行情况的一种基本电气计算，常规潮流计算的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。潮流计算的结果是电力系统稳定计算和故障分析的基础。通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。

具体表现在以下方面：

(1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。

(2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。

(3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。

(4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。

总结为在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。此外，电力系统潮流计算也是计算系统动态稳定和静态稳定的 - 1 -

基础。所以潮流计算是研究电力系统的一种很重要和基础的计算。

1.2潮流计算的要求

电力系统运行必须满足一定技术和经济上的要求。这些要求够成了潮流问题中某些变量的约束条件，常用的约束条件如下：

1.节点电压应满足

Uimin?Ui?Uimax(i?1,2,?n)

从保证电能质量和供电安全的要求来看，电力系统的所有电气设备都必须运行在额定电压附近。PU节点电压幅值必须按上述条件给定。因此，这一约束条件对PQ节点而言。

2.节点的有功功率和无功功率应满足

PGimin?PGi?PGimax?

Q?

Gimin?QGi?QGimax?

PQ节点的有功功率和无功功率，以及PU节点的有功功率，在给定是就必须满足上述条件，因此，对平衡节点的P和Q以及PU节点的Q应按上述条件进行检验。

3.节点之间电压的相位差应满足

|?ij|?|?i??j|?|?i??j|max

为了保证系统运行的稳定性，要求某些输电线路两端的电压相位不超过一定的数值。这一约束的主要意义就在于此。

因此，潮流计算可以归结为求解一组非线性方程组，并使其解答满足一定的约束条件。常用的方法是迭代法和牛顿法，在计算过程中，或得出结果之后用约束条件进行检验。如果不能满足要求，则应修改某些变量的给定值，甚至修改系统的运行方式，重新进行计算。

1.3潮流计算的优势

电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。 - 2 -

各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用于电力系统稳态研究，安全估计或最优潮流等，实际电力系统的潮流计算主要采用牛顿-拉夫逊法。借由MATLAB可以轻松实现计算复杂的电力系统潮流分布。

1.4潮流计算的用途 流计算是电力系统非常重要的分析计算，用以研究系统规划和运行中提出的各种问题。对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求；对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件(线路、变压器等)是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。

潮流计算是电力系统分析最基本的计算。除它自身的重要作用之外，在《电力系统分析综合程序》(PSASP)中，潮流计算还是网损计算、静态安全分析、暂态稳定计算、小干扰静态稳定计算、短路计算、静态和动态等值计算的基础

- 3 -

2 MATLAB简介

2.1MATLAB概述

MATLAB的含义是矩阵实验室（Matrix Laboratory），是美国MathWork公司于1982推出的一套高性能的数值计算可视化软件，，包括MATLAB主程序、SIMULINK 动态系统仿真包和各种专业工具箱它集数值分析，矩阵计算，信号处理和图形显示于一体，构成一个方便的，界面友好的用户环境，具有极强大的计算功能和极高的编程效率，特别适合于科学计算、数值分析、系统仿真和信号处理等任务。MATLAB程学设计语言结构完整，且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可以高效率的解决工业计算问题，特别是关于矩阵和矢量的计算。通过MATLAB语言，可以用类似数学公式的方式来编写算法，大大降低了程序需要的难度别难并节省了时间，从而可把主要的经历集中在算法的构思而不是编程上。学习运用MATLAB计算电力系统潮流分布是本次课程设计的重点，可以说，作为工科学生，会运用MATLAB来解决工程问题已成为必须。到目前为止，MATLAB已发展成为国际上最优秀的科技应用软件之一。它的功能十分强大，不仅仅可以实现计算潮流分布，还可以模拟仿真各式各样的数值系统，工程。这里将借助MATLAB来完成用直角牛顿-拉夫逊法计算电力系统潮流分布。

- 4 -

3牛顿-拉夫逊法概述

3.1牛顿-拉夫逊基本原理

潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采用牛顿-拉夫逊法。

牛顿--拉夫逊法(简称牛顿法)在数学上是求解非线性代数方程式的有效方法。其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程。即通常所称的逐次线性化过程。

对于非线性代数方程组：

f(x)?0 即 fi(x1,x2,?,xn)?0 (i?1,2,?,n) (3-1-1) 在待求量x的某一个初始估计值x(0)附近，将上式展开成泰勒级数并略去二阶及以上的高阶项，得到如下的经线性化的方程组：

f(x(0))?f'(x(0))?x(0)?0 (3-1-2) 上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量 ?x(0)??[f'(x(0))]?1f(x(0)) (3-1-3) 将?x(0)和x(0)相加，得到变量的第一次改进值x(1)。接着就从x(1)出发，重复上述计算过程。因此从一定的初值x(0)出发，应用牛顿法求解的迭代格式为：

f'(x(k))?x(k)??f(x(k)) (3-1-4) x(k?1)?x(k)??x(k) (3-1-5)

上两式中：f'(x)是函数f(x)对于变量x的一阶偏导数矩阵，即雅可比矩阵J;k为迭代次数。

- 5 -

有上式可见，牛顿法的核心便是反复形式并求解修正方程式。牛顿法当初始估计值x(0)和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。

牛顿潮流算法突出的优点是收敛速度快，若选择到一个较好的初值，算法将具有平方收敛特性，一般迭代4~5次便可以收敛到一个非常精确的解。而且其迭代次数与所计算网络的规模基本无关。牛顿法也具有良好的收敛可靠性，对于对以节点导纳矩阵为基础的高斯法呈病态的系统，牛顿法也能可靠收敛。牛顿法所需的内存量及每次迭代所需时间均较高斯法多。

牛顿法的可靠收敛取决于有一个良好的启动初值。如果初值选择不当，算法有可能根本不收敛或收敛到一个无法运行的节点上。对于正常运行的系统，各节点电压一般均在额定值附近，偏移不会太大，并且各节点间的相位角差也不大，所以对各节点可以采用统一的电压初值(也称为平直电压)，如假定：

Ui(0)?1 ?i(0)?0 或 ei(0)?1 fi(0)?0 (i?q1,2,?,n;i?s) (3-1-6) 这样一般能得到满意的结果。但若系统因无功紧张或其它原因导致电压质量很差或有重载线路而节点间角差很大时，仍用上述初始电压就有可能出现问题。解决这个问题的办法可以用高斯法迭代1~2次，以此迭代结果作为牛顿法的初值。也可以先用直流法潮流求解一次以求得一个较好的角度初值，然后转入牛顿法迭代。

3.2直角牛顿-拉夫逊法潮流计算求解过程

以下讨论的是用直角坐标形式的牛顿—拉夫逊法潮流的求解过程。当采用直角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量e,f,e,f...e,f1122nn由于平衡节点的电压向量是给定的，因此待求两共2(n?1)

需要2(n-1)个方程式。事实上，除了平衡节点的功率方程式在迭代过程中没有约束作用以外，其余每个节点都可以列出两个方程式。对PQ节点来说，Pis和Qis是给定的，因而可以写出

- 6 -

?Pi?

i

p?e?(Ge?Bf

is

ij?i

ij

j

ij

is

ij?i

ij

j

ij

j

)?

j

f?(Gf

jj?i

ij

j?i

j

?Q?Q?

f?(Ge?Bf)?ej?(Gij

f

?Bijej)?0?

?

（3-2-1） ?

?Bijej)?0?j

?

对PV节点来说，给定量是Pis和Vis，因此可以列出

?Pi?Pis?ei?(Gijej?Bij

j?i2

f

j

)?

f?(Gf

ij?i

ij

j

?Vi?Vis?(ei?

22

f

2i

)?0

?Bijej)?0?

?

? （3-2-2）

??

求解过程大致可以分为以下步骤：

（1）形成节点导纳矩阵Y （2）将各节点电压设初值U，

（3）将节点初值代入式(2-2-1)或式(2-2-2)，求出修正方程式的常数项向

量

（4）将节点电压初值代入求式，求出雅可比矩阵元素 （5）求解修正方程，求修正向量 （6）求取节点电压的新值

（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步

重新开始进行狭义次迭代，否则转入下一步

（8）计算支路功率分布，PV节点无功功率和平衡节点柱入功率。

以直角坐标系形式表示 ①. 迭代推算式

采用直角坐标时,节点电压相量及复数导纳可表示为:

??e?jfViiiYij?Gij?jBij

(3-2-3)

将以上二关系式代入上式中,展开并分开实部和虚部;假定系统中的第1,2,?,m号为P—Q节点,第m+1,m+2,?,n-1为P—V节点,根据节点性质的不同,得到如下迭代推算式:

⑴对于PQ节点

- 7 -

nn

?PG?

i?Pi?ei?(ijej?Bijfj)?fi?(Gijfj?Bijej)j?1j?1?

?

nn

? (3-2-4) ?Qi?Qi?fi?(Gijej?Bijfj)?ei?(Gijfj?Bijej)?

j?1j?1??

i?1,2,?,m

⑵对于PV节点

n

n

?P?

i?Pi?ei?(Gijej?Bijfj)?fij?1?(Gijfj?Bijej)j?1?

? (3-2-5)

?V2?f2?I?V2i?(e2ii)?

i?m?1,m?2,?,n?1

⑶对于平衡节点

平衡节点只设一个,电压为已知,不参见迭代,其电压为:

Vn?en?jfn (3-2-6)

②. 修正方程式迭代式共包括2(n-1)个方程.选定电压初值及变量修正量符号,代入方程并按泰勒级数展开,略去?ei,?fi二次方程及以后各项,得到一组线性方程组或线性化了的方程组，常称修正方程组:

?W??J?U (3-2-7)

???P1

???e1??

?Q?

?1????f?1?????

??P?

????m?

??em?W????Qm?

???f??P?

?U??m??m?1?

??e?m?1??U2m?1

?

?? ???

??fm?1???

?????P???e?n?1?

??n?1??U2

?n?1??

????fn?1??

- 8 -

???P1??e

1

?

???Q1??e

1????

???Pm??e1?

???Qm??e1J??

??Pm?1??e

1

?

???U2m?1??e

1????

???Pn?1??e1?2???Un?1??e

1?

??P1

?f1??Q1?f1???Pm?f1??Qm?f1?Pm?1?f1

??????

??P1?em??Q1?em???Pm?em??Qm?em?Pm?1?em

??P1?fm??Q1?fm???Pm?fm??Qm?fm?Pm?1?fm??U2m?1

?fm

???Pn?1?fm??U2n?1

?fm

??P1?em?1??Q1?em?1???Pm?em?1??Qm?em?1?Pm?1?em?1??U2m?1?em?1

???Pn?1?em?1??U2n?1?em?1

??P1?fm?1??Q1?fm?1???Pm?fm?1??Qm?fm?1?Pm?1?fm?1

??????

??P1?en?1??Q1?en?1???Pm?en?1??Qm?en?1?Pm?1?en?1

??U2m?1??U2m?1

?

?f1?em

???Pn?1?f1

??

???Pn?1?em

??U2m?1??U2m?1

?

?fm?1?en?1

???Pn?1?fm?1

??

???Pn?1?en?1

??U2n?1??U2n?1

?

?f1?em??U2n?1??U2n?1

?

?fm?1?en?1

??P?1

?fn?1?

???Q1??fn?1?

??????Pm??fn?1?

???Qm??fn?1?

??Pm?1??fn?1?

?2

??Um?1??fn?1??

?????Pn?1??fn?1?

?

??U2n?1??fn?1??

(3-2-8) ③.雅可比矩阵各元素的算式

式(3-2-8)中, 雅可比矩阵中的各元素可通过对式(3-2-4)和(3-2-5)进行偏导而求得.当j?i时, 雅可比矩阵中非对角元素为

???Pi??Qi

????(Gijei?Bijfi)???ej??fj?

???Pi??Qi?

??Bijei?Gijfi? (3-2-9) ??fj??ej

????U2??U2

???0

?ej?fj??

当j?i时,雅可比矩阵中对角元素为:

- 9 -

??Pn

i

?e???(G?B?

ijejijfj)?Giiei?Biifi?ij?1?

??Pn

i

?f???(G?

ijfj?Bijej)?Giifi?Biiei?

jj?1?

??Qn

i?

?e?

i?(Gijfj?Bijej)?Giifi?Biiei?j?1??

??Qn? (3-2-10)

i

?f???(Gij?ej?Bijfj)?Giiei?Biifi?

jj?1?

??U2?

i

?e??2e?

i

j??

??U2

i?

?f??2fi?

i??

由式(2-2-9)和(2-2-10)看出,雅可比矩阵的特点:

⒈矩阵中各元素是节点电压的函数,在迭代过程中,这些元素随着节点电压的变化而变化;

⒉导纳矩阵中的某些非对角元素为零时,雅可比矩阵中对应的元素也是为零.若Yij?0,则必有Jij?0;

⒊雅可比矩阵不是对称矩阵;(i?q1,2,?,n;i?s)

雅可比矩阵各元素的表示如下式(2-2-11):

?(Gijei

H??P??Bijfi)

i

?e??(j?i)

ij??

j??(Ge ??ijj?Bijfj)?Giiei?Biifi(j?i)

j?i

N??P?Bijei?Gijfi)

i?(j?i)

ij??f??

j??(Gf?Be)?Be?Gf??ijjijjiiiiii(j?i) j?i

?Bijei?Gijfi)

M??Qi

?e???(j?i)

ij?

j?(Gf?Be)?Be

??ijjijjiii?Giifi(j?i)

j?i

- 10 -

Gijei?Bijfi)?(j?i)??Qi?Lij??? ?(Ge?Bf)?Ge?Bf?fjijjiiiiii(j?i)??ijj?j?i

??U2

i?0(j?i) Rij????2e(j?i)?ej?i

??U2

i?0(j?i) Sij????2f?fji(j?i)?

- 11 -

4设计过程

4.1题目D

网络接线如图，各支路导纳和各节点功率已标么值标于图中。其中节点1是按给定功率发电的发电厂。设节点5电压保持为定值1.06不变，试分析该网络的潮流分布。方法不限，求解精度为10e-5。

手算过程见最后

- 12 -

4.2牛顿拉夫逊法程序流程图

6

- 13 -

4.3设计程序

%开始

clc

clear

disp('节点总数为：');

N=5

disp('平衡节点为：');

5

disp('PQ节点为：');

JD=[1,2,3,4]

Y=[10.834-32.5j -1.667+5j -1.667+5j -2.5+7.5j -5+15j; -1.667+5j 12.917-38.75j -10+30j 0 -1.25+3.75j; -1.667+5j -10+30j 12.917-38.75j -1.25+3.75j 0; -2.5+7.5j 0 -1.25+3.75j 3.75-11.25j 0;

-5+15j -1.25+3.75j 0 0 6.25-18.75j];

disp('导纳矩阵为:Y=');%求出导纳矩阵

ei=[1;1;1;1;1.06];fi=[0;0;0;0;0];P=[0.2;-0.45;-0.40;-0.60;0];Q=[0.20;-0.15;-0.05;-0.10;0]; disp(Y);wucha=0.0001;k=0;

G=real(Y);B=imag(Y);

whilewucha>0.00001

for m=1:4

Pi=0;Qi=0;

for n=1:5

Pi=Pi+ei(m)*(G(m,n)*ei(n)-B(m,n)*fi(n))+fi(m)*(G(m,n)*fi(n)+B(m,n)*ei(n)); % Pi（0） Qi=Qi+fi(m)*(G(m,n)*ei(n)-B(m,n)*fi(n))-ei(m)*(G(m,n)*fi(n)+B(m,n)*ei(n)); I(m)= (Pi-Qi*1i)/conj(ei(m)+fi(m)*1i);%节点注入电流Ii（0）,conj:求共轭 end

DP(m)=P(m)-Pi; %有功功率的不平衡量

DQ(m)=Q(m)-Qi; %无功功率的不平衡量

H(m,m)=-B(m,m)*ei(m)+G(m,m)*fi(m)+imag(I(m));

N(m,m)=G(m,m)*ei(m)+B(m,m)*fi(m)+real(I(m));

J(m,m)=-G(m,m)*ei(m)-B(m,m)*fi(m)+real(I(m));

L(m,m)=-B(m,m)*ei(m)+G(m,m)*fi(m)-imag(I(m));

JJ(2*m-1,2*m-1)=H(m,m);JJ(2*m-1,2*m)=N(m,m);JJ(2*m,2*m-1)=J(m,m);JJ(2*m,2*m)=L(m,m);

end

for m=1:4

for n=1:4

if n==m

else

- 14 -

H(m,n)=-B(m,n)*ei(m)+G(m,n)*fi(m);

N(m,n)=G(m,n)*ei(m)+B(m,n)*fi(m);

J(m,n)=-B(m,n)*fi(m)-G(m,n)*ei(m);

L(m,n)=G(m,n)*fi(m)-B(m,n)*ei(m);

JJ(2*m-1,2*n-1)=H(m,n);JJ(2*m-1,2*n)=N(m,n);JJ(2*m,2*n-1)=J(m,n);JJ(2*m,2*n)=L(m,n);%公式P130

end

end

end

disp('雅克比矩阵JJ:');

disp(JJ);

for m=1:4

DW(2*m-1)=DP(m);DW(2*m)=DQ(m);%形成DP和DQ不平衡量的矩阵 end

DY=JJ\DW';%解修正方程式，A\B=A逆*B

wucha=max(abs(DY));%abs绝对值

disp('第M次修正方程的解DY:');%

disp(DY);

for n=1:4

fi(n)=fi(n) +DY(2*n-1); ei(n)=ei(n) +DY(2*n); %计算各节点电压新值

end

u=ei+fi*i;

disp('节点电压的第C(k)次近似值:');

disp(u);

a=abs(u);%a的模

disp('各点的电压模值a为(节点号从小到大排列):');

disp(a);

b=angle(u);%a的幅角

disp('各点的电压幅角b为(节点号从小到大排列):');

disp(b);

k=k+1;

disp('迭代次数:');

disp(k);

end

for n=1:5 %平衡功率

I(n)=Y(5,n)*u(n);

end

disp('平衡节点的功率S5;');

S5=u(5)*sum(conj(I))%计算平衡节点的功率(公式P139)

for m=1:5

for n=1:5

S(m,n)=u(m)*(conj(u(m))-conj(u(n)))*conj(-Y(m,n));%计算m-n各支路功率

S2(n,m)=u(n)*(conj(u(n))-conj(u(m)))*conj(-Y(n,m));%计算n-m各支路功率conj共轭 DS(m,n)=S(m,n)+S2(n,m);%计算各支路功率的损耗(公式P135 4-51a)

- 15 -

end

end

c=sum(DS(:))/2

disp('各支路m-n功率');disp(S)

%disp('各支路n-m功率');disp(S2)

disp('各支路功率损耗');disp(DS)

disp('总损耗');disp(c)

4.4运行结果

节点总数为：

N =

5

平衡节点为：

ans =

5

PQ节点为：

JD =

1 2 3 4

导纳矩阵为:Y=

10.8340 -32.5000i -1.6670 + 5.0000i -1.6670 + 5.0000i -2.5000 + 7.5000i -5.0000

+15.0000i

-1.6670 + 5.0000i 12.9170 -38.7500i -10.0000 +30.0000i 0 -1.2500 +

3.7500i

-1.6670 + 5.0000i -10.0000 +30.0000i 12.9170 -38.7500i -1.2500 + 3.7500i 0 -2.5000 + 7.5000i 0 -1.2500 + 3.7500i 3.7500 -11.2500i 0 -5.0000 +15.0000i -1.2500 + 3.7500i 0 0

6.2500 -18.7500i

雅克比矩阵JJ:

- 16 -

33.4000 10.5340 -5.0000 -1.6670 -5.0000 -1.6670 -7.5000 -2.5000 -11.1340 31.6000 1.6670 -5.0000 1.6670 -5.0000 2.5000 -7.5000 -5.0000 -1.6670 38.9750 12.8420 -30.0000 -10.0000 0 0

1.6670 -5.0000 -12.9920 38.5250 10.0000 -30.0000 0 0 -5.0000 -1.6670 -30.0000 -10.0000 38.7500 12.9170 -3.7500 -1.2500

1.6670 -5.0000 10.0000 -30.0000 -12.9170 38.7500 1.2500 -3.7500 -7.5000 -2.5000 0 0 -3.7500 -1.2500 11.2500 3.7500

2.5000 -7.5000 0 0 1.2500 -3.7500 -3.7500 11.2500

第M次修正方程的解DY:

-0.0473

0.0430

-0.0863

0.0154

-0.0922

0.0141

-0.1076

0.0093

节点电压的第C(k)次近似值:

1.0430 - 0.0473i

1.0154 - 0.0863i

1.0141 - 0.0922i

1.0093 - 0.1076i

1.0600

各点的电压模值a为(节点号从小到大排列):

1.0440

1.0191

1.0183

1.0151

1.0600

各点的电压幅角b为(节点号从小到大排列):

-0.0453

-0.0848

-0.0907

-0.1062

迭代次数:

1

雅克比矩阵JJ:

- 17 -

33.1594 13.0920 -5.1360 -1.9751 -5.1360 -1.9751 -7.7040 -2.9621 -12.5811 33.6083 1.9751 -5.1360 1.9751 -5.1360 2.9621 -7.7040 -4.9331 -2.1241 38.3848 16.0302 -29.5988 -12.7427 0 0

2.1241 -4.9331 -16.8890 38.0788 12.7427 -29.5988 0 0 -4.9168 -2.1516 -29.5009 -12.9078 38.1553 16.2729 -3.6876 -1.6135

2.1516 -4.9168 12.9078 -29.5009 -17.0729 38.0553 1.6135 -3.6876 -7.3011 -3.3304 0 0 -3.6505 -1.6652 11.0516 4.3956

3.3304 -7.3011 0 0 1.6652 -3.6505 -5.5956 10.8516

第M次修正方程的解DY:

-0.0004

-0.0075

0.0017

-0.0101

0.0021

-0.0108

0.0032

-0.0131

节点电压的第C(k)次近似值:

1.0355 - 0.0477i

1.0053 - 0.0846i

1.0033 - 0.0901i

0.9963 - 0.1044i

1.0600

各点的电压模值a为(节点号从小到大排列):

1.0366

1.0089

1.0074

1.0017

1.0600

各点的电压幅角b为(节点号从小到大排列):

-0.0461

-0.0839

-0.0896

-0.1044

迭代次数:

2

雅克比矩阵JJ:

- 18 -

32.9334 12.9537 -5.0977 -1.9648 -5.0977 -1.9648 -7.6466 -2.9466 -12.5850 33.3371 1.9648 -5.0977 1.9648 -5.0977 2.9466 -7.6466 -4.8857 -2.0986 38.0494 15.8301 -29.3142 -12.5897 0 0 2.0986 -4.8857 -16.6941 37.6790 12.5897 -29.3142 0 0 -4.8664 -2.1233 -29.1988 -12.7379 37.7997 16.0624 -3.6499 -1.5922

2.1233 -4.8664 12.7379 -29.1988 -16.8446 37.6305 1.5922 -3.6499 -7.2110 -3.2738 0 0 -3.6055 -1.6369 10.9774 4.3253

3.2738 -7.2110 0 0 1.6369 -3.6055 -5.4960 10.6556

第M次修正方程的解DY:

1.0e-003 *

-0.0025

-0.0892

0.0098

-0.1231

0.0128

-0.1327

0.0241

-0.1731

节点电压的第C(k)次近似值:

1.0354 - 0.0477i

1.0052 - 0.0845i

1.0032 - 0.0901i

0.9961 - 0.1044i

1.0600

各点的电压模值a为(节点号从小到大排列):

1.0365

1.0087

1.0073

1.0016

1.0600

各点的电压幅角b为(节点号从小到大排列):

-0.0461

-0.0839

-0.0896

-0.1044

迭代次数:

3

- 19 -

雅克比矩阵JJ:

32.9307 12.9524 -5.0973 -1.9646 -5.0973 -1.9646 -7.6459 -2.9464 -12.5846 33.3340 1.9646 -5.0973 1.9646 -5.0973 2.9464 -7.6459 -4.8851 -2.0984 38.0451 15.8281 -29.3106 -12.5882 0 0 2.0984 -4.8851 -16.6922 37.6740 12.5882 -29.3106 0 0 -4.8658 -2.1230 -29.1950 -12.7362 37.7951 16.0602 -3.6494 -1.5920

2.1230 -4.8658 12.7362 -29.1950 -16.8424 37.6252 1.5920 -3.6494 -7.2098 -3.2731 0 0 -3.6049 -1.6366 10.9764 4.3243

3.2731 -7.2098 0 0 1.6366

第M次修正方程的解DY:

1.0e-007 *

-0.0013

-0.1352

0.0054

-0.1876

0.0087

-0.2038

0.0232

-0.2889

节点电压的第C(k)次近似值:

1.0354 - 0.0477i

1.0052 - 0.0845i

1.0032 - 0.0901i

0.9961 - 0.1044i

1.0600

各点的电压模值a为(节点号从小到大排列):

1.0365

1.0087

1.0073

1.0016

1.0600

各点的电压幅角b为(节点号从小到大排列):

-0.0461

-0.0839

-0.0896

-0.1044

- 20 - -3.6049 -5.4951 10.6529

迭代次数: 4

平衡节点的功率S5; S5 =

1.2982 + 0.2445i c =

0.0482 + 0.1445i

各支路m-n功率

0 0.2469 + 0.0815i 0.2793 + 0.0806i 0.5489 + 0.1333i 0.0954i

-0.2431 - 0.0701i 0 0.1891 - 0.0121i 0 - 0.0677i

-0.2746 - 0.0664i -0.1887 + 0.0132i 0 0.0633 + 0.0033i -0.5370 - 0.0977i 0 -0.0630 - 0.0023i 0 0.8895 + 0.1387i 0.4087 + 0.1058i 0 0

各支路功率损耗

0 0.0038 + 0.0113i 0.0047 + 0.0142i 0.0119 + 0.0356i 0.0433i

0.0038 + 0.0113i 0 0.0004 + 0.0011i 0.0127 + 0.0381i

0.0047 + 0.0142i 0.0004 + 0.0011i 0 0.0003 + 0.0010i 0.0119 + 0.0356i 0 0.0003 + 0.0010i 0

0.0144 + 0.0433i 0.0127 + 0.0381i 0 0

总损耗

0.0482 + 0.1445i

- 21 -

-0.8751 - -0.3960 0 0 0 0.0144 + 0 0 0 0

心得体会

通过此次课程设计，使我更加扎实的掌握了有关电力网络潮流算法方面的知识，在设计过程中虽然遇到了一些问题，但经过一次又一次的思考，一遍又一遍的检查终于找出了原因所在，也暴露出了前期我在这方面的知识欠缺和经验不足。实践出真知，通过亲自动手计算，使我们掌握的知识不再是纸上谈兵。 过而能改，善莫大焉。在课程设计过程中，我和我的团队不断发现自身错误，不断改正，不断领悟，最终的获得成功。这次课程设计终于顺利完成了，在设计中遇到了很多问题，最后在老师的指导下，终于迎刃而解。

课程设计虽然是一门专业课，但是却给我很多专业知识以及专业技能上的提升，使我对抽象的理论有了具体的认识。

回顾起此课程设计，至今我仍感慨颇多，从理论到实践，在这段日子里，可以说过得很充实，不仅巩固了以前所学过的知识，而且学到了很多在书本上所没有学到过的知识。通过这次课程设计使我懂得了理论与实际相结合是很重要的，只有理论知识是远远不够的，只有把所学的理论知识与实践相结合起来，从理论中得出结论，才能真正为社会服务。

在此次课程设计过程中，也对团队精神的进行了考察，让我们在合作起来更加默契，在成功后一起体会喜悦的心情。果然是团结就是力量，只有互相之间默契融洽的配合才能换来最终完美的结果。

此次设计也让我明白了自学的重要，有什么不懂不明白的地方要及时请教或上网查询，只要认真钻研，动脑思考，动手实践，就没有弄不懂的知识，收获颇丰。

再次感谢我的指导老师和团队中每一个成员的无私帮助

- 22 -

参考文献

[1].陈珩.《电力系统稳态分析》[M],中国电力出版社，2007，第三版

[2]. 陶海英. 基于MATLAB的电力系统稳态仿真分析[J], 江西电力职业技术学院学报, 20xx年 03期

- 23 -