用高斯赛德尔法求节点电压

1、武汉理工大学电力系统分析课程设计说明书目录1绪论 11.1 设计背景11.2 设计目标12设计原理.32.1 高斯-赛德尔迭代法32.2 Matlab简介42.2.1 Matlab概述42.2.2 matlab GUI 简介52.2.3 GUI 设计模板及设计窗口52.2.4 GUI 设计的基本操作62.2.5MATPOWER62.3 节点导纳矩阵92.3.1 节点导纳矩阵的性质：92.3.2节点导纳矩阵的意义102.3.3节点导纳矩阵的特点112.4 潮流计算122.4.1概述122.4.2 潮流计算的节点类型122.4.3 潮流计算基本方程132.4.4潮流计算的约束条件153计算过程及步

2、骤.163.1计算要求163.2详细计算173.2.1第一次迭代173.2.2第二次迭代183.2.3第三次迭代193.2.4第四次迭代194 MATPOWER中建模解决问题214.1 MATPOWER的技术规则214.1.1数据文件格式214.1.2控制选项224.2 MATPOWER实际建模234.2.1 建模中的实际问题234.2.2 建模结果显示24总结29致谢31参考文献32附录33用高斯赛德尔法求节点电压1绪论1.1 设计背景电力是衡量一个国家经济发展的主要指标，也是反映人民生活水平的重要标志，它已成为现代工农业生产、交通运输以及城乡生活等许多方面不可或缺的能源和动力。电力系统是由

3、发电、输电、变电、配电和用电等环节组成的电能生产与消费系统。它的功能是将自然界的一次能源通过发电动力装置转化成电能，再经输电、变电和配电将电能供应到各用户。为实现这一功能，电力系统在各个环节和不同层次还具有相应的信息与控制系统，对电能的生产过程进行测量、调节、控制、保护、通信和调度，以保证用户获得安全、经济、优质的电能。电力系统的出现，使电能得到广泛应用，推动了社会生产各个领域的变化，开创了电力时代，出现了近代史上的第二次技术革命。20世纪以来，电力系统的发展使动力资源得到更充分的开发，工业布局也更为合理，使电能的应用不仅深刻地影响着社会物质生产的各个侧面，也越来越广地渗透到人类日常生活的各个

4、层面。电力系统的发展程度和技术水准已成为各国经济发展水平的标志之一。潮流计算是在给定电力系统网络结构、参数和决定系统运行状态的边界条件的情况下确定系统稳态运行状态的一种基本方法，是电力系统规划和运营中不可缺少的一个重要组成部分。可以说，它是电力系统分析中最基本、最重要的计算，是系统安全、经济分析和实时控制与调度的基础。是电力系统研究人员长期研究的一个课题。MATLAB自1980年问世以来，它的强大的矩阵处理功能给电力系统的分析、计算带来许多方便。在处理潮流计算时，其计算机软件的速度已无法满足大电网模拟和实时控制的仿真要求，而高效的潮流问题相关软件的研究已成为大规模电力系统仿真计算的关键。随着计

5、算机技术的不断发展和成熟，对MATLAB潮流计算的研究为快速、详细地解决大电网的计算问题开辟了新思路。1.2 设计目标 在已知的节点导纳矩阵下，告诉部分节点类型及数据，要通过高斯-赛德尔法求出未知的节点电压，通过分析，找到要掌握的问题有以下几点:高斯-赛德尔法到底是什么，通过手算还是软件。根据题目要求掌握高斯塞德尔法，能够求解出题目的节点电压，熟练掌握所学知识，并能适当拓展，学会使用各种工具软件，提高自身能力，开拓视野。2设计原理2.1 高斯-赛德尔迭代法潮流计算高斯赛德尔迭代法(Gauss一Seidel method)是求解电力系统潮流的方法。潮流计算高斯赛德尔迭代 法又分导纳矩阵

6、迭代法和阻抗矩阵迭代法两种。前者 是以节点导纳矩阵为基础建立的赛德尔迭代格式;后者是以节点阻扰矩阵为基础建立的赛德尔迭代格式。 高斯赛德尔迭代法这是数学上求解线性或非 线性方程组的一种常用的迭代方法。要掌握著名德国数学家高斯和数学家赛德尔研究出的这种非线性代数方程组的迭代解法，还需要学习德国数学家雅可比研究出的雅可比迭代法，其推导过程如下；设线性方程组 的系数矩阵A可逆且主对角线元素均不为零，令并将A分解成 从而（1）可写成 令 其中 以为迭代矩阵的迭代法（公式） 称为雅可比迭代法，用向量的分量来表示， 其中为初始向量.由雅可比迭代公式可知，在迭代的每一步计算过程中是用的全部分量来计

