实验目的
根据所给的电力系统，编制潮流计算程序，通过计算机进行调试，最后完成一个切实可行的电力系统计算应用程序。通过自己设计电力系统计算程序使同学对电力系统分析有进一步理解，同时加强计算机实际应用能力的训练。 程序计算原理 1、概述
应用计算机进行电力系统计算，首先要掌握电力系统相应计算的数学模型；其次是运用合理的计算方法；第三则是选择合适的计算机语言编制计算程序。
建立电力系统计算的相关数学模型，就是建立用于描述电力系统相应计算的有关参数间的相互关系的数学方程式。该数学模型的建立往往要突出问题的主要方面，即考虑影响问题的主要因素，而忽略一些次要因素，使数学模型既能正确地反映实际问题，又使计算不过于复杂。
运用合理的计算方法，就是要求所选用的计算方法能快速准确地得出正确结果，同时还应要求在解算过程中占用内存少，以利提高计算机的解题规模。
选择合适的语言编写程序，就是首先确定用什么计算机语言来编制程序；其次是作出计算的流程图；第三根据流程图用选择的语言编写计算程序。然后上机调试，直到语法上无错误。本程序采用C语言进行编程。
所编制的程序难免存在逻辑错误，因此先用一个已知结果的系统作为例题进行计算。用程序计算的结果和已知结果相比较，如果结果相差甚远就要逐步分析程序的计算步骤，查出问题的出处；如果结果比较接近，则逐步分析误差来源；直到结果正确为止。 2、电力系统潮流计算的程序算法
潮流计算是电力系统分析中的一种最基本的计算，它的任务是对给定的运行条件确定系统的运行状态，如母线上的电压（幅值及相角）、网络中的功率分布及功率损耗等。
目前计算机潮流计算的方法主要有牛顿-拉夫逊算法和PQ分解法。牛顿-拉夫逊算法是数学上求解非线形方程组的有效方法，具有较好的收敛性，曾经是潮流计算中应用比较普遍的方法。PQ快速分解法是从牛顿-拉夫逊算法演变而来的，是将纯数学的牛顿-拉夫逊算法及电力系统具体特点相结合并进行简化及改进而得出的。PQ快速分解法比牛顿-拉夫逊算法大大提高了计算速度和节省了内存，故而本程序以PQ快速分解法进行潮流计算。
1）形成节点导纳矩阵 （1）自导纳的形成
对节点i其自导纳Yii是节点i以外的所有节点都接地时节点i对地的总导纳。显然，Yii应等于及节点i相接的各支路导纳之和，即
式中，yi0为节点i及零电位节点之间的支路导纳；yij为节点i及节点
j之间的支路导纳。
（2）互导纳的形成
对节点i及节点k之间的互导纳是节点i、k之间的支路导纳的负值，即Yik??yik
不难理解Yki?Yik。若节点i和k没有支路直接相连时，便有Yik=0 （3）含变压器支路的处理
若节点p、q间接有变压器，如下图所示，则可作出其∏型等值电路为：
图1 变压器∏型等值电路
则p、q的自导纳和节点间的互导纳分别为
2）计算不平衡功率△P、△Q并形成修正方程式
对每一个PQ节点或每一个PV节点都可以根据下列公式计算出有功功率增量△P
Pi?Pis?Pi?Pis?Vi?Vj(Gijcos?ij?Bijsin?ij)?????(i?1,2,j?1nn,n?1)
而对于每一个PQ节点还可以根据下面的公式计算出无功功率增量△Q
Qi?Qis?Qi?Qis?Vi?Vj(Gijsin?ij?Bijcos?ij)????(i?1,2,j?1,m)
在有功功率增量和无功功率增量不满足如下约束条件时
利用PQ分解法则可以形成如下修正方程
?P1??V??B11?????????????B1,n?1?1?????V1?1?B?????B?????BP?2?222,n?1???21?V??V??????????????????????????????????22?
2?????????Bn?1,1??Bn?1,2??Bn?1,n?1???V??????n?1n?1?P?n?1???????Vn?1?Q1??B11?????????????B1,m?V1?????V1?B?????B?????BQ2?222,m???21?VV2?????????????????????????????????2?
???????Bm,1?????Bm,2???Bm,m???V????m?Qm?????Vm?3）利用因子表法求解修正方程
???????????在电网计算中经常遇到这样的问题，对方程组需要反复多次求解，而
每次求解仅改变常数项F，系数矩阵保持不变。按照一般的高斯消去法，对每一改变的常数项，形成包括常数项及系数矩阵在内的增广矩阵，然后消去回代求出其解。可以看出，每次对增广矩阵中A矩阵元素的消元都是重复的，为了避免这种重复，我们把对相同的系数矩阵重复进行的消去及对不同的常数项进行的消去分开进行，因此对系数矩阵的消去只需进行一次，并在消去的过程中将对常数项进行消去运算的运算因子保存下来，形成所谓因子表，这就是因子表法。因为因子表记录了高斯消去法对常数项进行消去的全部信息，利用它便可对不同常数项进行消去，形成上三角矩阵，最后求出全部未知数。
在使用PQ分解法时，其系数矩阵是在迭代过程中保持不变的，所以为了节省内存和缩短运算时间我们采取了因子表法。同时由于电网的节点导纳矩阵矩阵是稀疏阵和对称阵，于是我们可以采取只保存系数矩阵的上三角阵来使运算更为简化。
若线性方程组一般形式如下：
?a11??a12??a13???a1n??x1??f1????????a??a??a?????22232n??x2???f2???????????????a33??a3n???x3???f3??????? ?????????????????????????????????????????????????????a??x??f?nn??n??n???a11??a12??a13???a1n????????a??a??a?22232n??T其中A???????????????a33??a3n?称为系数矩阵，X??x1??x2??x3??xn?称为未知数
???????????????????????????????????????????????????a?nn??向量， F??f1??f2??f3??理：
fn?称为常数项向量。将矩阵A的元素进行如下处
T(k)(k?1)(k?1)(k)aij?aij?aikakj????(k?1,2,,i?1;??j?k?1,k?2,,n)
,n)
a?(i)ij(i?1)aija(i?1)ii????(j?i?1,i?2,
电力系统潮流计算C语言程序及说明
