还剩18页未读,继续阅读
本资源只提供10页预览,全部文档请下载后查看!喜欢就下载吧,查找使用更方便
文本内容:
第一章引言导纳矩阵计算是电力系统分析最基本的计算除它自身的重要作用之外,在《电力系统分析综合程序》PSASP中,还是网损计算、静态安全分析、暂态稳定计算、小干扰静态稳定计算、短路计算、静态和动态等值计算的基础对规划中的电力系统,通过计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求;对运行中的电力系统,通过计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全,系统中所有母线的电压是否在允许的范围以内,系统中各种元件线路、变压器等是否会出现过负荷,以及可能出现过负荷时应事先采取哪些预防措施等
1.1课程设计目的与要求
1.
1.1设计目的
1.掌握电力系统潮流计算导纳矩阵的基本原理;
2.掌握并能熟练运用一门计算机语言(MATLAB语言或C语言或C++语言);
3.采用计算机语言对潮流计算导纳矩阵进行计算机编程
1.
1.2设计要求
1.程序源代码;
2.给定题目的输入,输出文件;
3.程序说明;
4.给定系统的程序计算过程;
5.给定系统的手算过程(至少迭代2次)
1.2设计内容
1.根据电力系统网络推导电力网络数学模型,写出节点导纳矩阵;
2.赋予各节点电压变量(直角坐标系形式)初值后,求解不平衡量;
3.形成雅可比矩阵;
4.求解修正量后,重新修改初值,从2开始重新循环计算;
5.求解的电压变量达到所要求的精度时,再计算各支路功率分布、功率损耗和平衡节点功率;
6.上机编程调试;
7.计算分析给定系统导纳矩阵分析并与手工计算结果做比较分析;
8.书写课程设计说明书,准备答辩第二章导纳矩阵计算概述
2.1导纳矩阵计算简介电力系统潮流导纳矩阵计算是研究电力系统稳态运行情况的一种计算,它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态各母线的电压各元件中流过的功率,系统的功率损耗等等在电力系统规划的设计和现有电力系统运行方式的研究中,都需要利用潮流计算来定量的分析比较供电方案或运行方式的合理性可靠性和经济性此外,电力系统的潮流计算也是计算机系统动态稳定和静态稳定的基础,所以潮流计算是研究电力系统的一种和重要和基础的计算电力系统潮流计算也分为离线计算和在线计算两种,前者主要用于系统规划设计和安排系统的运行方式,后者则用于正在运行系统的经常监视及实时控制利用电子数字计算机进行潮流计算从50年代中期就已经开始了在这20年内,潮流计算曾采用了各种不同的方法,这些方法的发展主要围绕着对潮流计算的一些基本要求进行的,对潮流计算的要求可以归纳为以下几点
(1)计算方法的可靠性或收敛性;
(2)对计算机内存量的要求;
(3)计算速度;
(4)计算的方便性和灵活性
2.2导纳矩阵计算的意义1在电网规划阶段通过潮流计算合理规划电源容量及接入点合理规划网架选择无功补偿方案满足规划水平的大、小方式下潮流交换控制、调峰、调相、调压的要求2在编制年运行方式时在预计负荷增长及新设备投运基础上选择典型方式进行潮流计算发现电网中薄弱环节供调度员日常调度控制参考并对规划、基建部门提出改进网架结构加快基建进度的建议3正常检修及特殊运行方式下的潮流计算用于日运行方式的编制指导发电厂开机方式有功、无功调整方案及负荷调整方案满足线路、变压器热稳定要求及电压质量要求4预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案第三章导纳矩阵设计题目
3.1导纳矩阵计算题目表3-1在100MVA的基础上6条母线的系统图支路母线节点电阻电抗比率
1120.
0000.
3001.
0252140.
0970.
4073160.
1230.
5184250.
2820.
6405350.
7231.
0506430.
0000.
1331.
1007460.
0800.370图3-1系统图该系统图是一个6节点,7支路的电力网络,通过手工计算和计算机计算得到该系统图的导纳矩阵,计算结果进行比较,观察两次运算结果是否一样
3.2对课题的分析及求解思路此电力系统是一个6节点,7支路的电力网络综合比较牛顿拉夫逊法(直角坐标、极坐标)、PQ分解法等多种求解方法的特点,最后确定采用牛顿拉夫逊法(极坐标)因为此方法所需解的方程组最少通过对本课程设计的理解分析,将整个系统的简化图分析如图3-2图3-2系统简化图第4章导纳矩阵计算算法及手工计算
4.1节点电压方程在电路中我们学过利用节点电压方程求解某几条支路的电流,现以下图4-1与图4-2为例推导节点电压方程组图4-1节点电压法为例图4-2用电流源代替电压源为例图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的自导纳,称为相应节点之间的互导纳因此,在一般情况下,在电力网络中有n个节点,则可以按式(4-3)的形式列出n个节点方程式,也可用矩阵的形式表示其中4-4)分别为节点注入电流列向量及节点电压列向量;4-5)为节点导纳矩阵,其中对角元素为节点i的自导纳,非对角线为节点i与节点j之间的互导纳
4.2节点导纳矩阵节点导纳矩阵既可根据自导纳和互导纳的定义直接求取,也可根据电路知识中找出改网络的关联矩阵,在节点电压方程的矩阵形式进行求解本章节我们主要讨论的是直接求解导纳矩阵根据节点电压方程章节我们知道,在利用电子数字计算机计算电力系统运行情况是,多采用形式的节点方程式其中阶数等于电力网络的节点数从而可以得到n个节点时的节点导纳矩阵方程组(4-6)如下4-6)由此可以得到n个节点导纳矩阵它反映了网络的参数及接线情况,因此导纳矩阵可以看成是对电力网络电气特性的一种数学抽象由导纳短阵所联系的节点方程式是电力网络广泛应用的一种数学模型通过上面的讨论,可以看出节点导纳矩阵的有以下特点
(1)导纳矩阵的元素很容易根据网络接线图和支路参数直观地求得,形成节点导纳矩阵的程序比较简单
(2)导纳矩阵为对称矩阵由网络的互易特性易知
(3)导纳矩阵是稀疏矩阵它的对角线元素一般不为零,但在非对角线元素中则存在不少零元素在电力系统的接线图中,一般每个节点与平均不超过3~4个其他节点有直接的支路连接因此,在导纳矩阵的非对角线元素中每行仅有3~4个非零元素,其余的都是零元素,而且网络的规模越大,这种现象越显著导纳矩阵的对称性和稀疏性对于应用计算机求解电力系统问题有很大的影响如果能充分地利用这两个特点,如在程序设计中储存导纳矩阵的对角元素和上三角元素(或下三角元素),排除零元素的储存和运算,就可以大大地节省储存单元和提高计算速度节点导纳矩阵的形式可归纳如下
(1)导纳矩阵的阶数等于电力网络的节点数
(2)导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连得不接地支路数
(3)导纳矩阵各对角元素,即节点的自导纳等于相应节点之间的支路导纳之和
(4)导纳矩阵非对角元素,即节点之间的互导纳等于相应节点之间的支路导纳的负值而在电力系统中进行潮流计算时,往往要计算不同接线下的运行状况,例如,改变变压器主抽头时,潮流分布也随之变化,以及改变其他设备参数进行计算潮流分布,此时就需要导出变化时的导纳矩阵就需要对所设计的程序进行参数设定,而不需要重复上述步骤去导出所求的导纳矩阵
4.3变压器的∏型等值电路在电力系统潮流计算中,往往要计算节点导纳矩阵,而我们计算节点导纳矩阵采用节点电压法来实现,如在变压器构成的电力系统中,需要将变压器模型转变成变压器∏型等值电路,在利用电路知识列节点电压方程,从而导出所需的导纳矩阵图4-3双绕组变压器的∏型等值电路(i,j为节点)而在电力系统潮流计算中一般采用标幺值进行计算,标幺值公式如下所以,变压器转变成∏型等值电路时,我们采用标幺值计算,使所求参数为变压器变比k的函数而在一个已经归算好的电力系统网中,若改变变压器的分接头来进行调压,这时变压器的等值电路参数也会相应得改变,此时采用∏型等值电路进行折算就显得较为方便下面是变压器的∏型等值电路分析过程如不计如果采用标么值计算,元件参数都应归算到同一基准值时得标么值,才能在同一个等值电路上分析和计算励磁支路的影响,双绕组变压器可用其阻抗与一个理想变压器串联的电路表示,如图所示理想变压器只有一个参数,那就是变比现以变压器阻抗按实际变比归算到低压侧的情况为例,推导出双绕组变压器的∏型等值电路流入和流出理想变压器的功率相等、分别为变压器高、低绕组的实际电压)(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-4),变压器的∏型等值电路(见图4-5),变压器的∏型等值电路
1.1:11212图4-4变压器的∏型等值电路所求变压器参数数据如下
1211.025:12图4-5变压器的∏型等值电路所求变压器参数数据如下所求电力系统网络的各个节点的相关参数如下(手工计算)通过手工计算得到了导纳矩阵,可以知道导纳矩阵的特点,手工计算非常繁琐,当节点很多时,将会更加的麻烦,所以要利用计算机来计算下章将介绍用计算机来得到导纳矩阵,并加以比较第五章Matlab概述
5.1Matlab简介MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平 MATLAB和Mathematica、Maple并称为三大数学软件它在数学类科技应用软件中在数值计算方面首屈一指MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等,主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域
5.2Matlab的应用MATLAB的基本数据单位是矩阵,它的指令表达式与数学、工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的事情简捷得多,并且mathwork也吸收了像Maple等软件的优点使MATLAB成为一个强大的数学软件在新的版本中也加入了对C,FORTRAN,C++,JAVA的支持可以直接调用用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用,此外许多的MATLAB爱好者都编写了一些经典的程序,用户可以直接进行下载就可以用MATLAB产品族可以用来进行以下各种工作 ●数值分析 ●数值和符号计算 ●工程与科学绘图 ●控制系统的设计与仿真 ●数字图像处理技术 ●数字信号处理技术 ●通讯系统设计与仿真●财务与金融工程MATLAB的应用范围非常广,包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域附加的工具箱(单独提供的专用MATLAB函数集)扩展了MATLAB环境,以解决这些应用领域内特定类型的问题
5.
3.Matlab中的一些命令1)一般MATLAB命令格式为[输出参数1,输出参数2,……]=(命令名)(输入参数1,输入参数2,……)输出参数用方括号,输入参数用圆括号如果输出参数只有一个可不使用括号2)可用↑、↓键来重现已输入的数据或命令用←、→键来移动光标进行修改3)所有MATLAB命令都用小写字母大写字母和小写字母分别表示不同的变量4)常用数有特定的名字,如pi(=
3.141596)、Inf(=∞)、NaN则表示不定型求得的结果(如0/0)5)矩阵的输入要一行一行的进行,每行各元素用空格或(,)分开,每行用(;)分开6)MATLAB书写格式为A=[123;456;789]在MATLAB中运行如下程序可得到A矩阵a=[123;456;789]a=1234567897)需要显示命令的计算结果时,则语句后面不加“;”号,否则要加“;”号运行下面两种格式可以看出他们的区别a=[123;456;789]a=[123;456;789];a=(不显示计算结果)1234568)当输入语句过长需要换行时,应加上“…”后再回车,则可连续输入9)diary命令使用该命令可以在窗口中以ASCII码形式记录所有的输入和输出但这个命令不是存储数据,而是存储输入与屏幕上输出的内容它可以记录下工作的过程在每个工作过程之前使用该命令,工作结束后使用diaryoff则能将整个工作过程记录下来10)save命令该命令存储定义的变量或演算结果,也可以用来存储指定的变量11)[d1d2d
3..]=sizea求矩阵的大小,对m*n二维矩阵,第一个为行数m第二个为列数n如果输入calearabc,则表示清除工作空间中指定变量abc;如果仅仅输入calear命令,则清除整个工作空间与此同时,MATLAB具有强大的矩阵运算功能但由于我们在求节点导纳矩阵时用的不多因此这里我们只作简单介绍1在MATLAB中表示一个矢量要用方括号,而列矢量的输入只需在行矢量输入格式基础上加转置符‘即可如x=[123;456]x=123456而x=[123;456]加转置符x=142536注意上面两式的区别2下面三条命令可以产生一个行矢量1a=linspacexyn2a=logspacexyn3a=[x:n:y]第一条命令可以在线性空间产生一个值在10x至10y之间间隔点数为n的行矢量一组数据第二条命令可以在对数空间产生一个值在x至y之间等间隔的行矢量一组数据其行矢量的起始值是x终值为y点数为n第三条命令产生X至y步长为n的行矢量3矩阵的加减乘除等和其它语言书写一样但要注意的是在运算符前面加有.则表示是元素对元素的操作.第六章导纳矩阵计算流程图及源程序
6.1导纳矩阵计算流程图程序编程框图如下软件启动输入节点数n支路数n1数组Wk=1s=Wk1g=Wk2t=complexWi3Wi4Ysg=Ysg-
1./t*Wi6;Ygs=Ysg;Ygg=Ygg+
1./t*Wi6^2+Wi
5./2;k=k+1Yss=Yss+
1./t+Wi
5./2;假k=n1真图6-1潮流计算流程图
6.2导纳矩阵计算源程序及结果节点导纳矩阵程序如下:n=input请输入节点数n=;n1=input请输入支路数n1=;W=input请输入由支路参数形成的矩阵W=;M=input请输入由节点号及其对地阻抗形成的矩阵M=;Y=zerosn;fori=1:nifMi2~=0;s=Mi1;Yss=
1./Mi2;endendfori=1:n1ifWi7==0s=Wi1;g=Wi2;elses=Wi2;g=Wi1;endt=complexWi3Wi4;Ysg=Ysg-
1./t*Wi6;Ygs=Ysg;Ygg=Ygg+
1./t*Wi6^2+Wi
5./2;Yss=Yss+
1./t+Wi
5./2;enddisp导纳矩阵Y=;dispY按照程序调试的如下结果请输入节点数n=6请输入支路数n1=7请输入由支路参数形成的矩阵W=[
120.
0000.
30001.0251;
140.
0970.
40701.0000;
160.
1230.
51801.0000;
250.
2820.
64001.0000;
350.
7231.
05001.0000;
430.
0000.
13301.1001;
460.
0800.
37001.0000]请输入由节点号及其对地阻抗形成的矩阵M=[10;20;30;40;50;60]导纳矩阵总结这次课程设计的主要任务是用计算机编写程序求出潮流计算中要用到的导纳矩阵为了确定结果是否正确,先做了一个手工计算,得到手算结果,然后比较两种方法(手算法与计算机算法)所得到的结果,看程序是否正确通过MATLAB仿真,由程序得到的结果与手算的结果一致,不同之处可能是由于手算法与计算机算法,在取计算结果时,保留的精度不一样所致因此,可认为两种方法所得结果是一样的在保证手算结果正确的基础上,也就确定了程序的准确性同时,通过此次课程设计也让我明白了要使电力系统运行的稳定,必须经过精密的设计和计算在课程设计过程中,学会了导纳矩阵的修改及其应用,加深了电力系统潮流计算的理解在这次课程设计中,我们尽量按照老师的要求做,但在具体的操作过程中,还是出现了很多的问题搞完这个课程设计让我感觉电力系统分析和MATLAB都是非常有用的课程由于对MATLAB的知识了解的比较少,在编程的时候出现了卡壳的情况这次课程设计后,我一定要重新对MATLAB这门课程做进一步的了解为在以后的毕业课程设计中能多些方案也为我子以后走上工作岗位,提升自己的专业技能,打下扎实的基础还有就是在十几天的课程设计中,使我养成了很好的学习习惯,和对学习知识的严谨的态度,同时也养成了积极查阅相关资料的好习惯,好习惯的养成是来之不易的,我相信在以后的学习和工作中,我将继续保持这些良好的习惯,并积极努力的学习让自己更上一层楼同时在此也感谢一直指导我的老师,此次课程设计的完成与老师的指导师分不开的,终在我们的一起努力下,完成了这门课程设计在此对指导老师致以衷心的感谢还要感谢学院图书馆为我提供丰富的参考资料,也感谢班上同学给了我很多宝贵的意见和参考,使我获益很多参考文献
[1]孟祥萍.《电力系统分析》.高等教育出版社.2004
[2]李维波.《MATLAB在电气工程中应用》.中国电力出版社.2007
[3]吴天明.赵新力.刘建存《MATLAB电力系统设计与分析》.国防工业出版社.2007
[4]邱晓燕.刘天琪《电力系统分析的计算机算法》.中国电力出版社.2009
[5]王沫然.MATLAB
6.0与科学计算.北京:电子工业出版社2001
[6]华智明,张瑞林.电力系统.重庆重庆大学出版社,1997
[7]方富淇电力系统分析.北京水利电力出版社,1989
[8]刘军.MATLAB在电力系统分析中的应用[D]电力系统及自动化学报.
2000.4
[9]刘卫国.陈昭平.张颖.MATLAB程序设计与应用[M].北京高等教育出版社.
2002.
655.0+j
13.
050.0+j
5.
030.0+j
18.0P5=
50.
16125431.1:
11.025:1输出导纳矩阵。