电力系统潮流计算的MATLAB辅助程序设计潮流计算程序

1、电力系统潮流计算的MATLAB辅助程序设计潮流计算，通常指负荷潮流，是电力系统分析和设计的主要组成部分，对系 统规划、安全运行、经济调度和电力公司的功率交换非常重要。此外，潮流计算 还是其它电力系统分析的基础，比如暂态稳定，突发事件处理等。现代电力系统 潮流计算的方法主要:高斯法、牛顿法、快速解耦法和MATLAB的M语言编写的 MATPOWER4.1，这里主要介绍高斯法、牛顿法和快速解耦法。高斯法的程序是 lfgauss，其与lfybus、busout和lineflow程序联合使用求解潮流功率lfybus、 busout和lineflow程序也可与牛顿法的 lfnewton程序和快速解耦法的

2、decouple程序联合使用。（读者可以到MATPOWER主页下载MATPOWER4.1,然后将 其解压到MATLAB目录下，即可使用该软件进行潮流计算）一、高斯-赛德尔法潮流计算使用的程序：高斯-赛德法的具体使用方法读者可参考后面的实例，这里仅介绍各程序的 编写格式：lfgauss：该程序是用高斯法对实际电力系统进行潮流计算，需要用到 busdata和linedata两个文件。程序设计为输入负荷和发电机的有功MW和无功 Mvar，以及节点电压标幺值和相角的角度值。根据所选复功率为基准值将负荷和 发电机的功率转换为标幺值。对于PV节点，如发电机节点，要提供一个无功功 率限定值。当给定电压过高或

3、过低时，无功功率可能超出功率限定值。在几次迭 代之后（高斯-塞德尔迭代为10次），需要检查一次发电机节点的无功出力，如果 接近限定值，电压幅值进行上下5%的调整，使得无功保持在限定值内。lfybus：这个程序需要输入线路参数、变压器参数以及变压器分接头参数。 并将这些参数放在名为linedata的文件中。这个程序将阻抗转换为导纳，并得 到节点导纳矩阵。busout：该程序以表格形式输出结果，节点输出包括电压幅值和相角，发电 机和负荷的有功和无功功率，以及并联电容器或电抗器的有功和无功功率。lineflow:该程序输出线路的相关数据，程序设计输出流入线路终端的有功 和无功的功率、线损以及节点功率

4、，还包含整个系统的有功和无功损耗。lfnewton是牛顿-拉夫逊法对实际电力系统潮流计算开发的程序，数据准备 和程序格式和高斯-赛德尔法一样，包括程序lfybus，busout和lineflowodecouple是快速解耦法对实际电力系统潮流计算开发的程序，同高斯法和 牛顿法一样需要用到三个程序：lfybus、busout、lineflowo二、数据准备为了在MATLAB环境下用高斯法进行潮流计算，必须定义下列变量：基准功 率，功率允许误差，加速因子和最大迭代次数。上述变量命名（小写字母）为： basemva、accuracy、accel 和 maxiter，一般规 定为： basemva=1

5、00 ； accuracy=0.001； accel=1.6； maxiter=80；输入文件准备的第一步是给节点编 号，节点号码必须是连续的，但节点数据输入不一定按顺序来编写。此外，还需 要下列数据文件：节点数据文件busdata :节点信息输入格式为单行输入，输入的数据形成一个 矩阵，叫做busdata矩阵。第一列为节点号；第二列为节点类型；第三列和第四 列分别为节点电压幅值（标幺值）和相角（单位为度）；第五列和第六列分别为负荷 的有功功率和无功功率；第七列到十列分别为发电机的有功功率、无功功率、最 小无功出力和最大无功出力；最后一列为并联电容器注入无功功率。第二列的编 码用0、1、2来区

