目录
1任务书 (2)
2.模型简介及等值电路 (3)
3.设计原理 (4)
4.修正方程的建立 (6)
5.程序流程图及MATLAB程序编写 (8)
6.结果分析 (14)
7.设计总结 (17)
8.参考文献 (17)
《电力系统分析》
课程设计任务书
题目极坐标表示的牛顿拉夫逊法潮流计算程序设计
学生姓名学号专业班
级
电气工程及其
自动化
设计内容与要求1. 设计要求
掌握MATLAB语言编程方法；理解和掌握运用计算机进行潮流计算的基本算法原理；针对某一具体电网，进行潮流计算程序设计。
其目的在于加深学生对电力系统稳态分析中课程中基本概念和计算方法的理解，培养学生运用所学知识分析和解决问题的能力。
2. 内容
1）学习并掌握MATLAB语言。
2）掌握变压器非标准变比概念及非标准变比变压器的等值电路。掌握节点导纳矩阵的概念及导纳矩阵的形成和修改方法。
3）掌握电力系统功率方程、变量和节点分类。
4）掌握利用极坐标表示的牛-拉法进行潮流计算的方法和步骤。5）选择一个某一具体电网，编制程序流程框图。
6）利用MATLAB语言编写该模型的潮流计算程序，并上机调试程序，对计算结果进行分析。
7）整理课程设计论文。
起止时间2013 年7 月 4 日至2013 年7月10日指导教师签名年月日
系（教研室）主任
签名
年月日学生签名年月日
2 模型简介及等值电路 2.1 课程设计模型：模型3
电力网络接线如下图所示，各支路阻抗标幺值参数如下：Z12=0.02+j0.06，Z13=0.08+j0.24， Z23=0.06+j0.18， Z24=0.06+j0.12， Z25=0.04+j0.12， Z34=0.01+j0.03， Z45=0.08+j0.24， k=1.1。该系统中，节点1为平衡节点，保
持1
1.060V j =+&为定值；节点2、3、4都是PQ 节点，节点5为PV 节点，给定的注入功率分别为:
20.200.20
S j =+，
3-0.45-0.15
S j =，
40.400.05
S j =--，
50.500.00S j =-+，5 1.10V =&。各节点电压（初值）标幺值参数如下： 节点
1 2 3 4 5
Ui （0）=ei （0）+jfi （0）
1.06+j0.0
1.0+j0.0 1.0+j0.0
1.0+j0.0
1.1+j0.0
计算该系统的潮流分布。计算精度要求各节点电压修正量不大于10-5。
图2-1
2.2 模型分析
节点1是平衡节点，节点2、3、4是PQ节点，节点5是PV节点。由题可得等值电路模型中各节点之间的导纳：y12=5.000-j15.000，y13=1.2500-j3.7500，y22=0.2750-j0.8250，y23=1.667-j5.000，y24=3.333-j6.667，y25=2.7500-j8.2500，y34=10.0000-j30.0000，y55=-0.25+j0.75
2.3 等值电路模型
在图2-2中，将图2-1中的编号重新编排，节点2、3、4、5、1替换为1、2、3、4、5。则各节点之间的导纳变为y12=1.667-j5，y13=3.333-j6.667，y14=2.75-j8.25，y15=5-j15，y52=1.25-j3.75，y23=10-j30，y34=1.25-j3.75，y11=0.275-j0.825，y44=-0.25+j0.75。
3 设计原理
本题采用了题目要求的牛顿－拉夫逊潮流计算的方法。牛顿法是数学中解决非线性方程式的典型方法，有较好的收敛性。
3.1潮流计算的定解条件
题中所给图表示一个五节点的简单电力系统，n 个节点电力系统的潮流方程的一般形式是
.
.
1
n i i
j ij j i
P jQ V Y V =-=∑ (1,2,3,...)i n =
或
.
.
1
n j i i i ij j V P jQ V Y =+=∑ (1,2,3,...)i n =
按变量的不同，一般将节点分为三种类型。 （1） PQ 节点
这类节点的有功功率和无功功率是给定的，节点(,)V δ是待求量。通常变电所都是这一类型节点。由于没有发电设备，故其发电功率为零。属于这一类节点的有按给定有功、无功功率发电的发电厂母线和没有其他电源的变电所母线。 （2） PV 节点
这类节点有功功率P 和电压幅值V 是给定的，节点的无功功率Q 和电压的相位δ是待求量。这类节点必须有足够的可调无功容量，用以维持给定的电压幅值，因此又称为电压控制节点。一般选择有一定无功储备的发电厂和具有可调无功电源设备的变电所作为PV 节点。 （3） 平衡节点
在潮流分布算出以前，网络中的功率损耗是未知的，因此，网络中至少有一个节点的有功功率P 不能给定，这个节点承担了系统系统的有功功率平衡，故称之为平衡节点。潮流计算时，一般只设一个平衡节点。
3.2潮流计算的约束条件
(1)所有节点电压必须满足
min max i i V V V ≤≤
(2)所有电源节点的有功功率和无功功率必须满足
min max
min max
}Gi Gi Gi Gi Gi Gi P P P Q Q Q ≤≤≤≤
(3)某些节点之间的电压应满足
max
i j i j
δδδδ-<-
3.3牛顿-拉夫逊的基本原理
设欲求解的非线性代数方程为
()0f x =
设方程的真实解为x ，则必有()0f x =。用牛顿-拉夫逊法求方程真实解x 的步骤如下：
首先选取余割合适的初始估值(0)
x 作为方程
(0)()0f x =的解，若恰巧有(0)()0f x =，则方程的真实解即为(0)x x =若(0)()0f x ≠，则做下一步。取
(1)(0)(0)x x x =+? 则(1)(0)(0)
()()f x f x x =+?
其中(0)x ?为初始估值的增量，即(0)
(1)(0)x
x x ?=-。设函数()f x 具有任意阶导数
(1)(0)(0)(0)'(0)(0)''(0)(0)2()()()()[()]/2...f x f x x f x f x x f x x =+?=+?+?+ 若所取的
(0)
x ?足够小，则含
(0)2()x ?的项及其余的一切高阶项均可略去，并使其等于零，即：
(1)(0)(0)(0)'(0)(0)
()()()()0f x f x x f x f x x =+?≈+?= 故得
(0)(0)
'(0)
()
()f x x
f x ?=-
可见，只要'(0)()f x ≠0，即可根据上式求出第一次的修正估值(1)x ，若恰巧有
(1)()f x =0，则方程的真实解即为(1)x x =。若(1)
()0f x ≠，则用上述方法由(1)x 再
确定第二次的修正估值。如此反复叠代下去，直到求得真实解x 为止。设第K 次的估值为第（K+1）次的修正估值，则有
()(1)
()
'()
()()k k k k f x x
x
f x +=-
迭代过程的收敛数据为
1
()()k f x ε<或2
()k x ε?<
其中，1
ε
，
2
ε
为预先给定的小正数。
4 修正方程的建立
极坐标表示的牛拉法修正方程为：
??????????????????????????
????????????????????????????????????
???
?=??????????????????????????????????????????????δδδδn P
nn np n n n n pn pp P P P P n p
n p n p
n p n P U U U U H H N H N H H H N H N H J J L J L J H H N H N H J J L J L J H H N H N H P P Q P Q P 222111
22112211
22222221212222222121111212111111121211112211//
式中留出了（n-m ）行空格和（n-m ）列空列。式中的有功，无功功率不平衡量Q P i
i
??.
分别由式（3-1a ）,式（3-1b ）可得为
∑==+-=?n
j j ij ij ij ij i i B G U U P P 1j i sin cos ）
（δδ ∑==--=?n j j ij ij ij ij i i B G U U Q Q 1
j i cos sin ）
（δδ （3-1a ，b ）
而式中雅可比矩阵的各个元素则分别为 j i P H δ??=
ij ；j j i ij U U P N ??=；j i Q J δ??=ij ; j j
i
ij U U Q L ??= （3-2） 修正方程中将i U ?改为/i i U U ?只是为使公式（3-2）中个偏导数的表示形式上更相似，为求取这些偏导数，可将i P 、i Q 分别展开如下
()
()
2121cos sin sin cos j n
i i
ii i j ij ij ij ij j j i j n
i i
ii i j ij ij ij ij j j i
P U G U U G B Q U B U U G B δδδδ==≠==≠=++=-+-∑∑ （3-3a,b ）
计及
cos cos()
sin()sin ()sin sin()
cos()cos ()cos cos()
sin()sin ()
sin sin()cos()cos ()ij i j i j ij
j i j ij i j i j ij j i j ij i j i j ij i i j ij i j i j ij i i j δδδδδδδδδδδδδδδδδδδδδδδδδδδδδδδδδδδδ??-?=
=-=?-?-??-==--=-?-?-?
??-==--=-??-??-==-=??-?
?
??
??
??
?
??
（3-4）
j i ≠时，由于对特定的j ，只有该特定节点的j δ，从而特定的ij i j δδδ=-是
变量，由式（3-2）到式（3-4）可得
()()sin cos cos sin i ij i j ij ij ij ij j i ij i j ij ij ij ij i P H U U G B Q
J U U G B δδδδδδ??
=
=-???
???
==-+???
（3-5a ）
相似的，由于对特定的j ，只有该特定节点的j U 是变量，可得
()()cos sin sin cos i ij i j ij ij ij ij i i
ij i j ij ij ij ij j
P N U U G B U Q L U U G B U δδδδ??
=
=+???
???
==-??? （3-5b ） j=i 时，由于i δ是变量，所有ij i j δδδ=-都是变量，可得
()()1
1
sin cos cos sin j n
i
ii i j ij ij ij ij j i j i j n
i
ii i j ij ij ij ij j i j i P H U U G B Q J U U G B δδδδδδ==≠==≠??==--?
??????
==+????
∑∑ （3-5c ） 相似的，由于i U 是变量，可得
()()21
21
cos sin 2sin cos 2j n
i
ii i j ij ij ij ij i ii j i j i j n
i
ii i j ij ij ij ij i ii j i j i P N U U G B U G U Q L U U G B U B U δδδδ==≠==≠??==++?
??????
==--????∑∑ （3-5d ）
5 设计流程图及程序的编写
输入数据
形成节点导纳矩阵
设定节点起始计算电压δ U
置迭代次数
应用公式计算不平衡量△P△Q
雅克比矩阵J是否已经全部
形成，i>n?
按公式计算雅克比矩阵元素J
增大节点号，i=i+1
解修正方程，由△P,△Q和J计算电压修正量△u和△δ
求出△u和△δ
迭代是否收敛，△u，△δ﹤ε
计算平衡节点功率S和线路的功率S
结束
K=k+1
△u=U+△u △δ=δ+△δ
置节点号i=1
程序中用到的符号所代表的意义：
Y 代表导纳矩阵
JJ 代表雅克比矩阵
?
pp代表有功功率的不平衡量P i
?
qq 代表无功功率的不平衡量Q i
?U
uu 代表各节点电压和相角的不平衡量δ?
U 代表各节点的电压
S 代表线路的功率
Q4 代表平衡节点的注入无功功率
K 代表迭代次数
N1 代表PQ节点和PV节点的总数
m,n代表系统中的节点总数，把平衡节点标为最大号
S5 平衡节点功率
d 变压器对地导纳
%The following program for load calculation is based on MATLAB6.5 %以下部分为输入原始数据（到标示‘///’标志为止）。
%求取节点导纳矩阵
g(1,1)=0.275;b(1,1)=0.825;
g(1,2)=1.667;b(1,2)=-5;
g(1,3)=3.333;b(1,3)=-6.667;
g(1,4)=2.75;b(1,4)=-8.25;
g(1,5)=5;b(1,5)=-15;
g(2,1)=1.667;b(2,1)=-5;
g(2,2)=-0;b(2,2)=0;
g(2,3)=10;b(2,3)=-30;
g(2,4)=0;b(2,4)=0;
g(2,5)=1.25;b(2,5)=-3.75;
g(3,1)=3.33;b(3,1)=-6.667;
g(3,2)=10;b(3,2)=-30;
g(3,3)=0;b(3,3)=0;
g(3,4)=1.25;b(3,4)=-3.75;
g(3,5)=0;b(3,5)=0;
g(4,1)=2.75;b(4,1)=-8.25;
g(4,2)=0;b(4,2)=0;
g(4,3)=1.24;b(4,3)=-3.75;
g(4,4)=-0.25;b(4,4)=0.75;
g(4,5)=0;b(4,5)=0;
g(5,1)=5;b(5,1)=-15;
g(5,2)=1.25;b(5,2)=-3.75;
g(5,3)=0;b(5,3)=0;
g(5,4)=0;b(5,4)=0;
课设报告 福建工程学院软件学院 题目：汇编计算器 班级： 1301 姓名 学号： 指导老师： 日期：
目录 1、设计目的 (3) 2、概要设计 (3) 2.1 系统总体分析 (3) 2.2 主模块框图及说明 (3) 3、详细设计 (4) 3.1 主模块及子模块概述 (4) 3.2各模块详运算 (4) 4、程序调试 (7) 4.1 运行界面分析 (7) 算法分析 (7) 4.2 调试过程与分析 (9) 5、心得体会 (11) 5.1 设计体会 (11) 5.2 系统改进 (11) 附录： (11)
1、设计目的 本课程设计是一次程序设计方法及技能的基本训练，通过实际程序的开发及调试，巩固课堂上学到的关于程序设计的基本知识和基本方法，进一步熟悉汇编语言的结构特点和使用，达到能独立阅读、设计编写和调试具有一定规模的汇编程序的水平。 2、概要设计 用8086汇编语言编写一个能实现四则混合运算、带括号功能的整数计算器程序。程序能实现键盘十进制运算表达式的输入和显示（例如输入：“1+2*(3-4)”），按“=”后输出十进制表示的运算结果。 2.1 系统总体分析 在8086的操作环境下，该计算器分成输入，数据存储，运算功能，输出几个大模块，实现了使用者使用该计算器时输入一个算式，能让系统进行计算。此计算器的实现功能是基本的数学的四则运算，结果范围在0~65535。 2.2 主模块框图及说明 此流程图简要的表现出了所要实现的功能以及一些功能的大概算法，同时也是我编写的一个总体的框架。 程序流程图说明：通过流程图，可以看出程序运行时，首先输出提示语气，当用户输入后，程序根据所输入内容进行判断，通过判断的结果来决定调用哪个功能模块，首要先要要判断的是否为0-9，“+”“-”“*”“/”这些字符，若不是就会报错，实则根据运算符号调用其功能模块完成运算。最后将运算的结果显示在主频幕上，返回主程序，使用户可以重新输入。
摘要 本文，首先简单介绍了基于在MALAB中行潮流计算的原理、意义，然后用具体的实例，简单介绍了如何利用MALAB去进行电力系统中的潮流计算。 众所周知，电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各线的电压、各元件中流过的功率、系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性、可靠性和经济性。 此外，在进行电力系统静态及暂态稳定计算时，要利用潮流计算的结果作为其计算的基础；一些故障分析以及优化计算也需要有相应的潮流计算作配合；潮流计算往往成为上述计算程序的一个重要组成部分。以上这些，主要是在系统规划设计及运行方式安排中的应用，属于离线计算范畴。 牛顿－拉夫逊法在电力系统潮流计算的常用算法之一，它收敛性好，迭代次数少。本文介绍了电力系统潮流计算机辅助分析的基本知识及潮流计算牛顿－拉夫逊法，最后介绍了利用MTALAB程序运行的结果。 关键词：电力系统潮流计算，牛顿－拉夫逊法，MATLAB
ABSTRACT This article first introduces the flow calculation based on the principle of MALAB Bank of China, meaning, and then use specific examples, a brief introduction, how to use MALAB to the flow calculation in power systems. As we all know, is the study of power flow calculation of power system steady-state operation of a calculation, which according to the given operating conditions and system wiring the entire power system to determine the operational status of each part: the bus voltage flowing through the components power, system power loss and so on. In power system planning power system design and operation mode of the current study, are required to quantitatively calculated using the trend analysis and comparison of the program or run mode power supply reasonable, reliability and economy. In addition, during the power system static and transient stability calculation, the results of calculation to take advantage of the trend as its basis of calculation; number of fault analysis and optimization also requires a corresponding flow calculation for cooperation; power flow calculation program often become the an important part. These, mainly in the way of system design and operation arrangements in the application areas are off-line calculation. Newton - Raphson power flow calculation in power system is one commonly used method, it is good convergence of the iteration number of small, introduce the trend of computer-aided power system analysis of the basic knowledge and power flow Newton - Raphson method, introduced by the last matlab run results. Keywords：power system flow calculation, Newton – Raphson method, matlab
课程设计报告 学生姓名:学号： 学院:电气工程学院 班级: 题目: 电力系统潮流计算 职称: 副教授 指导教师：李翠萍职称: 副教授 2014年 01月10日
1 潮流计算的目的与意义 潮流计算的目的：已知电网的接线方式与参数及运行条件，计算电力系统稳态运行各母线电压、个支路电流与功率及网损。对于正在运行的电力系统，通过潮流计算可以判断电网母线电压、支路电流和功率是否越限，如果有越限，就应采取措施，调整运行方式。对于正在规划的电力系统，通过潮流计算，可以为选择电网供电方案和电气设备提供依据。潮流计算还可以为继电保护和自动装置定整计算、电力系统故障计算和稳定计算等提供原始数据。 潮流计算的意义： (1)在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求。 (2)在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。 (3)正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。 (4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。 2 潮流计算数学模型 1.变压器的数学模型： 变压器忽略对地支路等值电路：
2.输电线的数学模型： π型等值电路： 3 数值方法与计算流程 利用牛顿拉夫逊法进行求解，用MATLAB 软件编程，可以求解系统潮流分 布根据题目的不同要求对参数进行调整，通过调节变压器变比和发电厂的电压，求解出合理的潮流分布，最后用matpower 进行潮流分析，将两者进行比较。 牛顿—拉夫逊法 1、牛顿—拉夫逊法概要 首先对一般的牛顿—拉夫逊法作一简单的说明。已知一个变量X 函数为： 0)(=X f 到此方程时，由适当的近似值) 0(X 出发，根据： ,......)2,1() ()() ()() () 1(='-=+n X f X f X X n n n n 反复进行计算，当) (n X 满足适当的收敛条件就是上面方程的根。这样的方 法就是所谓的牛顿—拉夫逊法。 这一方法还可以做下面的解释，设第n 次迭代得到的解语真值之差，即) (n X 的误差为ε时，则： 0)()(=+εn X f 把)() (ε+n X f 在) (n X 附近对ε用泰勒级数展开 0......)(! 2)()()()(2 )() () (=+''+ '+=+n n n n X f X f X f X f εεε 上式省略去2ε以后部分 0)()()()(≈'+n n X f X f ε
前言 电力系统潮流计算是研究电力系统稳态运行情况的一种计算,它根据给定的运行件及系统接线情况确定整个电力系统各部分的运行状态。在电力系统规划设计和现有电力系统运行方式的研究中,都需要利用潮流计算来定量分析、比较供电方案或运行方式的合理性、可靠性和经济性。本次课程设计任务是闭式网络的潮流计算，用到的方法为牛顿拉夫逊极坐标法潮流计算。 牛顿法是数学中解决非线性方程式的典型方法，有较好的收敛性。解决电力系统潮流计算问题是以导纳距阵为基础的，因此，只要在迭代过程中尽可能保持方程式系数距阵的稀疏性，就可以大大提高牛顿法潮流程序的放率。自从20 世纪60 年代中期利用了最佳顺序消去法以后，牛顿法在收敛性、内存要求、速度方面都超过了阻抗法，成为直到目前仍在广泛采用的优秀方法。
目录 1任务书 (2) 2.模型简介及等值电路 (3) 3.设计原理 (4) 4.修正方程的建立 (7) 5.程序流程图及MATLAB程序编写 (9) 6.结果分析 (15) 7.设计总结 (25) 8.参考文献 (26)
《电力系统分析》 课程设计任务书 题目极坐标表示的牛顿拉夫逊法潮流计算程序设计学生姓名学号专业班级
设计内容与要求1. 设计要求 掌握MATLAB语言编程方法；理解和掌握运用计算机进行潮流计算的基本算法原理；针对某一具体电网，进行潮流计算程序设计。 其目的在于加深学生对电力系统稳态分析中课程中基本概念和计算方法的理解，培养学生运用所学知识分析和解决问题的能力。 2. 内容 1）学习并掌握MATLAB语言。 2）掌握变压器非标准变比概念及非标准变比变压器的等值电路。掌握节点导纳矩阵的概念及导纳矩阵的形成和修改方法。 3）掌握电力系统功率方程、变量和节点分类。 4）掌握利用极坐标表示的牛-拉法进行潮流计算的方法和步骤。 5）选择一个某一具体电网，编制程序流程框图。 6）利用MATLAB语言编写该模型的潮流计算程序，并上机调试程序，对计算结果进行分析。 7）整理课程设计论文。 起止时间2013 年7 月 4 日至2013 年7月10日指导教师签名年月日 系（教研室）主任签 名 年月日学生签名年月日 2 模型简介及等值电路 2.1课程设计模型：模型3
课程设计报告 学生:学号： 学院: 班级: 题目: 电力系统潮流计算课程设计
课设题目及要求 一 .题目原始资料 1、系统图：两个发电厂分别通过变压器和输电线路与四个变电所相连。 2、发电厂资料： 母线1和2为发电厂高压母线，发电厂一总装机容量为（ 300MW ），母线3为机压母线，机压母线上装机容量为（ 100MW ），最大负荷和最小负荷分别为50MW 和20MW ；发电厂二总装机容量为（ 200MW ）。 3、变电所资料： (一) 变电所1、2、3、4低压母线的电压等级分别为：35KV 10KV 35KV 10KV (二) 变电所的负荷分别为： 60MW 40MW 40MW 50MW (三) 每个变电所的功率因数均为cos φ=0.85； 变电所1 变电所母线 电厂一 电厂二
(四) 变电所1和变电所3分别配有两台容量为75MVA 的变压器，短路损 耗414KW ，短路电压（%）=16.7；变电所2和变电所4分别配有两台容 量为63MVA 的变压器，短路损耗为245KW ，短路电压（%）=10.5； 4、输电线路资料： 发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为Ω17.0，单位长度的电抗为Ω0.402，单位长度的电纳为S -610*2.78。 二、 课程设计基本容： 1. 对给定的网络查找潮流计算所需的各元件等值参数，画出等值电路图。 2. 输入各支路数据，各节点数据利用给定的程序进行在变电所在某一负荷 情况下的潮流计算，并对计算结果进行分析。 3. 跟随变电所负荷按一定比例发生变化，进行潮流计算分析。 1) 4个变电所的负荷同时以2%的比例增大； 2) 4个变电所的负荷同时以2%的比例下降 3) 1和4号变电所的负荷同时以2%的比例下降，而2和3号变电所的 负荷同时以2%的比例上升； 4. 在不同的负荷情况下，分析潮流计算的结果，如果各母线电压不满足要 求，进行电压的调整。（变电所低压母线电压10KV 要求调整围在9.5-10.5 之间；电压35KV 要求调整围在35-36之间） 5. 轮流断开支路双回线中的一条，分析潮流的分布。（几条支路断几次） 6. 利用DDRTS 软件，进行绘制系统图进行上述各种情况潮流的分析，并进 行结果的比较。 7. 最终形成课程设计成品说明书。 三、课程设计成品基本要求： 1. 在读懂程序的基础上画出潮流计算基本流程图 2. 通过输入数据，进行潮流计算输出结果 3. 对不同的负荷变化，分析潮流分布，写出分析说明。 4. 对不同的负荷变化，进行潮流的调节控制，并说明调节控制的方法，并 列表表示调节控制的参数变化。 5. 打印利用DDRTS 进行潮流分析绘制的系统图，以及潮流分布图。
目录 一、设计任务和性能指标 (1) 1.1设计任务 (2) 1.2性能指标 (2) 二、设计方案 (2) 3 3 4 5 5 6 6 7 7 20 20 20 20 21 参考文献 (21) 附录1、系统硬件电路图 (22) 附录2、硬件实物图 (23) 附录3、器件清单 (24)
一、设计任务和性能指标 1.1设计任务 自制一个单片机最小系统，包括复位电路，采用外部小键盘输入数据，能够实现加法、乘法及一个科学计算，计算结果显示在四位一体的数码管上。 要求用Protel 画出系统的电路原理图（要求以最少组件，实现系统设计所要 显 位 监测模块采用二极管和扬声器（实验室用二极管代替）组成电路。 键盘电路采用4*4矩阵键盘电路。 显示模块采用4枚共阳极数码管和74ls273锁存芯片构成等器件构成。 整个单片机的接口电路： P0用于显示输出； P1用于键扫描输入； P2用于数码管位选控制； P3用于键盘扩展（部分运算符输入）；
三.系统硬件设计 3.1单片机最小系统 单片机最小系统就是支持主芯片正常工作的最小电路部分，包括主控芯片、复位电路和晶振电路。 主控芯片选取STC89C52RC芯片，因其具有良好的性能及稳定性，价格便宜应用方便。 扩展键：“log”，“ln”，“x^2”“小数点”，“开方” 共计25个按键，采用4*4矩阵键盘，键盘的行和列之间都有公共端相连，四行和四列的8个公共端分别接P1.0~P1.7，这样扫描P1口就可以完成对矩阵键盘的扫描，通过对16个按键进行编码，从而得到键盘的口地址，对比P1口德扫描结果和各按键的地址，我们就可以得到是哪个键按下，从而完成键盘的功能。 以下为键盘接口电路的硬件电路图
东北电力大学课程设计改革试用任务书: 电力系统潮流计算课程设计任务书 设计名称：电力系统潮流计算课程设计 设计性质：理论计算，计算机仿真与验证 计划学时：两周 一、设计目的 1.培养学生独立分析问题、解决问题的能力； 2.培养学生的工程意识，灵活运用所学知识分析工程问题的能力 3.编制程序或利用电力系统分析计算软件进行电力系统潮流分析。 二、原始资料 1、系统图：IEEE14节点。 2、原始资料：见IEEE14节点标准数据库 三、课程设计基本内容： 1.采用PSAT仿真工具中的潮流计算软件计算系统潮流； 1)熟悉PSAT仿真工具的功能； 2)掌握IEEE标准数据格式内容； 3)将IEEE标准数据转化为PSAT计算数据； 2.分别采用NR法和PQ分解法计算潮流，观察NR法计算潮流中雅可比矩阵的变化情况， 分析两种方法计算潮流的优缺点； 3.分析系统潮流情况，包括电压幅值、相角，线路过载情况以及全网有功损耗情况。
4.选择以下内容之一进行分析： 1)找出系统中有功损耗最大的一条线路，给出减小该线路损耗的措施，比较各种措施 的特点，并仿真验证； 2)找出系统中电压最低的节点，给出调压措施，比较各种措施的特点，并仿真验证； 3)找出系统中流过有功功率最大的一条线路，给出减小该线路有功功率的措施，比较 各种措施的特点，并仿真验证； 5.任选以下内容之一作为深入研究：（不做要求） 1)找出系统中有功功率损耗最大的一条线路，改变发电机有功出力，分析对该线路有 功功率损耗灵敏度最大的发电机有功功率，并进行有效调整，减小该线路的损耗； 2)找出系统中有功功率损耗最大的一条线路，进行无功功率补偿，分析对该线路有功 功率损耗灵敏度最大的负荷无功功率，并进行有效调整，减小该线路的损耗； 3)找出系统中电压最低的节点，分析对该节点电压幅值灵敏度最大的发电机端电压， 并有效调整发电机端电压，提高该节点电压水平； 四、课程设计成品基本要求： 1.绘制系统潮流图，潮流图应包括： 1)系统网络参数 2)节点电压幅值及相角 3)线路和变压器的首末端有功功率和无功功率 2.撰写设计报告，报告内容应包括以下几点： 1)本次设计的目的和设计的任务； 2)电力系统潮流计算的计算机方法原理，分析NR法和PQ分解法计算潮流的特点； 3)对潮流计算结果进行分析，评价该潮流断面的运行方式安全性和经济性； 4)找出系统中运行的薄弱环节，如电压较低点或负载较大线路，给出调整措施； 5)分析各种调整措施的特点并比较它们之间的差异； 6)结论部分以及设计心得； 五、考核形式 1.纪律考核：学生组织出勤情况和工作态度等； 2.书面考核：设计成品的完成质量、撰写水平等； 3.答辩考核：参照设计成品，对计算机方法进行电力系统潮流计算的相关问题等进行答辩； 4.采用五级评分制：优、良、中、及格、不及格五个等级。
科学计算器课程设计报告C课程设计 集团标准化工作小组 [Q8QX9QT-X8QQB8Q8-NQ8QJ8-M8QMN]
计算机科学与技术学部 C++课程设计 题目科学计算器 学部计算机科学与技术 班级计科1103 指导教师李军 姓名刘明 学号 2012年6月27日
摘要 计算器的产生和发展是建立在电子计算机基础之上的。硬件方面，自1946年第一台电子计算机诞生以来，计算机技术的发展可谓日新月异，从庞大的只能在实验室里供研究使用的计算机到如今能适应不同环境满足不同需求的各种各样的计算机；运算速度从每秒几千次到每秒几百亿次；处理器从焊有上百万个电子管的大的惊人的电子板到只有指甲大小的集成电路；现在计算机在硬件方面的发展已达到了每三个月更新换代一次的惊人速度。软件方面，也已从机器语言、汇编语言、高级语言发展到现如今的第四代语言——非结构化、面向对象、可视化的语言。 在这个计算器里面，我们实现了简单的四则运算以及更高功能的科学计算，它的外观简洁美观，使人们能快捷简单的操作。能准确的得到计算结果，大大减少了数字计算所需要的时间，为人们的生活带来便利。此系统在Windows 7环境下，使用VC++ 进行编写。 简单计算器包括双目运算和单目运算功能，双目运算符包含基本的四则运算及乘幂功能，单目运算符包含正余弦，对数，开方，阶乘，倒数，进制转换等运算。可对其输入任意操作数，包括小数和整数及正数和负数进行以上的所有运算并能连续运算。并且包含清除，退格功能等。我们所做的计算器其功能较Windows 7下的计算器还是很不够多，没有其菜单的实现功能项，没有其小巧的标准计算器。 关键词：计算器；运算；VC++等
>> %本程序的功能是用牛顿拉夫逊法进行潮流计算n=input('请输入节点数:n='); nl=input('请输入支路数:nl='); isb=input('请输入平衡母线节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n); e=zeros(1,n);f=zeros(1,n);V=zeros(1,n); O=zeros(1,n);S1=zeros(nl); for i=1:nl if B1(i,6)==0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2; Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; end %求导纳矩阵 disp('导纳矩阵Y='); disp(Y); G=real(Y);B=imag(Y); for i=1:n e(i)=real(B2(i,3)); f(i)=imag(B2(i,3)); V(i)=B2(i,4); end for i=1:n S(i)=B2(i,1)-B2(i,2); B(i,i)=B(i,i)+B2(i,5); end P=real(S);Q=imag(S); ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; while IT2~=0 IT2=0;a=a+1; for i=1:n if i~=isb C(i)=0; D(i)=0; for j1=1:n C(i)= C(i)+G(i,j1)*e(j1)-B(i,j1)*f(j1); D(i)= D(i)+G(i,j1)*f(j1)+B(i,j1)*e(j1); end
目录 一、设计任务 (1) 1.1 课程设计要求 (1) 1.2 课程设计题目 (1) 1.3 课程设计基本容 (2) 二、问题分析 (3) 2.1 节点设置及分类 (3) 2.2 参数求取 (3) 2.3 计算方法 (4) 三、问题求解 (7) 3.1 等值电路的计算 (7) 3.2画出系统等值电路图： (7) 3.3 潮流计算 (8) 四、误差分析 (29) 五、心得体会及总结 (38) 附录: (39) 参考文献 (39) 程序 (39)
电力系统潮流计算课程设计 一、设计任务 1.1 课程设计要求 1、在读懂程序的基础上画出潮流计算基本流程图 2、通过输入数据，进行潮流计算输出结果 3、对不同的负荷变化，分析潮流分布，写出分析说明。 4、对不同的负荷变化，进行潮流的调节控制，并说明调节控制的方法，并 列表表示调节控制的参数变化。 5、打印利用DDRTS进行潮流分析绘制的系统图，以及潮流分布图。 1.2 课程设计题目 系统图：两个发电厂分别通过变压器和输电线路与四个变电所相连。 变电所1 变电所2 母线电厂一电厂二
发电厂资料： 母线1和2为发电厂高压母线，发电厂一总装机容量为（ 300MW ），母线3为机压母线，机压母线上装机容量为（ 100MW ），最大负荷和最小负荷分别为50MW 和20MW ；发电厂 二总装机容量为（ 200MW ）。 变电所资料： (一)变电所1、2、3、4低压母线的电压等级分别为：35KV 10KV 35KV 10KV (二)变电所的负荷分别为：60MW 40MW 70MW 50MW (三)每个变电所的功率因数均为cos φ=0.85； (四)变电所1和变电所3分别配有两台容量为75MVA 的变压器，短路损耗414KW ，短路 电压（%）=16.7；变电所2和变电所4分别配有两台容量为63MVA 的变压器，短路损耗为245KW ，短路电压（%）=10.5； 输电线路资料： 发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为 Ω17.0，单位长度的电抗为Ω0.402，单位长度的电纳为S -610*2.78。 1.3 课程设计基本容 1. 对给定的网络查找潮流计算所需的各元件等值参数，画出等值电路图。 2. 输入各支路数据，各节点数据利用给定的程序进行在变电所在某一负荷情况下的潮 流计算，并对计算结果进行分析。 3. 跟随变电所负荷按一定比例发生变化，进行潮流计算分析。 1) 4个变电所的负荷同时以2%的比例增大； 2) 4个变电所的负荷同时以2%的比例下降 3) 1和4号变电所的负荷同时以2%的比例下降，而2和3号变电所的负荷同时以 2%的比例上升； 4. 在不同的负荷情况下，分析潮流计算的结果，如果各母线电压不满足要求，进行电 压的调整。（变电所低压母线电压10KV 要求调整围在9.5-10.5之间；电压35KV 要求调整围在35-36之间） 5. 轮流断开支路双回线中的一条，分析潮流的分布。（几条支路断几次） 6. 利用DDRTS 软件，进行绘制系统图进行上述各种情况潮流的分析，并进行结果的比 较。 7. 最终形成课程设计成品说明书。
高级语言程序（JAVA）课程设计报告 系部名称：商学系专业班级：营销*** 学生姓名：墨璇 墨兰学号： ********** ********** 指导教师：王芬教师职称：讲师 2014年06月26日
目录 一、课程设计目的及意义 .................................... 错误!未定义书签。 二、课程设计任务 .......................................... 错误!未定义书签。 2.1程序设计要求....................................... 错误!未定义书签。 三、课程设计时间 .......................................... 错误!未定义书签。 四、课程设计地点 .......................................... 错误!未定义书签。 五、课程设计内容 .......................................... 错误!未定义书签。 5.1开发工具与平台..................................... 错误!未定义书签。 (1).开发工具 ...................................... 错误!未定义书签。 (2).开发平台 ...................................... 错误!未定义书签。 5.2设计思路........................................... 错误!未定义书签。 5.3 程序测试 .......................................... 错误!未定义书签。 5.4实验总结........................................... 错误!未定义书签。 六、课程设计感想 .......................................... 错误!未定义书签。 七、附录（程序代码） ...................................... 错误!未定义书签。
课程设计说明书 题目电力系统分析系(部) 专业(班级) 姓名 学号 指导教师 起止日期
电力系统分析课程设计任务书系（部）：专业：指导教师：
目录 一、潮流计算基本原理 1.1潮流方程的基本模型 1.2潮流方程的讨论和节点类型的划分 1.3、潮流计算的意义 二、牛顿－拉夫逊法 2.1牛顿-拉夫逊法基本原理 2.2节点功率方程 2.3修正方程 2.4牛顿法潮流计算主要流程 三、收敛性分析 四、算例分析 总结 参考文献
电力系统分析潮流计算 一、潮流计算基本原理 1.1潮流方程的基本模型 电力系统是由发电机、变压器、输电线路及负荷等组成，其中发电机及负荷是非线性元件，但在进行潮流计算时，一般可以用接在相应节点上的一个电流注入量来代表。因此潮流计算所用的电力网络系由变压器、输电线路、电容器、电抗器等静止线性元件所构成，并用集中参数表示的串联或并联等值支路来模拟。结合电力系统的特点，对这样的线性网络进行分析，普通采用的是节点法，节点电压与节点电流之间的关系 V Y I = （1－1） 其展开式为 j n j ij i V Y I ∑==1 ),,3,2,1 (n i = （1－2） 在工程实际中，已经的节点注入量往往不是节点电流而是节点功率，为此必须应用联 系节点电流和节点功率的关系式 i i i i V jQ P I * -= ),,3,2,1(n i = （1－3） 将式（1－3）代入式（1－2）得到 j n j ij i i i V Y V jQ P ∑=* =-1 ),,3,2,1(n i = （1－4） 交流电力系统中的复数电压变量可以用两种极坐标来表示 i j i i e V V θ= （1－5） 或 i i i jf e V += （1－6） 而复数导纳为
目录 摘要................................................. - 1 - 1.设计意义与要求..................................... - 2 - 1.1设计意义 ...................................... - 2 - 1.2设计要求（具体题目）........................... - 2 - 2.题目解析........................................... - 3 - 2.1设计思路 ...................................... - 3 - 2.2详细设计 ...................................... - 4 - 2.2.1节点类型.................................. - 4 - 2.2.2待求量 ................................... - 4 - 2.2.3导纳矩阵.................................. - 4 - 2.2.4潮流方程.................................. - 5 - 2.2.5牛顿—拉夫逊算法.......................... - 6 - 2.2.5.1牛顿算法数学原理：................... - 6 - 2.2.5.2修正方程............................. - 7 - 2.2.5.3收敛条件............................. - 9 - 3.结果分析.......................................... - 10 - 4.小结.............................................. - 11 - 参考文献............................................ - 12 -
计算机技术综合课程设计 设计题目锅炉液位控制系统学生姓名史婷艳 专业班级自动化1302班学号20134460203 指导老师洪镇南 2017年1 月3日
目录 前言 (2) 1 锅炉汽包水位控制对象与控制指标 (4) 1.1锅炉汽包水位的特征 (4) 1.2汽包水位动态特性 (4) 1.2.1汽包水位在给水流量W作用下的动态特性 (4) 1.2.2汽包水位在蒸汽流量D扰动下的动态特性 (5) 1.2.3燃料量B扰动下汽包水位的动态特性 (6) 2. 汽包水位控制方案 (7) 2.1单冲量控制方式 (7) 2.2 双冲量控制方式 (8) 2.3 三冲量控制方式 (9) 3. 三冲量串级PID控制 (11) 3.1 串级PID控制 (11) 3.2 智能整定PID控制 (12) 4 汽包水位模糊控制器设计及仿真 (12) 4.1 输入输出变量 (12) 4.2 隶属度函数 (15) 4.3基于MATLAB/Simulink 环境建立的系统仿真分析 (16) 4.3.1 基于MATLAB/Simulink 的系统模型 (16)
4.3.2 仿真结果分析 (18) 总结与体会 (18) 参考文献 (20) 前言 锅炉是典型的复杂热工系统，目前，中国各种类型的锅炉有几十万台，由于设备分散、管理不善或技术原因，使大多数锅炉难以处于良好工况，增加了锅炉的燃料消耗，降低了效率。同时，锅炉工作过程中各项指标的调节难以建立数学模型，具有非线性、不稳定性、时滞等特点，所以如何改善对锅炉的控制，保证其正常工作，提高效率一直是人们关注的焦点。而汽包液位是锅炉安全、稳定运行的重要指标，保证液位在给定范围内，对于高蒸汽品质、减少设备损耗和运行损耗、确保整个网络安全运行具有要意义。 现代锅炉的特点之一就是蒸发量显著提高，汽包容积相对变小，水位变化速度很快，稍不注意就容易造成汽包满水或者烧成干锅，这都对汽包液位控制系统提出了更高的要求。汽包液位过高，会影响汽包内汽液分离效果，使汽包出口的饱和蒸汽带水增多，蒸汽带水会使汽轮机产生水冲击，引起轴封破损、叶片断裂等事故。同时会使饱和蒸汽中含盐量增高，降低过热蒸汽品质，增加在过热器管壁和汽轮机叶片上的结垢。水位过低，则可能破坏自然循环锅炉汽水循环系统中某些薄弱环节，以致局部水冷管壁被烧坏，严重时会造成爆炸事故。 目前，对汽包液位位控制大多采用常规PID控制方式，从控制方式来看，它们要么系统结构简单成本低，不能有效的控制锅炉汽包“虚假水位”现象，要么能够在一定程度上控制“虚假现象”，系统却过于复杂，成本投入过大。常用的蒸汽锅炉液位调节系统有三种基本结构：单冲量调节系统结构、双冲量调节系统结构、串级三冲量调节系统结
0引言 潮流是配电网络分析的基础，用于电网调度、运行分析、操作模拟和设计规划，同时也是电压优化和网络接线变化所要参考的内容。潮流计算通过数值仿真的方法把电力系统的详细运行情况呈现给工作人员，从而便于研究系统在给定条件下的稳态运行特点。随着市场经济的发展，经济利益是企业十分看重的，而线损却是现阶段阻碍企业提高效益的一大因素。及时、准确的潮流计算结果，可以给出配电网的潮流分布、理论线损及其在网络中的分布，从而为配电网的安全经济运行提供参考。从数学的角度来看，牛顿-拉夫逊法能有效进行非线性代数方程组的计算且具有二次收敛的特点，具有收敛快、精度高的特点，在输电网中得到广泛应用。随着现代计算机技术的发展，利用编程和相关软件，可以更好、更快地实现配电网功能，本文就是结合牛顿-拉夫逊法的基本原理，利用C++程序进行潮流计算，计算结果表明该方法具有良好的收敛性、可靠性及正确性。 1牛顿-拉夫逊法基本介绍 1.1潮流方程 对于N 个节点的电力网络（地作为参考节点不包括在内），如果网络结构和元件参数已知，则网络方程可表示为： =&& YV I （1-1） 式中，Y 为N*N 阶节点导纳矩阵；&V 为N*1维节点电压列向量；&I 为N*1维节点注入电流列向量。如果不计网络元件的非线性，也不考虑移相变压器，则Y 为对称矩阵。 电力系统计算中，给定的运行变量是节点注入功率，而不是节点注入电流，这两者之间有如下关系： ??=&&&EI S （1-2） 式中，&S 为节点的注入复功率，是N*1维列矢量；? &S 为&S 的共轭；
??i diag ??=???? &&E V 是由节点电压的共轭组成的N*N 阶对角线矩阵。 由（1-1）和（1-2），可得： ??=&&&S EYV 上式就是潮流方程的复数形式，是N 维的非线性复数代数方程组。将其展开，有： ?i i i ij j j i P jQ V Y V ∈-=∑&& j=1,2，….，N （1-3） 式中，j i ∈表示所有和i 相连的节点j ，包括j i =。 将节点电压用极坐标表示，即令i i i V V θ=∠&，代入式（1-3）中则有： ()i i i i ij ij j j j i P jQ V G jB V θθ∈-=∠-+∠∑ ()()cos sin i j ij ij ij ij j i V V G jB j θθ∈=+-∑ 故有： () ()cos sin sin cos i i j ij ij ij ij j i i i j ij ij ij ij j i P V V G B Q V V G B θθθθ∈∈?=+?? =-?? ∑∑ i=1,2,…,N （1-4） 式（1-4）是用极坐标表示的潮流方程。 而节点功率误差： (cos sin )θθ∈?=-+∑SP i i i j ij ij ij ij j i P P V V G B (1-5) (cos sin )θθ∈?=--∑SP i i i j ij ij ij ij j i Q Q V V G B (1-6) 式中：SP i P ，SP i Q 为节点i 给定的有功功率及无功功率。 1.2牛顿-拉夫逊法基本原理 1.2.1牛拉法的一般描述 牛拉法是把非线性方程式的求解过程变成反复对相应的线性方程式的求解过程，即非线性问题通过线性化逐步近似，这就是牛拉法的核心。下面以非线性方程式的求解过程来进行说明。 设电力网络的节点功率方程一般形式如下：
课程设计论文 基于MATLAB的电力系统潮流计算 学院:电气工程学院 专业:电气工程及自动化 班级:电自0710班 学号:0703110304 姓名: 马银莎
内容摘要 潮流计算是电力系统最基本最常用的计算。根据系统给定的运行条件，网络接线及元件参数，通过潮流计算可以确定各母线的电压（幅值和相角），各支路流过的功率，整个系统的功率损耗。潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节。因此，潮流计算在电力系统的规划计算，生产运行，调度管理及科学计算中都有着广泛的应用。 潮流计算在数学上是多元非线性方程组的求解问题，牛顿—拉夫逊Newton-Raphson法是数学上解非线性方程组的有效方法，有较好的收敛性。运用电子计算机计算一般要完成以下几个步骤：建立数学模型，确定解算方法，制订计算流程，编制计算程序。 关键词 牛顿-拉夫逊法（Newton-Raphson）变压器及非标准变比无功调节 高斯消去法潮流计算Mtlab
一 .电力系统潮流计算的概述 在电力系统的正常运行中，随着用电负荷的变化和系统运行方式的改变，网络中的损耗也将发生变化。要严格保证所有的用户在任何时刻都有额定的电压是不可能的，因此系统运行中个节点出现电压的偏移是不可避免的。为了保证电力系统的稳定运行，要进行潮流调节。 随着电力系统及在线应用的发展，计算机网络已经形成，为电力系统的潮流计算提供了物质基础。电力系统潮流计算是电力系统分析计算中最基本的内容，也是电力系统运行及设计中必不可少的工具。根据系统给定的运行条件、网络接线及元件参数，通过潮流计算可以确定各母线电压的幅值及相角、各元件中流过的功率、整个系统的功率损耗等。潮流计算是实现电力系统安全经济发供电的必要手段和重要工作环节，因此潮流计算在电力系统的规划设计、生产运行、调度管理及科学研究中都有着广泛的应用。它的发展主要围绕这样几个方面：计算方法的收敛性、可靠性；计算速度的快速性；对计算机存储容量的要求以及计算的方便、灵活等。 常规的电力系统潮流计算中一般具有三种类型的节点：PQ 、PV 及平衡节点。一个节点有四个变量，即注入有功功率、注入无功功率，电压大小及相角。常规的潮流计算一般给定其中的二个变量：PQ 节点（注入有功功率及无功功率），PV 节点（注入有功功率及电压的大小），平衡节点（电压的大小及相角）。 1、变量的分类： 负荷消耗的有功、无功功率——1L P 、1L Q 、2L P 、2L Q 电源发出的有功、无功功率——1G P 、1G Q 、2G P 、2G Q 母线或节点的电压大小和相位——1U 、2U 、1δ、2δ 在这十二个变量中，负荷消耗的有功和无功功率无法控制，因它们取决于用户，它们就称为不可控变量或是扰动变量。电源发出的有功无功功率是可以控制的自变量，因此它们就称为控制变量。母线或节点电压的大小和相位角——是受控制变量控制的因变量。其中, 1U 、2U 主要受1G Q 、2G Q 的控制, 1δ、2δ主要受 1G P 、2G P 的控制。这四个变量就是简单系统的状态变量。 为了保证系统的正常运行必须满足以下的约束条件： 对控制变量 max min max min ;Gi Gi Gi Gi Gi Gi Q Q Q P P P <<<< 对没有电源的节点则为 0;0==Gi Gi Q P 对状态变量i U 的约束条件则是 m a x m i n i i i U U U <<
潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 一、设计原理 1．牛顿-拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不
平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y。 （2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 二、设计内容 1.设计流程图