7、算的所有分量，显然在计算第i个分量时，已经计算出的最新分量, 没有被利用，从直观的角度看，最新的分量可能比旧的分量要好些，因此，对这些最新计算出来的第K+1次近似的分量加以利用，就得到所谓解方程组的高斯赛德尔迭代法。具体如下:把矩阵A分解成 其中，分别为A 的主对角元除外的下三角和上三角部分，于是，方程组便可写成 即 其中 ， 以为迭代矩阵构成的迭代法（公式） 称为高斯-赛德尔迭代法，用向量表示的形式为 收敛判据：复数模型：实数模型：，2.2 Matlab简介2.2.1 Matlab概述MATLAB (Matrix Laboratory)为美国Mathworks公司1983年首次推出的一套高性

8、能的数值分析和计算软件，其功能不断扩充，版本不断升级。 MATLAB将矩阵运算、数值分析、图形处理、编程技术结合在一起，为用户提供了一个强有力的科学及工程问题的分析计算和程序设计工具，它还提供了专业水平的符号计算、文字处理、可视化建模仿真和实时控制等功能，是具有全部语言功能和特征的新一代软件开发平台。MATLAB具有编程效率高、用户使用方便、扩充能力强、语句简单，内涵丰富、高效方便的矩阵和数组运算、方便的绘图功能等特点，给用户带来了极大的方便。MATLAB 已发展成为适合众多学科，多种工作平台、功能强大的大型软件。在欧美等国家的高校，MATLAB已成为线性代数、自动控制理论、数理统计、数字信号

9、处理、时间序列分析、动态系统仿真等高级课程的基本教学工具。成为攻读学位的本科、硕士、博士生必须掌握的基本技能。在设计研究单位和工业开发部门，MATLAB被广泛的应用于研究和解决各种具体问题。在中国，MATLAB也已日益受到重视，短时间内就将盛行起来，因为无论哪个学科或工程领域都可以从MATLAB中找到合适的功能。2.2.2 matlab GUI 简介 图形用户界面（GUI）是用户与计算机程序之间的交互方式，是用户与计算机进行信息交流的方式。计算机在屏幕显示图形和文本，若有扬声器还可产生 声音。用户通过输入设备，如：键盘、鼠标、跟踪球、绘制板或麦克风，与计算机通讯。用户界面设定了如何观看和如何感

10、知计算机、操作系统或应用程序。通常， 多是根据悦目的结构和用户界面功能的有效性来选择计算机或程序。图形用户界面或GUI是包含图形对象，如：窗口、图标、菜单和文本的用户界面。以某种方式 选择或激活这些对象，通常引起动作或发生变化。最常见的激活方法是用鼠标或其它点击设备去控制屏幕上的鼠标指针的运动。按下鼠标按钮，标志着对象的选择或 其它动作。Matlab作为强大的数学计算软件，同样也提供了图像用户界面设计的功能。在matlab中，基本的图形用户界面对象包含3类：用户控件对象（uicontrol）、下拉式菜单对象（uimenu）、和快捷菜单对象（uicontexmenu）。根据这些对象可以设计出界面

11、友好、操作方便的图形用户界面。2.2.3 GUI 设计模板及设计窗口Matlab为GUI设计准本了四个模板，分别是Blank GUI(默认)、GUI with Uicontronl（带控件对象的GUI）、GUI with Axes and Menu（带坐标轴与菜单的GUI）、Modal Question Dialog（带模式问话对话框的GUI模板），GUI 设计模板如图3-1所示。当用户选择不同模板时，在GUI设计模板界面的右边就会显示与该模板对应的GUI图形。图2-1 GUI设计模板选择设计模板后就进如GUI 设计窗口，GUI设计窗口由菜单栏、工具栏、控件工具栏以及图形对象设计区组成。在GU

12、I设计窗口的工具栏上有位置调整器、菜单编辑器、tab顺序编辑器、属性查看器等可视化设计工具。控件工具栏包括Push Button、Check Box、Edit Box、Popup Menu、Axes、table等控件对象，他们是构成GUI 的基本元素。2.2.4 GUI 设计的基本操作为了添加对象控件，可以从GUI 设计窗口的控件工具栏中选择一个对象，然后以拖曳方式在对象设计区建立该对象，其对象创建方式方便、简单。在GUI 设计窗口创建对象后，通过双击该对象，就会显示该对象的属性查看器，通过它可以设计该对象的属性值。在选中对象的前提下，单击鼠标右键，会弹出一个快捷菜单，可以从中某个子菜单进行相