6、分PQ节点、平衡节点和PV节点：0表示PQ节点，输入正的有功功率（MW）和无功功率（Mvar），并且要设定 节点电压初始估计值，一般幅值和相角分别设为1和0，若已经给定初始值，则 用其给定值来代替1和0。1表示平衡节点，且已知该节点的电压幅值和相角。2表示PV节点，要设定该节点的节点电压幅值和发电机的有功功率（MW）， 并设定发电机的无功最小出力和最大出力（Mvar）。线路数据文件linedata线路数据用节点对的方法来确定，数据包含在称为 linedata的矩阵中。第一列和第二列为节点号码，第三列到第五列为线路电阻、 电抗及该线路电纳值的一半，以标幺值表示。最后一列为变压器分接头设定值， 对

7、线路来说，需要输入1。线路输入为无输入顺序，对变压器来说，左侧的节点 号设为分接头端。zdata是线路数据输入变量，包括四项，前两项是节点编号，后两项是线路电 阻和电抗，均以标幺值表示，函数返回节点导纳矩阵。三、潮流计算的MATLAB程序清单1. lfgauss.m程序清单% Power flow solution by Gauss-Seidel method Vm=0; delta=0; yload=0; deltad =0;nbus = length（busdata（:,1）;kb=;Vm=; delta=; Pd=; Qd=; Pg=; Qg=; Qmin=; Qmax=;Pk=; P=

8、; Qk=; Q=; S=; V=;for k=1:nbusn=busdata(k,1);kb(n)=busdata(k,2); Vm(n)=busdata(k,3); delta(n)=busdata(k, 4);Pd(n)=busdata(k,5); Qd(n)=busdata(k,6); Pg(n)=busdata(k,7); Qg(n) = busdata(k,8);Qmin(n)=busdata(k, 9); Qmax(n)=busdata(k, 10);Qsh(n)=busdata(k, 11);if Vm(n) = accuracy & iter = maxiter iter=i

