贵州大学电气工程学院
电力系统分析课程设计报告书
题目: 高压输电网潮流的计算机算法程序设计
专业：电气工程及其自动化
班级：电自093
学号：0908040179
学生姓名：文绵丽
指导教师：韩松
2012年7月4日
目录
第一章课程设计概述 (3)
1.1 设计目的 (3)
1.2 设计要求 (3)
1.3 设计题目 (3)
1.4 设计内容 (3)
1.5设计时间 (4)
第二章（针对该课程设计题目的）设计思路 (4)
2.1 潮流计算题目 (4)
2.2 对课题的分析及求解思路 (5)
第三章电力系统潮流计算概述 (6)
3.1电力系统叙述 (6)
3.2潮流计算简介 (6)
3.3潮流计算的意义及其发展 (7)
第四章导纳矩阵的原理及计算方法 (8)
4.1自导纳和互导纳的确定方法 (8)
4.2潮流计算的基本方程 (13)
4.3电力系统节点分类 (16)
4.4潮流计算的约束条件 (17)
第五章程序流程图及程序代码 (18)
5.1 潮流计算流程图 (18)
5.2 潮流计算程序代码 (19)
第六章运行结果分析 (27)
6.1 正常运行结果分析 (27)
6.1.1 MATLAB运行结果 (27)
6.1.2 PowerWorld运行结果 (33)
6.1.3 综合分析 (34)
6.2 非正常运行结果分析 (35)
6.2.1 若某发电站发电量减半时对系统断电影响 (35)
6.2.2 若某发电站发电量减为零时对系统断电影响 (36)
第七章总结 (37)
第一章课程设计概述
1.1 设计目的
1. 掌握电力系统潮流计算的基本原理和电力系统运行方式的变化；
2. 掌握并能熟练运用一门计算机语言（MATLAB语言或C语言或C++语言）；
3. 采用计算机语言对潮流计算进行计算机编程计算。
1.2 设计要求
1. 程序源代码；
2.选定算例的输入，输出文件；
3. 程序说明；
4. 选定算例的程序计算过程；
5. 选定算例的手算过程（至少迭代2次）（可选）。
1.3 设计题目
高压输电网潮流的计算机算法程序设计（PQ分解法、牛顿-拉夫逊法）
或中压配电网潮流的计算机算法程序设计（前推后代法、同伦延拓法等）或电力系统短路故障的计算机算法程序设计（要求不限）
1.4 设计内容
1.根据电力系统网络推导电力网络数学模型，写出节点导纳矩阵；
2.赋予各节点电压变量（直角坐标系形式）初值后，求解不平衡量；
3.形成雅可比矩阵；
4.求解修正量后，重新修改初值，从2开始重新循环计算；
5.求解的电压变量达到所要求的精度时，再计算各支路功率分布、功率损耗和平衡节
点功率;
6.上机编程调试；连调；
7.计算分析给定系统潮流分析并与手工计算结果作比较分析。
8.准备计算机演示答辩，书写该课程设计说明书（必须计算机打印）。1.5设计时间
2012年春季第17周至第18周
第二章（针对该课程设计题目的）设计思路2.1 潮流计算题目
○5○4○3
50MW25MW 10Mvar 180MW 100Mvar ○1○2
50MW 30Mvar
如图所示，变压器参数为：Z t1=0.006198+j0.22727, Z t2=0.0016529+j0.04752,变比分别为231/110kV,231/121kV;220kV线路参数为：
Z=0.01219+j0.06508;110kV线路：Z=0.13429+j0.20661; 23段，Z=0.13429+j0.20661；45段，Z=0.01219+j0.06508，所有阻抗都已经归算至220kV侧。
2.2 对课题的分析及求解思路
此电力系统是一个5节点，5支路的电力网络。综合比较牛顿拉夫逊法（直角坐标、极坐标）、PQ分解法等多种求解方法的特点，最后确定采用牛顿拉夫逊法（极坐标）。因为此方法所需解的方程组最少。
该系统中，母线5为平衡节点，保持定值电压ù=233.2+j0kV,4节点为PV 节点，注入有功功率为50MW，其余节点均为PQ节点，注入功率分别为?1=40+j30MVA，?2=-50-j30MVA,?3=-180-j100MVA.
首先计算节点数据矩阵如下，
% (bus#)(volt) (ang) (p) (q) (bus type)
bus=[
1 1.00 0.00 -0.10 0.60 1;
2 1.00 0.00 -0.50 -0.30 1;
3 1.00 0.00 -1.80 -1.00 1;
4 1.0
5 0.00 0.15 0.00 2;
5 1.0
6 0.00 0.00 0.00 3];
然后计算线路数据矩阵如下，
% b#1 b#2 ( R ) ( X ) (G) ( B ) ( K )
line = [
1 2 0.13429 0.20661 0.0 0.0 0;
2 3 0.13429 0.20661 0.0 0.0 0;
4 5 0.01219 0.06508 0.0 0.0 0;
4 3 0.0016
5 0.04752 0.0 0.0 0.954545;
5 1 0.00620 0.22727 0.0 0.0 1.05];
这就形成了计算潮流的原始数据，把它代入程序中就可以计算潮流分布。
第三章电力系统潮流计算概述
3.1电力系统叙述
电力工业发展初期，电能是直接在用户附近的发电站（或称发电厂）中生产的，各发电站孤立运行。随着工农业生产和城市的发展，电能的需要量迅速增加，而热能资源和水能资源丰富的地区又往往远离用电比较集中的城市和工矿区，为了解决这个矛盾，就需要在动力资源丰富的地区建立大型发电站，然后将电能远距离输送给电力用户。同时，为了提高供电的可靠性以及资源利用的综合经济性，又把许多分散的各种形式的发电站，通过送电线路和变电所联系起来。这种由发电机、升压和降压变电所，送电线路以及用电设备有机连接起来的整体，即称为电力系统。
现代电力系统提出了“灵活交流输电和新型直流输电”的概念。灵活交流输电技术是指运用固态电子器件与现代自动控制技术对交流电网的电压、相位角、阻抗、功率以及电路的通断进行实时闭环控制，从而提高高压输电线路的诉讼能力和电力系统的稳态水平。新型直流输电技术是指应用现电力电子技术的最新成果，改善和简化变流站的造价等。
运营方式管理中，潮流是确定电网运行方式的基本出发点：在规划领域，需要进行潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了电网在预想操作预想下的电网的潮流分布以及校验运行的可靠性。在电力系统调度运行的多个领域都涉及到电网潮流计算。潮流是确定电力网咯运行状态的基本因素，潮流问题是研究电力系统稳态问题的基础和前提。
3.2潮流计算简介
电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压。各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量的分析比较供电方案或运行方式的合理性。可靠性和经济性。此外，电力系统的潮流计算也是计算机系统动态稳定和静态稳定的基础，所以潮流计算是研究电力系统的一种和重要和基础的计算。
电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。
利用电子数字计算机进行潮流计算从50年代中期就已经开始了。在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的，对潮流计算的要求可以归纳为以下几点：
（1）计算方法的可靠性或收敛性；
（2）对计算机内存量的要求；
（3）计算速度；
（4）计算的方便性和灵活性。
3.3潮流计算的意义及其发展
电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算，即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有的电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和稳态分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采用牛顿—拉夫逊法。
运行方式管理中，潮流是确定电网运行方式的基本出发点；在规划领域，需要进行潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了多个在预想操作情况下电网的潮流分布以及校验运行可靠性。在电力系统调度运行的多个领域问题是研究电力系统稳态问题的基础和前提。
在用数字见算计算机解电力系统潮流问题的开始阶段，普遍采取以节点导纳矩阵为基础的逐次代入法。这个方法的原理比较简单，要求的数字计算机内存量比较差下，适应50年代电子计算机制造水平和当时电力系统理论水平，但它的收敛性较差，当系统规模变大时，迭代次数急剧上升，在计算中往往出现迭代不收敛的情况。这就迫使电力系统的计算人员转向以阻抗矩阵为基础的逐次代入法。阻抗法改善了系统潮流计算问题的收敛性，解决了导纳无法求解的一些系统的潮流计算，在60年代获得了广泛的应用，阻抗法德主要缺点是占用计算机内
存大，每次迭代的计算量大。当系统不断扩大时，这些缺点就更加突出，为了克服这些缺点，60年代中期发展了以阻抗矩阵为基础的分块阻抗法。这个方法把一个大系统分割为几个小的地区系统，在计算机内只需要存储各个地区系统的阻抗矩阵及它们之间联络的阻抗，这样不仅大幅度的节省了内存容量，同时也提高了计算速度。
克服阻抗法缺点是另一个途径是采用牛顿-拉夫逊法。这是数学中解决非线性方程式的典型方法，有较好的收敛性。在解决电力系统潮流计算问题时，是以导纳矩阵为基础的，因此，只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿法潮流程序的效率。自从60年代中期，牛顿法中利用了最佳顺序消去法以后，牛顿法在收敛性。内存要求。速度方面都超过了阻抗法，成为了60年代末期以后广泛采用的优秀方法。
第四章 导纳矩阵的原理及计算方法
4.1自导纳和互导纳的确定方法
电力网络的节点电压方程： B
B B I Y U =
B I 为节点注入电流列向量，注入电流有正有负，注入网络的电流为正，流出网络的电流为负。根据这一规定，电源节点的注入电流为正，负荷节点为负。既无电源又无负荷的联络节点为零，带有地方负荷的电源节点为二者代数之和。
B U 为节点电压列向量，由于节点电压是对称于参考节点而言的，因而需先选定参考节点。在电力系统中一般以地为参考节点。如整个网络无接地支路，则需要选定某一节点为参考。设网络中节点数为（不含参考节点），则B I ，B U 均为n*n 列向量。B Y 为n*n 阶节点导纳矩阵。
节电导纳矩阵的节点电压方程： B B B I Y U =
展开为：
：11121311
1212223222313233333
1
2
3n n n n n n nn n n Y Y Y Y I U Y Y Y Y I U Y Y Y Y I
U Y Y Y Y I U ????????????????????????=??????????
??????????????
B Y 是一个n*n 阶节点导纳矩阵，其阶数就等于网络中除参考节点外的节点数。 节点导纳矩阵的对角元素ii Y (i=1，2， n)成为自导纳。自导纳数ii Y 值上就等于在i 节点施加单位电压，其他节点全部接地时，经节点i 注入网络的电流，
因此，它可以定义为：
/(0,)ii i i j
Y I U U j i ==≠
节点i 的自导纳ii Y 数值上就等于与节点直接连接的所有支路导纳的总和。 节点导纳矩阵的非对角元素ij Y (j=1，2，…，n;i=1，2，…。，n;j=i)称互导纳，
由此可得互导纳ij Y 数值上就等于在节点i 施加单位电压，其他节点全部接地时，经节点j 注入网络的电流，因此可定义为：
/(0,)ji ji i j
Y I U U j i =≠≠ 节点j ，i 之间的互导纳ij Y 数值上就等于连接节点j ，i 支路到导纳的负值。显然，恒ij Y 等于ji Y 。互导纳的这些性质决定了节点导纳矩阵是一个对称稀疏矩阵。而且，由于每个节点所连接的支路数总有一个限度，随着网络中节点数的增
加非
零元素相对愈来愈少，节点导纳矩阵的稀疏度，即零元素数与总元素的比值就愈来愈高。
2.1.2节点导纳矩阵的性质及意义 节点导纳矩阵的性质：
（1）B Y 为对称矩阵，ij Y =ji Y 。如网络中含有源元件，如移相变压器，
则对称性不再成立。 （2）
B
Y 对无接地支路的节点，其所在行列的元素之和均为零，即
,,1
1
0,0n
n
i j
j i j i Y
Y ====∑∑。对于有接地支路的节点，其所在行列的元素之和等于
该点接地支路的导纳。利用这一性质，可以检验所形成节点导纳矩阵的正确性。 （3）
B Y 具有强对角性：对角元素的值不小于同一行或同一列中任一元素
（4）B Y
为稀疏矩阵，因节点i ，j 之间无支路直接相连时ij Y =0，这种情况在实际电力系统中非常普遍。矩阵的稀疏性用稀疏度表示，其定义为矩阵中的零元
素与全部元素之比，即 2
/S Z n =， 式中Z 为B Y 中的零元素。S 随节点数n 的
增加而增加：n=50，S 可达92%；n=100，S 可达90%；n=500，S 可达99%，充分利用节点导纳矩阵的稀疏性可节省计算机内存，加快计算速度，这种技巧称为稀疏技术。
节点导纳矩阵的意义：
B Y 是n*n 阶方阵，其对角元素 ii Y (i=1，2，----n)称为自导纳，非对角元素ij Y (i ，j=1，2， n ， i j ≠)称为互导纳。将节点电压方程B B B I Y U =展开为
1111211221222122n n n n n nn I Y Y Y U I Y Y Y U I Y Y Y Un ??????????????????=??????????????????
???? 可见， /(0,,1,2,,,)ii
i i j
Y I U U i j n i j ===≠ 表明，自导纳ii Y 在数值上等于仅在节点i 施加单位电压而其余节点电压均为零（即其余节点全部接地）时，经节点i 注入网络的电流。其显然等于与节点i 直接相连的所有
支路的导纳之和。同时可见/(0,,1,2,,)ij i j i
Y I U U i j n j i ===≠ 。表明，互导纳在数值上等于仅在节点j 施加单位电压而其余节点电压均为零时，经节点i 注入网络的电流，其显然等于(ij y -)即ij Y =ij y -。ij y 为支路的导纳，负号表示该电流流出网络。如节点ij 之间无支路直接相连，则该电流为0，从而ij Y =0。
注意字母几种不写法的不同意义：粗体黑字表示导纳矩阵，大写字母ij Y 代矩阵B Y 中的第i 行第j 列元素，即节点i 和节点j 之间的互导纳。小写字母i ，j 支路的导纳等于支路阻抗的倒数数，
1/ij ij y Z =。
根据定义直接求取节点导纳矩阵时，注意以下几点：
1)。节点导纳矩阵是方阵，其阶数就等于网络中除去参考节点外的节点数。参考节点一般取大地，编号为零。
2)。节点导纳矩阵是稀疏矩阵，其各行非零非对角元素就等于与该行相对应节点所连接的不接地支路数。
3)。节点导纳矩阵的对角元素就等于各该节点所连接导纳的总和。因此，与没有接地支路的节点对应的行或列中，对角元素为非对角元素之和的负值。
4)。节点导纳矩阵的非对角元素等于连接节点i ，j 支路导纳的负值。因此，一般情况下，节点导纳矩阵的对角元素往往大于非对角元素的负值。
5)。节点导纳矩阵一般是对称矩阵，这是网络的互易特性所决定的。从而，一般只要求求取这个矩阵的上三角或下三角部分。
2.1.3非标准变比变压器等值电路
变压器型等值电路更便于计算机反复计算,更适宜于复杂网络的潮流计算.双绕组变压器可用阻抗与一个理想变压器串联的电路表示.理想变压器只是一个参数,那就是变比1U /U2K =。现在变压器阻抗按实际变比归算到低压侧为例,推导出变压器型等值电路.
a 双绕组变压器原理图
b 变压器阻抗归算到低压侧等值模型
流入和流出理想变压器的功率相等
11
12/U I U I K =
1
2
/I I K =
式中, 1U /U2K =是理想变压器的变比,1U 和 2U 分别为变压器高,低绕组的实际电压.从图b 直接可得：
122T
U K
U I Z =+
从而可得:
12T 1T 2122
T T U U Y U Y U I Z Z =-=-K K K K
12T 12T 2T T U U Y U I Y U Z Z =-=-K K
式中T T Y 1/Z =,又因节点电流方程应具有如下形式:
1111122
I Y U +Y U = 2211222-I Y U +Y U =
将式（1-8）与（1-9）比较，得：
2
11T Y =Y /K 12T Y
=-Y /K 21T Y =-Y /K 22T Y =Y
因此可得各支路导纳为:
12
12T 2121T 101112T 2202221T Y =-Y Y /Y =-Y Y /1Y Y Y Y 1Y Y Y Y =K
??=K
??-K
?
=-=K ??
K -=-=?
K ?
由此可得用导纳表示的变压器型等值电路:
图 c
4.2潮流计算的基本方程
在潮流问题中，任何复杂的电力系统都可以归纳为以下元件（参数）组成。（1）发电机（注入电流或功率）
（2）负荷（注入负的电流或功率）
（3）输电线支路（电阻，电抗）
（4）变压器支路（电阻，电抗，变比）
（5）母线上的对地支路（阻抗和导纳）
（6）线路上的对地支路（一般为线路充电点容导纳）
集中了以上各类型的元件的简单网络如图
(a) 潮流计算用的电网结构图
(b) 潮流计算等值网络
采用导纳矩阵时，节点注入电流和节点电压构成以下线性方程组
I=YU
其中
1
2
n
I
I
I=
I
??
??
??
??
??
??
??
1
2
U
U
U=
Un
??
??
??
??
??
??
可展开如下形式 n
i ij j j 1
I Y U (i=1,2,n )==∑
由于实际电网中测量的节点注入量一般不是电流而是功率，因此必须将式中
的注入电流用节点注入功率来表示。 节点功率与节点电流之间的关系为 i S =i i i i P jQ U I -=
式中i
Gi LDi P P P =-，i Gi LDi Q Q Q =-
因此用导纳矩阵时，PQ 节点可以表示为i S /i i
i i i
P jQ I U U =-==
把这个关系代入式中 ，得
1(1,2,)n
i i
ij j j i
P jQ Y U i n U =-==∑
就是电力系统潮流计算的数学模型-----潮流方程。它具有如下特点：
（1）它是一组代数方程，因而表征的是电力系统的稳定运行特性。 （2）它是一组非线性方程，因而只能用迭代方法求其数值解。
（3）由于方程中的电压和导纳既可以表为直角坐标，又可表为极坐标，因而潮流方程有多种表达形式---极坐标形式，直角坐标形式和混合坐标形式。
a 。取 i
i i U U θ=∠ ，||ij ij ij Y y β=∠，得到潮流方程的极坐标形式：
1
n
i i i i ij j i j P jQ U Y U θθ=-=∠∠∑
b 。 取 i
i i U e jf =+ ， ij ij ij Y G jB =+，得到潮流方程的直角坐标形式：
1111()()()()n
n
i i ij j ij j i ij j ij j j j n n i i ij j ij j i ij j ij j j j P e G e B f f G f B e Q f G e B f e G f B e ====?=-++?
?
??=--+??
∑∑∑∑
c 。取， i
i i U U θ=∠ ij ij ij Y G jB =+，得到潮流方程的混合坐标形式：
11(cos sin )(sin cos )n
i i j ij ij ij ij j n i i j ij ij ij ij j P U U G B Q U U G B θθθθ==?=+??
??=-??
∑∑
不同坐标形式的潮流方程适用于不同的迭代解法。例如：利用牛顿---拉夫逊迭代法求解，以直角坐标和混合坐标形式的潮流方程为方便；而P-Q 解耦法是在混合坐标形式的基础上发展而成，故当然采用混合坐标形式。 （4）它是一组n 个复数方程，因而实数方程数为2n 个但方程中共含4n 个变量：P ，Q ，U 和θ，i=1，2， ，n ，故必须先指定2n 个变量才能求解。
4.3电力系统节点分类
用一般的电路理论求解网络方程，目的是给出电压源(或电流源)研究网络内的电流(或电压)分布，作为基础的方程式，一般用线性代数方程式表示。然而在电力系统中，给出发电机或负荷连接母线上电压或电流(都是向量)的情况是很少的，一般是给出发电机母线上发电机的有功功率(P)和母线电压的幅值(U)，给出负荷母线上负荷消耗的有功功率(P)和无功功率(Q)。主要目的是由这些已知量去求电力系统内的各种电气量。所以，根据电力系统中各节点性质的不同，很自然地把节点分成三类：
① PQ 节点
对这一类点，事先给定的是节点功率(P ，Q)，待求的未知量是节点电压向量(U ，θ)，所以叫PQ 节点。通常变电所母线都是PQ 节点，当某些发电机的输出功率P 。Q 给定时，也作为PQ 节点。PQ 节点上的发电机称之为PQ 机(或PQ 给定型发电机)。在潮流计算中，系统大部分节点属于PQ 节点。 ② PU 节点
这类节点给出的参数是该节点的有功功率P 及电压幅值U ，待求量为该节点的无功功率Q 及电压向量的相角θ。这类节点在运行中往往要有一定可调节的无功电源。用以维持给定的电压值。通常选择有一定无功功率储备的发电机母线或者变电所有无功补偿设备的母线做PU 节点处理。PU 节点上的发电机称为PU 机(或PU 给定型发电机) ③ 平衡节点
在潮流计算中，这类节点一般只设一个。对该节点，给定其电压值，并在计算中取该节点电压向量的方向作为参考轴，相当于给定该点电压向量的角度为零。也就是说，对平衡节点给定的运行参数是U 和θ，因此有城为U θ节点，而待求量是该节点的P 。Q ，整个系统的功率平衡由这一节点承担。
关于平衡节点的选择，一般选择系统中担任调频调压的某一发电厂(或发电机)
，有时也可能按其他原则选择，例如，为提高计算的收敛性。可以选择出线 数多或者靠近电网中心的发电厂母线作平衡节点。 以上三类节点4个运行参数P 。Q 。U 。θ中，已知量都是两个，待求量也是两个，
只是类型不同而已。
4.4潮流计算的约束条件
电力系统运行必须满足一定技术和经济上的要求。这些要求够成了潮流问
题中某些变量的约束条件，常用的约束条件如下： 1. 节点电压应满足 2.
min
max (1,2,)i i i U U U i n ≤≤=
从保证电能质量和供电安全的要求来看，电力系统的所有电气设备都必须运行在额定电压附近。PU 节点电压幅值必须按上述条件给定。因此，这一约束条件对PQ 节点而言。
3. 节点的有功功率和无功功率应满足
4.
min max min max Gi Gi Gi Gi Gi Gi P P P Q
Q Q ≤≤?
?≤≤?
PQ 节点的有功功率和无功功率，以及PU 节点的有功功率，在给定是就必须满足上述条件，因此，对平衡节点的P 和Q 以及PU 节点的Q 应按上述条件进行检验。
5. 节点之间电压的相位差应满足
max ||||||ij i j i j θθθθθ=-<-
为了保证系统运行的稳定性，要求某些输电线路两端的电压相位不超过一定的数值。这一约束的主要意义就在于此。
因此，潮流计算可以归结为求解一组非线性方程组，并使其解答满足一定的约束条件。常用的方法是迭代法和牛顿法，在计算过程中，或得出结果之后用约束条件进行检验。如果不能满足要求，则应修改某些变量的给定值，甚至修改系 统的运行方式，重新进行计算。
第五章程序流程图及程序代码5.1 潮流计算流程图
牛顿---拉夫逊法计算潮流的基本步骤有如下几步：
①形成节点导纳矩阵Y B;
②设各节点电压初值；
③由初值可以求出修正方程中的不平衡量；
④求修正方程式的系数矩阵--雅可比矩阵的各个元素；
⑤解修正方程式，求各节点电压的变化量；
⑥计算各节点电压的新值，即修正后的值；
⑦运用各节点的新值自第三步开始下一次迭代；
⑧计算平衡节点功率和线路功率。
流程图如下：
输入原始数据
形成节点导纳矩阵设电压初值e(0)、f(0)
设迭代次数k=0
计算误差向量
?P(k)、?Q(k)、?U2(k)
收敛否
求雅可比矩阵元素
解修正方程，求解?e(k)、?f(k)
修正节点电压e(k+1)=e(k)-?e(k) f(k+1)=f(k)-?f(k)
K=K+1
K>Kmax
不收敛停机
求PU节点无功功率，
求平衡节点功率
求支路功率分布和损耗
停机
是
否
5.2 潮流计算程序代码
主程序：
clc;
clear all
format long;
global nSW nPQ nPV; global nb; global nl; global bus;
global line;
global Y;
global nodenum;
global lPQ;
global myf;
openfile;
change;
ybus;
NR;
PQ;
flow;
ret;
output;
打开bus line的数据文件的子程序：
function openfile
global nb nl;
global line;
global bus;
global myf;
[dfile,pathname]=uigetfile('*.m','Select Data File'); if pathname == 0
error(' you must select a valid data file')
else
lfile =length(dfile);
% strip off .m
eval(dfile(1:lfile-2));
end
[nl,ml]=size(line);
[nb,mb]=size(bus);
outfile='';
for I=1:length(dfile)
if dfile(I)=='.'
break;
end
outfile=strcat(outfile,dfile(I));
end
outfile=strcat(outfile,'.txt');
myf=fopen(outfile,'w');
节点重新编号：
function change
global nb;
global nl;
global nPQ;
基础工业工程课程设计 ————模特排时法
一、取题背景 (2) 台钻厂简介 (2) 产品介绍 (2) 产品零件来源 (2) 课题选择及课程设计目的 (3) 课程设计内容及要求 (3) 二、应用的基本原理 (3) 模特法的基本原理 (3) 模特法的21种基本动作 (4) 三、现行方案 (5) 现行方案的生产现场布置分布图 (5) 运用MOD法对台钻的装配流程进行记录、分析 (6) 四、问题分析 (20) 五、改进方案设计 (21) 改善后现场布置图 (22) 改善后的各工位动作因素表 (23) 六、计算正常时间及标准时间 (36) 计算正常时间 (36) 确定宽放率 (36) 计算标准时间 (36) 七、对比分析及改善评价 (37) 八、设计总结与心得体会 (38)
一、取题背景 本次课程设计收集的数据和资料取自于福州台钻厂。 （一）台钻厂简介 福州台钻厂创建于1958年，是福建工程学院的校办工厂，其前身是国家机械电子工业部台式钻床定点生产企业，国家渔业机械仪器行业定点企业，经过多年的发展以及福建工程学院雄厚的技术和科研支持，如今已成为具有上千万资产的小型企业，其出产的台钻床由于品质卓越、性能稳定连续多次被评为部优产品，增氧机获得了农业部“农业机械推广许可证书”。正因为福建工程学院科研与技术的支持和台钻厂自身坚持质量第一的宗旨，福州台钻厂获得了多项国家专利，并被评为一级信用企业、福州工业利税“百强企业”、福建省小型工业企业“百佳明星”企业和出口创汇先进单位等。企业的主导产品销往全国各地，并远销欧美等地。 多年来，台钻厂以经营为主体，同时又承担学校学生的认知实习、课程设计以及毕业设计等教学实习任务。 （二）产品介绍 企业主导产品有：“武夷山”牌台式钻床系列和水产养殖设备。 1.“武夷山”牌台式钻床，主要包括Z4112、4012、Z4116、Z4016、ZM406、Z4416、Z4416A、Z4020、Z4025、Z4032、S406、S4112等产品。其中Z4112台式钻床于1980年、1986年和1991年被国家机械部评为部优产品。 2．“武夷山”牌水产养殖机械，包括YC0.37、YC0.75、YC1.5、YC4.5等水车式增氧机和YL1.5叶轮式增氧机。企业建有专用的增氧机试验场所，定期开展整机测试工作，增氧能力和动力效率等指标均高于国家行业标准。 （三）产品零件来源 1. 自制件 该企业主要从事五大零件的自制，包括台钻的主轴、主轴套筒、主轴箱、底座、立柱（Z4112、Z4116系列还包括中间工作台的制造）。 2. 外购件 该企业外购零件主要有电机和轴承等标准件。 3. 外协件 该企业选择固定的三家认证厂家购买外协件。
网页制作课程设计报告 学院： 专业班级： 姓名： 学号： 成绩： 阅卷教师：
目录 1.设计目的 (1) 2.设计思想 (1) 2.1网站整体结构规划思想 (1) 2.2 主页设计思想 (1) 2.3子页的设计思想 (1) 3网页详细设计分析 (1) 4结论 (2)
1.设计目的 阐述该个人网站的设计意图和创意，简单介绍自己的个人网站。 2.设计思想 阐述网站的整体设计思想，包括： 2.1网站整体结构规划思想 要求阐述网站整体结构的选择、设计的思想，绘制网站结构草图。 2.2 主页设计思想 要求对主页的布局思路进行阐述和分析。 2.3子页的设计思想 要求对子页的设计以及网页对象的选取思路进行阐述和分析。 3网页详细设计分析 要求选取一张网页，对网页的设计实现过程进行阐述和分析，详细说明制作该网页的步骤，所使用的网页对象以及该网页对象的操作方法。
4结论 对整个设计报告做归纳性总结，并分析设计过程中的困难及如何解决的，最后提出展望。 一、设计目的 本课程的设计目的是通过实践使同学们经历Dreamweaver cs3开发的全过程和受到一次综合训练，以便能较全面地理解、掌握和综合运用所学的知识。结合具体的开发案例，理解并初步掌握运用Dreamweaver cs3可视化开发工具进行网页开发的方法；了解网页设计制作过程。通过设计达到掌握网页设计、制作的技巧。了解和熟悉网页设计的基础知识和实现技巧。根据题目的要求，给出网页设计方案，可以按要求，利用合适图文素材设计制作符合要求的网页设计作品。熟练掌握Photoshop cs3、Dreamweaver cs3等软件的的操作和应用。增强动手实践能力，进一步加强自身综合素
电压源与电流源的等效变换 一、实验目的 1、加深理解电压源、电流源的概念。 2、掌握电源外特性的测试方法。 二、原理及说明 1、电压源是有源元件，可分为理想电压源与实际电压源。理想电压源在一定的电流 范围内，具有很小的电阻，它的输出电压不因负载而改变。而实际电压源的端电压随着电流变化而变化，即它具有一定的内阻值。理想电压源与实际电压源以及它们的伏安特性如图4-1所示(参阅实验一内容)。 2、电流源也分为理想电流源和实际电流源。 理想电流源的电流是恒定的，不因外电路不同而改变。实际电流源的电流与所联接的电路有关。当其端电压增高时，通过外电路的电流要降低，端压越低通过外电路的电 并联来表示。图4-2为两种电流越大。实际电流源可以用一个理想电流源和一个内阻R S 流源的伏安特性。
3、电源的等效变换 一个实际电源，尤其外部特性来讲，可以看成为一个电压源，也可看成为一个电流源。两者是等效的，其中I S=U S/R S或 U S=I S R S 图4-3为等效变换电路，由式中可以看出它可以很方便地把一个参数为U s 和R s 的 电压源变换为一个参数为I s 和R S 的等效电流源。同时可知理想电压源与理想电流源两者 之间不存在等效变换的条件。 三、仪器设备 电工实验装置： DG011、 DG053 、 DY04 、 DYO31 四、实验内容 1、理想电流源的伏安特性 1)按图4-4(a)接线，毫安表接线使用电流插孔，R L 使用1KΩ电位器。 2)调节恒流源输出，使I S 为10mA。， 3)按表4-1调整R L 值,观察并记录电流表、电压表读数变化。将测试结果填入表4-1中。 2、实际电流源的伏安特性 按照图4-4(b)接线，按表4-1调整R L 值，将测试的结果填入表4-1中。
山东建筑大学本科学生课程设计任务书 课程设计题目*****工程项目成本指标分析 学院管理工程学院专业房地产经营与管理年级2011级 已知参数和设计要求： 本课程设计所要达到的目的：基于收集到的工程项目概预算书，分析某项目的成本指标，并在此基础上提出管理建议。使学生进一步深入掌握成本管理的理论知识，并将所学专业知识应用于工程项目管理实践，培养学生的实际工作能力。 主要任务： 从房地产公司、施工单位、咨询公司或从网上下载完成的招投标文书，每个同学的工程项目不同，需要从建设单位角度和施工单位角度制定该项目的成本控制方案。 设计成果： 1、网络上下载的估算报告或概预算书原文（电子版和纸板）； 2、成本指标分析课程设计文本（内容详见附表）。 主要要求： 1、课程设计要求每人独立完成，设计的结果由设计人员每人单独提交。 2、课程设计方案的深度：各种管理建议要落实到可以操作的层面上，比如某 一具体的措施能够明确到某个人或部门能够执行；至少要完成附表中的相关指标分析。 3、严格按照学校课程设计要求进行。 学生应完成的工作： ●成本规划与成本控制的理论； ●设计和施工阶段各成本控制要点； ●结合实际工程项目的成本构成分析成本管理的具体措施； ●结合实际工程项目分析成本各子项控制与项目总成本之间的关系。
目前资料收集情况（含指定参考资料）： 1、《建设工程成本计划与控制》（任宏主编，中国高等教育出版社）； 2、《工程项目管理》（丛培经主编，中国建筑工业出版社）； 3、《建设工程成本规划与控制》； 4、课程设计任务书； 5、工程管理相关网站； 6、其它相关资料。 课程设计的工作计划： 1、明确设计内容和要求，0.5天； 2、设计任务计划的制定，0.5天； 3、完成项目成本控制方案并提交成果，5天。 任务下达日期 2013 年 12 月 19 日完成日期 2013 年 12 月 27 日
Web应用开发技术 实验报告 专业：计算机科学与技术 班级： 学号： 姓名：
一、设计题目 个人网站 二、目的 １、本次设计是学生在学完ASP动态网站开发课程后的一次实践性很强的课程设计，是对ASP进行动态网站开发所学知识的综合运用。 2、掌握使用ASP技术进行网站开发设计。 3、通过本次实习，使学生加深所学知识内容的理解，并能积极地调动学生的学习兴趣，结合实际应用操作环境，真正做到理论与实际相结合。 三、功能需求描述 此网站可以对主人留言，来发表自己的心情，也可以把自己的联系方式写入其中，达到和睦相处、心灵的驿站的目的等。 四、总体设计
五、详细设计 （一）、我的主页 此页面为网站的主页，通过发布新心情，点击通讯录可以查看通讯录好友信息，点击留言板可以查看好友留言。 主要代码：