13、应的操作。在对象设计区右击鼠标，会显示与图形窗口有关的快捷菜单。2.2.5MATPOWER潮流计算经历了一个由手工利用交直流计算到应用数字电子计算机的发展过程，目前常见的潮流计算仿真软件主要有电科院的PASAP，美国的BPA，PSS/E等，但源代码都不公开，所以最佳方案还是基于MATLAB语言编写的电力潮流计算的软件MATPOWER。MATPOWER是一个用MATLAB的M文件编写的软件包，最新版本是MATPOWER4.1，浙江大学的甘德强教授在康奈尔工作期间是主要研发人员。 首先在MATPOWER主页下载相关的压缩文件，解压文件，将解压后文件放到MATLAB的搜索路径下。实际操作中注意完成以

14、上步骤后要在MATLAB中的SET PATH中将MATPOWER添加到MATLAB路径中，这步必不可少开始输入数据，定义数组给定PQ节点电压初值给定PV节点电压实部（或虚部）置迭代计数b=0计算PQ节点电压实部和虚部置迭代计数b=0先计算PV节点无功功率再用其计算PV节点电压实部和虚部计算PQ节点电压实部和虚部先计算PV节点无功功率再用其计算PV节点电压实部和虚部计算平衡节点的有功和无功求=+判断所有|是否<0.000001b=b+1N结果输出Y图2-2高斯-塞德尔法潮流计算的流程框图2.3 节点导纳矩阵2.3.1 节点导纳矩阵的性质：（1）为对称矩阵，=。如网络中含有源元件，如移相变压

15、器，则对称性不再成立。（2）对无接地支路的节点，其所在行列的元素之和均为零，即 。对于有接地支路的节点，其所在行列的元素之和等于该点接地支路的导纳。利用这一性质，可以检验所形成节点导纳矩阵的正确性。（3）具有强对角性：对角元素的值不小于同一行或同一列中任一元素。（4）为稀疏矩阵，因节点i ，j 之间无支路直接相连时=0，这种情况在实际电力系统中非常普遍。矩阵的稀疏性用稀疏度表示，其定义为矩阵中的零元素与全部元素之比，即 ， 式中Z 为中的零元素。S 随节点数n 的增加而增加：n=50，S可达92%；n=100，S 可达90%；n=500，S可达99%，充分利用节点导纳矩阵的稀疏性可节省计算机内

16、存，加快计算速度，这种技巧称为稀疏技术。节点导纳矩阵以导纳形式描述电力网络节点注入电流和节点电压关系的矩阵。它给出了电力网络连接关系和元件特性的全部信息。根据基尔霍夫电流定律可写出电力网络中的n个节点方程式： 可用矩阵的形式表示I=YV。其中， 对角元素为节点 i的自导纳，非对角线为节点 i与节点j之间的互导纳。节点导纳矩阵反映了网络的参数及接线情况，因此，节点导纳矩阵可以看成是对电力网络电气特性的一种数学抽象。由导纳矩阵所联系的节点方程式是电力网络广泛应用的一种数学模型。2.3.2节点导纳矩阵的意义是n*n阶方阵，其对角元素 (i=1，2，-n)称为自导纳，非对角元素(i，j=1，2，n，

17、)称为互导纳。将节点电压方程展开为： （2-11）可见 (2-12)表明，自导纳在数值上等于仅在节点i施加单位电压而其余节点电压均为零（即其余节点全部接地）时，经节点i注入网络的电流。其显然等于与节点i直接相连的所有支路的导纳之和。同时可见。表明，互导纳在数值上等于仅在节点j施加单位电压而其余节点电压均为零时，经节点i注入网络的电流，其显然等于()即=。为支路的导纳，负号表示该电流流出网络。如节点ij之间无支路直接相连，则该电流为0，从而=0。注意字母几种不写法的不同意义：粗体黑字表示导纳矩阵，大写字母代矩阵中的第i行第j列元素，即节点i和节点j之间的互导纳。小写字母i，j支路的导纳等于支路阻

