高斯赛德尔潮流计算c语言编程,高斯赛德尔法潮流计算.doc

高斯——赛德尔法潮流计算

? 潮流计算斯赛德尔迭代法(Gauss一Seidel method)求解电力系统潮流的方法。潮流计算高斯赛德尔迭代?法又分导纳矩阵迭代法和阻抗矩阵迭代法两种。前者 是以节导纳矩阵为基础建立的赛德尔迭代格式;后者是以节点阻扰矩阵为基础建立的赛德尔迭代格式。 高斯赛德尔迭代法这是数学上求解线性或非?线性方程组的一种常用的迭代方法。

N

Y

［1］系统节点的分类

根据给定的控制变量和状态变量的不同分类如下

P、Q节点(负荷节点)，给定Pi、Qi求Vi、Si，所求数量最多；

负荷节点，变电站节点(联络节点、浮游节点)，给定PGi、QGi的发电机节点，给定QGi的无功电源节点；

PV节点(调节节点、电压控制节点)，给定Pi、Qi求Qn、Sn，所求数量少，可以无有功储备的发电机节点和可调节的无功电源节点；

平衡节点(松弛节点、参考节点(基准相角)、S节点、VS节点、缓冲节点)，给定Vi，δi=0，求Pn、Qn(Vs、δs、Ps、Qs)。

［2］潮流计算的数学模型

线性的节点电压方程 YV=I

根据S=V可得非线性的节点电压方程(为I的共轭) YV=I==

节点功率与节点电流的关系:

在国外,对于复数变量不打点,其模要加绝对值符号;在国内,对于复数变量,在S、V、I上要打点，Y、Z上不打点，其模不加绝对值符号。

3) 式2—5

对于发电机Pi、Qi为正，对负荷来说Pi、Qi为负

4)展开YV=I得

上式代入式2—5得n维的非线性复数电压方程组

式2—6

该式为潮流计算的基本方程

［3］高斯—赛德尔法潮流计算

1)高斯法潮流计算

① 将式2—6展开成电压方程 式2—7

假设系统节点数是n，PQ节点数为m，m+1及之后的节点是PV节点，第n个节点是平衡节点。

展开式2—7得高斯法潮流计算的基本方程

式2—8

② 考虑到i=1时matlab中for语句的使用可写成

③ 由于平衡节点的电压和相角给定,不用计算,只要计算i=1—n-1节点的电压,但平衡节点的参数和变量要用于其他节点的电压计算.

式2—8的计算过程中有 i=1、2、···n-1

④ 特点：在计算i节点的k+1次电压时，所用的i节点前后(包括i节点)的电压都是k次迭代的结果。

2)高斯—赛德尔法潮流计算

① 在高斯法潮流计算中引入赛德尔法迭代方式即为高斯—赛德尔法潮流计算

② 对应式2—8的高斯—赛德尔法潮流计算的方程为

式2—9

在式2—9的计算中有

③ 特点：在计算i节点的k+1次电压时，1~i-1节点的电压用的是k+1次时的电压，而i~n-1节点的电压用的是k次时的电压，即在迭代过程中每个被求的电压新值立即被带入到下一个电压新值的计算中。

3)基于导纳矩阵的直角坐标高斯—赛德尔法潮流计算

① 设 展开式2—6并将实、虚部分列

式2—10

式2—11

② 令 式2—12

注：、中不包括j=i的参数和变量；

、中分别 有k+1次和k次的变量；

在、中没有单独列出。()

③

④ 将2—12代入式2—10和2—11得

式2—13

式2—14

⑤ 将式2—9展开，实、虚部分列，再将式2—12代入，得节点电压的实部、虚部

式2—15

式2—16

⑥ 对P、V节点，根据常数

式2—17

4)部分求解方程

对于P、Q节点：用式2—15求，用式2—16求

对于P、V节点：用式2—14求用式2—15求，式2—16求

5)为了加速收敛，引入加速因子α，α=1~1.8之间，复数电压: 式2—18

6)实数模型： 式2—19

) 式2—20

、是式2—15～式2—17计算出的值，、是考虑到α修正后的值，、是上一次用于迭代的实际值(不一定是式2—15~式2—17计算出的值)

7)三种加速过程

①每次求出的、立即用于求解下一个电压新值；

②每次求出的、同时立即用α进行修正，得到的、同时用于求解下一个电压新值；

③每次求出的、分别用α进行修正，得到的、分别用于求解下一个电压新值。

注：三种加速过程中，速度又快到慢依次为③②①。

8)收敛判据：

复数模型：

实数模型：，

9)三种收敛判据情况：

①用前后两次经α修正后的电压值；

②用前后两次式2—15~式2