导纳矩阵课程设计Word文档格式.docx

1、 电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控制。利用电子数字计算机进行潮流计算从50年代中期就已经开始了。在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的，对潮流计算的要求可以归纳为以下几点：|（1） 计算方法的可靠性或收敛性；（2） 对计算机内存量的要求；（3） 计算速度；（4） 计算的方便性和灵活性。导纳矩阵计算的意义（1）在电网规划阶段,通过潮流计算,合理规划电源容量及接入点,合理规划网架,选择无功补偿方案,满足规划水平的大、小方式下潮流交换控制、

2、调峰、调相、调压的要求。（2）在编制年运行方式时,在预计负荷增长及新设备投运基础上,选择典型方式进行潮流计算,发现电网中薄弱环节,供调度员日常调度控制参考,并对规划、基建部门提出改进网架结构,加快基建进度的建议。（3）正常检修及特殊运行方式下的潮流计算,用于日运行方式的编制,指导发电厂开机方式,有功、无功调整方案及负荷调整方案,满足线路、变压器热稳定要求及电压质量要求。（4）预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。第三章 导纳矩阵设计题目 导纳矩阵计算题目表3-1 在100MVA的基础上6条母线的系统图支路母线节点电阻电抗比率1243657图3-1 系统图该系统

3、图是一个6节点，7支路的电力网络，通过手工计算和计算机计算得到该系统图的导纳矩阵，计算结果进行比较，观察两次运算结果是否一样。 对课题的分析及求解思路 此电力系统是一个6节点，7支路的电力网络。综合比较牛顿拉夫逊法（直角坐标、极坐标）、PQ分解法等多种求解方法的特点，最后确定采用牛顿拉夫逊法（极坐标）。因为此方法所需解的方程组最少。通过对本课程设计的理解分析，将整个系统的简化图分析如图3-2：图3-2 系统简化图第四章导纳矩阵计算算法及手工计算 节点电压方程在电路中我们学过利用节点电压方程求解某几条支路的电流，现以下图4-1与图4-2为例推导节点电压方程组。 图4-1节点电压法为例 图4-2用

4、电流源代替电压源为例图4-1表示了一个具有两个电源和一个等值负荷的系统。为电源电势，为电源的内部导纳，为负荷的等值导纳，为各支路的导纳。如果以地为电压参考点，设节点1、2、3的电压为，根据基尔霍夫电流KCL法对节点1、2、3列节点电流方程得式（4-1）： （4-1）上式中左端为节点1、2、3流出的电流，右端为注入个节点的电流。由上式可以得到一个等效的等值电路图4-1。图4-2中利用了电流源代替的电压源。在图4-2中可知的式（4-2）： （4-2） 为等值电流源向网络注入的电流。将与式（4-1）联立得式（4-3）： （4-3）上式中称为节点1、2、3的自导纳，称为相应节点之间的互导纳。因此，在一

5、般情况下，在电力网络中有n个节点，则可以按式（4-3）的形式列出n个节点方程式，也可用矩阵的形式表示。其中 （4-4）分别为节点注入电流列向量及节点电压列向量； （4-5） 为节点导纳矩阵，其中对角元素为节点 i的自导纳，非对角线为节点 i与节点j之间的互导纳。 节点导纳矩阵节点导纳矩阵既可根据自导纳和互导纳的定义直接求取，也可根据电路知识中找出改网络的关联矩阵，在节点电压方程的矩阵形式进行求解。本章节我们主要讨论的是直接求解导纳矩阵。根据节点电压方程章节我们知道，在利用电子数字计算机计算电力系统运行情况是，多采用形式的节点方程式。其中阶数等于电力网络的节点数。从而可以得到n个节点时的节点导纳

