MATLAB在电力系统潮流计算中的应用.docxVIP

MATLAB在电力系统潮流计算中的应用 电气与自动化工程学院 电气工程及其自动化专业 2012级 蒋志涛 背景 电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算。它的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。电力系统潮流计算的结果是电力系统稳定计算和故障分析的基础。在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。 电力系统潮流计算属于稳态分析范畴，不涉及系统元件的动态特性和过渡过程。因此其数学模型不包含微分方程，是一组高阶非线性方程。非线性代数方程组的解法离不开迭代，因此，潮流计算方法首先要求它是能可靠的收敛，并给出正确答案。随着电力系统规模的不断扩大，潮流问题的方程式阶数越来越高，目前已达到几千阶甚至上万阶，对这样规模的方程式并不是采用任何数学方法都能保证给出正确答案的。这种情况促使电力系统的研究人员不断寻求新的更可靠的计算方法。 MATLAB相关功能和模块介绍 本程序以m文件的形式编写，由外部文件（excel文件）将所需数据输入，由程序读入，进行计算并进行迭代，当数据精度符合要求的时候，停止迭代，输出结果将显示在下方窗口中，并输出到文件夹中的外部文件（excel）中，计算过程结束。 程序设计 程序代码： %%Data Input N=xlsread(DATA,1,B2); R=xlsread(DATA,1,B3:Z3); X=xlsread(DATA,1,B4:Z4); B=xlsread(DATA,1,B5:Z5); P=xlsread(DATA,1,B6:Z6); Q=xlsread(DATA,1,B7:Z7); KT=xlsread(DATA,1,B8:Z8); P0=xlsread(DATA,1,B9:Z9); Q0=xlsread(DATA,1,B10:Z10); Un=xlsread(DATA,1,B11); U1=xlsread(DATA,1,B12); %%Define Starting Value S=zeros(1,N+1); S2=zeros(1,N); S1=zeros(1,N); U=zeros(1,N+1); Ux=zeros(1,N); Uy=zeros(1,N); UH=zeros(1,N+1); %voltage that convert to high voltage side Z=zeros(1,N); SLD=zeros(1,N+1); S0=zeros(1,N); for n=1:1:N Z(n)=R(n)+1i*X(n); end for n=1:1:N+1 SLD(n)=P(n)+1i*Q(n); end for n=1:1:N S0(n)=P0(n)+1i*Q0(n); end %%Calculate Point Load S(N+1)=SLD(N+1); for n=2:1:N S(n)=SLD(n)-1i*(B(n-1)+B(n))*(Un^2)/2; end %%Calculate back and forth to get more accurate Un Delt=1; S2(N)=SLD(N+1); K=0; U(1)=U1; while Delt0.001 K=K+1; for n=N:-1:2 %calculate back to get S(1) S1(n)=S2(n)+(real(S2(n))^2+imag(S2(n))^2)*Z(n)/(Un^2); S2(n-1)=S(n)+S0(n)+S1(n); end S1(1)=S2(2)+(real(S2(2))^2+imag(S2(2))^2)*Z(2)/(Un^2); S(1)=S1(1); for n=1:1:N %calculate forth to get new Un Uy(n)=(real(S1(n))*real(Z(n))+imag(S1(n))*imag(Z(n)))/U(n); Ux(n)=(real(S1(n))*imag(Z(n))-imag(S1(n))*real(Z(n)))/U(n);