电力系统稳态分析大课后复习根据高斯赛德尔法潮流计算

,.电力系统稳态分析姓名:学号：学院(系):自动化学院专业:电气工程基于Matlab的高斯和高斯—赛德尔法的潮流计题 目: 算指导老师：,.2014年12月摘要电力系统潮流计算是电力系统稳态运行分析中最基本和最重要的计算之一，是电力系统其他分析计算的基础，也是电力网规划、运行研究分析的一种方法，在电力系统中具有举足轻重的作用。经典算法有高斯法，高斯－赛德尔迭代法及牛顿法等，近年来学者们开始应用非线性规划法及智能算法等优化方法求解潮流问题，提高了收敛的可靠性。感谢阅读高斯－赛德尔迭代法开始于上世纪50年代，是一种直接迭代求解方程的算法，既可以解线性方程组，可以解非线性方程组。高斯法求解节点电压的特点是:在计算节点i第k+1次的迭代电压时，前后所用的电压都是第k次迭代的结果，整个一轮潮流迭代完成后，把所有计算出的电压新值用于下一轮电压新值的计算过程中。该计算方法简单，占用计算机内存小，能直接利用迭代求解节点电压方程，对电压初值的选取要求不是很严格。但它的收敛性能较差，系统规模增大时，迭代次数急剧上升。谢谢阅读本文首先对高斯—赛德尔算法进行了综述，然后推导了该算法的计算过程，通过MATLAB软件计算了该算法的实例。感谢阅读关键字：潮流计算 高斯法 高斯－赛德尔法 迭代,.AbstractPowerflowcalculationistheoneofthemostbasicandthemostimportantcalculationinthesteadystateanalysisofpowersystem.Itisthefoundationofotheranalyticalcalculationofpowersystem,amethodofanalysisandplanning,operationofpowernetwork.Soitplaysadecisiveroleinthepowersystem.TheclassicalalgorithmistheGaussmethod,Gauss-SeideliterativemethodandNewton'smethod,inrecentyears.Scholarsbegantoapplicatenonlinearprogrammingmethodandintelligentalgorithmoptimizationmethodforsolvingpowerflowproblem,enhancesthereliabilityofconvergence.感谢阅读Gauss-Seideliterativemethodbeganinthe50'soflastcentury,isa谢谢阅读,.directiterationequationalgorithm,whichcansolvethelinearequationandnonlinearequations.CharacteristicsofGauss'smethodtocalculatethenodevoltageis:intheiterativecalculationofnodei’sK+1-timesvoltage,thevoltageisusedtheresultsofK-timesiterative.Aftercompletingthewholeroundofpowerflowiteration,allvoltagevalueisusedtocalculatethenextroundofnewvoltagevalueof.Themethodissimpleandcapturessmallmemory.Italsocandirectlyusetheiterativesolutionofthenodevoltageequation.theselectionofinitialvaluesarenotverystrict.Butithaspoorconvergenceperformance.Thesystemscaleincreases,whenthenumberofiterationsrise.精品文档放心下载ThispapergivesanoverviewoftheGaussSeidelalgorithmatthefirst.ThenitshowthecalculationprocessofthisalgorithmthroughtheMATLABsoftware.谢谢阅读Keywords:GaussGauss-Seideliterativemethodthemethodofpowerflowcalculation谢谢阅读目录1高斯迭代法和高斯—赛德尔迭代法概述 52节点导纳矩阵 6,.2.1不定导纳矩阵....................................................................2.2导纳矩阵.......................................................................高斯迭代法........................................................................谢谢阅读高斯-赛德尔迭代法..................................................................谢谢阅读4.1高斯-赛德尔法的原理............................................................4.2关于高斯法和高斯-赛德尔法的讨论................................................5实例验证..........................................................................5.1案例描述....................................................................5.2模型的建立...................................................................105.3案例程序流程图................................................................115.4案例程序....................................................................135.5程序运行步骤和结果............................................................176结果分析..........................................................................207总结..............................................................................217参考文献..........................................................................22,.一高斯迭代法和高斯—赛德尔迭代法概述电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算。它的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。电力系统潮流计算的结果是电力系统稳定计算和故障分析的基础。给定电力系统的网络结构，参数和决定系统运行状况的边界条件，电力系统的稳态运行状态便随之确定。潮流计算就是要通过数值仿真的方法把电力系统的详细运行状态呈现给运行和工作人员，以便研究系统在给定条件下的稳定运行特点。潮流计算是电力系统分析中最基本、最重要的计算，是电力系统运行、规划以及安全性、可靠性分析和优化的基础，也是各种电磁暂态和机电暂态分析的基础和出发点。谢谢阅读20世纪50年代中期，随着电子计算机的发展，人们开始在计算机上用数学模拟的方法进行潮流计算。最初在计算机上实现的潮流计算方法是以导纳矩阵为基础的高斯迭代法（Gauss法）。这种方法内存需求小，但收敛性差。后来在高斯迭代法上进行改进，这就是高斯——赛德尔迭代法(Gauss一Seidelmethod)，潮流计算高斯—赛德尔迭代法，分为导纳矩阵迭代法和阻抗矩阵迭代法两种。前者是以节点导纳矩阵为基础建立的赛德尔迭代格式，后者是以节点阻扰矩阵为基础建立的赛德尔迭代格式。高斯——赛德尔迭代法这是数学上求解线性或非线性方程组的一种常用的迭代方法。感谢阅读,.牛顿-拉夫逊方法是解非线性代数方程组的一种基本方法，在潮流计算中也得到了应用。20世纪60年代中后期，系数矩阵技术和编号优化技术的提出使牛顿-拉夫逊的解题规模和计算效率进一步提高，至今仍是潮流计算中的广泛采用的优秀算法。谢谢阅读20世纪70年代中期，Stott在大量计算实践的基础上提出了潮流计算的快速分解法，是潮流计算的速度大大提高，可以应用于在线，但是直至20世纪谢谢阅读年代末期才对快速分解法潮流的收敛性给出了比较满意的解释。精品文档放心下载由于潮流计算在电力系统中的特殊地位和作用，对其计算方法有如下较高精品文档放心下载的要求：要有可靠的收敛性，对不同的系统及不同的运行条件都能收敛；精品文档放心下载占用内存小、计算速度快；调整和修改容易，使用灵活方便。本文使用的高斯法和高斯－赛德尔迭代法，开始于上世纪50年代，是一种谢谢阅读直接迭代求解方程的算法，既可以解线性方程组，可以解非线性方程组。高斯法求解节点电压的特点是:在计算节点i第k+1次的迭代电压时，前后所用的电压都是第k次迭代的结果，整个一轮潮流迭代完成后，把所有计算出的电压新值用于下一轮电压新值的计算过程中。高斯-赛德尔法是刚刚计算出的x值在下次迭代中被立即使用。两种方法都计算方法简单，占用计算机内存小，能直接利用迭代求解节点电压方程，对电压初值的选取要求不是很严格，但收敛性能较差，系统规模增大时，迭代次数急剧上升。谢谢阅读二 节点导纳矩阵不定导纳矩阵,.令连通的电力网络的节点数为N，大地作为节点未包括在内。网络中有b条支路，包括接地支路。如果把地节点增广进来，电网的（N+1）×b阶节点支路的关联矩阵A0，b阶支路导纳矩阵是yb，定义（N+1）×（N+1）阶节点导精品文档放心下载纳矩阵Y0为YAyAT(2-1)00b0并有网络方程..(2-2)YVI00导纳矩阵选地节点为电压参考点，将它排在第N+1位，令参考点点位为零，则可将节点不定导纳矩阵表示的网络方程(2-1)写成分块的形式精品文档放心下载Yy0VI][yoTy00][0]=[I0（2-3）展开后有YV=I（2-4）和yTV=I0(2-5)0式（2-5）中Y为N×N阶矩阵，V和I分别为N维节点电压和电流列矢量，精品文档放心下载I0为流入地节点的电流。三 高斯迭代法高斯迭代法是最早在计算机上实现的潮流计算方法。这种方法编程简单，在精品文档放心下载某些应用领域，如配电网计算潮流计算中还有应用。另外，也用于为牛顿-拉夫谢谢阅读逊法提供初值。,.考察基于节点导纳矩阵的高斯迭代法。在网络方程（2-4）中，将平衡点s排感谢阅读在最后，并将导纳矩阵写成分块的形式，取出前n个方程有谢谢阅读(3-1)YnVn+YsVs=In平衡节点s的电压V给定，n个节点的注入电流矢量I已知，则有sn(3-2)YnVn=In−YsVs实际电力系统给定量是n个节点的注入功率。注入电流和注入功率之间的关系是谢谢阅读Ii=̃i=1,2,3….n(3-3) i i i iIn=[](3-4)i精品文档放心下载n可以得到0Yn=L+D+UY11其中L=D=Y⋱,Y,⋮0][Y[Y…YY44]0Y…YU=[Y]0代入式（2-2），经过整理可得到V=D−1{I−YV−LV−UV}（3-5）nnssnn考虑到电流和功率的关系式，（3-5）可以写成为V（k+1）=1{̃iiVi