6、矩阵方程组（4-6）如下： （4-6） 由此可以得到n个节点导纳矩阵：它反映了网络的参数及接线情况，因此导纳矩阵可以看成是对电力网络电气特性的一种数学抽象。由导纳短阵所联系的节点方程式是电力网络广泛应用的一种数学模型。 通过上面的讨论，可以看出节点导纳矩阵的有以下特点： （1）导纳矩阵的元素很容易根据网络接线图和支路参数直观地求得，形成节点导纳矩阵的程序比较简单。 （2）导纳矩阵为对称矩阵。由网络的互易特性易知（3）导纳矩阵是稀疏矩阵。它的对角线元素一般不为零，但在非对角线元素中则存在不少零元素。在电力系统的接线图中，一般每个节点与平均不超过34个其他节点有直接的支路连接。因此，在导纳矩阵的非

7、对角线元素中每行仅有34个非零元素，其余的都是零元素，而且网络的规模越大，这种现象越显著。导纳矩阵的对称性和稀疏性对于应用计算机求解电力系统问题有很大的影响。如果能充分地利用这两个特点，如在程序设计中储存导纳矩阵的对角元素和上三角元素（或下三角元素），排除零元素的储存和运算，就可以大大地节省储存单元和提高计算速度。节点导纳矩阵的形式可归纳如下：（1）导纳矩阵的阶数等于电力网络的节点数。（2）导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连得不接地支路数。（3）导纳矩阵各对角元素，即节点的自导纳等于相应节点之间的支路导纳之和。（4）导纳矩阵非对角元素，即节点之间的互导纳等于相应节点之间的支

8、路导纳的负值。而在电力系统中进行潮流计算时，往往要计算不同接线下的运行状况，例如，改变变压器主抽头时，潮流分布也随之变化，以及改变其他设备参数进行计算潮流分布，此时就需要导出变化时的导纳矩阵就需要对所设计的程序进行参数设定，而不需要重复上述步骤去导出所求的导纳矩阵。变压器的型等值电路在电力系统潮流计算中，往往要计算节点导纳矩阵，而我们计算节点导纳矩阵采用节点电压法来实现，如在变压器构成的电力系统中，需要将变压器模型转变成变压器型等值电路，在利用电路知识列节点电压方程，从而导出所需的导纳矩阵。图 4-3双绕组变压器的型等值电路（i，j为节点）而在电力系统潮流计算中一般采用标幺值进行计算，标幺值公

9、式如下：所以，变压器转变成型等值电路时，我们采用标幺值计算，使所求参数为变压器变比k的函数。而在一个已经归算好的电力系统网中，若改变变压器的分接头来进行调压，这时变压器的等值电路参数也会相应得改变，此时采用型等值电路进行折算就显得较为方便。下面是变压器的型等值电路分析过程： 如不计如果采用标么值计算，元件参数都应归算到同一基准值时得标么值，才能在同一个等值电路上分析和计算。励磁支路的影响，双绕组变压器可用其阻抗与一个理想变压器串联的电路表示，如图所示。理想变压器只有一个参数，那就是变比现以变压器阻抗按实际变比归算到低压侧的情况为例，推导出双绕组变压器的型等值电路。流入和流出理想变压器的功率相等

10、： （分别为变压器高、低绕组的实际电压） （4-7） （4-8）联立（4-7）、（4-8）两个公式解得： （4-9） （4-10）根据电路原理节点1、2的节点电流方程具有如下形式： （4-11）将式（4-9）、（4-10）与式（4-11）比较得（4-12）: （4-12） 因此可以的得到各支路导纳为 （4-13） 导纳矩阵的手工计算在计算电力系统网络的潮流分布时，我们需要把变压器转化成变压器的型等值电路来进行计算器等效导纳，以下是变压器的型等值电路（见图4-4），变压器的型等值电路（见图4-5），变压器的型等值电路 :1 1 2 1 2 图4-4变压器的型等值电路 所求变压器参数数据如下： 1

11、 2 1 :1 2 图4-5变压器所求电力系统网络的各个节点的相关参数如下（手工计算）：通过手工计算得到了导纳矩阵，可以知道导纳矩阵的特点，手工计算非常繁琐，当节点很多时，将会更加的麻烦，所以要利用计算机来计算。下章将介绍用计算机来得到导纳矩阵，并加以比较。第五章 Matlab概述 Matlab简介MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决

