1、电力系统稳态分析课程设计doc内蒙古科技大学本科生课程设计说明书题 目基于MATLAB的电力系统复杂潮流分析学生姓名：XXX学 号：专 业：电气工程及其自动化班 级：电气2011 1班指导教师：XXX摘要 1Abstract 2内蒙古科技大学课程设计任务书 1概述- 3-51.1待设计电气设备系统图51.2电力系统概述51.3潮流计算的概述61.4电力系统稳态分析潮流计算总结 71.5 MATLAB的相关知识与学习 - 82程序设计142.1设计思路 142.2设计流程思路 152.3设计程序 163设计结果204总结-38参考文献基于MATLAB的电力系统复杂潮流分析摘要电力系统潮流计算是电
2、网分析的基础应用，是对复杂电力系统正常和故障条 件下稳态运行状态的计算。给定电力系统的网络结构、参数和决定电力系统运行 状况的边界条件，确定电力系统运行的方法之一是朝流计算。牛顿-拉夫逊法（简称牛顿法）在数学上是求解非线性代数方程式的有效方 法。其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求 解的过程。即通常所称的逐次线性化过程。MATLAB是一种交互式、面向对象的程序设计语言，广泛应用于工业界与 学术界，主要用于矩阵运算.采用迭代法，通过建立矩阵的修正方程来依次迭代， 逐步逼近真值来计算出电力网的电压，功率分布。使用MATLAB软件进行编程，在很大程度上节省了内存，减少了
3、计算量。 通过对本题计算我们了解了一些工程计算和解决工程问题的方法。关键词：潮流计算，稳态运行，牛顿 一拉夫逊法，MATLABBased on MATLAB electrical power system complex tidalcurre nt an alysisAbstractPower Flow Analysis Grid computing is the basis of applications, the complex power system under normal and fault conditions for the calculation of steady stat
4、e operation. Given the power system network structure, parameters and decisions operation of the power system boundary conditions, to determine the method of operati on of the power system is one of North Korea flow calculati on.PQ decompositi on method is the form of polar coord in ates Newt on - t
5、he widening trend of a simplified calculation method. P-Q decomposition method adopted by the specific characteristics of the power system ariysis, Newt on s Law ofthe Jacobian matrix formula has effectively simplified and improved. As a result of these simplified formula that invoIves only the coef
6、ficient matrix, the balance of power has not cha nged node equati ons and the conv erge nce criteri on, becausethe results will not reduce the accuracy.MATLAB is an in teractive, object-orie nted program ming Ian guage, widely used in industry and academia, mainly for matrix calculation. Using itera
7、tion, the ame ndme nt through the establishme nt of matrix iterative equatio n to tur n, gradually moving towards a true value to calculate the voltage electricity grid, power distributio n.Use MATLAB software program ming, sav ing memory to a large exte nt, reduce the amount of computation. By this
8、 calculation we understand that a number of engin eeri ng calculatio n and solve engin eeri ng problemsKeywords : The trend, the PQ decomposition method, MATLAB内蒙古科技大学课程设计任务书课程名称电力系统稳态分析课程设计设计题目两机五节点网络潮流计算一PQ节点法指导教师刘景霞 时间1周一、 教学要求电力系统稳态分析课程设计以设计和优化电力系统的潮流分析为重点， 提高学生综合能力为目标，尽可能结合实际工程进行。设计内容的安排要充分考虑学校
9、现有的设备，设计时 间及工程实际需要，并使学生初步学会运用所学知识解决工程中的实际问题。二、 设计资料及参数（一） 设计原始资料1、 待设计电气设备系统图2、 电力系统网络各元件参数3、 电力系统电气元件的使用规范4、 电力工程电气设计手册（二） 设计参考资料1、 电力系统稳态分析，陈珩，中国电力出版社，2007,第三版2、 电力系统分析，韩祯祥，浙江大学出版社，2005,第三版3、 电力系统分析课程实际设计与综合实验，祝书萍，中国电力出版社，2007，第一版三、 设计要求及成果1.根据给定的参数或工程具体要求，收集和查阅资料；学习相关软件（软件自选） 。2.在给定的电力网络上画出等值电路图。
10、3.运用计算机进行潮流计算。4.编写设计说明书。基本要求：1.编写潮流计算程序；2.在计算机上调试通过（？）；3.运行程序并计算出正确结果（？）；4.写出课程设计报告（包括以下内容）（1份）（1）程序框图；（2）源程序；（3）符号说明表；（4）算例及计算结果5.编写计算说明书（1份）。四、 进度安排根据给定的参数或工程具体要求，收集和查阅资料（半天）学习软件（MATLA或C语言等）（一天半）编程计算复杂系统潮流计算（三天）编写计算设计书（一天）五、 评分标准课程设计成绩评定依据包括以下几点：1） 工作态度（占10%；2） 基本技能的掌握程度（占20%；3） 程序编写是否合理是否有运行结果（40
11、%）;4） 课程设计说明书编写水平（占30%）=5） 分为优、良、中、合格、不合格五个等级。考核方式：设计期间教师现场检查；评阅设计报告1概述1.1待设计电气设备系统图0.45+j0.150.4+j0.0550.04+j0.12；：.G 1(0.2+j0.2)0.6+j0.1系统接线图（其中节点1为平衡节点，节点2、3、4、5为PQ节点。）1.2 电力系统概述电力工业发展初期，电能是直接在用户附近的发电站 （或称发电厂）中生产的，各发电站孤立运行。随着工农业生产和城市的发展，电能的需要量迅速增加， 而热能资源（如煤田）和水能资源丰富的地区又往往远离用电比较集中的城市和 工矿区，为了解决这个矛盾
12、，就需要在动力资源丰富的地区建立大型发电站， 然 后将电能远距离输送给电力用户。同时，为了提高供电可靠性以及资源利用的综 合经济性，又把许多分散的各种形式的发电站，通过送电线路和变电所联系起来。 这种由发电机、升压和降压变电所，送电线路以及用电设备有机连接起来的整体， 即称为电力系统。电力系统加上发电机的原动机（如汽轮机、水轮机），原动机的力能部分（如热 力锅炉、水库、原子能电站的反应堆）、供热和用热设备，则称为动力系统。现代电力系统提出了 “灵活交流输电与新型直流输电” 的概念。灵活交流输 电技术是指运用固态电子器件与现代自动控制技术对交流电网的电压、相位角、 阻抗、功率以及电路的通断进行实
13、时闭环控制， 从而提高高压输电线路的输送能 力和电力系统的稳定水平。新型直流输电技术是指应用现电力电子技术的最新成 果，改善和简化变流站的造价等。运行方式管理中，潮流是确定电网运行方式的基本出发点； 在规划领域，需 要进行潮流分析验证规划方案的合理性； 在实时运行环境，调度员潮流提供了电 网在预想操作情况下电网的潮流分布以校验运行可靠性。 在电力系统调度运行的多个领域都涉及到电网潮流计算。潮流是确定电力网络运行状态的基本因素，潮 流问题是研究电力系统稳态问题的基础和前提。1.3潮流计算的概述电力系统潮流计算是研究电力系统稳态运行情况的一种计算， 它根据给定的运行条件及系统接线情况确定整个电力系
14、统各部分的运行状态：各母线的电压， 各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力 系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行 方式的合理性。可靠性和经济性。此外，电力系统潮流计算也是计算系统动态稳 定和静态稳定的基础。所以潮流计算是研究电力系统的一种很重要和基础的计 算。电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统 规划设计和安排系统的运行方式，后者则用于正在运行系统的经常监视及实时控 制。利用电子数字计算机进行电力系统潮流计算从 50年代中期就已经开 始。在这20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围
15、 绕着对潮流计算的一些基本要求进行的。对潮流计算的要求可以归纳为下面几 占：八、(1)计算方法的可靠性或收敛性；(2)对计算机内存量的要求；(3)计算速度；(4)计算的方便性和灵活性。电力系统潮流基本步骤1、 导纳矩阵的形成2、 节点类型(1)PV点：注入有功功率P为给定值，电压也保持在给定数值。(2)PQ节点：注入有功功率和无功功率是给定的(3)平衡节点：用来平衡全电网的功率。选一容量足够大的发电机担任平衡 全电网功率的职责。平衡节点的电压大小与相位是给定的， 通常以它的相角为参量，即取其电压 相角为0个独立的电力网中只设一个平衡点。3、 高斯迭代法电力系统潮流计算问题在数学上是一组多元非线
16、性方程式求解问题， 其解法都离不开迭代。因此，对潮流计算方法，首先要求它能可靠地收敛，并给出正确答案。由于电力系统结构及参数的一些特点，并且随着电力系统不断扩大， 潮流计算的方程式阶数也越来越高，对这样的方程式并不是任何数学方法都能保 证给出正确答案的。这种情况成为促使电力系统计算人员不断寻求新的更可靠方 法的重要因素。1.4电力系统稳态分析潮流计算总结潮流计算针对电力革统各正常运行方式，而静态安全分析则要研究各 种运行方式下个别系统元件退出运行后系统的状况。其目的是校验系统是 否能安全运行，即是否有过负荷的元件或电压过低的母线等。原则上讲，静态安全分析也可 U用潮流计算来代替。但是一般静态安
17、全分析需要校验的状态数非常多，用严格的潮流计算来分析这些状态往往计算量过大，因 此不得不寻求一些特殊的算法以满足要求。利用电子数字计算机进行电力系统潮流计算从 20世纪50年代中期就己开始，此后，潮流计算曾采用了各种不同的方法，这些方法的发展主要 是围绕着对潮流计算的一些基本要求进行的，对潮流计算的要求可以归纳 为下面几点：(1)计算方法的可靠性或收敛性(2)对计算速度和内存量的要求(3)计算的方便性和灵活性现在主要的研究都是围绕着牛顿拉夫逊法来进行的，牛顿法是数学中解决非线性方程式的典型方法，有较好的收敛性。解决电力系统潮流计算 问题是以导纳距阵为基础的，因此，只要在迭代过程中尽可能保持方程
18、式 系数距阵的稀疏性，就可以大大提高牛顿法潮流程序的放率。自从 20世纪60年代中期利用了最佳顺序消去法以后，牛顿法在收敛性、内存要求、速 度方面都超过了阻抗法，成为直到目前仍在广泛采用的优秀方法。同样优 秀的方法还有很多，fast decoupled, PQ 分解法，等等。1.5 MATLAB的相关知识与学习MATLAB言及其应用一门基础工具课。MATLAB! MathWorks公司开发的一套高性能的数值计算和可视化的仿真软件， 包括MATLA主程序、SIMULINK动态系统仿真包和各种专业工具箱。它基于矩阵运算，集数值计算、图形显示、图形 界面设计等功能和模糊逻辑等 30多个专业工具箱于一
19、体，具有极强大的计算功 能和极高的编程效率，特别适合于科学计算、数值分析、系统仿真和信号处理等 任务。MATLAB! 一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界， 主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动态分析、 绘图等方面也有强大的功能。MATLAE程学设计语言结构完整，且具有优良的移植性，它的基本数据元素 是不需要定义的数组。它可以高效率的解决工业计算问题，特别是关于矩阵和矢 量的计算。通过MATLAB!言，可以用类似数学公式的方式来编写算法，大大降 低了程序需要的难度别难并节省了时间， 从而可把主要的经历集中在算法的构思 而不是编程上。表3-1管
20、理命令和函数help在线帮助文件doc装入超文本说明whatM MAT MEX文件的目录列表type列出M文件lookfor通过help条目搜索关键字which定位函数和文件Demo运行演示程序Path控制MATLA的搜索路径管理变量和工作空间Who列出当前变量Whos列出当前变量（长表）Load从磁盘文件中恢复变量Save保存工作空间变量Clear从内存中清除变量和函数Pack整理工作空间内存Size矩阵的尺寸Len gth向量的长度disp显示矩阵或与文件和操作系统有关的命令cd改变当前工作目录Dir目录列表Delete删除文件Gete nv获取环境变量值!执行DOS操作系统命令Unix执
21、行UNIX操作系统命令并返回结果Diary保存MATLA任务控制命令窗口Cedit设置命令行编辑Clc清命令窗口Home光标置左上角Format设置输出格式Echo底稿文件内使用的回显命令more在命令窗口中控制分页输出启动和退出MATLABQuit退出MATLABStartup引用MATLA时所执行的M文件Matlabrc主启动M文件般信息InfoMATLA系统信息及Mathworks公司信息Subscribe成为MATLA的订购用户hostidMATLA主服务程序的识别代号Whats new在说明书中未包含的新信息Ver版本信息操作符和特殊字符+加减*矩阵乘法*数组乘法A矩阵幕A数组幕左除
22、或反斜杠/右除或斜杠./数组除KronKronecker张量积:冒号()圆括号方括号.小数点父目录继续J逗号J分号%注释!感叹号转置或引用=赋值=相等关系操作符&逻辑与|逻辑或逻辑非xor逻辑异或逻辑函数Exist检查变量或函数是否存在Any向量的任一元为真，则其值为真All向量的所有元为真，则其值为真Fi nd找出非零元素的索引号三角函数Sin正弦Si nh双曲正弦Asi n反正弦As in h0.00001;dt(1)=0;u(1)=1.06;for m=2:5for n=1:5 pt(n )=u(m)*u( n)*(G(m, n)*cos(dt(m)-dt( n)+B(m, n)*si n(dt(m)-dt(n);enddisp(各节点有功率不平衡量为：)dp(m)=p(m)-sum(pt)endfor m=2:5u仁diag(u);u2=i nv(u1(2:5,2:5);u3=flipud(rot90(u(2:5);u4=u3(1:4,:);dp仁flipud(rot90(dp(2:5);dp2=
