高斯赛德尔法潮流计算

高斯——赛德尔法潮流计算潮流计算高斯——赛德尔迭代法(Gauss一Seidel method)是求解电力系统潮流的方法。

潮流计算高斯——赛德尔迭代法又分导纳矩阵迭代法和阻抗矩阵迭代法两种。

前者是以节点导纳矩阵为基础建立的赛德尔迭代格式;后者是以节点阻扰矩阵为基础建立的赛德尔迭代格式。

高斯——赛德尔迭代法这是数学上求解线性或非线性方程组的一种常用的迭代方法。

本实验通过对电力网数学模型形成的计算机程序的编制与调试，获得形成电力网数学模型：高斯---赛德尔法的计算机程序，使数学模型能够由计算机自行形成，即根据已知的电力网的接线图及各支路参数由计算程序运行形成该电力网的节点导纳矩阵和各节点电压、功率。

通过实验教学加深学生对高斯---赛德尔法概念的理解，学会运用数学知识建立电力系统的数学模型，掌握数学模型的形成过程及其特点，熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。

高斯---赛德尔法潮流计算框图［1］系统节点的分类根据给定的控制变量和状态变量的不同分类如下①P、Q节点（负荷节点），给定Pi、Qi求Vi、Si，所求数量最多；②负荷节点，变电站节点（联络节点、浮游节点），给定P Gi、Q Gi的发电机节点，给定Q Gi的无功电源节点；③PV节点（调节节点、电压控制节点），给定P i、Q i求Q n、S n，所求数量少，可以无有功储备的发电机节点和可调节的无功电源节点；④平衡节点（松弛节点、参考节点（基准相角）、S节点、VS节点、缓冲节点），给定V i，δi=0，求P n、Q n(V s、δs、P s、Q s)。

［2］潮流计算的数学模型1)线性的节点电压方程YV=I根据S=V错误!未找到引用源。

可得非线性的节点电压方程(错误!未找到引用源。

为I的共轭) YV=I=错误!未找到引用源。

=错误!未找到引用源。

节点功率与节点电流的关系:错误!未找到引用源。

2)在国外,对于复数变量不打点,其模要加绝对值符号;在国,对于复数变量,在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—174）部分求解方程对于P、Q节点：用式2—15求错误!未找到引用源。

，用式2—16求错误!未找到引用源。

对于P、V节点：用式2—14求错误!未找到引用源。

用式2—15求错误!未找到引用源。

，式2—16求错误!未找到引用源。

5）为了加速收敛，引入加速因子α，α=1~1.8之间，复数电压: 错误!未找到引用源。

式2—186)实数模型：错误!未找到引用源。

式2—19错误!未找到引用源。

）式2—20错误!未找到引用源。

、错误!未找到引用源。

是式2—15～式2—17计算出的值，错误!未找到引用源。

、错误!未找到引用源。

是考虑到α修正后的值，错误!未找到引用源。

、错误!未找到引用源。

是上一次用于迭代的实际值（不一定是式2—15~式2—17计算出的值）7）三种加速过程①每次求出的错误!未找到引用源。

、错误!未找到引用源。

立即用于求解下一个电压新值；②每次求出的错误!未找到引用源。

、错误!未找到引用源。

同时立即用α进行修正，得到的错误!未找到引用源。

、错误!未找到引用源。

同时用于求解下一个电压新值；③每次求出的错误!未找到引用源。

、错误!未找到引用源。

分别用α进行修正，得到的错误!未找到引用源。

、错误!未找到引用源。

分别用于求解下一个电压新值。

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

8）收敛判据：复数模型：错误!未找到引用源。

实数模型：错误!未找到引用源。

，错误!未找到引用源。

9）三种收敛判据情况：①用前后两次经α修正后的电压值；②用前后两次式2—15~式2—17计算出来的值；③前一次用α修正的值，后一次用式2—15~式2—17计算出的值。

10）高斯—赛德尔法是用前后两次迭代的最大电压误差作收敛判据，ε取10-5~10-6,牛顿法是用最大功率误差为收敛判据，ε取10-3~10-5，所以后者为好。