12、方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 Matlab的应用MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情

13、简捷得多，并且mathwork也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C+ ，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。MATLAB 产品族可以用来进行以下各种工作： 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真 数字图像处理技术 数字信号处理技术 通讯系统设计与仿真 财务与金融工程 MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和

14、测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用 MATLAB 函数集）扩展了 MATLAB 环境，以解决这些应用领域内特定类型的问题。. Matlab中的一些命令1）一般MATLAB命令格式为 输出参数1，输出参数2，=（命令名）（输入参数1，输入参数2，） 输出参数用方括号，输入参数用圆括号如果输出参数只有一个可不使用括号。2）可用、键来重现已输入的数据或命令。用、键来移动光标进行修改。3）所有MATLAB命令都用小写字母。大写字母和小写字母分别表示不同的变量。4）常用数有特定的名字，如pi（=） 、Inf（=）、NaN则表示不定型求得的结果（如0/0）。5）

15、矩阵的输入要一行一行的进行，每行各元素用空格或（，）分开，每行用（；）分开。 6）MATLAB书写格式为A=1 2 3 ；4 5 6 ；7 8 9 在MATLAB中运行如下程序可得到A矩阵 a=1 2 3;4 5 6;7 8 9 a = 1 2 3 4 5 6 7 8 97）需要显示命令的计算结果时，则语句后面不加“；”号，否则要加“；”号。 运行下面两种格式可以看出他们的区别7 8 9 a=1 2 3;7 8 9; a = （不显示计算结果） 1 2 38）当输入语句过长需要换行时，应加上“”后再回车，则可连续输入。9）diary命令使用该命令可以在窗口中以ASCII码形式记录所有的输入和输

16、出。但这个命令不是存储数据，而是存储输入与屏幕上输出的内容。它可以记录下工作的过程。在每个工作过程之前使用该命令，工作结束后使用diary off 则能将整个工作过程记录下来。10）save 命令该命令存储定义的变量或演算结果，也可以用来存储指定的变量。11）d1,d2,d3,.=size（a） 求矩阵的大小，对m*n二维矩阵，第一个为行数m,第二个为列数n。如果输入 calear a b c ，则表示清除工作空间中指定变量a,b,c；如果仅仅输入calear命令，则清除整个工作空间。 与此同时，MATLAB具有强大的矩阵运算功能, 但由于我们在求节点导纳矩阵时用的不多, 因此这里我们只作简单

17、介绍。1） 在MATLAB中表示一个矢量要用方括号， 而列矢量的输入只需在行矢量输入格式基础上加转置符（）即可。 如 x=1 2 3;4 5 6 x = 而 x=1 2 3;4 5 6 （加转置符） x = 1 4 2 5 3 6注意上面两式的区别。2） 下面三条命令可以产生一个行矢量1a=linspace（x,y,n） 2a=logspace（x,y,n） 3a=x:n:y第一条命令可以在线性空间产生一个值在10x至10y之间间隔点数为n的行矢量（一组数据）。第二条命令可以在对数空间产生一个值在x至y之间等间隔的行矢量（一组数据）。其行矢量的起始值是x, 终值为y, 点数为n。第三条命令产生

18、X至y步长为n的行矢量。3） 矩阵的加, 减, 乘, 除等, 和其它语言书写一样。但要注意的是在运算符前面加有（.）则表示是元素对元素的操作.第六章 导纳矩阵计算流程图及源程序 导纳矩阵计算流程图程序编程框图如下：软件启动输入节点数n支路数n1数组Wk=1s=W（k,1）g=W（k,2）t=complex（W（i,3）,W（i,4）Y（s,g）=Y（s,g）-1./（t*W（i,6）;Y（g,s）=Y（s,g）; Y（g,g）=Y（g,g）+1./（t*（W（i,6）2）+W（i,5）./2; k=k+1Y（s,s）=Y（s,s）+1./t+W（i,5）./2; 假k=n1真图 6-1 潮流计

