1、1高斯赛德尔法潮流计算潮流计算高斯赛德尔迭代法(Gauss 一 Seidel method)是求解电力系统潮流的方法。潮流计算高斯赛德尔迭代 法又分导纳矩阵迭代法和阻抗矩阵迭代法两种。前者 是以节点导纳矩阵为基础建立的赛德尔迭代格式;后者是以节点阻扰矩阵为基础建立的赛德尔迭代格式。 高斯赛德尔迭代法这是数学上求解线性或非 线性方程组的一种常用的迭代方法。本实验通过对电力网数学模型形成的计算机程序的编制与调试，获得形成电力网数学模型：高斯-赛德尔法的计算机程序，使数学模型能够由计算机自行形成，即根据已知的电力网的接线图及各支路参数由计算程序运行形成该电力网的节点导纳矩阵和各节点电压、功率。通过实
2、验教学加深学生对高斯-赛德尔法概念的理解，学会运用数学知识建立电力系统的数学模型，掌握数学模型的形成过程及其特点，熟悉各种常用应用软件，熟悉硬件设备的使用方法，加强编制调试计算机程序的能力，提高工程计算的能力，学习如何将理论知识和实际工程问题结合起来。高斯- 赛德尔法潮流计算框图N开始输入数据，定义数组给定 PQ 节点电压初值给定 PV 节点电压实部（或虚部）置迭代计数 b=0计算 PQ 节点电压实部和虚部先计算 PV 节点无功功率再用其计算 PV 节点电压实部和虚部计算平衡节点的有功和无功判断所有 | |是否0Y(2*i-1,j)=Y(2*i-1,j)-r/(r2+x2);Y(2*i,j)=
3、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;endif i*jmVm=0;for i=1:n-1j=1;A1=0;A2=0;if ij for 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*
4、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);if ihg=Y(2*i-1,i);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
5、);9A(i)=atan(F(i)/E(i);A(i)=A(i)*180/pi;continueendE(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)-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 VcVmVm=Vc;endV(i)=v;A(i)=atan(F(i)/E(i);A(i)=A(i)*180/pi;endk=k+1;
6、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);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 03.9002 -66.9808 3.1120 63.4921 0-0.7547 -0.8299 1.5846 0 02.6415 3.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.3333k =1011P =-1.6000 -2.0000 -3.7000 5.0000 0.5238Q =-0.8000 -1.0000 -1.3000 1.3885 0.5238V =0.8885 1.0817 1.0579 1.0500 1.0500A =-11.6107 -0.4133 1.1798 0.0028 0
