电力系统分析课程设计

2009-2010学年度下学期电力系统分析课程设计

电力系统短路故障的计算机

算法程序设计

2010年 6 月 1 日

第1章 第2章 第3章 第4章 第5章 第6章 第7章 第8章

目 录

课程设计目的.........................................................2 课程设计内容........................................................2 课程设计原理........................................................2 选择语言理由........................................................3 程序主框图及主要数据变量说明....................5 程序及说明.............................................................9 程序结果及分析....................................................12 参考文献..................................................................17

第1章 课程设计目的

根据所给的电力系统，编制短路电流计算程序，通过计算机进行调试，最后完成一个切实可行的电力系统计算应用程序。通过自己设计电力系统计算程序使同学们对电力系统分析有进一步理解，同时加强计算机实际应用能力的训练。

第2章 课程设计内容

电力系统故障的计算程序设计及编制和调试。

第3章 课程设计原理

1．计算机计算原理

应用计算机进行电力系统计算，首先要掌握电力系统相应计算的数学模型；其次是运用合理的计算方法；第三则是选择合适的计算机语言编制计算程序。

建立电力系统计算的相关数学模型，就是建立用于描述电力系统相应计算的有关参数间的相互关系的数学方程式。该数学模型的建立往往要突出问题的主要方面，即考虑影响问题的主要因素，而忽略一些次要因素，使数学模型既能正确地反映实际问题，又使计算不过于复杂。

运用合理的计算方法，就是要求所选用的计算方法能快速准确地得出正确结果，同时还应要求在解算过程中占用内存少，以利提高计算机的解题规模。

选择合适的语言编写程序，就是首先确定用什么计算机语言来编制程序；其次是作出计算的流程图；第三根据流程图用选择的语言编写计算程序。然后上机调试，直到语法上无错误。

所编制的程序难免存在逻辑错误，因此先用一个已知结果的系统作为例题进行计算。用程序计算的结果和已知结果相比较，如果结果相差甚远就要逐步分析程序的计算步骤，查出问题的出处；如果结果比较接近，则逐步分析误差来源；直到结果正确为止。

最后将所编制出的正确计算程序，用于电力系统的实际计算。 2．电力系统短路计算计算机算法

一般在电力系统短路计算中，多数情况下只要计算短路电流、电压的周期分量起始值。因此，电力系统短路电流计算的数学模型，可归结为求解短路故障初始状态下的等值电路稳态解的问题。对于三相对称短路，可建立一相等值电路的计算模型，对于不对称短路，则可应用对称分量法建立系统的正、负、零序网络，

从而建立故障计算的序网模型。

（1）对称短路计算

有关数学模型和原理框图以及已知结果的例题，参见《电力系统分析》第六章。常用的计算方法为节点导纳矩阵法或节点阻抗矩阵法，其形成方法分别参见《电力系统分析》第四章。

（2）简单不对称短路计算

有关数学模型、原理框图以及已知结果的例题参见《电力系统分析》第八章8-4节。

第4章 选择语言理由

我们使用的是第四代计算机语言的MATLAB，利用其丰富的函数资源，它的优点如下：

1）语言简洁紧凑，使用方便灵活，库函数极其丰富。MATLAB程序书写形式自由，利用起丰富的库函数避开繁杂的子程序编程任务，压缩了一切不必要的编程工作。由于库函数都由本领域的专家编写，用户不必担心函数的可靠性。可以说，用MATLAB进行科技开发是站在专家的肩膀上。

具有FORTRAN和C等高级语言知识的读者可能已经注意到，如果用FORTRAN或C语言去编写程序，尤其当涉及矩阵运算和画图时，编程会很麻烦。例如，如果用户想求解一个线性代数方程，就得编写一个程序块读入数据，然后再使用一种求解线性方程的算法（例如追赶法）编写一个程序块来求解方程，最后再输出计算结果。在求解过程中，最麻烦的要算第二部分。解线性方程的麻烦在于要对矩阵的元素作循环，选择稳定的算法以及代码的调试动不容易。即使有部分源代码，用户也会感到麻烦，且不能保证运算的稳定性。解线性方程的程序用FORTRAN和C这样的高级语言编写，至少需要四百多行，调试这种几百行的计算程序可以说很困难。以下用MATLAB编写以上两个小程序的具体过程。

MATLAB求解下列方程，并求解矩阵A的特征值。 Ax=b,其中：

A= 32 13 45 67 23 79 85 12 43 23 54 65