19、算流程图导纳矩阵计算源程序及结果节点导纳矩阵程序如下:n=input（请输入节点数：n=）;n1=input（请输入支路数：n1=W=input（请输入由支路参数形成的矩阵：W=M=input（请输入由节点号及其对地阻抗形成的矩阵：M=Y=zeros（n）;for i=1:n if M（i,2）=0; s=M（i,1）; Y（s,s）=1./M（i,2）; endn1 if W（i,7）=0 s=W（i,1）;g=W（i,2）; else s=W（i,2）;g=W（i,1）; t=complex（W（i,3）,W（i,4）; Y（s,g）=Y（s,g）-1./（t*W（i,6）; Y（g,s）

20、=Y（s,g）; Y（s,s）=Y（s,s）+1./t+W（i,5）./2; disp（导纳矩阵Y= disp（Y）按照程序调试的如下结果：n=6n1=7W=1 2 0 1; 1 4 0 0; 1 6 0 0; 2 5 0 0; 3 5 0 0; 4 3 0 1; 4 6 0 0M=1 0;2 0;3 0;4 0;5 0;6 0导纳矩阵总结这次课程设计的主要任务是用计算机编写程序求出潮流计算中要用到的导纳矩阵。为了确定结果是否正确，先做了一个手工计算，得到手算结果，然后比较两种方法（手算法与计算机算法）所得到的结果，看程序是否正确。通过MATLAB仿真，由程序得到的结果与手算的结果一致，不同之

21、处可能是由于手算法与计算机算法，在取计算结果时，保留的精度不一样所致。因此，可认为两种方法所得结果是一样的。在保证手算结果正确的基础上，也就确定了程序的准确性。同时，通过此次课程设计也让我明白了要使电力系统运行的稳定，必须经过精密的设计和计算。在课程设计过程中，学会了导纳矩阵的修改及其应用，加深了电力系统潮流计算的理解。 在这次课程设计中，我们尽量按照老师的要求做，但在具体的操作过程中，还是出现了很多的问题。搞完这个课程设计让我感觉电力系统分析和MATLAB都是非常有用的课程。由于对MATLAB的知识了解的比较少，在编程的时候出现了卡壳的情况。这次课程设计后，我一定要重新对MATLAB这门课程

22、做进一步的了解。为在以后的毕业课程设计中能多些方案。也为我子以后走上工作岗位，提升自己的专业技能，打下扎实的基础。 还有就是在十几天的课程设计中，使我养成了很好的学习习惯，和对学习知识的严谨的态度，同时也养成了积极查阅相关资料的好习惯，好习惯的养成是来之不易的，我相信在以后的学习和工作中，我将继续保持这些良好的习惯，并积极努力的学习。让自己更上一层楼。同时在此也感谢一直指导我的老师，此次课程设计的完成与老师的指导师分不开的，终在我们的一起努力下，完成了这门课程设计。在此对指导老师致以衷心的感谢。还要感谢学院图书馆为我提供丰富的参考资料，也感谢班上同学给了我很多宝贵的意见和参考，使我获益很多。参

23、考文献1 孟祥萍. 电力系统分析高等教育出版社.20042 李维波. MATLAB在电气工程中应用中国电力出版社.2007 3 吴天明.赵新力.刘建存 MATLAB电力系统设计与分析国防工业出版社.20074 邱晓燕.刘天琪 电力系统分析的计算机算法中国电力出版社.20095 王沫然 .与科学计算 .北京:电子工业出版社,20016 华智明，张瑞林 .电力系统 .重庆：重庆大学出版社，19977 方富淇 电力系统分析 .北京：水利电力出版社，19898 刘军.MATLAB 在电力系统分析中的应用D：电力系统及自动化学报.9 刘卫国.陈昭平.张颖.MATLAB 程序设计与应用M.北京：高等教育出版社.