18、抗的倒数数，。根据定义直接求取节点导纳矩阵时，注意以下几点：1) 节点导纳矩阵是方阵，其阶数就等于网络中除去参考节点外的节点数。参考节点一般取大地，编号为零。2) 节点导纳矩阵是稀疏矩阵，其各行非零非对角元素就等于与该行相对应节点所连接的不接地支路数。3) 节点导纳矩阵的对角元素就等于各该节点所连接导纳的总和。因此，与没有接地支路的节点对应的行或列中，对角元素为非对角元素之和的负值。4) 节点导纳矩阵的非对角元素等于连接节点i，j支路导纳的负值。因此，一般情况下，节点导纳矩阵的对角元素往往大于非对角元素的负值。5) 节点导纳矩阵一般是对称矩阵，这是网络的互易特性所决定的。从而，一般只要求求取这

19、个矩阵的上三角或下三角部分。2.3.3节点导纳矩阵的特点 (1）导纳矩阵是稀疏矩阵。它的对角线元素一般不为零，但在非对角线元素中则存在不少零元素。在电力系统的接线图中，一般每个节点与平均不超过34个其他节点有直接的支路连接。因此，在导纳矩阵的非对角线元素中每行仅有34个非零元素，其余的都是零元素，而且网络的规模越大，这种现象越显著。 (2) 导纳矩阵为对称矩阵。由网络的互易特性易知。导纳矩阵的对称性和稀疏性对于应用计算机求解电力系统问题有很大的影响。如果能充分地利用这两个特点，如在程序设计中储存导纳矩阵的对角元素和上三角元素（或下三角元素），排除零元素的储存和运算，就可以大大地节省储存单元和提

20、高计算速度。2.4 潮流计算2.4.1概述 电力系统潮流计算的结果是计算出各节点电压、相位、支路功率、网络损耗等。通过计算可以得到节点电压，可以评估电压是否满足电能质量指标的要求；节点电压相位是评估系统稳定性的重要参数。潮流计算结果是各种控制技术措施、调度方案、规划方案评估的依据。电力系统潮流计算属于稳态分析范畴，不涉及系统元件的动态特性和过渡过程。因此其数学模型不包含微分方程，是一组高阶非线性方程。非线性代数方程组的解法离不开迭代，因此，潮流计算方法首先要求它是能可靠的收敛，并给出正确答案。在用数字计算机求解电力系统潮流问题的开始阶段，人们普遍采用以节点导纳矩阵为基础的高斯-赛德尔迭代法（一

21、下简称导纳法）。这个方法的原理比较简单，要求的数字计算机的内存量也比较小，适应当时的计算机制作水平和电力系统理论水平，于是电力系统计算人员转向以阻抗矩阵为主的逐次代入法。2.4.2 潮流计算的节点类型用一般的电路理论求解网络方程，目的是给出电压源(或电流源)研究网络内的电流(或电压)分布，作为基础的方程式，一般用线性代数方程式表示。然而在电力系统中，给出发电机或负荷连接母线上电压或电流(都是向量)的情况是很少的，一般是给出发电机母线上发电机的有功功率(P)和母线电压的幅值(U)，给出负荷母线上负荷消耗的有功功率(P)和无功功率(Q)。主要目的是由这些已知量去求电力系统内的各种电气量。所以，根据

22、电力系统中各节点性质的不同，很自然地把节点分成三类：PQ节点对这一类点，事先给定的是节点功率(P，Q)，待求的未知量是节点电压向量(U，)，所以叫PQ节点。通常变电所母线都是PQ节点，当某些发电机的输出功率P。Q给定时，也作为PQ节点。PQ节点上的发电机称之为PQ机(或PQ给定型发电机)。在潮流计算中，系统大部分节点属于PQ节点。PU节点这类节点给出的参数是该节点的有功功率P及电压幅值U，待求量为该节点的无功功率Q及电压向量的相角。这类节点在运行中往往要有一定可调节的无功电源。用以维持给定的电压值。通常选择有一定无功功率储备的发电机母线或者变电所有无功补偿设备的母线做PU节点处理。PU节点上的

23、发电机称为PU机(或PU给定型发电机)平衡节点在潮流计算中，这类节点一般只设一个。对该节点，给定其电压值，并在计算中取该节点电压向量的方向作为参考轴，相当于给定该点电压向量的角度为零。也就是说，对平衡节点给定的运行参数是U和，因此有城为U节点，而待求量是该节点的P。Q，整个系统的功率平衡由这一节点承担。关于平衡节点的选择，一般选择系统中担任调频调压的某一发电厂(或发电机)，有时也可能按其他原则选择，例如，为提高计算的收敛性。可以选择出线数多或者靠近电网中心的发电厂母线作平衡节点。以上三类节点4个运行参数P、Q、U、中，已知量都是两个，待求量也是两个，只是类型不同而已。2.4.3 潮流计算基本方

24、程电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。采用导纳矩阵时，节点注入电流和节点电压构成如式(2-7)所示线性方程组可展开如下形式： (2-18) 由于实际电网中测量的节点注入量一般不是电流而是功率，因此必须将式中的注入电流用节点注入功率来表示。节点功率与节点电流

25、之间的关系为： (2-19)式中，因此用导纳矩阵时，PQ节点可以表示为把这个关系代入式中 ，得 （2-20）式（2-20）就是电力系统潮流计算的数学模型-潮流方程。它具有如下特点：1：它是一组代数方程，因而表征的是电力系统的稳定运行特性。2：它是一组非线性方程，因而只能用迭代方法求其数值解。3：由于方程中的电压和导纳既可以表为直角坐标，又可表为极坐标，因而潮流方程有多种表达形式-极坐标形式，直角坐标形式和混合坐标形式。（1）取 ，得到潮流方程的极坐标形式： (2-21)(2) 取 ， ，得到潮流方程的直角坐标形式： (2-22)(3) 取 ，得到潮流方程的混合坐标形式： (2-23)不同坐标形

26、式的潮流方程适用于不同的迭代解法。例如：利用牛顿-拉夫逊迭代法求解，以直角坐标和混合坐标形式的潮流方程为方便；而P-Q解耦法是在混合坐标形式的基础上发展而成，故当然采用混合坐标形式。4: 它是一组n个复数方程，因而实数方程数为2n个但方程中共含4n个变量：P，Q，U和，i=1，2，n，故必须先指定2n个变量才能求解。2.4.4潮流计算的约束条件电力系统运行必须满足一定的技术和经济上的要求。这些要求构成了潮流问题中某些变量的约束条件，常用的约束条件如下：节点电压应满足小于节点最大额定电压并大于最小额定电压，即： (2-24)从保证电能质量和供电安全的要求来看，电力系统的所有电气设备都必须运行在额

27、定电压附近。PV节点电压幅值必须按上述条件给定。因此，这一约束条件对PQ节点而言。节点的有功功率和无功功率应满足小于节点最大额定功率并大于最小额定功率，即： (2-25) PQ节点的有功功率和无功功率，以及PV节点的有功功率，在给定时就必须满足上述条件，因此，对平衡节点的P和Q以及PV节点的Q应按上述条件进行检验。节点之间电压的相位差应满足小于最小额定相角差，即： (2-26)为了保证系统运行的稳定性，要求某些输电线路两端的电压相位不超过一定的数值。这一约束的主要意义就在于此。因此，潮流计算可以归结为求解一组非线性方程组，并使其解答满足一定的约束条件。常用的方法是迭代法和牛顿法，在计算过程中，

28、或得出结果之后用约束条件进行检验。如果不能满足要求，则应修改某些变量的给定值，甚至修改系统的运行方式，重新进行计算。3计算过程及步骤3.1计算要求如图所示的简单电力系统，节点导纳矩阵为：3-j9 -2+j6 -1+j3 0-2+j6 3.666-j11 -0.666+j12 -1+j3 -1+j3 -0.666+j12 3.666-j11 -2+j60 -1+3j -2+j6 3-j9节点2、3的注入功率已知，节点1为平衡节点，节点4为PV节点，用高斯-塞德尔法求节点2、3、4的电压，只迭代一次，取节点2、3的电压初值均为10。(图中各值均为标么值)。图3.1 电力系统模型图设 ， 由高斯-赛

29、德尔迭代式 得： 另外，我们可以通过多次迭代求得，其详细迭代过程如下：3.2详细计算3.2.1第一次迭代 由题目所给条件得，根据式3-2可计算第一次迭代后所得到的U2即（4-1）同理求得节点三的电压即（4-2）由于节点4的无功功率未知所以先求计算节点4的无功功率如下式求得无功功率后即可求得节点4第一次迭代后的电压值保持电压幅值不变，修正后= 3.2.2第二次迭代保持电压幅值不变，修正后3.2.3第三次迭代由前此结果得保持电压幅值不变，修正后3.2.4第四次迭代保持电压幅值不变，修正后结果如下表所示：表3-1 迭代计算结果次数 电压U2U3U4第一次迭代1.8509-j0.2662.6593-j

30、0.82651.041-j0.1377第二次迭代3.286-j1.566-4.506-j0.1670.1845-j0.336第三次迭代-3.928+j0.793-3.3330-j1.323-3.760-j0.6882第四次迭代-1.855+j2.5860.054+j3.3290.4646+j0.9242综上所述利用高斯塞德尔法第三次迭代后得到的电压值分别为U2=-1.855+j2.586，U3=0.054+j3.329,U4 =0.4646+j0.9242结果如下表所示：4 MATPOWER中建模解决问题4.1 MATPOWER的技术规则4.1.1数据文件格式 在潮流计算时，把电网的各种参数写

31、成MATPOWER的数据文件格式，一般为M文件，MATPOWER4.0提供两种数据文件格式：version.1和version.2,我采用version.2格式。每个电网用变量名为“mpc”的结构体来定义，结构体mpc的不同字段用baseMVA、bus、branch、gen等来定义和返回电网的具体参数。在这些字段中，除baseMVA是标量外，其他都是矩阵。列的数据类似于标准的IEEE和PTI列的数据格式。MATPOWER中提供了大量的IEEE制定的典型模型，规范细节可以再caseformat.m中看到。mpc不同字段的简要介绍如下：（1）字段baseMVA是标量，设置基准电压，一般为100MV

32、A。 （2）字段bus是一个矩阵，用来设置电网母线参数。矩阵每一行都对应一个单一的母线，列的数据格式为bus_i、type、Pd、Qd、Gs、Bs、area、Vm、Va、baseKV、zone、Vmax、Vmin。bus_i用来设置母线编号；type用来设置母线类型，1为PQ节点母线，2为PV节点母线，3为平衡节点母线，4为孤立节点母线；Pd、Qd用来设置母线注入负荷的有功功率和无功功率；Gs、Bs用来设置该母线并联电导和电纳；baseKV用来设置母线基准电压；Vm和Va用来设置母线电压幅值和相位初值；Vmax、Vmin用来设置母线工作时最高、最低电压幅值；area、zone用来设置电网断面号

33、和分区号，一般为1。（3）字段gen 为一个矩阵，用来设置接入电网的发电机参数。矩阵的每一行对应一个发电机，列的数据格式为bus、Pg、Qg、Qmax、Qmin、Vg、mBase、status、Pmax、Pmin。bus用来设置接入发电机的母线编号；Pg、Qg用来设置接入发电机的有功功率和无功功率；Pmax、Pmin用来设置接入发电机的有功功率最大和最小允许值；Qmax、Qmin用来设置接入发电机的无功功率最大和最小允许值；Vg用来设置接入发电机的工作电压；mBase用来设置接入发电机的基准功率；status用来设置发电机的工作状态，1表示投入，0退出。（4）字段branch用来设置电网中各支

34、路参数。矩阵的每一行都对应一个单一支路，列的数据格式为fbus、tbus、r、x、b、rateA、rateB、rateC、ratio、angle、status、angmin、angmax。Fbus、tbus用来设置该支路由起始节点编号和终止节点编号；r、x、b用来设置该支路的电阻、电抗、电纳；rateA、rateB、rateC用来设置该支路长期、短期、紧急允许功率；ratio用来设置该支路的变比，导线时为0，变压器时为变比；angle为支路的相位角度，导线时为0；status用来设置支路工作状态，1为投入，0退出运行；angmin、angmax用来设置支路相位角度的最大最小差值。4.1.2控制

35、选项MATPOWER软件不但能够进行交流潮流计算，还能进行直流潮流、最优潮流等计算，进行计算时还可以选用不同的算法及输出格式。为实现以上功能，MATPOWER采用一个选项向量“mpoption”来达到对选项的控制。在变量名PF_ALG中通过改变值来达到改变算法，1对应牛顿法，2对应快速解耦算法（XB版本），3对应快速解耦算法（BX版本），4为高斯-赛德尔法，系统默认为1，即牛顿法。在PF_MAX_IT中设置牛顿法的最大迭代次数，PF_MAX_IT_GS设置高斯-赛德尔法的最大迭代次数，系统默认为1000次。为实现系统算法从牛顿法到高斯-赛德尔法转换，使用方法如下所示：在MATPOWER中文手册

36、还有对MATPOWER的使用详细介绍（在MATPOWER/docs中有manual.pdf）,在MATPOWER中有很多M文件很重要，比如runpf.m为定义一个运算潮流程序，Ext2int.m为将数据矩阵从外部节点编号转换为内部节点编号，loadcase.m为将数据从文件或者结构体重导入到数据矩阵中，还有很多，这些M文件使我们的建模得以实现，同时也是我们很好找到我们模型中问题。在我的建模过程中遇到了很多这样的问题，比如error('loadcase: syntax error or undefined data matrix(ices) in the filen%s', er

37、r5); mpc = loadcase(casedata);通过不断分析，加深理解，这些错误都是因为搭建的模型不足以有结果，所以系统的这些功能方便了我在编程时查找错误，这是这是它的一大优势。另外，MATPOWER还提供了很多实例，如case_ieee30.m 就是 IEEE30 节点系统，这些实例给了初学者很多启发，同时也起到了规范作用。4.2 MATPOWER实际建模4.2.1 建模中的实际问题 首先要将bus模型建立起来，因为建模方便的需要，将题目的1243点编号变成1432，所以点1为平衡节点，2、4点是PQ节点，点3是PV节点，根据bus的中的type类型的定义，1为3,2、4为1,3

38、为2。Bus的矩阵如下：% bus data% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vminmpc.bus = 1 3 0 0 0 0 1 1.04 0 100 1 1.1 0.94; 2 1 100 -50 0 0 1 1 0 100 1 1.1 0.94; 3 2 0 0 0 0 1 1.05 0 100 1 1.1 0.94; 4 1 -50 20 0 0 1 1 0 100 1 1.1 0.94; ; 接着要搭建支路模型，一开始我错误的判断支路是3条，其实后来仔细分析是5条，这都是因为基本功不扎实导致的，branch的矩

39、阵如下：% branch data% fbus tbus r x b rateA rateB rateC ratio angle status angmin angmaxmpc.branch = 1 4 0.316 0.949 0 0 0 0 0 0 1 -360 360; 3 4 0.0554 0.998 0 0 0 0 0 0 1 -360 360; 2 4 0.316 0.949 0 0 0 0 0 0 1 -360 360; 1 2 0.316 0.949 0 0 0 0 0 0 1 -360 360; 2 3 0.316 0.949 0 0 0 0 0 0 1 -360 360; ;

40、然后其他一些语句的编写形成了我的M文件，然后试着运行，但是却一次次失败，通过对IEEE提供的范例，然后我试着对IEEE的case与文件对比，在它们的差异中，对IEEE的case逆推，找出关键原因在于发电机模型必不可少，这可能是MATPOWER的开发中作者的意图在于对电力系统的完整分析，数据来源上充沛，而局部的节点导纳矩阵的建模确由于loadcase.m、runpf.m等函数文件中严格要求而受到限制。由于这次的目标仅仅是求出节点电压，所以我试着加入发电机的模型，是它的潮流计算能得出节点电压，发电机的模型建立如下：% generator data% bus Pg Qg Qmax Qmin Vg m

41、Base status Pmax Pmin mpc.gen = 1 0 0 300 -300 1 100 1 250 10; 3 163 0 300 -300 1 100 1 300 10; ;4.2.2 建模结果显示在MATLAB的命令窗口输入：计算机输出的结果如下： 图4.1 高斯-赛德尔法MATPOWER结果显示图 系统默认高斯-赛德尔法为1000次迭代，显示结果中会给出四个节点的电压幅值相位，而且会给出最大最小值，会给出各支路的损耗，节点功率。但这仅仅是发电机模型随机选取的结果，随机时会对节点电压产生较大的影响比如，会有下面的结果:或者出现这次会发现电压发生很大的变化，说明发电机模型会

42、有很大的影响，所以发电机模型要认真建立，通过分析，我把发电机模型设为如下: % generator data% bus Pg Qg Qmax Qmin Vg mBase status Pmax Pmin mpc.gen = 1 0 0 300 -300 1 100 1 250 10; 3 30 0 300 -300 1 100 1 300 10; ;运行结果如下: 以上便是运行结果。 总结拿到课设题目之后的第一件事情就是查找相关资料以熟悉电力系统潮流计算的相关理论知识，无疑这一步是至关重要的一步。经过深入学习电力系统潮流计算的理论知识，是我真确的明白潮流计算的重要性及潮流计算的计算流程，这为我

43、后面编写潮流计算程序打下了坚实的基础。接着不得不做的一件事就是深入学习MATLAB，MATLAB是完成毕业设计的工具。MATLAB 有着强大的数值计算的功能，它的编程语言同样简短而强大。通过查阅相关资料，我了解了高斯-赛德尔法的具体含义，还有MATLAB是很好的分析电力系统潮流的软件，所以我有了大致的思路，运用MATPOWER软件解决问题十分方便，而且了解这个电力系统广泛应用的软件对今后的学习工作十分有益。 我主要参考了去年买的一本MATPOWER/Simulink电力系统建模与仿真，在cornell主页下载了MATPOWER软件，按照相关步骤安装，然后我按照题目进行建模，在实际操作中遇到了相

44、当多的问题，首先是runpf功能多次运行不成功，每次报错都是同一原因，即loadcase里的250行对错误5的定义，由于半路出家，我对它的意思不能理解，通过对范例的逆推，我总结出这是由于建模的不完整导致的，所以，我通过加入发电机模型，使这个模型得以完成。同时也有像mpc = loadcase(casedata);这个错误我经常遇到，一时间也不能理解，直到现在我也仅仅理解成负荷点建模不正确，到底什么意思，估计还需要更深入地研究。在整个探究过程中，我不断遇到问题，不断上网找资料，从雅可比、高斯-赛德尔、到MATPOWER中文手册，在这些资料中我受益匪浅，同时我也想很多同学请教，他们热心帮助给了我很

45、多启发，特别是袁辉同学热心帮助我，给了我很大鼓励。这几天时间里有时一整天都没什么进展，有时会因为一个小小的启发而欢喜，喜怒哀乐各种滋味也许只有认认真真去探究才能得到。 MATPOWER是康奈尔大学在90年代开发出来的，令我震撼的是90年代这所遐迩闻名的常春藤高校就做出这么有前瞻性、创新性的工作，这使得MATLAB、Simulink在电力系统中应用更为广泛。现在的浙江大学甘德强教授也是当年的主要研发人员。随着电气、电子、信息、计算机的学科交叉性越来越明显，如何把握好学科的发展方向，如何使我们的电气专业更好的向前发展，美国高校给了我们很大的启发，美国的高校电气一级学科全称电气与计算机，这也说明了计

46、算机融入到电气工程领域已是必然的趋势，只有掌握好计算机知识，才能把握好电气发展的生命线。另外，随着电力系统的发展，新能源如光伏发电、风能等领域的发展，电力系统肯定会越来越复杂，比如风能大规模的接入电力系统，对电力系统稳定性可靠性的影响，以及替代等功率的火电厂时经济性，这些问题只会使我们的建模越来越困难，现有的比如牛顿法、PQ法很有可能满足不了实际要求，如何是我们数学学科应用到电力系统分析中是我们今后发展的背景。对于我个人而言，这次课设的启发在于思维的拓展，至少在学科结合上，以后我会在学好主干学科时，花一些时间在数学方法、计算机软件上。记得大二期末时我们参加基础强化训练时，第一次接触到MATLA

47、B感觉很难，很生僻，还有点抵触情绪，现在想想真的很幼稚，选择了这条路那只有勇往直前，不要畏惧什么，学习时不要被动，主动的学习会发现另一片天，生活就是这样，主动地生活，每天都是希望，只要有希望在，我们每天都会动力十足，信心满满。虽然课程设计过程中经历了很多，付出了很多，但由于能力和时间的有限仍有许多值得改进的地方。例如登陆界面及潮流计算主界面的美观性有待提高，软件与用户的互动交流还不够好，这些问题都有待改善。总之，通过这次毕业设计，使我学会了如何更好的学习，如何更好的查资料，如何更好的找信息。感谢老师给予我们这样一个动手设计的机会，让我着实提高了自己的动手设计能力。致谢通过这次的课程设计，我不仅

48、做到了学以致用，还锻炼了我的独立思考、学会分析问题和遇到困难时能勇于尝试和探索寻找方法解决问题的能力。本次近两个星期的课程设计能够顺利完成离不开宋老师的尽责辅导和同学们的热心帮助。由于宋老师在课上对于潮流计算的讲解非常仔细，再加上对课本进行了温习，因此成功的完成了本次设计。在设计过程中，无论是在构思步骤，参数计算，还是建模上，宋老师在课上的讲解都给与了我极大的帮助，在我设计遇到困难无法进行时，同学们也给予了帮助和引导。在此我要向宋老师以及寻求过帮助的同学们致以最衷心的感谢和深深的敬意。最后，在设计的完成之际，我要感谢宋老师能从百忙之中抽出时间来批评指正我的设计，提出宝贵的意见，其他同学也给予我很多的帮助，在这里一并感谢他们。参考文献1何仰赞，汪增银