−YV−∑i−1YV(k)−∑nYV(k)}i=1,2,…nissj=1ijjj=i+1ijj(3-6),.给定，i=1，2,…n，代入上式中可得电压新值，逐次迭代直到前后两次迭精品文档放心下载代求得的电压值的差小于某一收敛精度为止。这是高斯迭代法的基本解算步骤。谢谢阅读四高斯—赛德尔迭代法高斯-赛德尔法的原理每次迭代要从节点1扫描到节点n。在计算V（k+1）时，V（k+1），j=感谢阅读i j1,2,…,i−1已经求出，若若迭代是一个收敛过程，它们应比Vj（k），j=1,2,…,i−谢谢阅读谢谢阅读j i高斯—赛德尔迭代的思想，即一旦求出电压新值，在最后的迭代中立即使用。精品文档放心下载这种方法的迭代公式是V（k+1）=1i Yii

̃−YisVsi−1V(k+1){S−∑Yijj=1Vi

n−∑YijV(k)}i=1,2,…nj=i+1j高斯—赛德尔法比高斯迭代法的收敛性好。关于高斯法和高斯-赛德尔法的讨论对于形如 f(x)=0的非线性代数方程组，总可以写成x=φ(x)的形式，于是，有如下的高斯迭代公式：{x（0）=x0x（k+1）=φ（x（k））高斯迭代法的收敛性主要由ϕ（x∗）≝∂φ(x)|∂xT x=x∗

（4-1）（4-2）（4-3）（4-4）的谱半径决定。x∗是x的解点。当ϕ（x∗）的谱半径小于1时，高斯迭代法感谢阅读可以收敛，ϕ（x∗）的谱半径越小高斯迭代法的收敛性越好。感谢阅读求解式（4-3）有高斯法和高斯-赛德尔法。高斯法的迭代过程为感谢阅读,.x(k)=φ(x(k),x(k),…x(k))i=1,2…n(4-5)i高斯-赛德尔法的迭代公式是x(k1)=φ(k1),x(k1),…x(k1)，(k)(k))i=1,2…n(4-6)i即刚刚计算出的x值在下次迭代中被立即使用，当max(|x(k1) x(k)|)<，时，迭代收敛。感谢阅读 对于连通的电力网络，各节点的电压是相关的，而不管两个节点之间是否有谢谢阅读支路直接相连。由于Y矩阵是高度稀疏的，由高斯迭代法的公式（4-5）可见，谢谢阅读计算节点i的电压时，只有和节点i有支路直接相连的节点j的电压对V有贡献。感谢阅读这种方法在迭代修正时利用的信息较少，收敛性较差，其优点是内存需求较少。感谢阅读五 实例验证案例描述（有变压器支路的情况）如图所示的一个三母线电力系统，在母线①和母线精品文档放心下载③之间的输电线的母线③端连接着一个纵向串联加压器，可在同一电压等级改变谢谢阅读电压幅值。该系统的网络元件用图（2）所示的等值电路表示，串联之路用电阻感谢阅读和电抗表示，并联支路用电纳表示。支路（1,3）用一个变比可调的等值变压器谢谢阅读之路表示，非标准变比t=1.05，在节点①侧。试形成该网络的节点导纳矩阵。精品文档放心下载假定节点①的注入功率S1=2.0 j1.0，节点②的注入功率是S2=0.5 j0.415，感谢阅读节点③是Vθ节点，V3=1.0。试基于节点导纳矩阵的高斯法和高斯-赛德尔法计精品文档放心下载算潮流。,.图1三母线电力系统 图2等值电路图3电路的相关参数模型建立首先对支路编号并规定串联之路的正方向如图（3）所示，则可得到广义节点支路关联矩阵A。A中行与节点对应，列与支路对应。A矩阵为谢谢阅读−1−1/t0100A=[10−1010]011001支路导纳矩阵为Yn=diag[y12y13y23y10y20y30]感谢阅读0.2494−4.9875i0000000.9901−9.9010i0000=000.49505−4.9505i0000000.01i00[00000.03i0000000.02i],.建立节点导纳矩阵如下：Y=AybAT根据节点导纳矩阵可写出式（2-5）的表达式V1]=D−1{In−YsVs−LVn−UVn}V=[V2Y0Y00Yns3s0Ỹ000,Y,L=[Y210]U=[00]当使用高斯迭代法进行计算时，其基本的迭代公式如下：V(k+1)=(1.1474−13.958i)−1{−2+1i+(0.2494−4.9875i)V(k)+0.9430−9.430i}VV(k+1)=(0.74445−9.908i)−1{0.5−0.415i+(0.2494−4.9875i)V(k)+0.49505−4.9505i}V2i=1,2,…n将上式写成简单迭代法的高斯法迭代格式为V(k+1)=f(V(k),V(k))V=f(V(k),V(k))2212当使用高斯—赛德尔迭代法进行计算时，其基本的迭代公式如下谢谢阅读V(k+1)=(1.1474−13.958i)−1{−2+1i+(0.2494−4.9875i)V(k)+0.9430−9.430i}V1V(k+1)=(0.74445−9.908i)−1{0.5−0.415i+(0.2494−4.9875i)V(k+1)+0.49505V2−4.9505i},.i=1,2,…n将上式写成简单迭代法的高斯法迭代格式为V=f((k),(k))VV高斯—赛德尔迭代法与高斯法的不同点在与，算出一个V(k+1)后，立即在下感谢阅读一步运算中应用。案例程序流程图报告基于matlab软件编程实现上述电力系统潮流计算，采用高斯迭代法感谢阅读和高斯-赛德尔迭代法。程序流程图如下所示：开始输入原始数据——节点个数N、节点之路关联矩阵A，支路导纳矩阵形成节点导纳矩阵Y给定节点注入功率S、节点电压初值V0、平衡节点电Vs给定判断迭代结束值e,.高斯法V(k1)=f(V(k),V(k))Fact=|max({Vi(k1)−Vi(k)}|感谢阅读NoFact<高斯—赛德尔法V(k1)=f(V(k1),V(k))Fact1=|max({Vi(k1)−Vi(k)}|感谢阅读Fact1<No结果输出（迭代次数，迭代过程）结束案例程序程序如下：clear;,.formatshort;N=input('请输入节点个数：N=');A=input('请输入广义节点支路关联矩阵：A=\n');A%显示关联矩阵感谢阅读Yb=input('请输入各支路导纳:\n');谢谢阅读n=2*N;YY=zeros(n);%计算节点导纳矩阵fori=1:(n*n)m=floor(i/n);if(i==(n*m+m+1))YY(i)=Yb(m+1);endYY(n*n)=Yb(n);endY1=A*YY;disp('节点导纳矩阵为:');Yn=Y1*(A.')DD=diag(Yn);D=zeros(N);fori=1:N*Nm=floor(i/N);,.if(i==(m*N+m+1))D(i)=DD(m+1);endD(N*N)=DD(N);%对角矩阵DendL=tril(Yn)-D;%下三角矩阵LU=triu(Yn)-D;%上三角矩阵ULL=zeros(N-1);DD=zeros(N-1);UU=zeros(N-1);LL=L(1:(N-1),1:(N-1));DD=D(1:(N-1),1:(N-1));UU=U(1:(N-1),1:(N-1));%取矩阵n-1*n-1感谢阅读%输入相关的信息S=input('请输入节点1到节点n-1的注入功率S=\n');V0=input('请给定节点电压初值V0=\n');精品文档放心下载Vs=input('请输入平衡节点N的电压Vs=\n');谢谢阅读e=input('请输入收敛判断值:e=\n');精品文档放心下载Vk=zeros(N-1,1);Vk1=zeros(N-1,1);Vk=V0;D1=inv(DD);,.Ys=Yn(N,1:N-1);Ys0=(Ys*Vs).';I=zeros(1,2);Vl=zeros(N-1,1);Vl1=zeros(N-1,1);Vl1=V0;Vll1(1,1:2)=(V0);%记录每次迭代过程中各节点电压值精品文档放心下载Vll2(1,1)=0;j=1;l=1;%高斯迭代法while(j==1)%循环计算l=l+1;Vl=Vl1;I=(S')./(Vl');Lk=(LL)*(Vl.');Uk=(UU)*(Vl.');Vl11=D1*(I-Ys0-Lk-Uk);%根据高斯迭代法感谢阅读Vl1=Vl11.';FVl=Vl1-Vl;F=max(FVl);Fact=sqrt((real(F))^2+(imag(F))^2);谢谢阅读,.Vll1(l,1:2)=(Vl1);Vll2(l,1)=Fact;P=[Vll1,Vll2];ifFact<e %判断误差值是否满足条件j=0; %若满足条件，设置标志位，停止迭代endendk=1;V1(1,1:2)=(V0);%记录每次迭代过程中各节点电压值感谢阅读V2(1,1)=0;jj=1;FVk=zeros(1,2);%高斯赛德尔迭代while(jj==1)%循环计算k=k+1;Vk1=Vk;I1=conj(S(1))/conj(Vk(1));谢谢阅读Lk1=LL(1,1)*Vk(1)+LL(1,2)*Vk(2);谢谢阅读Uk1=UU(1,1)*Vk(1)+UU(1,2)*Vk(2);感谢阅读Vk(1)=D1(1,1)*(I1-Ys0(1)-Lk1-Uk1);精品文档放心下载I2=conj(S(2))/conj(Vk(2));感谢阅读Lk2=LL(2,1)*Vk(1)+LL(2,2)*Vk(2);谢谢阅读,.Uk2=UU(2,1)*Vk(1)+UU(2,2)*Vk(2);谢谢阅读Vk(2)=D1(2,2)*(I2-Ys0(2)-Lk2-Uk2);谢谢阅读FVk(1)=Vk(1)-Vk1(1);FVk(2)=Vk(2)-Vk1(2);FVk=max(FVk.');Fact1=sqrt((real(FVk))^2+(imag(FVk))^2);感谢阅读V1(k,1:2)=(Vk);V2(k,1)=Fact1;V=[V1,V2];ifFact1<e %判断误差值是否满足条件jj=0; %若满足条件，设置标志位，停止迭代endendfprintf('高斯迭代次数%d\n',(l-1));%输出结果精品文档放心下载A=[0:l-1].';P=[A,Vll1,Vll2];disp('整个迭代过程如下：');disp(' 迭代次数 V1(k)

V2(k)max(V(k+1)-V(k))');disp(P);fprintf('高斯赛德尔迭代次数%d\n',(k-1));%输出结果A=[0:k-1].';感谢阅读,.V=[A,V1,V2];disp('整个迭代过程如下：');disp(' 迭代次数 V1(k) V2(k)max(V(k+1)-V(k))');disp(V);%绘制精度曲线P=[1:l-1];plot(P,Vll2(2:l));holdon;P1=[1:k-1];plot(P1,V2(2:k),'-.');holdoff;title('高斯法和高斯赛德尔法的精度曲线')legend('高斯迭代法','高斯-赛德尔迭代法')感谢阅读xlabel('迭代次数')ylabel('精度')grid程序运行步骤和结果程序的运行步骤如下：输入节点个数：N=3，输入广义节点支路关联矩阵：[-1-1/1.050100;10-1010;011001]，精品文档放心下载和支路导纳矩阵,.[0.2494-4.9875i,0.9901-9.9010i,0.49505-4.9505i,0.01i,0.03i,0.02i]后运行结果如下谢谢阅读导纳矩阵的计算结果：输入节点1到节点2的注入功率[-2-1i,0.5+0.415i]。精品文档放心下载输入节点1和2电压初值[11]。输入平衡节点N的电压Vs=1。输入收敛判断值:e=0.00001后，运行的结果如下：精品文档放心下载输入的参数高斯迭代法的计算结果,.高斯赛德尔的计算结果,.再将计算出来每次迭代出来的精度制成精度曲线，曲线如下：精品文档放心下载高斯法和高斯赛德尔法的精度曲线0.16高斯迭代法0.14高斯-赛德尔迭代法0.120.1度0.08精0.060.040.02002468101214迭代次数六结果分析通过对上面的数据整理分析，可以看出，节点导纳的计算结果为精品文档放心下载1.1474-13.9580i -0.2494+4.9875i -0.9430+9.4295i精品文档放心下载-0.2494+4.9875i0.7445-9.9080i-0.4950+4.9505i-0.9430+9.4295i-0.4950+4.9505i1.4852-14.8315i与书上的例2.3结果相同。当使用高斯法计算书上的例7.1后，通过matlab的仿