资源描述：
《matlab在电力系统潮流计算中的应用》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、MATLAB在电力系统潮流计算中的应用电气与自动化工程学院电气工程及其自动化专业2012级蒋志涛背景电力系统潮流计算是研究电力系统稳态运行情况的i种基本电气计算。它的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。电力系统潮流计算的结果是电力系统稳定计算和故障分析的基础。在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮
2、流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控屮，则采用在线潮流计算。电力系统潮流计算属于稳态分析范畴,不涉及系统元件的动态特性和过渡过程。因此其数学模型不包含微分方程，是一组高阶非线性方程。非线性代数方程组的解法离不开迭代，因此，潮流计算方法首先要求它是能可靠的收敛，并给出正确答案。随着电力系统规模的不断扩大，潮流问题的方程式阶数越来越高，目前已达到几千阶甚至上万阶，对这样规模的方程式并不是采用任何
3、数学方法都能保证给出止确答案的。这种情况促使电力系统的研究人员不断寻求新的更可靠的计算方法。二.MATLAB相关功能和模块介绍本程序以m文件的形式编写，由外部文件（excel文件）将所需数据输入，由程序读入，进行计算并进行迭代，当数据精度符合要求的时候，停止迭代，输出结果将显示在下方窗口中，并输出到文件夹中的外部文件（excel）中，计算过程结朿。三.程序设计程序代码：%%DataInputN=xlsread('DATA',1,'B2');R二xlsread('DATA',1,'B3:Z3');X=xlsread
4、(，DATA',1,'B4:Z4');B^xlsreadCDATA',1,'B5:Z5');P二xlsread('DAT/V,1,，B6:Z6,);Q二xlsread('DATA',1,'B7:Z7');KT二xlsread('DATA',1,，B8:Z8，);PO=xlsreadCDATA',1,'B9:Z9');QO=xlsread(，DATA，,1,'B10:Z10');Un=xlsread(,DATA',1,'Bll');Ul=xlsread(，DATA',1,'B12');%%DefineStartin
5、gValueS=zeros(1,N+l);S2二zeros(1,N);Sl=zeros(1,N);U二zeros(1,N+l);Ux二zeros(1,N);Uy二zeros(1,N);UH二zeros(1,N+l);%voltagethatconverttohighvoltagesideZ二zeros(1,N);SLD=zeros(1,N+l);SO二zeros(1,N);forn=l:1:NZ(n)=R(n)+li*X(n);endforn=l:1:N+1SLD(n)=P(n)+li*Q(n);endforn=
6、l:1:NS0(n)=P0(n)+li*Q0(n);end%%CalculatePointLoadS(N+1)=SLD(N+1);forn=2:1:NS(n)=SLD(n)-li*(B(n-l)+B(n))*(UrT2)/2;end%%CaleulatebackandforthtogetmoreaccurateUnDelt=l;S2(N)=SLD(N+1);K二0;U⑴二Ul;whileDelt>0.001K二K+l;forn=N:-1:2%calculatebacktogetS(1)51(n)=S2(n)+(r
7、eal(S2(n))2+imag(S2(n))2)*Z(n)/(Un2);52(n-1)=S(n)+S0(n)+S1(n);endSl(l)=S2(2)+(rea1(S2(2)厂2+imag(S2(2))"2)*Z(2)/(Un"2);S(1)=S1(1);forn=l:1:N%calculateforthtogetnewUnUy(n)=(real(S1(n))*real(Z(n))+imag(SI(n))*imag(Z(n)))/U(n);Ux(n)=(real(SI(n))*imag(Z(n))-imag(S
8、I(n))*real(Z(n)))/U(n);UH(n+l)=sqrt((U(n)-Uy(n)厂2+Ux(n厂2);ifKT(n)==0U(n+l)=sqrt((U(n)-Uy(n)厂2+Ux(n厂2);elseU(n+l)=sqrt((U(n)-Uy(n))2+Ux(n)2)/KT(n);endendDelt=abs(UH(N+l)-Un);Un=UH(N+l);end
