高斯赛德尔法潮流计算

1、高斯赛德尔法潮流计算  潮流计算高斯赛德尔迭代法(gauss一seidel method)是求解电力系统潮流的方法。潮流计算高斯赛德尔迭代 法又分导纳矩阵迭代法和阻抗矩阵迭代法两种。前者 是以节点导纳矩阵为基础建立的赛德尔迭代格式;后者是以节点阻扰矩阵为基础建立的赛德尔迭代格式。 高斯赛德尔迭代法这是数学上求解线性或非 线性方程组的一种常用的迭代方法。本实验通过对电力网数学模型形成的计算机程序的编制与调试，获得形成电力网数学模型：高斯-赛德尔法的计算机程序，使数学模型能够由计算机自行形成，即根据已知的电力网的接线图及各支路参数由计算程序运行形成该电力网的节点导纳矩

2、阵和各节点电压、功率。通过实验教学加深学生对高斯-赛德尔法概念的理解，学会运用数学知识建立电力系统的数学模型，掌握数学模型的形成过程及其特点，熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。高斯-赛德尔法潮流计算框图置迭代计数b=0给定pq节点电压初值给定pv节点电压实部（或虚部）输入数据，定义数组开始计算pq节点电压实部和虚部先计算pv节点无功功率再用其计算pv节点电压实部和虚部计算平衡节点的有功和无功求=+ b=b+1判断所有|是否<0.000001 n y结果输出结束1系统节点的分类根据给定的

3、控制变量和状态变量的不同分类如下 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潮流计算的数学模型1) 线性的节点电压方程 yv=i根据s=v可得非线性的节点电压方程(为i的共轭) yv=i= 节点功率与节点电流的关

4、系: 2) 在国外,对于复数变量不打点,其模要加绝对值符号;在国内,对于复数变量,在s、v、i上要打点，y、z上不打点，其模不加绝对值符号。 3） 式25 对于发电机pi、qi为正，对负荷来说pi、qi为负4）展开yv=i得 上式代入式25得n维的非线性复数电压方程组 式26 该式为潮流计算的基本方程3高斯赛德尔法潮流计算 1)高斯法潮流计算 将式26展开成电压方程 式27 假设系统节点数是n，pq节点数为m，m+1及之后的节点是pv节点，第n个节点是平衡节点。 展开式27得高斯法潮流计算的基本方程 式28 考虑到i=1时matlab中for语句的使用可写成 由于平衡节点的电压和相角给定,不用

5、计算,只要计算i=1n-1节点的电压,但平衡节点的参数和变量要用于其他节点的电压计算. 式28的计算过程中有 i=1、2、···n-1 特点：在计算i节点的k+1次电压时，所用的i节点前后(包括i节点)的电压都是k次迭代的结果。 2）高斯赛德尔法潮流计算 在高斯法潮流计算中引入赛德尔法迭代方式即为高斯赛德尔法潮流计算 对应式28的高斯赛德尔法潮流计算的方程为 式29 在式29的计算中有 特点：在计算i节点的k+1次电压时，1i-1节点的电压用的是k+1次时的电压，而in-1节点的电压用的是k次时的电压，即在迭代过程中每个被求的电压新值立即被带入到下一个电压新值的计

6、算中。3）基于导纳矩阵的直角坐标高斯赛德尔法潮流计算 设 展开式26并将实、虚部分列 式210 式211 令 式212 注：、中不包括j=i的参数和变量； 、中分别 有k+1次和k次的变量； 在、中没有单独列出。（） 将212代入式210和211得 式213 式214 将式29展开，实、虚部分列，再将式212代入，得节点电压的实部、虚部 式215 式216 对p、v节点，根据常数 式2174）部分求解方程对于p、q节点：用式215求，用式216求对于p、v节点：用式214求用式215求，式216求5）为了加速收敛，引入加速因子，=11.8之间，复数电压: 式2186)实数模型： 式219 ）