9、ter+1;for n = 1:nbus; YV = 0+j*0;for L = 1:nbr;if (nl(L) = n & mline(L) = 1), k=nr(L); YV = YV + Ybus(n,k)*V(k);elseif (nr(L) = n & mline(L)=1), k=nl(L); YV = YV + Ybus(n,k)*V(k);end end Sc = conj(V(n)*(Ybus(n,n)*V(n) + YV); Sc = conj(Sc); DP(n) = P(n) - real(Sc); DQ(n) = Q(n) - imag(Sc); if kb(n) =

10、 1 S(n) =Sc; P(n) = real(Sc); Q(n) = imag(Sc); DP(n) =0; DQ(n)=0; Vc(n) = V(n); elseif kb(n) = 2 Q(n) = imag(Sc); S(n) = P(n) + j*Q(n);if Qmax(n) = 0 Qgc = Q(n)*basemva + Qd(n) - Qsh(n); if abs(DQ(n) = 10 if DV(n) = 0.045 if Qgc Qmax(n), Vm(n) = Vm(n) - 0.005; DV(n)=DV(n)+.005; end else, end else,en

11、d else,end end if kb(n) = 1 Vc(n) = (conj(S(n)/conj(V(n) - YV )/ Ybus(n,n); else, end if kb(n) = 0 V(n) = V(n) + accel*(Vc(n)-V(n); elseif kb(n) = 2 Vcl = imag(Vc(n); VcR = sqrt(Vm(n)”2 - Vc2);Vc(n) = VcR + j*VcI;V(n) = V(n) + accel*(Vc(n) -V(n); endendmaxerror=max( max(abs(real(DP), max(abs(imag(DQ

12、); if iter = maxiter & maxerror accuracyfprintf(nWARNING: Iterative solution did not converged after ) fprintf(%g, iter), fprintf( iterations.nn)fprintf(Press Enter to terminate the iterations and print the results n) converge = 0; pause, else, endendif converge = 1tech= (ITERATIVE SOLUTION DID NOT

13、CONVERGE); elsetech=(Power Flow Solution by Gauss-Seidel Method);endk=0;for n = 1:nbusVm(n) = abs(V(n); deltad(n) = angle(V(n)*180/pi;if kb(n) = 1S(n)=P(n)+j*Q(n);Pg(n) = P(n)*basemva + Pd(n);Qg(n) = Q(n)*basemva + Qd(n) - Qsh(n);k=k+1;Pgg(k)=Pg(n);elseif kb(n) =2k=k+1;Pgg(k)=Pg(n);S(n)=P(n)+j*Q(n);

14、Qg(n) = Q(n)*basemva + Qd(n) - Qsh(n); endyload(n) = (Pd(n)- j*Qd(n)+j*Qsh(n)/(basemva*Vm(n)”2); endPgt = sum(Pg); Qgt = sum(Qg); Pdt = sum(Pd); Qdt = sum(Qd); Qsht = sum(Qsh); busdata(:,3)=Vm; busdata(:,4)=deltad;clear AcurBus DP DQ DV L Sc Vc VcI VcR YV converge delta2.lfybus.m程序清单% This program o

15、btains the Bus Admittance Matrix for power flow solution j=sqrt(-1); i = sqrt(-1);nl = linedata(:,1); nr = linedata(:,2); R = linedata(:,3);X = linedata(:,4); Bc = j*linedata(:,5); a = linedata(:, 6); nbr=length(linedata(:,1); nbus = max(max(nl), max(nr);Z = R + j*X; y= ones(nbr,1)./Z;%支路导纳for n = 1

16、:nbrif a(n) =0 a(n) = 1; else endYbus=zeros(nbus,nbus); % 将Ybus初始化为0%非对角元素的数值Ybus(nl(k),nr(k)=Ybus(nl(k),nr(k)-y(k)/a(k); Ybus(nr(k),nl(k)=Ybus(nl(k),nr(k);end end% 对角元素的数值 for n=1:nbusfor k=1:nbr if nl(k)=n Ybus(n,n) = Ybus(n,n)+y(k)/(a(k)2) + Bc(k); elseif nr(k)=n Ybus(n,n) = Ybus(n,n)+y(k) +Bc(k)

17、; else, end end endclear Pgg3. busout.m程序清单% This program prints the power flow solution in a tabulated form % on the screen.disp(tech)fprintf(Maximum Power Mismatch = %g n, maxerror)fprintf(No. of Iterations = %g nn, iter)head = Bus Voltage Angle LoadGeneration Injected,No. Mag. Degree MW Mvar MW M

18、var Mvar;disp(head)for n=1:nbusfprintf( %5g, n), fprintf( %7.3f, Vm(n),fprintf( %8.3f, deltad(n), fprintf( %9.3f, Pd(n),fprintf( %9.3f, Qd(n), fprintf( %9.3f, Pg(n), fprintf( %9.3f , Qg(n), fprintf( %8.3fn, Qsh(n)endfprintf( n), fprintf(fprintf( %9.3f, Pdt), fprintf(fprintf( %9.3f, Pgt), fprintf()To

19、tal%9.3f, Qdt),%9.3f, Qgt), fprintf( %9.3fnn, Qsht)4.lineflow.m程序清单% This program is used in conjunction with lfgauss or lfNewton % for the computation of line flow and line losses.SLT = 0;fprintf(n)fprintf(Line Flow and Losses nn)fprintf( -Line-Power at bus & line flow-Line loss-Transformern()fprin

20、tf(from toMWMvarMVAMWMvartapn()for n = 1:nbusbusprt = 0;for L = 1:nbr;if busprt =0fprintf(n),fprintfC %6g, n),fprintf(%9.3f(, P(n)*basemva)fprintf(%9.3f, Q(n)*basemva), fprintf(%9.3fn, abs(S(n)*basemva)busprt = 1;else, endif nl(L)=n k = nr(L);In = (V(n) - a(L)*V(k)*y(L)/a(L)2 + Bc(L)/a(L)2*V(n);Ik =

21、 (V(k) - V(n)/a(L)*y(L) + Bc(L)*V(k);Snk = V(n)*conj(In)*basemva;Skn = V(k)*conj(Ik)*basemva;SL = Snk + Skn;SLT = SLT + SL;elseif nr(L)=n k = nl(L);In = (V(n) - V(k)/a(L)*y(L) + Bc(L)*V(n);Ik = (V(k) - a(L)*V(n)*y(L)/a(L)2 + Bc(L)/a(L)2*V(k);Snk = V(n)*conj(In)*basemva;Skn = V(k)*conj(Ik)*basemva;SL

22、 = Snk + Skn;SLT = SLT + SL;else, endif nl(L)=n | nr(L)=nfprintf(%12g, k), fprintf(%9.3f, real(Snk), fprintf(%9.3f, imag(Snk)fprintf(%9.3f, abs(Snk),fprintf(%9.3f, real(SL), if nl(L) =n & a(L) = 1 fprintf(%9.3f, imag(SL), fprintf(%9.3fn, a(L) else, fprintf(%9.3fn, imag(SL) end else, end end end SLT

23、= SLT/2;fprintf( n), fprintf( Total loss)fprintf(%9.3f, real(SLT), fprintf(%9.3fn, imag(SLT) clear Ik In SL SLT Skn Snk5.lfnewton.m程序清单% Power flow solution by Newton-Raphson methodns=0; ng=0; Vm=0; delta=0; yload=0; deltad=0;nbus = length(busdata(:,1);kb=;Vm=; delta=; Pd=; Qd=; Pg=; Qg=; Qmin=; Qma

24、x=;Pk=; P=; Qk=; Q=; S=; V=;for k=1:nbusn=busdata(k,1);kb(n)=busdata(k,2); Vm(n)=busdata(k,3); delta(n)=busdata(k, 4);Pd(n)=busdata(k,5); Qd(n)=busdata(k,6); Pg(n)=busdata(k,7); Qg(n) = busdata(k,8);Qmin(n)=busdata(k, 9); Qmax(n)=busdata(k, 10);Qsh(n)=busdata(k, 11);if Vm(n) = accuracy & iter = maxi

25、terfor ii=1:mfor k=1:mA(ii,k)=0;%初始化雅可比矩阵end, enditer = iter+1;for n=1:nbusnn=n-nss(n);lm=nbus+n-ngs(n)-nss(n)-ns;J11=0; J22=0; J33=0; J44=0;for ii=1:nbrif mline(ii)=1if nl(ii) = n | nr(ii) = nif nl(ii) = n , l = nr(ii); endif nr(ii) = n , l = nl(ii); endJ11=J11+ Vm(n)*Vm(l)*Ym(n,l)*sin(t(n,l)- delt

26、a(n) + delta(l);J33=J33+ Vm(n)*Vm(l)*Ym(n,l)*cos(t(n,l)- delta(n) + delta(l); if kb(n)=1J22=J22+ Vm(l)*Ym(n,l)*cos(t(n,l)- delta(n) + delta(l);J44=J44+ Vm(l)*Ym(n,l)*sin(t(n,l)- delta(n) + delta(l); else, end if kb(n) = 1 & kb(l) =1lk = nbus+l-ngs(l)-nss(l)-ns;ll = l -nss(l);% J1的非对角元素A(nn, ll) =-Vm

27、(n)*Vm(l)*Ym(n,l)*sin(t(n,l)- delta(n) + delta(l);if kb(l) = 0 % J2的非对角元素 A(nn, Ik) =Vm(n)*Ym(n,l)*cos(t(n,l)-delta(n) + delta(l); end if kb(n) = 0 % J3的非对角元素 A(lm, ll) =-Vm(n)*Vm(l)*Ym(n,l)*cos(t(n,l)- delta(n)+delta(l); endif kb(n) = 0 & kb(l) = 0 % J4的非对角元素 A(lm, lk) =-Vm(n)*Ym(n,l)*sin(t(n,l)- d

28、elta(n) + delta(l);end else end else , end else, end end Pk = Vm(n)”2*Ym(n,n)*cos(t(n,n)+J33; Qk = -Vm(n)”2*Ym(n,n)*sin(t(n,n)-J11;if kb(n) = 1 P(n)=Pk; Q(n) = Qk; end % Swing bus P if kb(n) = 2 Q(n)=Qk;if Qmax(n) = 0Qgc = Q(n)*basemva + Qd(n) - Qsh(n); if iter 2 if Qgc Qmax(n), Vm(n) = Vm(n) - 0.01

29、;end else, end else,end else,end end if kb(n) = 1 A(nn,nn) = J11; % J1对角元素 DC(nn) = P(n)-Pk; end if kb(n) = 0 A(nn,lm) = 2*Vm(n)*Ym(n,n)*cos(t(n,n)+J22; % J2对角元素 A(lm,nn)= J33;% J3对角元素A(lm,lm) =-2*Vm(n)*Ym(n,n)*sin(t(n,n)-J44; % J4对角元素 DC(lm) = Q(n)-Qk; end end DX=ADC; for n=1:nbus nn=n-nss(n); lm=n

30、bus+n-ngs(n)-nss(n)-ns;if kb(n) = 1delta(n) = delta(n)+DX(nn); endif kb(n) = 0Vm(n)=Vm(n)+DX(lm); endendmaxerror=max(abs(DC);if iter = maxiter & maxerror accuracyfprintf(nWARNING: Iterative solution did not converged after ) fprintf(%g, iter), fprintf( iterations.nn)fprintf(Press Enter to terminate

31、the iterations and print the results n) converge = 0; pause, else, endendif converge = 1tech= (ITERATIVE SOLUTION DID NOT CONVERGE); elsetech=(Power Flow Solution by Newton-Raphson Method);endV = Vm.*cos(delta)+j*Vm.*sin(delta); deltad=180/pi*delta; i=sqrt(-1);k=0;for n = 1:nbusif kb(n) = 1k=k+1;S(n

32、)= P(n)+j*Q(n);Pg(n) = P(n)*basemva + Pd(n);Qg(n) = Q(n)*basemva + Qd(n) - Qsh(n);Pgg(k)=Pg(n);Qgg(k)=Qg(n);elseif kb(n) =2k=k+1;S(n)=P(n)+j*Q(n);Qg(n) = Q(n)*basemva + Qd(n) - Qsh(n);Pgg(k)=Pg(n);Qgg(k)=Qg(n);endyload(n) = (Pd(n)- j*Qd(n)+j*Qsh(n)/(basemva*Vm(n)”2); end busdata(:,3)=Vm; busdata(:,4

33、)=deltad;Pgt = sum(Pg); Qgt = sum(Qg); Pdt = sum(Pd); Qdt = sum(Qd); Qsht = sum(Qsh);6.decouple.m程序清单% Fast Decoupled Power Flow Solutionns=0; Vm=0; delta=0; yload=0; deltad=0;nbus = length(busdata(:,1);kb=;Vm=; delta=; Pd=; Qd=; Pg=; Qg=; Qmin=; Qmax=;Pk=; P=; Qk=; Q=; S=; V=;for k=1:nbusn=busdata(

34、k,1);kb(n)=busdata(k,2); Vm(n)=busdata(k,3); delta(n)=busdata(k, 4);Pd(n)=busdata(k,5); Qd(n)=busdata(k,6); Pg(n)=busdata(k,7); Qg(n) = busdata(k,8);Qmin(n)=busdata(k, 9); Qmax(n)=busdata(k, 10);Qsh(n)=busdata(k, 11);if Vm(n) = accuracy & iter =maxiter %检验不平衡功率 iter = iter+1; id=0; iv=0;for n=1:nbus

35、 nn=n-nss(n); J11=0; J33=0;for ii=1:nbrif mline(ii)=1if nl(ii) = n | nr(ii) = nif nl(ii) = n,l = nr(ii); endif nr(ii) = n,l = nl(ii); endJ11=J11+ Vm(n)*Vm(l)*Ym(n,l)*sin(t(n,l)- delta(n) + delta(l);J33=J33+ Vm(n)*Vm(l)*Ym(n,l)*cos(t(n,l)- delta(n) + delta(l); else , end else, end endPk = Vm(n)”2*Ym(

36、n,n)*cos(t(n,n)+J33;Qk = -Vm(n)”2*Ym(n,n)*sin(t(n,n)-J11;if kb(n) = 1 P(n)=Pk; Q(n) = Qk; end % Swing bus P if kb(n) = 2 Q(n)=Qk;% Between the 1th & 6th iterations% the Mvar of generator buses are% tested. If not within limits Vm(n)% is changed in steps of 0.05 pu to% bring the generator Mvar within

37、Qgc = Q(n)*basemva + Qd(n) - Qsh(n); if Qmax(n) = 0 if iter = 10 if Qgc Qmax(n),Vm(n) = Vm(n) - 0.005;end % the specified limits. else, end else, end else,end end if kb(n) = 1 id = id+1;DP(id) = P(n)-Pk;DPV(id) = (P(n)-Pk)/Vm(n); end if kb(n) = 0 iv=iv+1;DQ(iv) = Q(n)-Qk;DQV(iv) = (Q(n)-Qk)/Vm(n);en

38、d endDd=-B1DPV;DV=-B2DQV;id=0;iv=0;for n=1:nbus if kb(n) = 1 id = id+1;delta(n) = delta(n)+Dd(id); end if kb(n) = 0 iv = iv+1;Vm(n)=Vm(n)+DV(iv); endendmaxerror=max(max(abs(DP),max(abs(DQ);if iter = maxiter & maxerror accuracyfprintf(nWARNING: Iterative solution did not converged after ) fprintf(%g,

39、 iter), fprintf( iterations.nn)fprintf(Press Enter to terminate the iterations and print the results n) converge = 0; pause, else, endendITERATIVE SOLUTION DID NOT CONVERGE); elseif converge = 1 tech=(tech=(Power Flow Solution by Fast Decoupled Method);endk=0;V = Vm.*cos(delta)+j*Vm.*sin(delta);delt

40、ad=180/pi*delta;clear A; clear DC; clear DX i=sqrt(-1);for n = 1:nbusif kb(n) = 1S(n)=P(n)+j*Q(n);Pg(n) = P(n)*basemva + Pd(n);Qg(n) = Q(n)*basemva + Qd(n) - Qsh(n);k=k+1;Pgg(k)=Pg(n);elseif kb(n) =2S(n)=P(n)+j*Q(n);Qg(n) = Q(n)*basemva + Qd(n) - Qsh(n);k=k+1;Pgg(k)=Pg(n);endyload(n) = (Pd(n)- j*Qd(

41、n)+j*Qsh(n)/(basemva*Vm(n)”2);endbusdata(:,3)=Vm; busdata(:,4)=deltad;Pgt = sum(Pg); Qgt = sum(Qg); Pdt = sum(Pd); Qdt = sum(Qd); Qsht = sum(Qsh); clear Pk Qk DP DQ J11 J33 B1 B1inv B2 B2inv DPV DQV Dd delta ib id ii iv jb jj四、30节点电力系统计算实例潮流计算时，必须将前面的六个程序保存在MATLAB目录下格式为.m的文件，然后在MATLAB的命令窗口输入如下命令：cle

42、arbasemva = 100; accuracy = 0.001; accel = 1.8; maxiter = 100;%30节点电力系统% 母线一母线电压相角 负载 发电机注入功率%编号节点幅值角度有功无功有功无功无功最小值无功最大值无功busdata=111.060.00.00.00.00.0000221.0430.021.7012.740.00.0-40500301.00.02.41.20.00.0000401.060.07.61.60.00.0000521.010.094.219.00.00.0-40400601.00.00.00.00.00.0000701.00.022.810.

43、90.00.0000821.010.030.030.00.00.0-30400901.00.00.00.00.00.00001001.00.05.82.00.00.0-624191121.0820.00.00.00.00.00001201.0011.27.5000001321.071000.000-6240140106.21.600000150108.22.500000160103.51.800000170109.05.800000180103.20.900000190109.53.400000200102.20.7000002101017.511.2000002201000.00000023

44、0103.21.600000240108.76.700004.32501000.000000260103.52.3000002701000.0000002801000.000000290102.40.900000300%线路数据%bus buslinedata=1 2344252661R0.01920.04520.05700.01320.04720.05810.01190X0.05750.18520.17370.03790.19830.17630.041410.61.91/2 B0.026400.020400.018400.004200.020900.018700.0045001 for111

45、111100lines00570.04600.11600.010201670.02670.08200.008501680.01200.04200.004501690.00.20800.00.9786100.556000.9699110.2080019100.1100014120.256000.93212130.1400011214.1231.2559011215.0662.1304011216.0945.1987011415.2210.1997011617.0824.1923011518.1073.2185011819.0639.1292011920.0340.0680011020.0936.

46、2090011017.0324.0845011021.0348.0749011022.0727.1499012122.0116.0236011523.1000.2020012224.1150.1790012324.1320.2700012425.1885.3292012526.2544.3800012527.1093.20870128270.396000.9682729.2198.4153012730.3202.6027012930.2399.453301828.0636.20000.02141628.0169.05990.0651最后运行程序输入以下命令：Ifybus%形成节点导纳矩阵Ifg

47、auss%高斯-赛德尔法潮流计算busout%屏幕显示潮流计算结果lineflow%计算并显示线路潮流和损耗将Ifgauss变为Ifnewton/decouple,即可使用牛顿-拉夫逊法/快速解耦法进行潮流计算，输入以上4个命令行后，即可得到潮流计算结果：Command WindowC:.UmGLMg 牵 ipijiMkDQCLuiEt% .电力幕统分析程序的 trt-m is a case-insens it it? natchinstaalYou can inprave the performance of your code by using cKacl nane nal ches an

48、d we th&refare recamend that you update yaui usas? accoidingly, Alt rrnat ivlyj you can disablt this warning uing warning (J offnMAILAB: dispatcher; InesaclCsseMatchJ ) Ibis farriinE will btcoiic an. error in fuiu.ce rLeases.Pov&r FIdv Salution by GsusESeidel Methad MaxLnuR Power Misnatch = 0. 00095

49、 6B4No. af lie rat ions - 31and will be used1TJL、Na.Mag,DeerreMMvatrKWBvar-UIU匚M-arI1 q X 3 d 5 目 7 9 9LQ LI L2 L3 A -d,060】 0d3,022-013.010-012 005-010 051.0di 032.057 071.0d30-000 -5a 196 -3a 002 -9-明 9 -14. 33 D -1L. 396 -13.149 -12 111 -14. 432 -16. D2J -14. 432 -15-301 -15. 30 D -16.1900, 0002L

50、 7002, 400 L如。94b 200 cl oao22 SOO so. oao o oao s sao o oao1】.200 o oao 土 2000.000 L2.7D0 lr200 K600 L9.000CJ.000 L0.900 3(1.000 0.000 2.ODO 0.000 7.5D0 0.ODO l.fiOO260.950d 0.0000.0000.000D.OQQ0.000D.OQQD.000D.OQQD.OOOD.OQQD.OOOD.OOOD.OOO-17.010册.网6 。 Q00 A ooo35. 995 oa OOO o. ooa30. 7590. DOO0.

51、DOO16. U30. DOO10. 4Q60. DOO0,000 cl oao 0, 000 cl oao 0. 000 cl oao 0. 000 CL 0-30 0 000 19 0d0 0 000 0. 000 0. 000 o. oaoJ 1 4 LCommand V/indoutf151617IS19202122232425627?829301. D331. D451. 039i. raa 1.D251. D291. D321. D331. D271.D22L DL91. 00L场I- 0L1L D050. 9四-16.276-15.379-16.1B7-熟 B-17.(149-1

52、6.3E-16.JGS-5.455-16.60-16.329-16-J23-16.335T5E13-12.0E6-1& 133-13.016E.200 3.50-3 5.000 3.200 9.50Q 2.200 17.500 d. ooa 3.200 B. 700 o. ooa 5.500 o. ooa o.ooa 2. ioa 10.600500 L. BOO EL BOO 0. 000iOO 0. 700IL. 200 0. ooo L. 600 6. 700 oa ooo 2. 300 oa ooo 0. ooo 0. 900 L9000. 000Q. 0000. 00 D0 00 D

53、CL 00 D0 00 D0. 00 00 00 D0. 00 05 OQDCL 0000, 000CL 0000, 000CL 0000, 000o. oao0. Odd0. 0000. 0000. 0000. 0000. 0000. 0000. 0000. MOos oaoo, aooos oao0.000os 0-300.000a.odo a.ooo 0.0D0 a. ooo a.odo a. ooo a.ooo a.ooo a.ooo 4.3D0 Ou 000 0.000 Ou 000 0.000 du 000 0.000lot: al293.400126. 200300, 95012

54、5.03923.300Varn.ins : Could ncd f ind an. tKact (cas_3en3it vt) jiatcb far * linef lav-C:Us0rsLisoweipingDocunaits电力系,骨析KAILAfl程用LLLwflaw. n is a case-insensitive nalch and will he used inst cad.Yuu Dan iiiprova the parfornaiice uf your code by using ezacrtJ g眼聂成t皿场3尸j匚usage sccDiditiely- Alternativ

55、ely, ycu can disable This vamms usihevarning(?Jdi spot cher: InexactCaeMat ch3 )-IhiE varninff will bacone an error in future releases.A.r* 1Cammand WindowLineFlow andLoosesLine-Power nbus & line fLowL Ltie1O Ir-diisfotneifron toMVMvarMWNvartapI260.950-17.010261. 50n bL177.743-22. L4QL7. 1175.油L0.51

56、7383.975 L2533. 35d2. a 077. 0792L3. 30036. L264CL 1971 -172. 2B232 667L75. 3505.46L0.517445.7022- 720d5. 7811. 106-0. 59632.99。 704as. ooaW.9963. 173d61.905-0.河at. m2.0172.对33-2.iOO-2002. 6S1r1-30.390,95490. 442-3077. 073d78. 034-3u oa?7B. 095Ou 771.3J51-7. 600-6007. 7672-X 233M.7131. 106-0.5193-77

57、. 2634 43277.弱d0.77】3我670.132-L7.6247E. 3130,605.3lL2蚓.3L14*27如.192o.aoo4. &我0. 328-9. 2 DO16 9959 瓦 7212-79. 955& 1713D. 2552.9958. 178Ai-N. 21010. 4S717. B49Q. 15-.637T* 1irrCommand V/inriowsQ.QOOa. ooo0. DdO2-59.3533. 229E9 9452. D472.2634-69.627ta.O572. 0260. 6051.181-37.537-1. 91537. 5360. 368-

58、0.593329.634-3.71229. 766Q. 103-D.653g27.6B7-7.31828. G3B0. DQOt-5930.978LG15.32S0. 6561 瓦 342-0. QOQ1.279a.飒2338.S40-g. 5752. L310. OSO-13.085T-22.a 00-LQ. 90025. 27254.36-L2.5418. 3140. 151-1.687a-37. 1701. 31737. L930. 3dS-0-593s-30.000a. 75930. 010-29.4313. 5429. 5990. 103-D.55329-0.570-2. 3都2.

59、4330a 00090.000(L 000o. oao-27.6073. 91L29. 0360. DOO1.S93Ll0,003-L5.65315653A 000MlL027-7316.74728, 640-0. oooOrSLl10-r80017,t)00L：.962旦-15.32(L 02315. SdD-0a 0001.27g9-27.731-5.936P0, 369-o, m0.3L12a柬GE3. 569g, G%oa D81Du isaAL75-347土 3936, 9200.0140.037M* ICommand Window2115.723g.18-551Oa 1300u23

60、fiJt227.5B24. 437乩SU0. 052D. 107ILQ.QDOL6 H316. LIS9-Ou 003L6, 1-1LilCL 0000.M112-k200-7.5001% 4794-44P13-9. 94 L45, 237Q. DOO4.686L3-0u02-LQ. 27-110. 271CL 000Ou 132LI7.3522. 42SS, 219a. 0740, 155137.a52土獭成：斜CL 27Dud 28167.2063,3707, 9550. 0530.112130-00010. 40610. 40612o.aziLQ. -106io. las0. DOO旗1