1、摘 要潮 流 计 算 是 电 力 系 统 非 常 重 要 的 分 析 计 算 ， 用 以 研 究 系 统 规 划 和 运 行 中提 出 的 各 种 问 题 。 对 规 划 中 的 电 力 系 统 ， 通 过 潮 流 计 算 可 以 检 验 所 提 出 的 电力 系 统 规 划 方 案 能 否 满 足 各 种 运 行 方 式 的 要 求 ； 对 运 行 中 的 电 力 系 统 ， 通 过潮 流 计 算 可 以 预 知 各 种 负 荷 变 化 和 网 络 结 构 的 改 变 会 不 会 危 及 系 统 的 安 全 ，系 统 中 所 有 母 线 的 电 压 是 否 在 允 许 的 范 围 以 内 ，
2、 系 统 中 各 种 元 件 (线 路 、 变压 器 等 )是 否 会 出 现 过 负 荷 ， 以 及 可 能 出 现 过 负 荷 时 应 事 先 采 取 哪 些 预 防 措施 等 。 潮 流 计 算 是 电 力 系 统 分 析 最 基 本 的 计 算 。 除 它 自 身 的 重 要 作 用 之 外 ， 在 电 力 系 统 分 析 综 合 程 序 (PSASP)中 ， 潮 流 计 算 还 是 网 损 计 算 、 静 态 安 全分 析 、 暂 态 稳 定 计 算 、 小 干 扰 静 态 稳 定 计 算 、 短 路 计 算 、 静 态 和 动 态 等 值 计算 的 基 础 。传统的潮流计算程序缺乏
3、图形用户界面，结果显示不直接难与其他分析功能集成。网络原始数据输入工作大量且易于出错。本文采用 MATLAB 语言运行WINDOWS 操作系统的潮流计算软件。而采用 MATLAB 界面直观，运行稳定，计算准确。关键词：电力系统潮流计算；牛顿拉夫逊法潮流计算；MATLAB目 录一、 概述1.1设计目的与要求.1.1.1 设计目的1.1.2 设计要求.1.2 设计题目1.3 设计内容.二、 电力系统潮流计算概述.2.1 电力系统简介2.2 潮流计算简介2.3 潮流计算的意义及其发展. 三、 潮流计算设计题目3.1 潮流计算题目3.2 对课题的分析及求解思路四、 潮流计算算法及手工计算.4.1 变压
4、器的型等值电路4.2 节点电压方程4.3 节点导纳矩阵.4.4 导纳矩阵在潮流计算中的应用.4.5 潮流计算的手工计算五、 Matlab 概述5.1 Matlab 简介5.2 Matlab 的应用5.3 矩阵的运算.5.3.1 与常数的运算.5.3.2 基本数学运算5.3.3 逻辑关系运算5.4 Matlab 中的一些命令.六、 潮流计算流程图及源程序6.1 潮流计算流程图6.2 潮流计算源程序图.6.3 运行计算结果.总结参考文献第一章1.1 设计目的与要求1.1.1 设计目的1. 掌握电力系统潮流计算的基本原理；2. 掌握并能熟练运用一门计算机语言（MATLAB 语言或 C 语言或 C+语
5、言） ；3. 采用计算机语言对潮流计算进行计算机编程。1.1.2 设计要求1. 程序源代码；2. 给定题目的输入，输出文件；3. 程序说明；4. 给定系统的程序计算过程；5. 给定系统的手算过程（至少迭代 2 次） 。1.2设计题目电力系统潮流计算（牛顿-拉夫逊法、P-Q 分解法）1.3设计内容1. 根据电力系统网络推导电力网络数学模型，写出节点导纳矩阵；2. 赋予各节点电压变量（直角坐标系形式）初值后，求解不平衡量；3. 形成雅可比矩阵；4. 求解修正量后，重新修改初值，从 2 开始重新循环计算；5. 求解的电压变量达到所要求的精度时，再计算各支路功率分布、功率损耗和平衡节点功率；6. 上机
6、编程调试；7. 计算分析给定系统潮流分析并与手工计算结果做比较分析；8. 书写课程设计说明书，准备答辩。第二章21 电力系统叙述电力工业发展初期，电能是直接在用户附近的发电站（或称发电厂）中生产的，各发电站孤立运行。随着工农业生产和城市的发展，电能的需要量迅速增加，而热能资源和水能资源丰富的地区又往往远离用电比较集中的城市和工矿区，为了解决这个矛盾，就需要在动力资源丰富的地区建立大型发电站，然后将电能远距离输送给电力用户。同时，为了提高供电的可靠性以及资源利用的综合经济性，又把许多分散的各种形式的发电站，通过送电线路和变电所联系起来。这种由发电机、升压和降压变电所，送电线路以及用电设备有机连接
7、起来的整体，即称为电力系统。现代电力系统提出了“灵活交流输电和新型直流输电”的概念。灵活交流输电技术是指运用固态电子器件与现代自动控制技术对交流电网的电压、相位角、阻抗、功率以及电路的通断进行实时闭环控制，从而提高高压输电线路的诉讼能力和电力系统的稳态水平。新型直流输电技术是指应用现电力电子技术的最新成果，改善和简化变流站的造价等。运营方式管理中，潮流是确定电网运行方式的基本出发点：在规划领域，需要进行潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了电网在预想操作预想下的电网的潮流分布以及校验运行的可靠性。在电力系统调度运行的多个领域都涉及到电网潮流计算。潮流是确定电力网咯运行状
8、态的基本因素，潮流问题是研究电力系统稳态问题的基础和前提。1.2 潮流计算简介电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压。各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量的分析比较供电方案或运行方式的合理性。可靠性和经济性。此外，电力系统的潮流计算也是计算机系统动态稳定和静态稳定的基础，所以潮流计算是研究电力系统的一种和重要和基础的计算。电力系统潮流计算也分为离线计算和在线计算两种，前者主要用于系统规划设计和安排系统的运行方式，后者则用
9、于正在运行系统的经常监视及实时控制。利用电子数字计算机进行潮流计算从 50 年代中期就已经开始了。在这 20年内，潮流计算曾采用了各种不同的方法，这些方法的发展主要围绕着对潮流计算的一些基本要求进行的，对潮流计算的要求可以归纳为以下几点：（1） 计算方法的可靠性或收敛性；（2） 对计算机内存量的要求；（3） 计算速度；（4） 计算的方便性和灵活性。2.3 潮流计算的意义及其发展电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算，即节点电压和功率分布，用以检查系统各元件是否过负荷。各点电压是否满
10、足要求，功率的分布和分配是否合理以及功率损耗等。对现有的电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和稳态分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实际电力系统的潮流技术那主要采用牛顿拉夫逊法。运行方式管理中，潮流是确定电网运行方式的基本出发点；在规划领域，需要进行潮流分析验证规划方案的合理性；在实时运行环境，调度员潮流提供了多个在预想操作情况下电网的潮流分布以及校验运行可靠性。在电力系统调度运行的多个领域问题是研究电力系统稳态问题的基础和前提。在用数字见算计算机解电力系统潮流问题的开始阶段，普
11、遍采取以节点导纳矩阵为基础的逐次代入法。这个方法的原理比较简单，要求的数字计算机内存量比较差下，适应 50 年代电子计算机制造水平和当时电力系统理论水平，但它的收敛性较差，当系统规模变大时，迭代次数急剧上升，在计算中往往出现迭代不收敛的情况。这就迫使电力系统的计算人员转向以阻抗矩阵为基础的逐次代入法。阻抗法改善了系统潮流计算问题的收敛性，解决了导纳无法求解的一些系统的潮流计算，在 60 年代获得了广泛的应用，阻抗法德主要缺点是占用计算机内存大，每次迭代的计算量大。当系统不断扩大时，这些缺点就更加突出，为了克服这些缺点，60 年代中期发展了以阻抗矩阵为基础的分块阻抗法。这个方法把一个大系统分割为
12、几个小的地区系统，在计算机内只需要存储各个地区系统的阻抗矩阵及它们之间联络的阻抗，这样不仅大幅度的节省了内存容量，同时也提高了计算速度。克服阻抗法缺点是另一个途径是采用牛顿-拉夫逊法。这是数学中解决非线性方程式的典型方法，有较好的收敛性。在解决电力系统潮流计算问题时，是以导纳矩阵为基础的，因此，只要我们能在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高牛顿法潮流程序的效率。自从 60 年代中期，牛顿法中利用了最佳顺序消去法以后，牛顿法在收敛性。内存要求。速度方面都超过了阻抗法，成为了 60 年代末期以后广泛采用的优秀方法。第三章3.1 潮流计算题目Table I Line Data
13、 of the 6-Bus System on 100 MVA BaseLineNumberBusFromNumberTo R XTapRatio1 1 2 0.000 0.300 1.0252 1 4 0.097 0.4073 1 6 0.123 0.5184 2 5 0.282 0.6405 3 5 0.723 1.0506 4 3 0.000 0.133 1.1007 4 6 0.080 0.3703.2 对课题的分析及求解思路此电力系统是一个 6 节点，7 支路的电力网络。综合比较牛顿拉夫逊法（直角坐标、极坐标） 、PQ 分解法等多种求解方法的特点，最后确定采用牛顿拉夫逊法（极坐标）
14、。因为此方法所需解的方程组最少。55.0+j13.050.0+j5.0 30.0+j18.0 P5=50.1612 54 31.1:11.025:1第四章4.1 变压器的型等值电路在电力系统潮流计算中，往往要计算节点导纳矩阵，而我们计算节点导纳矩阵采用节点电压法来实现，如在变压器构成的电力系统中，需要将变压器模型转变成变压器型等值电路（见图 1-1） ，在利用电路知识列节点电压方程，从而导出所需的导纳矩阵。图 3-1 双绕组变压器的型等值电路（i，j 为节点）而在电力系统潮流计算中一般采用标幺值进行计算，标幺值公式如下： 与值 单 与标 与值 值 值 单如果采用标么值计算，元件参数都应归算到同
15、一基准值时得标么值，才能在同一个等值电路上分析和计算。所以，变压器转变成型等值电路时，我们采用标幺值计算，使所求参数为变压器变比 k 的函数。而在一个已经归算好的电力系统网中，若改变变压器的分接头来进行调压，这时变压器的等值电路参数也会相应得改变，此时采用型等值电路进行折算就显得较为方便。下面是变压器的型等值电路分析过程：如不计励磁支路的影响，双绕组变压器可用其阻抗与一个理想变压器串联的电路表示，如图所示。理想变压器只有一个参数，那就是变比 k= 。现以12U变压器阻抗按实际变比归算到低压侧的情况为例，推导出双绕组变压器的型等值电路。流入和流出理想变压器的功率相等：( 、 分别为变压器高、低绕
16、组的实际电压） （3-12/kUI1U2YT:-1） （3-2）2/kI联立（3-1） 、 （3-2）两个公式解得：（3-. 12122TZkTYUI3）（3-.1122 2TkTI4）根据电路原理节点 1、2 的节点电流方程具有如下形式：（3-.112.2IYU5）将式（3-3） 、 （3-4）与式（3-5）比较得(3-6):(3-6)21T21TY/k/因此可以的得到各支路导纳为（3-121T202TT221TYy/k1kY/Yy7） 4.2 节点电压方程在电路中我们学过利用节点电压方程求解某几条支路的电流，现以下图 3-2-1 与图 3-2-2 为例推导节点电压方程组。图 3-2-1 节
17、点电压法为例 图 3-2-2 用电流源代替电压源为例图 3-2-1 表示了一个具有两个电源和你一个等值负荷的系统。 、 为电1e2源电势， 、 为电源的内部导纳， 为负荷的等值导纳， 、 、 为各1y2 3y4y56支路的导纳。如果以地为电压参考点，设节点 1、2、3 的电压为，根据基尔霍夫电流 KCL 法对节点 1、2、3 列节点电流方程得式（3-8）：（3-8）121311216215232.43153230yVyVyyVyVye上式中左端为节点 1、2、3 流出的电流，右端为注入个节点的电流。由上式可以得到一个等效的等值电路图 3-2-2。图 3-2-2 中利用了电流源代替的电压源。在图
18、 1-2-2 中可知的式（3-9）：（3-9） 1122Iye为等值电流源向网络注入的电流。将与式（3-8）联立得式（3-10）：(3-112131212232313233YVYVIYYI10）上式中 称为节点146256346, ,yyYy1、2、3 的自导纳， 称为相应16134235,Y y节点之间的互导纳。因此，在一般情况下，在电力网络中有 n 个节点，则可以按式（3-10）的形式列出 n 个节点方程式，也可用矩阵的形式表示 。其中IYV.1.2.nII .1.2.nV分别为节点注入电流列向量及节点电压列向量； 121212nnnYYYY为节点导纳矩阵，其中对角元素 为节点 i 的自导
19、纳，非对角线 为节点 i 与i ijY节点 j 之间的互导纳。4.3 节点导纳矩阵节点导纳矩阵既可根据自导纳和互导纳的定义直接求取，也可根据电路知识中找出改网络的关联矩阵，在节点电压方程的矩阵形式进行求解。本章节我们主要讨论的是直接求解导纳矩阵。根据节点电压方程章节我们知道，在利用电子数字计算机计算电力系统运行情况是，多采用 形式的节点方程IYV式。其中阶数等于电力网络的节点数。从而可以得到 n 个节点时的节点导纳矩阵方程组（3-11）如下：(3-11）11211 2 232 123.nnnnnYVYVIYVYVI由此可以得到 n 个节点导纳矩阵：11212 212nnnnYYYY 它反映了网
20、络的参数及接线情况，因此导纳矩阵可以看成是对电力网络电气特性的一种数学抽象。由导纳短阵所联系的节点方程式是电力网络广泛应用的一种数学模型。通过上面的讨论，可以看出节点导纳矩阵的有以下特点：（1）导纳矩阵的元素很容易根据网络接线图和支路参数直观地求得，形成节点导纳矩阵的程序比较简单。（2）导纳矩阵为对称矩阵。由网络的互易特性易知 。ijjiY（3）导纳矩阵是稀疏矩阵。它的对角线元素一般不为零，但在非对角线元素中则存在不少零元素。在电力系统的接线图中，一般每个节点与平均不超过34 个其他节点有直接的支路连接。因此，在导纳矩阵的非对角线元素中每行仅有 34 个非零元素，其余的都是零元素，而且网络的规
21、模越大，这种现象越显著。导纳矩阵的对称性和稀疏性对于应用计算机求解电力系统问题有很大的影响。如果能充分地利用这两个特点，如在程序设计中储存导纳矩阵的对角元素和上三角元素（或下三角元素） ，排除零元素的储存和运算，就可以大大地节省储存单元和提高计算速度。节点导纳矩阵的形式可归纳如下：（1）导纳矩阵的阶数等于电力网络的节点数。（2）导纳矩阵各行非对角元素中非零元素的个数等于对应节点所连得不接地支路数。（3）导纳矩阵各对角元素，即节点的自导纳等于相应节点之间的支路导纳之和。（4）导纳矩阵非对角元素，即节点之间的互导纳等于相应节点之间的支路导纳的负值。而在电力系统中进行潮流计算时，往往要计算不同接线下
22、的运行状况，例如，改变变压器主抽头时，潮流分布也随之变化，以及改变其他设备参数进行计算潮流分布，此时就需要导出变化时的导纳矩阵就需要对所设计的程序进行参数设定，而不需要重复上述步骤去导出所求的导纳矩阵。4.4 导纳矩阵在潮流计算中的应用导纳矩阵在潮流计算中的应用起到重要的作用，前面我们介绍了根据系统网络的接线盒参数形成节点导纳矩阵的方法。尽管形成节点导纳矩阵的原理是简单的，但如果采用手算的方法，即使节点数不多的系统也仍然有相当大的工作量。因此只有应用计算机才能快速而准确地完成这些计算任务。本章节我们介绍形成系统节点导纳矩阵的实用程序。为了形成节点导纳矩阵，必须知道电力系统的接线图。从前面的讨论
23、知道，网络接线由节点及连接两个节点的支路确定。实际上，只有输入了各支路两端的节点号，就相当输入了系统的接线图。在计算潮流分布时，我们必须先导出该网络的导纳矩阵，而进行潮流计算时解非线性的节点电压方程的有关方法中，高斯-塞德尔迭代法和牛顿-拉夫逊迭代法是计算机潮流计算中常用的基本方法。这两种方法既可用以解线性方程组，也课用以解非线性方程组。高斯-塞德尔迭代法由于其简单而在早期的潮流计算程序中得以采用。但嗣后就逐渐被牛顿型算法所取代。目前这种方法多半与牛顿型算法配合使用以弥补后者的不足。牛顿-拉夫逊法的收敛性较好，但对初值的要求比较严格，是当前广泛采用的计算机潮流算法。运用计算机进行潮流计算，一般
24、要完成以下几个步骤：建立数学模型、确立计算方法、制定计算机流程并编制程序、上机计算及对计算结果进行分析。因此我们可以知道导纳矩阵在潮流计算中是很重要的。本节只是对导纳矩阵在潮流分布的计算机算法一些简单的描述，我们将在下一章对其进行详细讲解。4.5 潮流计算的手工计算在计算电力系统网络的潮流分布时，我们需要把变压器转化成变压器的型等值电路来进行计算器等效导纳，以下是变压器 的 型等值电路（见图 3-1） ，1T变压器 的型等值电路（见图 3-2） ，其等效导纳的计算过程如下：2T变压器 的型等值电路1T图 3-1所求变压器 参数数据如下：1TTT1210T2220Tj0.31Yj3.sj.yj.
25、50sk1.05.Yj3.j.0793sk1.05y j.j.81z变压器 的型等值电路2T图 3-2所求变压器 参数数据如下：2TTT1230T2240Tj0.13Yj7.518sj.yj6.3sk1Yj7.518j.6214sk1.1y j.j835z 所求电力系统网络的各个节点的相关参数如下： 161401223141456621205Yyy.39j.827054j.32j0.79j3.250.987j.32sj5s.j.Yy.39.827sj 3425256334034550j13.56j.08.576j4.18sy.7j.8sY6.5j0.3.49j0.61.49j8.1sj.y 3
26、6434056625359j.41sY.82.50.j2.34j6.85j0.214.3j.120sy.j.s71.9j.1.j.9s 6164Y0.9j.820.53.820.4.2第五章5.1 Matlab 简介MATLAB 是 由 美 国 mathworks 公 司 发 布 的 主 要 面 对 科 学 计 算 、 可 视 化 以及 交 互 式 程 序 设 计 的 高 科 技 计 算 环 境 。 它 将 数 值 分 析 、 矩 阵 计 算 、 科 学 数 据可 视 化 以 及 非 线 性 动 态 系 统 的 建 模 和 仿 真 等 诸 多 强 大 功 能 集 成 在 一 个 易 于 使用
27、的 视 窗 环 境 中 ， 为 科 学 研 究 、 工 程 设 计 以 及 必 须 进 行 有 效 数 值 计 算 的 众 多科 学 领 域 提 供 了 一 种 全 面 的 解 决 方 案 ， 并 在 很 大 程 度 上 摆 脱 了 传 统 非 交 互 式程 序 设 计 语 言 （ 如 C、 Fortran） 的 编 辑 模 式 ， 代 表 了 当 今 国 际 科 学 计 算 软件 的 先 进 水 平 。MATLAB 和 Mathematica、 Maple 并 称 为 三 大 数 学 软 件 。 它 在 数 学 类 科 技应 用 软 件 中 在 数 值 计 算 方 面 首 屈 一 指 。 M
28、ATLAB 可 以 进 行 矩 阵 运 算 、 绘 制 函数 和 数 据 、 实 现 算 法 、 创 建 用 户 界 面 、 连 接 其 他 编 程 语 言 的 程 序 等 ， 主 要 应用 于 工 程 计 算 、 控 制 设 计 、 信 号 处 理 与 通 讯 、 图 像 处 理 、 信 号 检 测 、 金 融建 模 设 计 与 分 析 等 领 域 。 5.2 Matlab 的应用MATLAB 的 基 本 数 据 单 位 是 矩 阵 ， 它 的 指 令 表 达 式 与 数 学 、 工 程 中 常 用 的形 式 十 分 相 似 ， 故 用 MATLAB 来 解 算 问 题 要 比 用 C， F
29、ORTRAN 等 语 言 完 成 相同 的 事 情 简 捷 得 多 ， 并 且 mathwork 也 吸 收 了 像 Maple 等 软 件 的 优 点 ,使MATLAB 成 为 一 个 强 大 的 数 学 软 件 。 在 新 的 版 本 中 也 加 入 了 对C， FORTRAN， C+ ， JAVA 的 支 持 。 可 以 直 接 调 用 ,用 户 也 可 以 将 自 己 编 写 的实 用 程 序 导 入 到 MATLAB 函 数 库 中 方 便 自 己 以 后 调 用 ， 此 外 许 多 的 MATLAB爱 好 者 都 编 写 了 一 些 经 典 的 程 序 ， 用 户 可 以 直 接
30、进 行 下 载 就 可 以 用 。MATLAB 产 品 族 可 以 用 来 进 行 以 下 各 种 工 作 ： 数 值 分 析 数 值 和 符 号 计 算 工 程 与 科 学 绘 图 控 制 系 统 的 设 计 与 仿 真 数 字 图 像 处 理 技 术 数 字 信 号 处 理 技 术 通 讯 系 统 设 计 与 仿 真 财 务 与 金 融 工 程 MATLAB 的 应 用 范 围 非 常 广 ， 包 括 信 号 和 图 像 处 理 、 通 讯 、 控 制 系 统 设计 、 测 试 和 测 量 、 财 务 建 模 和 分 析 以 及 计 算 生 物 学 等 众 多 应 用 领 域 。 附 加
31、的工 具 箱 （ 单 独 提 供 的 专 用 MATLAB 函 数 集 ） 扩 展 了 MATLAB 环 境 ， 以 解 决这 些 应 用 领 域 内 特 定 类 型 的 问 题 。 5.3.1 与常数的运算常数与矩阵的运算即是同该矩阵的每一元素进行计算。但需注意进行数除时，常数通常只能做除量。基本函数运算中，矩阵的函数运算是矩阵预算中最实用的部分，常用的主要有以下几个：det(a) 求矩阵 a 的行列式eig(a) 求矩阵 a 的特征值inv(a)或 a(-1) 求矩阵 a 的逆矩阵rank(a) 求矩阵 a 的秩trace(a) 求矩阵 a 的迹（对角线元素之和）我们进行工程计算时常常遇到
32、矩阵对应元素之间的运算。这种运算不同于前面讲的数学运算，为有所区别，我们称之为数组运算。5.3.2 基本数学运算数组的加、减与矩阵的加、减运算完全相同。而乘除法运算有相当大的区别，数组的乘除法是指两同维数组对应元素之间的乘除法，它们的运算符为“.*”和“./”或“.。 ”前面讲过常数与矩阵的除法运算中常数只能做除数。在数组运算中有了“对应关系”的规定，数组与常数之间的除法运算没有任何限制。另外，矩阵的数组运算中还有幂运算（运算符.） 、指数运算（exp） 、对数运算(log)、和开方运算（ sqrt） 、等，有了“对应元素”的规定，数组的运算实质上就是针对数组内部的每个元素进行的。矩阵的幂运算
33、与数组的幂运算有很大的区别。5.3.4. 逻辑关系运算逻辑运算是 MATLAB 中数组运算所特有的一种运算形式，也是几乎所有的高级语言普遍适用的一种运算。5.4. Matlab 中的一些命令1）一般 MATLAB 命令格式为 输出参数 1，输出参数 2，=（命令名） （输入参数 1，输入参数2，）输出参数用方括号，输入参数用圆括号如果输出参数只有一个可不使用括号。 2）可用、键来重现已输入的数据或命令。用、键来移动光标进行修改。3）所有 MATLAB 命令都用小写字母。大写字母和小写字母分别表示不同的变量。4）常用数有特定的名字，如 pi（=3.141596） 、Inf（=） 、NaN 则表示
34、不定型求得的结果（如 0/0） 。5）矩阵的输入要一行一行的进行，每行各元素用空格或（， ）分开，每行用（；）分开。6）MATLAB 书写格式为 A=1 2 3 ；4 5 6 ；7 8 9 在 MATLAB 中运行如下程序可得到 A 矩阵a=1 2 3;4 5 6;7 8 9a = 1 2 34 5 67 8 97）需要显示命令的计算结果时，则语句后面不加“；”号，否则要加“；”号。运行下面两种格式可以看出他们的区别a=1 2 3;4 5 6;7 8 9 a=1 2 3;4 5 6;7 8 9;a = （不显示计算结果） 1 2 34 5 68）当输入语句过长需要换行时，应加上“”后再回车，则
35、可连续输入。9）diary 命令使用该命令可以在窗口中以 ASCII 码形式记录所有的输入和输出。但这个命令不是存储数据，而是存储输入与屏幕上输出的内容。它可以记录下工作的过程。在每个工作过程之前使用该命令，工作结束后使用 diary off 则能将整个工作过程记录下来。格式 diary(文件名)（扩展名）. diary off一般来说扩展名可取，m 这样就可在 MATLABBIN目录下存入该文件。 10）save 命令该命令存储定义的变量或演算结果，也可以用来存储指定的变量。 命令格式为save 文件名 . 扩展名 11）what 命令该命令可以在当目录下显示 MATLAB 文件和 MAT
36、数据文件12）dir 命令 显示当前目录下的所有文件.13）clear 命令 14）d1,d2,d3,=size(a) 求矩阵的大小，对 m*n 二维矩阵，第一个为行数 m,第二个为列数 n。如果输入 calear a b c ，则表示清除工作空间中指定变量 a,b,c；如果仅仅输入 calear 命令，则清除整个工作空间。 与此同时，MATLAB 具有强大的矩阵运算功能, 但由于我们在求节点导纳矩阵时用的不多, 因此这里我们只作简单介绍。1) 在 MATLAB 中表示一个矢量要用方括号， 而列矢量的输入只需在行矢量输入格式基础上加转置符()即可。如 x=1 2 3;4 5 6x =1 2 3
37、4 5 6而 x=1 2 3;4 5 6 (加转置符)x =1 42 53 6注意上面两式的区别。2) 下面三条命令可以产生一个行矢量1 a=linspace(x,y,n) 2 a=logspace(x,y,n) 3 a=x:n:y第一条命令可以在线性空间产生一个值在 10x至 10y之间间隔点数为 n 的行矢量(一组数据)。第二条命令可以在对数空间产生一个值在 x 至 y 之间等间隔的行矢量(一组数据)。其行矢量的起始值是 x, 终值为 y, 点数为 n。第三条命令产生 X 至 y 步长为 n 的行矢量。但是, 三个命令之间存在差别，下面的例子可以说明这一点。例一 x=logspace(0,
38、5,6)x =1 10 100 1000 10000 100000例二 x=linspace(0,10,11)x =0 1 2 3 4 5 6 7 8 9 10例三 x=0:1:10x =0 1 2 3 4 5 6 7 8 9 10通过上面三个例子可以看出例一, 例二中 n 代表选取的点数。而在例三 中 n 则表示步长. 我们应当注意它们的区别。3) 矩阵的加, 减, 乘, 除等 , 和其它语言书写一样。但要注意的是在运算符 前面加有(.)则表示是元素对元素的操作.4)以下是常用的运算命令运算命令名 功能Angle 求复数的角Min 求最小值Max 求最大值Sum 求和Roots 求多项式的根
39、Poly 由多项式的根求多项式的系数Polyval 求给定点多项式的值Polyder 多项式求导在进行潮流分布计算时，实际上是由多个简单系统构成的复杂系统，在求节点导纳矩阵时的要用到反馈的指令，所以在 MATLAB 中有下面几种命令可以解决两个系统间的连接问题。1）系统的并联parallel 命令可以实现两个系统的并联。示意图如下：u1 y1u + yu2 y2并联后的系统传递函数表示式为：其中 n1、d 1和 n2、d 2分别为 g1（s） 、12112()()ndgsgsg2（s）的传递函数分子、分母系数行矢量。命令格式： n，d = paralltl（n 1，d 1，n 2，d 2）a，
40、b，c，d = paralltl（a 1，b 1，c 1，d 1，a 2，b 2，c 2，d 2）2）系统的串联series 命令实现两个系统的串联，示意图如下：u1 y u2 y2串联后系统的传递函数为 系统 1系统 2g1（s） g2（s）1212()()ngssd命令格式： n，d = series（n 1，d 1，n 2，d 2） a，b，c，d = series（a 1，b 1，c 1，d 1，a 2，b 2，c 2，d 2）3）系统的反馈feedback 命令实现两个系统的反馈连接，示意图如下：u1 + y1y2 u2连接后系统的传递函数表示为： 1122()()gsnd命令格式：
41、 n，d = feedback（n 1，d 1，n 2，d 2）或： n，d = feedback（n 1，d 1，n 2，d 2，sign） a，b，c，d = feedback（a 1，b 1，c 1，d 1，a 2，b 2，c 2，d 2，sign）其中 sign 是指示 y2到 u1连接的符号，缺省时默认为负（即 sign = -1） 。4）系统的闭环cloop 命令可以将系统的输出反馈到系统的输入构成闭环系统，示意图如下：u y正、负反馈后闭环系统为： 11()gsnmd命令格式： n，d = cloop （n 1，d 1，sign） ac，bc，cc，dc = cloop（a，b，
42、c，d，sign）通过以上对 MATLAB 基本指令的了解，我们就可以对所求的电力系统网络的节点导纳矩阵进行画编程框架图。第六章g1（s）g2（s）g1（s）6.1 潮流计算流程图在上一章节，我们对 MATLAB 的基本指令有了初步的了解，我们就可以画出程序编程框图，按照程序编程框图，我们可以通过 MATLAB 软件进行编程，实现我们所求的节点导纳矩阵。程序编程框图如下：6.2 潮流计算源程序图节点导纳矩阵程序如下:n=input(请输入节点数：n=);n1=input(请输入支路数：n1=);B=input(请输入由支路参数形成的矩阵：B=);X=input(请输入由节点号及其对地阻抗形成的
43、矩阵：X=);Y=zeros(n);for i=1:nif X(i,2)=0;p=X(i,1);Y(p,p)=1./X(i,2);endend for i=1:n1if B(i,6)=0p=B(i,1);q=B(i,2);elsep=B(i,2);q=B(i,1);end Y(p,q)=Y(p,q)-1./(B(i,3)*B(i,5);Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+1./(B(i,3)*B(i,5)2)+B(i,4)./2;Y(p,p)=Y(p,p)+1./B(i,3)+B(i,4)./2;enddisp(导纳矩阵 Y=);disp(Y)按照程序调试的如下结果：请输入
44、节点数：n=6请输入支路数：n1=7请输入由支路参数形成的矩阵：B=1 2 0.300i 0 1.025 1;2 5 0.282+0.640i 0 1 0;5 3 0.723+1.050i 0 1 0;4 3 0.133i 0 1.100 1;4 6 0.080+0.370i 0 1 0;1 6 0.123+0.518i 0 1 0;1 4 0.097+0.407i 0 1 0请输入由节点号及其对地阻抗形成的矩阵：X=1 0;2 0;3 0;4 0;5 0;6 0导纳矩阵 Y=总结此次课程设计我在 MATLAB 编程，潮流计算， WORD 文档的编辑方面均有所提高，但也暴露出了一些问题；理论知识的储备还是不足，对 MATLAB的性能和特点还不能有一个全面的把握，对 WORD 软件也不是很熟练，相信通过以后的学习能弥补这些不足，从而达到一个新的层次。此次课程设计的完成，还要感谢老师的帮助以及同学的协力合作，团结才能出成绩。参考文献1 孟祥萍. 电力系统分析 高等教育出版社.20041 李维波. MATLAB 在电气工程中应用 中国电力出版社.20071 邱晓燕.刘天琪 电力系统分析的计算机算法 中国电力出版社 .20091 吴天明.赵新力.刘建存 MATLAB 电力系统设计与分析 国防工业出版社.2007