［4］编程程序步骤如下第一步：设定初值0max =∆V ，1=i 定义Z 矩阵,s 设定循环次数100=k第二步：用一判据（0)2,(==i Z ）先求PQ 节点用2-15式求)1(+k i e ，再代入2-16替代)(k i e 求)1(+k i f 。

则);(;;)1()1()1()()1()1()1()1()1(++++++++•=-=∆+=k i k i k i k i k i k i k i k i k e f arctg V V V jf e V i δif;;)1(max max )1(++∆=∆∆>∆k i k i V V V V 根据收敛判据5max 10-=<∆εV 输出代求量，即if;)7,(;)6,(10)1()1(6max ++-==<∆k i k i i Z V i Z V δ第三步：（0)!2,(=i Z ）求PV 节点用2-14求i Q再用2-16求)1(+k i f ，将其代入2-17，求)1(+k i e ，则);/(;;)1()1()1()()1()1()1()1()1(++++++++•=-=∆+=k i k i k i k i k i k i k i k i k e f arctg V V V jf e V i δif;)1(max max)1(++∆=∆∆>∆k i k i V V V V 根据收敛判据6max 10-=<∆εV 输出代求量，即if;)7,(;)3,(10)1()1(6max ++-==<∆k i k i i Z Q i Z V δ第四步：求平衡节点n利用式2-13和2-14式求i P 和i Q ,然后输出，即 ;)3,(;)2,(i i Q i Z P i Z ==最后输出Z 矩阵试验题目：用形成Y 阵的五节点系统，假定节点1、2、3为PQ 节点，节点4为PV 节点、节点5为平衡节点，试分别用高斯—赛德尔法潮流计算其潮流。

取 收敛判据为｜△m ax V ｜<610-。

给定：程序如下：clearclcI=[-2,-3,2,2,3];J=[4,5,3,1,1];R=[0,0,0.08,0.04,0.1];X=[0.015,0.03,0.3,0.25,0.35];K=[1.05,1.05,0.25,0.25,0];n=5;L=5;Y=zeros(2*n,n);for m=1:Li=I(m);j=J(m);r=R(m);x=X(m);k=K(m);if i*j==0Y(2*i-1,i)=Y(2*i-1,i)+r;Y(2*i,i)=Y(2*i,i)-x;endif i*j>0Y(2*i-1,j)=Y(2*i-1,j)-r/(r^2+x^2);Y(2*i,j)=Y(2*i,j)+x/(r^2+x^2);Y(2*j-1,i)=Y(2*i-1,j);Y(2*j,i)=Y(2*i,j);Y(2*i-1,i)=Y(2*i-1,i)+r/(r^2+x^2);Y(2*i,i)=Y(2*i,i)-x/(r^2+x^2)+k;Y(2*j-1,j)=Y(2*j-1,j)+r/(r^2+x^2);Y(2*j,j)=Y(2*j,j)-x/(r^2+x^2)+k;endif i*j<0i=-i;Y(2*i-1,j)=Y(2*i-1,j)-r/(r^2+x^2)/k;Y(2*i,j)=Y(2*i,j)+x/(r^2+x^2)/k;Y(2*j-1,i)=Y(2*i-1,j);Y(2*j,i)=Y(2*i,j);Y(2*i-1,i)=Y(2*i-1,i)+r/(r^2+x^2)/k^2;Y(2*i,i)=Y(2*i,i)-x/(r^2+x^2)/k^2;Y(2*j-1,j)=Y(2*j-1,j)+r/(r^2+x^2);Y(2*j,j)=Y(2*j,j)-x/(r^2+x^2);endendYP=[-1.6,-2.0,-3.7,5.0,0];Q=[-0.8,-1.0,-1.3,0,0];E=[1,1,1,1.05,1.05];F=[0,0,0,0,0];k=0;V=[1,1,1,1.05,1.05];A=[0,0,0,0,0];h=3;m=0.000001;Vm=1;while Vm>mVm=0;for i=1:n-1A1=0;A2=0;if i>jfor j=1:i-1g=Y(2*i-1,j);b=Y(2*i,j);e=E(j);f=F(j);A1=A1+g*e-b*f;A2=A2+g*f+b*e;endendfor j=i+1:ng=Y(2*i-1,j);b=Y(2*i,j);e=E(j);f=F(j);A1=A1+g*e-b*f;A2=A2+g*f+b*e;ende=E(i);f=F(i);p=P(i);q=Q(i);g=Y(2*i-1,i);b=Y(2*i,i);g=Y(2*i-1,i);b=Y(2*i,i);Q(i)=-b*(e^2+f^2)-e*A2+f*A1;q=Q(i);E(i)=g/(g^2+b^2)*((p*e+q*f)/(e^2+f^2)-A1)+b/(g^2+b^2)*((p*f-q*e)/( e^2+f^2)-A2);v=V(i);F(i)=sqrt(v^2-E(i)^2);A(i)=atan(F(i)/E(i));A(i)=A(i)*180/pi;continueendE(i)=g/(g^2+b^2)*((p*e+q*f)/(e^2+f^2)-A1)+b/(g^2+b^2)*((p*f-q*e)/( e^2+f^2)-A2);F(i)=g/(g^2+b^2)*((p*f-q*e)/(e^2+f^2)-A2)+b/(g^2+b^2)*((p*e+q*f)/( e^2+f^2)-A1);v=sqrt(E(i)^2+F(i)^2);Vc=v-V(i);Vc=abs(Vc);if Vc>VmVm=Vc;endV(i)=v;A(i)=atan(F(i)/E(i));A(i)=A(i)*180/pi;endk=k+1;endfor j=1:ne=E(j);f=F(j);g=Y(2*i-1,j);b=Y(2*i,j);P(n)=P(n)+E(n)*(g*e-b*f);Q(n)=Q(n)-E(n)*(g*f+b*f);endkPQVA运行结果：1.3787 -0.6240 -0.7547 0 0-6.2917 3.9002 2.6415 0 0 -0.6240 1.4539 -0.8299 0 03.9002 -66.9808 3.1120 63.4921 0-0.7547 -0.8299 1.5846 0 02.64153.1120 -35.7379 0 31.74600 0 0 0 00 63.4921 0 -66.6667 00 0 0 0 00 0 31.7460 0 -33.3333 k =11P =-1.6000 -2.0000 -3.7000 5.0000 0.5238-0.8000 -1.0000 -1.3000 1.3885 0.5238 V =0.8885 1.0817 1.0579 1.0500 1.0500A =-11.6107 -0.4133 1.1798 0.0028 0 >>。