7、式220 、是式215式217计算出的值，、是考虑到修正后的值，、是上一次用于迭代的实际值（不一定是式215式217计算出的值）7）三种加速过程每次求出的、立即用于求解下一个电压新值；每次求出的、同时立即用进行修正，得到的、同时用于求解下一个电压新值；每次求出的、分别用进行修正，得到的、分别用于求解下一个电压新值。注：三种加速过程中，速度又快到慢依次为。8）收敛判据：复数模型：实数模型：，9）三种收敛判据情况：用前后两次经修正后的电压值；用前后两次式215式217计算出来的值；前一次用修正的值，后一次用式215式217计算出的值。10）高斯赛德尔法是用前后两次迭代的最大电压误差作收敛判据，取1

8、0-510-6,牛顿法是用最大功率误差为收敛判据，取10-310-5，所以后者为好。4编程程序步骤如下第一步：设定初值，定义z矩阵,s设定循环次数第二步：用一判据（）先求pq节点用2-15式求，再代入2-16替代求。则 if 根据收敛判据输出代求量，即 if 第三步：（）求pv节点用2-14求再用2-16求，将其代入2-17，求，则 if 根据收敛判据输出代求量，即 if 第四步：求平衡节点n利用式2-13和2-14式求和,然后输出，即 最后输出z矩阵试验题目：用形成y阵的五节点系统，假定节点1、2、3为pq节点，节点4为pv节点、节点5为平衡节点，试分别用高斯赛德尔法潮流计算其潮流。取收敛判

9、据为<。给定：程序如下：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:l i=i(m);j=j(m);r=r(m);x=x(m);k=k(m); if i*j=0 y(2*i-1,i)=y(2*i-1,i)+r; y(2*i,i)=y(2*i,i)-x; end if i*j>0 y(2*i-1,j)=y(2*i-1,j)-r/(r2+x2); y(2*i,j

10、)=y(2*i,j)+x/(r2+x2); 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/(r2+x2); y(2*i,i)=y(2*i,i)-x/(r2+x2)+k; y(2*j-1,j)=y(2*j-1,j)+r/(r2+x2); y(2*j,j)=y(2*j,j)-x/(r2+x2)+k; end if i*j<0 i=-i; y(2*i-1,j)=y(2*i-1,j)-r/(r2+x2)/k; y(2*i,j)=y(2*i,j)+x/(r2+x2)/k; y(2*j-1,i)=y(2*i-1,j

11、); y(2*j,i)=y(2*i,j); y(2*i-1,i)=y(2*i-1,i)+r/(r2+x2)/k2; y(2*i,i)=y(2*i,i)-x/(r2+x2)/k2; y(2*j-1,j)=y(2*j-1,j)+r/(r2+x2); y(2*j,j)=y(2*j,j)-x/(r2+x2); 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>m

12、vm=0;for i=1:n-1 j=1; a1=0; a2=0; if i>j for j=1:i-1 g=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; end end for j=i+1:n g=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; end e=e(i);f=f(i);p=p(i);q=q(i);g=y(2*i-1,i);b=y(2*i,i); if i>h g=y(2*i-1,i);

13、 b=y(2*i,i); q(i)=-b*(e2+f2)-e*a2+f*a1; q=q(i); e(i)=g/(g2+b2)*(p*e+q*f)/(e2+f2)-a1)+b/(g2+b2)*(p*f-q*e)/(e2+f2)-a2); v=v(i); f(i)=sqrt(v2-e(i)2); a(i)=atan(f(i)/e(i); a(i)=a(i)*180/pi; continue end e(i)=g/(g2+b2)*(p*e+q*f)/(e2+f2)-a1)+b/(g2+b2)*(p*f-q*e)/(e2+f2)-a2); f(i)=g/(g2+b2)*(p*f-q*e)/(e2+f2

14、)-a2)+b/(g2+b2)*(p*e+q*f)/(e2+f2)-a1); v=sqrt(e(i)2+f(i)2); vc=v-v(i); vc=abs(vc); if vc>vm vm=vc; end v(i)=v; a(i)=atan(f(i)/e(i); a(i)=a(i)*180/pi;endk=k+1;endfor j=1:n e=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);endk pqva运行结果：y = 1.3787 -0.6240 -0.7547 0 0 -6.2917 3.9002 2.6415 0 0 -0.6240 1.4539 -0.8299 0 0 3.9002 -66.9808 3.1120 63.4921 0 -0.7547 -0.8299 1.5846 0 0 2.6415 3.1120 -35.7379 0 31.7460 0 0 0 0