电力系统分析中, 潮流计算的意义十分重大,它是在电网正常或故障情况下的稳定运行状态的计算，为了更加深入地理解潮流计算，本文基于牛顿-拉夫逊潮流计算原理，设计出一个8机28节点的电力系统网络结构，然后基于Matlab进行仿真运行分析。   各类配电网潮流算法性能通常从以下几个方面进行分析: 算法的收敛速度、稳定性、算法的复杂程度。潮流计算的稳定性对于维护电力系统稳定潮流计算。它根据给定的初始值和网络结构确定整个系统的运行状态, 确定各个母线上的电压幅值及相角、网络中的功率分布以及网络损耗等。潮流计算的基本方程是由电力系统的网络方程得到的。在电力系统中都存在一些静态装置, 例如变压器、输电线路、并联电容器和电抗器等, 它们可以由电阻、电感和电容等基础元件构成的等值电路来模拟。

  牛顿-拉夫逊法实质上就是切线法，是逐步线性化的方法。牛顿-拉夫逊计算法不仅用于求解单变量方程，还可以求解多变量非线性代数方程。

  设有单变量非线性方程   在求解此方程的时候，需要先给出解的近似值 x ( 0 ) x^{(0)} x(0),它的与真解的误差将满足下列方程式   将上述式子左边的函数在 x ( 0 ) x^{(0)} x(0)附近展开成泰勒级数得到：   如果差值 的值很小，差值的二次及以上阶次的各项均可以省去，于是可以简化成：   这是对于变量的修正量线性方程式（修正方程式），据此可以得到修正量   但是修正后的近似解 x ( 1 ) x^{(1)} x(1)与真解仍然存在误差，为了进一步逼近真解，可以一直迭代下去，最终可以得到迭代计算通式：   而迭代过程中的收敛判据为：   由此可见，牛顿-拉夫逊法实质上就是切线法，是逐步线性化的方法。牛顿-拉夫逊计算法不仅用于求解单变量方程，还可以求解多变量非线性代数方程。 设有n个联立的非线性代数方程   假定已经给出各变量的初值 x 1 ( 0 ) x1^{(0)} x1(0), x 2 ( 0 ) x2^{(0)} x2(0),…, x n ( 0 ) x n^{(0)} xn(0),令   分别为各个变量的修正量，使其满足方程组，即   将n个多元函数在初始值附近分别展开成泰勒级数，并略去二次及以上阶次的各项，可以得到牛顿-拉夫逊法的修正方程式   利用高斯消去法或者三角分解法可以解出修正量,然后对初始值近似解进行修正   如此反复迭代，在进行第k+1次迭代时，从而求出修正方程式：   得到修正量   并对各变量进行修正，最终可以表示为   其中，** X X X**是由n个变量组成的n维矩阵列向量， F （ X ） F（X） F（X）,是由n个多元函数组成的n维向量: J J J是n*n阶方阵，称为雅可比矩阵。迭代过程一直进行到满足收敛判据

  由于节点电压可以采用不同的坐标系表示，牛顿-拉夫逊潮流计算也将相应地采用不同的计算公式。   采用直角坐标时，节点电压可以表示为   导纳矩阵元素可以表示为   假定系统中的第1,2，…，m号节点为PQ节点，第i个节点的给定功率设为 P i s Pis Pis和 Q i s Qis Qis，对于该节点可以写出如下方程   假定系统中的第m+1，m+2，…，n-1号节点为PU节点，对于其中每一个节点可以写出如下方程   第n号节点为平衡节点，其中电压   是给定的，所以不参加迭代，对于PQ节点方程式和PU节点方程式来说，总共包含2(n-1)个方程式，带求的变量也是2（n-1）个，因此，我们可以得到如下的修正方程式   其中修正方程式还可以写成分块矩阵的形式   其中，   对于PQ节点   对于PU节点   雅可比矩阵各元素都是节点电压的函数，它们的数值将在迭代的过程中不断地改变，因此在用牛顿-拉夫逊法计算潮流首先要输入网络的原始数据以及各节点的给定值并且形成节点导纳矩阵，输入节点电压初值，设置迭代计数k=0,然后开始进入牛顿-拉夫逊算法的迭代过程。

（1） baseMVA   baseMVA是一个标量，用来设置基准容量。对于计算中采用有名值，可以根据实际情况设置，在本文设计中设置为100MVA。

（2） bus data矩阵   bus矩阵是关于电网母线节点参数设计的矩阵，通过设置母线的参数来达到计算的目的，如表3.1所示，bus data有如下参数

  其中，各项参数的含义为：

（3） generator data矩阵   generator data矩阵是用来设置电网中发电机的参数，各项参数名称如表3.2所示：

     其中，

bus用来设置接入发电机的母线编号

Pg用来设置接入发电机的有功功率，注意功率输入的是有名值

Qg用来设置接入发电机的无功功率

Qmax用来设置接入发电机的无功功率的最大允许值

Qmin用来设置接入发电机的无功功率的最小允许值

Vg用来设置接入发电机的工作电压，注意输入的是标幺值

mBase用来设置接入发电机的功率基准

status用来设置发电机的工作状态，1表示投入运行，2表示投出运行

Pmax用来设置接入发电机的无功功率的最大允许值

Pmin用来设置接入发电机的无功功率的最小允许值

其余的Pc1，Pc2，Qc1min，Qc1max，Qc2min，Qc2max，ramp_agc，ramp_10，ramp_30，ramp_q，apf均表示发电机其他特征量，在实际计算中可以设置为0

（4） branch data矩阵   brach data是一个矩阵，用来设置电网中各个支路之间的参数，矩阵参数名称如表3.3所示：

  其中，

（5） generator cost data矩阵   一般来说，generator cost data矩阵保持原有默认的设计，不作修改。

  基于Matpower的潮流计算矩阵的设计，需要改变如上矩阵的参数，使其满足收敛条件，最终得到相应的计算结果，因此熟悉Matpower矩阵的编写对于最终的仿真结果有着极其重要的意义。

  电力网络中的运行状态可以用节点方程或者回路方程进行描述，节点方程以母线电压作为待求量，节点方程的求解结果对于潮流计算应用广泛。   一般地，对于有n个独立节点的网络，可以列写出n个节点方程   可以写成矩阵形式   矩阵Y称为节点导纳矩阵，其对角线元素Y称为节点i的自导纳，等于接于节点i的所有支路导纳之和，非对角线元素成为节点ij之间的互导纳，等于直接连接于节点ij之间的支路导纳的负值，如果节点ij之间不存在直接支路，则导纳为0，因此，节点导纳矩阵是一个稀疏的对称矩阵。

  基于上述原理，编写出电力系统的节点导纳矩阵的计算代码