《[工学]电力系统稳态分析课程设计.doc》由会员分享，可在线阅读，更多相关《[工学]电力系统稳态分析课程设计.doc（40页珍藏版）》请在三一文库上搜索。
1、内蒙古科技大学电力系统稳态分析课程设计说明书 题 目：基于Matlab的两机五节点 网络潮流仿真计算牛拉法 姓 名： 学 号： 学 院：信息工程学院 专 业：电气工程及其自动化 班 级：10 电气 目录 目录II摘 要III内蒙古科技大学课程设计任务书IV第一章 电力系统潮流计算及MATLAB概述- 1 -第一节 潮流计算简介- 1 -第二节 MATLAB概述- 1 -第三节 矩阵的简单基本运算- 2 -第二章 潮流计算的数学模型- 3 -第一节 导纳矩阵的原理及计算方法- 3 -第二节 电力系统节点分类- 8 -第三节 潮流计算的约束条件- 9 -第三章 牛顿拉夫逊法概述- 10 -第一节
2、牛顿-拉夫逊法- 10 -第二节 用牛顿-拉夫逊法求解潮流计算的过程- 12 -第四章 牛顿拉夫逊法潮流具体计算- 17 -第一节 按照设计要求先列出手算过程- 17 -第五章 牛顿拉夫逊直角坐标潮流计算Matlab程序及运行结果- 23 -第一节 程序框图- 23 -第一节 Matlab程序- 24 -第二节 Matlab程序运行结果- 30 -第四节 本程序的符号说明- 34 -总结及感想- 35 -参考文献及资料- 36 -摘 要潮流计算的目的在于：确定电力系统的运行方式；检查系统中各元件是否过电压或过载；为电力系统继电保护的整定提供依据；为电力系统的稳定计算提供初值，为电力系统规划和经
3、济运行提供分析基础。牛顿迭代法（Newtons method）又称为牛顿-拉夫逊方法（Newton-Raphson method），牛顿-拉夫逊法(简称牛顿法)在数学上是求解非线性代数方程式的有效方法。其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程。即通常所称的逐次线性化过程。MATLAB是一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界，主要用于矩阵运算采用迭代法，通过建立矩阵的修正方程来依次迭代，逐步逼近真值来计算出电力网的电压，功率分布。本文采用牛顿-拉夫逊法解算电力稳态潮流，用手算和计算机算法对其进行设计。使用MATLAB软件进行编程，在很大程度
4、上节省了内存，减少了计算量。通过对本题计算我们了解了一些工程计算和解决工程问题的方法。 关键词： 电力系统潮流计算；牛顿拉夫逊法潮流计算；潮流计算的数学模型 MATLAB矩阵运算；程序； 内蒙古科技大学课程设计任务书课程名称电力系统稳态分析课程设计设计题目基于Matlab的两机五节点网络潮流仿真计算牛拉法指导教师刘景霞时间1周一、教学要求 电力系统稳态分析课程设计以设计和优化电力系统的潮流分析为重点，提高学生综合能力为目标，尽可能结合实际工程进行。设计内容的安排要充分考虑学校现有的设备，设计时间及工程实际需要，并使学生初步学会运用所学知识解决工程中的实际问题。二、设计资料及参数（一）设计原始资
5、料1、待设计电气设备系统图2、电力系统网络各元件参数3、电力系统电气元件的使用规范4、电力工程电气设计手册（二）设计参考资料1、电力系统稳态分析，陈珩，中国电力出版社，2007，第三版2、电力系统分析，韩祯祥，浙江大学出版社，2005，第三版3、电力系统分析课程实际设计与综合实验，祝书萍，中国电力出版社，2007，第一版三、设计要求及成果1.根据给定的参数或工程具体要求，收集和查阅资料；学习相关软件（软件自选）。2.在给定的电力网络上画出等值电路图。3.运用计算机进行潮流计算。4.编写设计说明书。基本要求：1.编写潮流计算程序；2.在计算机上调试通过（？）；3.运行程序并计算出正确结果（？）；
6、4.写出课程设计报告（包括以下内容）（1份）（1）程序框图；（2）源程序；（3）符号说明表；（4）算例及计算结果5.编写计算说明书（1份）。四、进度安排根据给定的参数或工程具体要求，收集和查阅资料（半天）学习软件（MATLAB或C语言等）（一天半）编程计算复杂系统潮流计算（三天）编写计算设计书（一天）五、评分标准课程设计成绩评定依据包括以下几点：1） 工作态度（占10%）；2） 基本技能的掌握程度（占20%）；3） 程序编写是否合理是否有运行结果(40%)；4） 课程设计说明书编写水平(占30%)。5） 分为优、良、中、合格、不合格五个等级。考核方式：设计期间教师现场检查；评阅设计报告系统接线
7、图其中节点1为平衡节点，节点2、3、4、5为PQ节点。VI 第一章 电力系统潮流计算及MATLAB概述第一节 潮流计算简介潮流计算是电力系统分析中的一种最基本的计算，它的任务是在给定的接线方式和运行条件下，确定系统的运行状态，如各母线上的电压（幅值和相角）、网络中的功率分布及功率损耗等，是电力系统的稳态计算。潮流计算是对电力系统正常运行状况的分析和计算，即电力系统中的电压、电流、功率的计算，即潮流计算；潮流计算方法很多：高斯塞德尔法、牛顿拉夫逊法、P-Q分解法、直流潮流法，以及由高斯塞德尔法、牛顿拉夫逊法演变的各种潮流计算方法。第二节 MATLAB概述目前电子计算机已广泛应用于电力系统的分析计
8、算，潮流计算是其基本应用软件之一。现有很多潮流计算方法。对潮流计算方法有五方面的要求：（1）计算速度快（2）内存需要少（3）计算结果有良好的可靠性和可信性（4）适应性好，亦即能处理变压器变比调整、系统元件的不同描述和与其它程序配合的能力强（5）简单实用。 MATLAB是一种交互式、面向对象的程序设计语言，广泛应用于工业界与学术界，主要用于矩阵运算，同时在数值分析、自动控制模拟、数字信号处理、动态分析、绘图等方面也具有强大的功能。 MATLAB程序设计语言结构完整，且具有优良的移植性，它的基本数据元素是不需要定义的数组。它可以高效率地解决工业计算问题，特别是关于矩阵和矢量的计算。MATLAB与C
9、语言和FORTRAN语言相比更容易被掌握。通过M语言，可以用类似数学公式的方式来编写算法，大大降低了程序所需的难度并节省了时间，从而可把主要的精力集中在算法的构思而不是编程上。 另外，MATLAB提供了一种特殊的工具：工具箱（TOOLBOXES）.这些工具箱主要包括：信号处理（SIGNAL PROCESSING）、控制系统（CONTROL SYSTEMS）、神经网络（NEURAL NETWORKS）、模糊逻辑(FUZZY LOGIC)、小波(WAVELETS)和模拟（SIMULATION）等等。不同领域、不同层次的用户通过相应工具的学习和应用，可以方便地进行计算、分析及设计工作。 MATLAB
10、设计中，原始数据的填写格式是很关键的一个环节，它与程序使用的方便性和灵活性有着直接的关系。原始数据输入格式的设计，主要应从使用的角度出发，原则是简单明了，便于修改。第三节 矩阵的简单基本运算矩阵是MATLAB数据存储的基本单元，而矩阵的运算是MATLAB语言的核心，在MATLAB语言系统中几乎一切运算均是以对矩阵的操作为基础的。矩阵的基本数学运算包括矩阵的四则运算、与常数的运算、逆运算、行列式运算、秩运算、特征值运算等基本函数运算，这里进行简单介绍。矩阵的加、减、乘运算符分别为“+，*” ，用法与数字运算几乎相同，但计算时要满足其数学要求 在MATLAB中矩阵的除法有两种形式：左除“”和右除“
11、/”。矩阵的转置，矩阵的求逆在MATLAB应用非常广泛。如：位阵的生成 eye(n) 、eye(m,n) ；零矩阵的生成 zeros(n) zeros(m,n) ；全1矩阵的生成 ones(n) ones(m,n) ；随机元素矩阵函数rand(n,m) 、rand(n) ；对角矩阵diag(V) % V=1 2 3 4; 伴随矩阵compan(P) % p=1,a1,a2,.,an ；上三角矩阵 下三角矩阵triu(B) tril(B) 矩阵的行列式det(A) ；矩阵求逆inv(A) ；矩阵的秩rank(A) ；矩阵三角分解L,U = lu(A) 矩阵奇异值分解cond(A) ；矩阵的范数N
12、 = norm(A,选项) ；矩阵的特征多项值与特征向量V,D=eig(A) ；矩阵的特征多项式、特征方程和特征根P = poly(A) 、V = roots(P) 第二章 潮流计算的数学模型第一节 导纳矩阵的原理及计算方法一 自导纳和互导纳的确定方法电力网络的节点电压方程： (2-1)为节点注入电流列向量，注入电流有正有负，注入网络的电流为正，流出网络的电流为负。根据这一规定，电源节点的注入电流为正，负荷节点为负。既无电源又无负荷的联络节点为零，带有地方负荷的电源节点为二者代数之和。为节点电压列向量，由于节点电压是对称于参考节点而言的，因而需先选定参考节点。在电力系统中一般以地为参考节点。如
13、整个网络无接地支路，则需要选定某一节点为参考。设网络中节点数为（不含参考节点），则，均为n*n列向量。为n*n阶节点导纳矩阵。节电导纳矩阵的节点电压方程： 展开为： ： (2-2)是一个n*n阶节点导纳矩阵，其阶数就等于网络中除参考节点外的节点数。 节点导纳矩阵的对角元素 (i=1，2，n)成为自导纳。自导纳数值上就等于在i节点施加单位电压，其他节点全部接地时，经节点i注入网络的电流，因此，它可以定义为： (2-3)节点i的自导纳数值上就等于与节点直接连接的所有支路导纳的总和。节点导纳矩阵的非对角元素 (j=1，2，n;i=1，2，。，n;j=i)称互导纳，由此可得互导纳数值上就等于在节点i施
14、加单位电压，其他节点全部接地时，经节点j注入网络的电流，因此可定义为： (2-4) 节点j，i之间的互导纳数值上就等于连接节点j，i支路到导纳的负值。显然，恒等于。互导纳的这些性质决定了节点导纳矩阵是一个对称稀疏矩阵。而且，由于每个节点所连接的支路数总有一个限度，随着网络中节点数的增加非零元素相对愈来愈少，节点导纳矩阵的稀疏度，即零元素数与总元素的比值就愈来愈高。二非标准变比变压器等值电路变压器型等值电路更便于计算机反复计算,更适宜于复杂网络的潮流计算.双绕组变压器可用阻抗与一个理想变压器串联的电路表示.理想变压器只是一个参数,那就是变比。现在变压器阻抗按实际变比归算到低压侧为例,推导出变压器
15、型等值电路. a 双绕组变压器原理图b变压器阻抗归算到低压侧等值模型流入和流出理想变压器的功率相等 (2-6)式中, 是理想变压器的变比,和 分别为变压器高,低绕组的实际电压.从图b直接可得： （2-7） 从而可得: （2-8）式中,又因节点电流方程应具有如下形式: （2-9）将式（1-8）与（1-9）比较，得： 因此可得各支路导纳为: （2-10）由此可得用导纳表示的变压器型等值电路:图 c 潮流计算的基本方程(a) 潮流计算用的电网结构图 在潮流问题中，任何复杂的电力系统都可以归纳为以下元件（参数）组成。（1）发电机（注入电流或功率）（2）负荷（注入负的电流或功率）（3）输电线支路（电阻，
16、电抗）（4）变压器支路（电阻，电抗，变比）（5）母线上的对地支路（阻抗和导纳）（6）线路上的对地支路（一般为线路充电点容导纳）集中了以上各类型的元件的简单网络如图(a).采用导纳矩阵时，节点注入电流和节点电压构成以下线性方程组 (2-11)其中 可展开如下形式 (2-12)由于实际电网中测量的节点注入量一般不是电流而是功率，因此必须将式中的注入电流用节点注入功率来表示。节点功率与节点电流之间的关系为 (2-13)式中，因此用导纳矩阵时，PQ节点可以表示为把这个关系代入式中 ，得（2-14）式（3-4 ）就是电力系统潮流计算的数学模型-潮流方程。它具有如下特点：（1）它是一组代数方程，因而表征的
17、是电力系统的稳定运行特性。（2）它是一组非线性方程，因而只能用迭代方法求其数值解。（3）由于方程中的电压和导纳既可以表为直角坐标，又可表为极坐标，因而潮流方程有多种表达形式-极坐标形式，直角坐标形式和混合坐标形式。 a。取 ，得到潮流方程的极坐标形式： (2-15)b。 取 ， ，得到潮流方程的直角坐标形式： (2-16)c。取， ，得到潮流方程的混合坐标形式： (2-17)（4）它是一组n个复数方程，因而实数方程数为2n个但方程中共含4n个变量：P，Q，U和，i=1，2，n，故必须先指定2n个变量才能求解。第二节 电力系统节点分类用一般的电路理论求解网络方程，目的是给出电压源(或电流源)研究
18、网络内的电流(或电压)分布，作为基础的方程式，一般用线性代数方程式表示。然而在电力系统中，给出发电机或负荷连接母线上电压或电流(都是向量)的情况是很少的，一般是给出发电机母线上发电机的有功功率(P)和母线电压的幅值(U)，给出负荷母线上负荷消耗的有功功率(P)和无功功率(Q)。主要目的是由这些已知量去求电力系统内的各种电气量。所以，根据电力系统中各节点性质的不同，很自然地把节点分成三类： PQ节点对这一类点，事先给定的是节点功率(P，Q)，待求的未知量是节点电压向量(U，)，所以叫PQ节点。通常变电所母线都是PQ节点，当某些发电机的输出功率P。Q给定时，也作为PQ节点。PQ节点上的发电机称之为
19、PQ机(或PQ给定型发电机)。在潮流计算中，系统大部分节点属于PQ节点。 PV节点这类节点给出的参数是该节点的有功功率P及电压幅值U，待求量为该节点的无功功率Q及电压向量的相角。这类节点在运行中往往要有一定可调节的无功电源。用以维持给定的电压值。通常选择有一定无功功率储备的发电机母线或者变电所有无功补偿设备的母线做PU节点处理。PU节点上的发电机称为PU机(或PV给定型发电机) 平衡节点在潮流计算中，这类节点一般只设一个。对该节点，给定其电压值，并在计算中取该节点电压向量的方向作为参考轴，相当于给定该点电压向量的角度为零。也就是说，对平衡节点给定的运行参数是U和，因此有城为U节点，而待求量是该
20、节点的P。Q，整个系统的功率平衡由这一节点承担。关于平衡节点的选择，一般选择系统中担任调频调压的某一发电厂(或发电机)，有时也可能按其他原则选择，例如，为提高计算的收敛性。可以选择出线数多或者靠近电网中心的发电厂母线作平衡节点。以上三类节点4个运行参数P、Q、U、中，已知量都是两个，待求量也是两个，只是类型不同而已。第三节 潮流计算的约束条件 电力系统运行必须满足一定技术和经济上的要求。这些要求够成了潮流问题中某些变量的约束条件，常用的约束条件如下：1. 节点电压应满足 (2-12)2.从保证电能质量和供电安全的要求来看，电力系统的所有电气设备都必须运行在额定电压附近。PU节点电压幅值必须按上
21、述条件给定。因此，这一约束条件对PQ节点而言。2. 节点的有功功率和无功功率应满足3. (2-13)PQ节点的有功功率和无功功率，以及PU节点的有功功率，在给定是就必须满足上述条件，因此，对平衡节点的P和Q以及PU节点的Q应按上述条件进行检验。4. 节点之间电压的相位差应满足 (2-14) 为了保证系统运行的稳定性，要求某些输电线路两端的电压相位不超过一定的数值。这一约束的主要意义就在于此。 因此，潮流计算可以归结为求解一组非线性方程组，并使其解答满足一定的约束条件。常用的方法是迭代法和牛顿法，在计算过程中，或得出结果之后用约束条件进行检验。如果不能满足要求，则应修改某些变量的给定值，甚至修改
22、系统的运行方式，重新进行计算。 第三章 牛顿拉夫逊法概述第一节 牛顿-拉夫逊法电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采用牛顿-拉夫逊法。牛顿
23、-拉夫逊法(简称牛顿法)在数学上是求解非线性代数方程式的有效方法。其要点是把非线性方程式的求解过程变成反复地对相应的线性方程式进行求解的过程。即通常所称的逐次线性化过程。对于非线性代数方程组： 即 (3-1)在待求量x的某一个初始估计值附近，将上式展开成泰勒级数并略去二阶及以上的高阶项，得到如下的经线性化的方程组： (3-2)上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量 (3-3)将和相加，得到变量的第一次改进值。接着就从出发，重复上述计算过程。因此从一定的初值出发，应用牛顿法求解的迭代格式为： (3-4) (3-5)上两式中：是函数对于变量x的一阶偏导数矩阵，即雅可比矩阵J;
24、k为迭代次数。有上式可见，牛顿法的核心便是反复形式并求解修正方程式。牛顿法当初始估计值和方程的精确解足够接近时，收敛速度非常快，具有平方收敛特性。牛顿潮流算法突出的优点是收敛速度快，若选择到一个较好的初值，算法将具有平方收敛特性，一般迭代45次便可以收敛到一个非常精确的解。而且其迭代次数与所计算网络的规模基本无关。牛顿法也具有良好的收敛可靠性，对于对以节点导纳矩阵为基础的高斯法呈病态的系统，牛顿法也能可靠收敛。牛顿法所需的内存量及每次迭代所需时间均较高斯法多。 牛顿法的可靠收敛取决于有一个良好的启动初值。如果初值选择不当，算法有可能根本不收敛或收敛到一个无法运行的节点上。对于正常运行的系统，各
25、节点电压一般均在额定值附近，偏移不会太大，并且各节点间的相位角差也不大，所以对各节点可以采用统一的电压初值(也称为平直电压)，如假定： 或 (3-6) 这样一般能得到满意的结果。但若系统因无功紧张或其它原因导致电压质量很差或有重载线路而节点间角差很大时，仍用上述初始电压就有可能出现问题。解决这个问题的办法可以用高斯法迭代12次，以此迭代结果作为牛顿法的初值。也可以先用直流法潮流求解一次以求得一个较好的角度初值，然后转入牛顿法迭代。第二节 用牛顿-拉夫逊法求解潮流计算的过程以下讨论的是用直角坐标形式的牛顿拉夫逊法潮流的求解过程。当采用直角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个分量由
26、于平衡节点的电压向量是给定的，因此待求共需要2(n-1)个方程式。事实上，除了平衡节点的功率方程式在迭代过程中没有约束作用以外，其余每个节点都可以列出两个方程式。 （3-7）对PQ节点来说，是给定的，因而可以写出 （3-8） （3-9）求解过程大致可以分为以下步骤：（1）形成节点导纳矩阵（2）将各节点电压设初值U，（3）将节点初值代入相关求式，求出修正方程式的常数项向量（4）将节点电压初值代入求式，求出雅可比矩阵元素（5）求解修正方程，求修正向量（6）求取节点电压的新值（7）检查是否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步（8）计算支路功率分布
27、，PV节点无功功率和平衡节点柱入功率。以直角坐标系形式表示. 迭代推算式 采用直角坐标时,节点电压相量及复数导纳可表示为: (3-10)将以上二关系式代入上式中,展开并分开实部和虚部;假定系统中的第1,2,m号为PQ节点,第m+1,m+2,n-1为PV节点,根据节点性质的不同,得到如下迭代推算式: 对于PQ节点 (3-11)对于PV节点 (3-12)对于平衡节点 平衡节点只设一个,电压为已知,不参见迭代,其电压为: (3-13). 修正方程选定电压初值及变量修正量符号之后代入式中,并将其按泰勒级数展开,略去二次方程及以后各项,得到修正方程如下: (3-14) .雅可比矩阵各元素的算式式(3-7
28、)中, 雅可比矩阵中的各元素可通过对式(3-11)和(3-12)进行偏导而求得.当时, 雅可比矩阵中非对角元素为 (3-15)当时,雅可比矩阵中对角元素为: (3-16) 第四章 牛顿拉夫逊法潮流具体计算 第一节 按照设计要求先列出手算过程手算过程：解：该系统中，节点1平衡节点，保持U1=1.06+j0.00为定值；其他节点为PQ节点，给定的注入功率分别为：， ，1.形成节点导纳矩阵Y= 2.计算各节点功率的不平衡量 取U1=1.06+j0.00计算各点功率得1.0*（-5.00*1.06-15.00*0.0）+0.0*（-5.00*0.0+15.00*1.06）+1.0*（10.834*1.
29、0+32.50*0.0）+0.0*（10.834*0.0-32.500*1.0）+1*（-1.667*1-5*0）+0*（-1.667*0+5*1）+1*（-1.667*1-5*0）+0*（-1.667*0+5*1）+1*（-2.5*1-7.5*0）+0*（-2.5*0-7.5*1） =-0.300*（-5*1.06-15*0）-15*0）-1*（-5*0+15*1.06）+0*(10.834*1+32.5*0)_1*(10.834*0-32.5*1)+0*（-1.667*1-5*1）-1*(-1.667*0+5*1）+0*（-1.667*1-5*1）-1*(-1.667*0+5*1）+0*(
30、-2.5*1-7.5*0)-1*（-2.5*1+7.5*1） =-0.90同理可得：-0.075，0，0 所以： 由得： 同理： 3.计算雅克比矩阵中各元素先计算各节点注入电流：由得： 同理： 然后计算雅各比矩阵各元素同理可计算雅各比矩阵中其他各元素如上 :雅可比矩阵J=4.解修正方程式各节点电压得第一次：，求得各节点的新值后，就可开始第二次迭代。每次迭代结果如下第二次：，第三次：，第四次：，经过四次迭代可满足的要求。5. 计算平衡节点功率和线路功率平衡节点功率 = -2.6697 - 8.0092i线路功率, ,， , ,， , ，上式没写出的均为0第五章 牛顿拉夫逊直角坐标潮流计算Matl
31、ab程序及运行结果第一节 程序框图第一节 Matlab程序GG0.08+j0.240.01+j0.030.06+j0.180.06+j0.180.04+j0.120.08+j0.240.02+j0.060.45+j0.150.4+j0.050.6+j0.1-(0.2+j0.2)51234clear;clc%重新编号，把原题中的节点1,2,3,4,5重新依次编号为5,1,2,3,4，其中1-4号为PQ节点，5号为平衡节点y=0;%输入原始数据，求节点导纳矩阵y (1,2)=1/(0.06+0.18i); y (1,3)=1/(0.06+0.18i); y (1,4)=1/(0.04+0.12i)
32、;y(1,5)=1/(0.02+0.06i);y(2,3)=1/(0.01+0.03i);y(2,5)=1/(0.08+0.24i);y(3,4)=1/(0.08+0.24i);y(4,5)=0;for i=1:5for j=i:5y(j,i)=y(i,j);endendY=0;%求互导纳for i=1:5for j=1:5if i=jY(i,j)=-y(i,j);endendend%求自导纳for i=1:5Y(i,i)=sum(y(i,:);endY %Y 为导纳矩阵G=real(Y);B=imag(Y);%原始节点功率S(1)=0.2+0.2i;S(2)=-0.45-0.15i;S(3)
33、=-0.4-0.05i;S(4)=-0.6-0.1i;S(5)=0;P=real(S);Q=imag(S);%赋初值U=ones(1,5);U(5)=1.06;e=zeros(1,5);ox=ones(8,1);fx=ones(8,1);count=0 %计算迭代次数while max(fx)1e-5for i=1:4 for j=1:4 H(i,j)=0;N(i,j)=0;M(i,j)=0;L(i,j)=0;oP(i)=0;oQ(i)=0; endendfor i=1:4for j=1:5 oP(i)=oP(i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j)+B(i,j)*
34、sin(e(i)-e(j);oQ(i)=oQ(i)-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j)-B(i,j)*cos(e(i)-e(j);endoP(i)=oP(i)+P(i); oQ(i)=oQ(i)+Q(i);endfx=oP,oQ;%求雅克比矩阵%当i=j时候求H,N,M,L 如下：for i=1:4for j=1:4if i=j H(i,j)=-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j)-B(i,j)*cos(e(i)-e(j);N(i,j)=-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j
35、);L(i,j)=H(i,j);M(i,j)=-N(i,j);endendendH,N,M,L%当i=j 时H,N,M,L如下：for i=1:4for j=1:5if i=jH(i,i)=H(i,i)+U(i)*U(j)*(G(i,j)*sin(e(i)-e(j)-B(i, j)*cos (e(i)-e(j); N(i,i)=N(i,i)-U(i)*U(j)*(G(i, j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j);M(i,i)=M(i,i)-U(i)*U(j)*(G(i,j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j); L(i,i)
36、=L(i,i)-U(i)*U(j)*(G(i,j)*sin(e(i)-e(j)-B(i,j)*cos(e(i)-e(j);endendN(i,i)=N(i,i)-2*(U(i)2*G(i,i);L(i,i)=L(i,i)+2*(U(i)2*B(i,i);endJ=H,N;M,L %J 为雅克比矩阵ox=-(inv(J)*fx);for i=1:4oe(i)=ox(i); oU(i)=ox(i+4)*U(i);endfor i=1:4e(i)=e(i)+oe(i); U(i)=U(i)+oU(i);endcount=count+1;endox,U,e,count%求节点注入的净功率i=5;for
37、 j=1:5P(i)=U(i)*U(j)*(G(i,j)*cos(e(i)-e(j)+B(i,j)*sin(e(i)-e(j)+P(i);Q(i)=U(i)*U(j)*(G(i,j)*sin(e(i)-e(j)-B(i,j)*cos(e(i)-e(j)+Q(i);endS(5)=P(5)+Q(5)*sqrt(-1);S%求节点注入电流I=Y*U第二节 Matlab程序运行结果(a) 运行结果Y值：迭代过程：电压值：平衡节点注入功率及电流：第四节 本程序的符号说明 本节主要介绍本设计任务书用到得指令。（1）Clc，clear清屏，清除以前显示结果；（2）disp 显示待输出的内容，如disp(节
38、点总数为：)，显示：节点总数为：；（3）while，end；，for，end；循环语句，用于要循环的地方；（4）/ 正除； 反除；如DY=JJDW（5）sum 求和；（6）conj 复数求共轭；（7） 矩阵的转置。（8）% 注释（9）： 自然数循环如1:5（10）abs 取绝对值如abs(DY)（11）max 取最大值如max(A)（12）for- end 循环语句（13）if-end语句 （14）if-else-end总结及感想说起做设计的感觉，起初的时候没有接触过MATLAB，所以编程感觉很吃力，但是通过自己的摸索，上网查看资料和借到同学的我们学院的MATLAB选修课的实验册，慢慢的入门了
39、，最后通过自己的努力，在同学们和老师的帮助下终于编写完程序并计算出符合要求的结果。在此谢谢奋发的同学朋友们和可爱的老师们！ 前一段时间我们还做了发电厂电气部分的课程设计，感觉课程设计真的很锻炼我们的能力，使我们能够认识到自己的不足，同时了解到我们所学知识的实际用途。这次我就深深感觉到自己的能力不足，对MATLAB软件的操作极不熟悉，今后一定要加强自己在动手方面的能力。参考文献及资料1.电力系统稳态分析，陈珩，中国电力出版社，2007，第三版2.信息工程学院基于MATLAB下的信号与系统实验指导书3.Matlab在电气工程中的应用，李维波，中国电力出版社，20074.精通Matlab7M. 朱仁峰，北京：清华大学出版社，20065.Matlab命令大全姚东等，人民邮电出版社，2000，第一版 6.电力系统分析课程实际设计与综合实验，祝书萍，中国电力出版社，第一版- 34 -