计算机辅助分析课程设计及matlab自动实现潮流调节与计算

重要声明：1、本文只是提供一定参考，有错勿怪。2.、本文采用的潮流调节及计算只针对本实验，其中有很多固化的数据与参数，以及特殊的公式、请根据要求的不同进行一定的调整（潮流计算的程序适用任何情况、不同的只是系统参数）。3、本文潮流调节为自动调节生成，调节变量有发电机2、3电压，三台变压器分接头，没有进行无功补偿，请根据具体情况加入或删减变量。4、自动调节潮流程序中，控制变量为各点电压首次迭代值、最终各点电压。最终选择标准为网损最低的一组数据。5、本实验报告为草版，有错误，但懒得改。6、鉴于个人所有文档免费得原则，本文为免费文档，请尽量下载。7、最后！都免费了，要求不要太高。要求及原始数据电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷.各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。在电力系统中可能发生的各种故障中，危害最大且发生概率较高的首推短路故障。产生短路故障的主要原因是电力设备绝缘损坏。短路故障分为三相短路、两相短路、单相接地短路及两相接地短路。其中三相短路时三相电流仍然对称，其余三类短路统成为不对称短路。短路故障大多数发生在架空输电线路。电力系统设计与运行时,要采取适当的措施降低短路故障的发生概率。短路计算可以为设备的选择提供原始数据。一、课程设计的目的与任务：该课程设计为电气工程及其自动化专业的主要专业课设计之一,是学习电力系统分析基础和电力系统计算机辅助分析课程后的一个重要的实践性教学环节。其目的在于巩固和加深对电力系统潮流和短路电流计算基本原理的理解，学习和掌握应用计算机进行电力系统设计和计算的方法，培养学生独立分析和解决问题的能力。二、课程设计的基本要求：掌握电力系统等值模型和参数计算，以及潮流和短路计算的基本原理，学会应用计算机计算系统潮流分布和短路电流的方法。三、课程设计选题原则：该课程设计是根据电力系统分析基础和电力系统计算机辅助分析课程内容，结合实际工程和科研的电力系统网络进行系统的潮流和短路电流计算。四、课程设计内容：1、课程设计地点：机房和教室2,设计内容（1）电力系统潮流计算（2）电力系统短路电流计算具体内容包括：1）系统等值模型和参数计算。2）利用现有的潮流和短路电流计算程序进行计算，或编制电力系统潮流和短路计算程序。3）计算结果分析。4）编写说明书。3、时间安排1、系统等值模型和参数计算（2天）2、读懂现有的算法程序或视需要编制算法程序并进行计算（4天）3、计算结果整理和分析（1天）4、编写说明书（2天）5、成绩考核（1天）五、课程设计主要参考资料：1、李光琦.电力系统暂态分析.北京：水利电力出版社.第二版2、陈拓.电力系统稳态分析.北京：水利电力出版社.第二版3、孟祥萍.基于MATLAB的电力系统设计与计算.长春：吉林人民出版社4、孟祥萍.电力系统分析.北京：高等教育出版社5、电力系统潮流和短路计算课程设计指导书.自编六、课程设计成绩考核：根据学生的平时表现，说明书，以及设计完成后的口试成绩等进行综合评定。课程结束时，要求学生有课程设计说明书。课程设计成绩以平时成绩（占50%）、设计说明书（占50%）进行综合评定。

原始数据与思考题目标电网单线图系统参数表1.线路参数表缴编号(km)线犍阻{Q/km}型ITT删九{fi/km}维循此半{SAm}4-5LGJ-240/301040.0470.41.78X10-64-6LGJ-120/701100.0741.47X10-65-7LGJ-120/25980.0791.60X10-66-9LGJ-95/551630.0921.80XI。/7-8LGJ-240/30880.0471.78XI。/8-9LGJ-240/301220.0471.78XI。/说明：线路零序电抗为正序由抗3倍。表2.变压器参数表

线辍号变压耀号变压敏:匕（kV）短路电压百分数（％）2-7SSPL-220000242±3X2.5%/20103-9SSPL-120000242±3X2.5%/1561-4SSPL-240000242±3X2.5%/17.511说明：变压器零序电抗与正序电抗相等，且均为A/Y。接法。表3.发电机参数表1发电机额定功率{MW}额定电压{kV}额定功率因数120016.50.852180180.85310013.80.85表4.发电机参数表2发电机母线名X(/｛。｝x“｛。｝Go(S)X“｛。｝Tq。（S）小S)110.14600.06088.960.09690.0969047.28220.89580.11986.000.86450.19690.53512.80331.31250.18138.591.25780.25000.6006.02表5.负荷数据表节点号有功负荷（MW）无功负荷（MVA）51255069030810035思考题1、牛顿拉夫逊法与PQ分解法有哪些联系？有哪些区别？二者的计算性能如何？2、为什么在用计算机对某网络初次进行潮流计算时往往是要调潮流，而并非任何情况下只一次送入初始值算出结果就行呢？要考虑什么条件？各变量是如何划分的？哪些可调？哪些不可调？3、发电机节点的注入无功为负值说明了什么？4、线路有功潮流最有可能的流向是？线路无功潮流最有可能的流向是？5、负荷功率因数对系统潮流有什么影响？6、潮流控制的主要手段有哪些？7、如何降低系统网损？8、如何校验主设备的短路容量？9、利用节点阻抗矩阵计算短路电流和利用导纳矩阵计算短路电流的算法有何异同？10、减小短路电流的措施有哪些?设计主要内容：1、完成调试潮流计算与短路计算程序调试；计算程序以《电力系统计算机辅助分析》所用教材中的例程为基础，要求掌握算法原理，建议进行适当改写以利于本次设计应用。2、根据设计资料准备计算数据；整理原始数据，使其符合计算程序的输入条件。要求比较采用变压器n型模型前后两种情况的差异。3、为指定电网安排合理的稳态运行方式；合理设置目标电网中的控制变量，以1号机为平衡机，将潮流分布调整为可行的运行方式（要求所有节点的电压不得低于L0p.U.,也不能高于L05P.U.）绘制潮流分布图。分析该运行方式的特点，提出进一步降低网损的改进的措施，并利用计算结果简单证明之。具体内容参考设计资料所附思考题。4、列表分析节点7发生各种短路故障时的短路电流，以及运行方式对故障电流的影响；5、完成设计说明书。要求结构清晰，即从电网模型建立，到初态潮流分布计算，到问题分析，再到调整后的方案计算。故障电流计算也类似。不必逐一回答所有思考题。参数计算、潮流调节、与潮流计算matlab源程序%参数计算的clc;clear;L=[104,110,98,163,88,122];r=[0.047,0.074,0.079,0.092,0.047,0.047];b=[1,78e-006,1,47e-006,1,60e-006,1,80e-006,1,78e-006,1,78e-006];disp(线路电抗X=');X1=100/220竺*0.4*L;%线路各段电抗(均为标幺值,取Ub=220KV,SB=1OOMVA,下同)disp(Xl);fori=l:6R(i)=100/220A2*L(i)*r(i);%加分号避免出现过程项，下同end %线路各段电阻disp(线路电阻R=); disp(R);%输出Rfori=l:6Z(i)=R(i)+j*Xl(i);Z0(i)=R(i)+3*j*Xl(i);enddisp('线路正序、负序阻抗Z1=Z2=);disp(Z);fori=l:6Bl(i)=220A2/100*b(i)*L(i);B(i)=O+j*Bl(i);end %线路各容纳之半disp(线路容纳之半B/2=');disp(B); %输出Bdisp(线路零序阻抗Z0:');disp(ZO);U0=[l1,10,6];SN=[240,220,120];UN=f242,242,242];fori=l:3;XT(i)=j*U0(i)*242A2/SN(i)/220A2;XT0(i)=XT(i);end;disp,变压器正、负序电抗值Xt：');disp(XT);dispC变压器零序电抗值XtO：');disp(XT0);x=[0.0608,0.1198,0.1813]; %发电机暂态电抗有名值UGN=[16.5,18,13.8];fori=l:3K(i)=242/UGN(i);endfori=l:3Xdg(i)=j*100/220A2*x(i)*K(i)*K(i);end %发电机暂态电抗dispC发电机暂态电抗Xdg=');disp(Xdg);%输出XdgPL=[125,90,100];QL=[5O,3O,35];fori=l:3ZL(i)=100/(PL(i)A2+QL(i)A2)*(PL(i)+j*QL(i));end %负荷节点阻抗dispC负荷节点阻抗ZIR; disp(ZL);%输出ZL%支数参数形成的矩阵Bl、支数参数形成的矩阵B2、节点号及其对地阻抗形成的矩阵X的形成B1=r1,4,XT(1),0,1,0;2,7,XT(2),0,1,0;3,9,XT(3),0,1,0;4,5,R(1)+X1(1)*j,B(1),1,0;4,6,R⑵+Xl(2)*j,B(2),l,0;5,7,R(3)+Xl(3)*j,B(3),l,0;6,9,R(4)+Xl(4)*j,B(4),l,0;7,8,R(5)+Xl(5)*j,B(5),l,0;8,9,R(6)+Xl(6)*j,B(6),l,0;];dispC支数参数形成的矩阵：Bl=');disp(Bl);%输出支数参数形成的矩阵B1P=[200,180,100];W=(P/lOO)+((P/O.85)*sqrt(l-O.85*O.85)/lOO)*j;B2=[W(1),0,105,1,0,1； W(2),0,l,l,0,3; W(3),0,l,1,0,3；0,0,l,0,0,2;0,1.25+0.5i,l,0,0,2;0,0.9+0.3i,1,0,0,2;0,0,1,0,0,2;0,1+0.35i,1,0,0,2;0,0,1,0,0,2;];dispC支数参数形成的矩阵：B2=');disp(B2);%输出支数参数形成的矩阵B2%潮流自动调节的与自动计算的n=9;nl=9;isb=l;pr=0.0001;mm=0;X=[1,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0;9,0];dispC*******************************************************************、♦disp(，只考虑各点电压水平，得到的满足要求的数据如下：)；foraa=0:5B2(2,4)=l+0.01*aa;forab=0:5B2(3,4)=l+0.01*ab;forba=0:5Bl(l,5)=0.95+0.025*ba;forbb=0:5Bl(2,5)=0.95+0.025*bb;forbc=0:5Bl(3,5)=0.95+0.025*bc;Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=zeros(1,n);S1=zeros(nl);fori=l:nifX(i,2)~=0;p=X(i,l);Y(p,p)=l./X(i,2);endendfori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2);q=Bl(i,l);endY(p,q)=Y(p,q)-l./(B1(i,3)*Bl(i,5));Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+l./(B1(i,3)*Bl(i,5)A2)+Bl(i,4),/2;Y(p,p)=Y(p,p)+1./B1(i,3)+Bl(i,4)./2;endG=real(Y);B=imag(Y);fori=l:ne(i)=real(B2(i,3))；f(i)=imag(B2(i,3))；V(i)=B2(i,4);endfori=l:nS(i)=B2(i,l)-B2(i,2);B(i,i)=B(i,i)+B2(i,5);endP=real(S);Q=imag(S);ICTl=0;IT2=l;N0=2*n;N=N0+l;a=0;whileIT2~=0IT2=0;a=a+l;fori=l:nifi~=isbC(i)=o；D(i)=O;forjl=l:nC(i)=C(i)+G(i,jl)*eGl)-B(i,jl)*f(jl);D(i)=D(i)+G(i,jl)*f(jl)+B(ijl)*e(j1);endPl=C(i)*e(i)+f(i)*D(i);Ql=f(i)*C(i)-D(i)*e(i);V2=e(i)A2+f(i)A2;ifB2(i,6)~=3DP=P(i)-Pl;DQ=Q(i)-Ql;forjl=l:nifjl~=isb&jl~=iXl=-G(i,jl)*e(i)-B(ijl)*f(i);X2=B(ijl)*e(i)-G(ijl)*f(i);X3=X2;X4=-X1;p=2*i-l;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X4;J(m,q)=X2;elseifj1==i&j1~=isbXl=-C(i)-G(i,i)*e(i)-B(i,i)*f(i)；X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i);X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);p=2*i-l;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2;endendelseDP=P(i)-Pl;DV=V(i)A2-V2;forjl=l:nifjl~=isb&jl~=iXl=-G(ijl)*e(i)-B(ijl)*f(i);X2=B(ijl)*e(i)-G(ijl)*f(i);X5=0;X6=0;p=2*i-l;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+l;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2;elseifjl==i&jl~=isbXl=-C(i)-G(i,i)*e(i)-B(i,i)*f(i)；X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);X5=-2*e(i);X6=-2*f(i);p=2*i-l;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+1;J(m,q)=Xl;J(m,N)=DP;q=q+I;J(p,q)=X6;J(m,q)=X2;endendendendendfork=3:N0kl=k+l;Nl=N;fork2=kl:NlJ(k,k2)=J(k,k2)./J(k,k);endJ(k,k)=l;ifk~=3k4=k-l;fork3=3:k4fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endifk==NO,break;endfork3=kl:N0fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endelsefork3=kl:N0fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endendendfork=3:2:N0-lL=(k+l)./2;e(L)=e(L)-J(k,N);kl=k+l;f(L)=f(L)-J(kl,N);endfork=3:N0DET=abs(J(k,N));ifDET>=prIT2=IT2+1;endendICT2(a)=IT2;ICT1=ICT1+1;fork=l:ndy(k)=sqrt(e(k)A2+f(k)A2);endfori=l:nDy(ICTl,i)=dy(i);endforvb=l:nifDy(l,vb)<=1.05&Dy(l,vb)>=1.00;zy(vb)=l;elsezy(vb)=0;endendend%用高斯消去法解"w=・J*V”fork=l:nV(k)=sqrt(e(k)A2+f(k)A2);O(k)=atan(f(k)./e(k))*180./pi;endE=e+f*j;forva=l:nifV(va)<=1.05&V(va)>=1.000;yj(va)=l;elseyj(va)=0;endendifyj(D*yj(2)*yj(3)*yj(4)*yj(5)*yj(6)*yj(7)*yj(8)*yj(9)==lifzy(1)*zy(2)*zy(3)*zy(4)*zy(5)*zy(6)*zy(7)*zy(8)*zy(9)==l;dispC发电机2出口电压：)；disp(B2(2,4))；dispC发电机3出口电压：上disp(B2(3,4))；dispC变压器1接头电压：,);disp(Bl(l,5));dispC变压器2接头电压：上disp(B1(2,5));dispC变压器3接头电压：)disp(B1(3,5));dispC各节点的电压大小V为(节点号从小到大排列)：)；disp(V);forp=l:nC(p)=0;forq=l:nC(p)=C(p)+conj(Y(p,q))*conj(E(q));endS(p)=E(p)*C(p);endfori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2)；q=Bl(i,l);endSi(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))-conj(E(q)))*conj(1./(Bl(i,3)*Bl(i,5))));endfori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2);q=Bl(i,l);endSj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))-conj(E(p)))*conj(1./(Bl(i,3)*Bl(i,5))));endZDS=0;fori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2)；q=Bl(i,l);endDS(i)=Si(p,q)+Sj(q,p);ZDS=DS(i)+ZDS;endmm=mm+l;ZDSS(mm)=sqrt(real(ZDS)A2+imag(ZDS)A2);dispC******************************************************************，＞fori=l:nZB224(mm)=B2(2,4)；ZB234(mm)=B2(3,4)；ZB115(mm)=Bl(l,5);ZB125(mm)=B1(2,5);ZB135(mm)=B1(3,5);endendendendendendendenddisp('******************************************************************)；'dispC满足各点电压要求的数据组共有：)；disp(mm);LOL=min(ZDSS);fordota=l:mmifZDSS(dota)==LOLmm=dota;endenddispC当加入线路损耗时，得到其中线路损耗最小的一组组号及关键数据为：);disp(mm);disp('发电机2出口电压：);disp(ZB224(mm));dispC发电机3出口电压：');disp(ZB234(mm));dispC变压器1接头电压：1);disp(ZB115(mm));dispf变压器2接头电压：1);disp(ZB125(mm));dispC变压器3接头电压：，)；disp(ZBl15(mm));B2(2,4)=ZB224(mm);B2(3,4)=ZB234(mm);Bl(l,5)=ZB115(mm);Bl(2,5)=ZB125(mm);Bl(3,5)=ZB135(mm);dispC用此组数据算的的潮流结果为：上disp('节点数：');disp(n);dispC节点支路数：);disp(nl);disp('平衡节点：上disp(isb);disp('精度：)；disp(pr);disp('Bl=');disp(Bl);disp('B2=');disp(B2);disp('X=');disp(X);Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=zeros(1,n);Sl=zeros(nl);dispC节点数：上disp(n);dispC节点支路数：,);disp(nl);disp('平衡节点：上disp(isb);disp('精度：上disp(pr);mm=O;X=[1,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0;9,0];fori=l:nifX(i,2)~=0;p=X(i,l);Y(p,p)=l./X(i,2);endendfori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2);q=Bl(i,l);endY(p,q)=Y(p,q)-l./(Bl(i,3)*Bl(i,5));Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+l./(Bl(i,3)*Bl(i,5)A2)+Bl(i,4)./2;Y(p,p)=Y(p,p)+l./Bl(i,3)+Bl(i,4)./2;endG=real(Y);B=imag(Y);fori=l:ne(i)=real(B2(i,3));f(i)=imag(B2(i,3));V(i)=B2(i,4)；endfori=l:nS(i)=B2(i,l)-B2(i,2);B(i,i)=B(i,i)+B2(i,5);endP=real(S);Q=imag(S);ICT1=O;IT2=1;N0=2*n;N=N0+l;a=0;whileIT2~=0IT2=0;a=a+l;fori=l:nifi~=isbC(i)=0;D(i)=0;forjl=l:nC(i)=C(i)+G(i,jl)*e(jl)-B(i,jl)*f(jl);D(i)=D(i)+G(i,jl)*f01)+B(i,jl)*e(jl);endPl=C(i)*e(i)+f(i)*D(i);Ql=f(i)*C(i)-D(i)*e(i);V2=e(i)A2+f(i)A2;ifB2(i,6)~=3DP=P(i)-Pl;DQ=Q(i)-Ql;forjl=l:nifjl~=isb&jl~=iXl=-G(ijl)*e(i)-B(ijl)*f(i)；X2=B(ijl)*e(i)-G(ijl)*f(i);X3=X2;X4=-X1;p=2*i-1;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2;elseifj1==i&j1~=isbXl=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i);X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);p=2*i-l;q=2*j1-1;J(p,q)=X3;J(p,N)=DQ;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X4;J(m,q)=X2;endendelseDP=P(i)-Pl;DV=V(i)A2-V2;forjl=l:nifjl~=isb&jl~=iXl=-G(ijl)*e(i)-B(ijl)*f(i);X2=B(ijl)*e(i)-G(ijl)*f(i);X5=0;X6=0;p=2*i-l;q=2*jl-l;J(p,q)=X5;J(p,N)=DV;m=p+1;J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X6;J(m,q)=X2;elseifj1==i&j1~=isbXl=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);X5=-2*e(i);X6=-2*f(i);p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+1;J(m,q)=X1;J(m,N)=DP;q=q+1;J(p,q)=X6;J(m,q)=X2;endendendendendfork=3:N0kl=k+l;Nl=N;fork2=kl:NlJ(k,k2)=J(k,k2),/J(k,k);endJ(k,k)=l;ifk~=3k4=k-l;fork3=3:k4fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endifk==NO,break;endfork3=kl:N0fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endelsefork3=kl:N0fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endendendfork=3:2:N0-lL=(k+1).⑵e(L)=e(L)-J(k,N);kl=k+l;f(L)=f(L)-J(kl,N);endfork=3:N0DET=abs(J(k,N));ifDET>=prIT2=IT2+1;endendICT2(a)=IT2;ICT1=ICT1+1;fork=l:ndy(k)=sqrt(e(k)A2+f(k)A2);endfori=l:nDy(ICTl,i)=dy(i);endend%用高斯消去法解"w=-J*V"disp('迭代次数上disp(ICTl);dispC没有达到精度要求的个数工disp(ICT2);fork=l:nV(k)=sqrt(e(k)A2+f(k)A2);O(k)=atan(f(k)./e(k))*18O./pi;endE=e+f*j;dispC各节点的实际电压标幺值E为(节点号从小到大排列)：1disp(E);disp('各节点的电压大小V为(节点号从小到大排列)：,)；disp(V);dispC各节点的电压角。为(节点号从小到大排列)：)；disp(O);forp=l:nC(p)=0；forq=l:nC(p)=C(p)+conj(Y(p,q))*conj(E(q));endS(p)=E(p)*C(p);enddisp。各节点的功率S为(节点号从小到大排列)：,);disp(S);dispC各条支路的首段功率Si为(顺序同您输入Bl时一样)：)；fori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2);q=Bl(i,l);endSi(p,q)=E(p)*(conj(E(p))*conj(B1(i,4)./2)+(conj(E(p)*B1(i,5))-conj(E(q)))*conj(1./(Bl(i,3)*Bl(i,5))));disp(Si(p,q));enddispC各条支路的末段功率Sj为(顺序同您输入Bl时一样)：)；fori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2);q=Bl(i,l);endSj(q,p)=E(q)*(conj(E(q))*conj(B1(i,4)./2)+(conj(E(q)./B1(i,5))-conj(E(p)))*conj(1./(Bl(i,3)*Bl(i,5))));disp(Sj(q,p));enddispC各条支路的功率损耗DS为(顺序同您输入B1时一样)：,);fori=l:nlifBl(i,6)==0p=Bl(i,l);q=Bl(i,2);elsep=Bl(i,2);q=Bl(i,l);endDS(i)=Si(p,q)+Sj(q,p);disp(DS(i));endfori=l:ICTlCs(i)=i;enddispC以下是各节点每次迭代后的电压值(如图所示),);plot(Cs,Dy),xlabel。迭代次数)ylabelC电压)title。电压迭代次数曲线,)；最后是课程设计报告二刃理N火弯课程设计课程名称： 电力系统分析 设计题目：基于Matlab计算程序的电力系统运行分析学 院： 电力工程学院 专业： 电气工程自动化 年 级： 2009 学生姓名易军（你可以称我易大神，可惜dota、lol都是菜鸟顶级水平）指导教师： 宋琪、谢静、单节杉 日期： 2013T-13 教务处制电力系统潮流计算是电力系统分析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件是否过负荷.各点电压是否满足要求，功率的分布和分配是否合理以及功率损耗等。对现有用力系统的运行和扩建，对新的电力系统进行规划设计以及对用力系统进行静态和暂态稳定分析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。在电力系统中可能发生的各种故障中，危害最大且发生概率较高的首推短路故障。产生短路故障的主要原因是电力设备绝缘损坏。短路故障分为三相短路、两相短路、单相接地短路及两相接地短路。其中三相短路时三相电流仍然对称，其余三类短路统成为不对称短路。短路故障大多数发生在架空输电线路。电力系统设计与运行时，要采取适当的措施降低短路故障的发生概率。短路计算可以为设备的选择提供原始数据。H标电网接线图

一、电网模型的建立设计中，选取Sb=100MVA,Ub=220KV,将所有支路的参数都折算到220KV电压等级侧本文利用matlab计算参数其源程序为：clc;clear;L=[104,110,98,163,88,122];r=[0.047,0.074,0.079,0.092,0.047,0.047];b=[1.78e-006,1.47e-006,1.60e-006,1.80e-006,1.78e-006,1.78e-006];dispC线路电抗X=*);Xl=l00/220A2*0.4*L; %线路各段电抗(均为标幺值，disp(Xl);fori=l:6R(i)=100/220A2*L(i)*r(i);enddispC线路电阻R=');disp(R);fori=l:6Z(i)=R(i)+j*Xl(i);Z0(i)=R(i)+3*j*Xl(i);enddispC线路正序、负序阻抗Z1=Z2=);disp(Z);fori=l:6Bl(i)=220A2/100*b(i)*L(i);B(i)=0+j*Bl(i);enddispC线路容纳之半B/2=');disp(B);dispC线路零序阻抗Z0:»);disp(ZO);U0=[l1,10,6];SN=[240,220J20J;UN=[242,242,242];fori=l:3;XT(i)=j*U0(i)*242A2/SN(i)/220A2;XT0(i)=XT(i);end;dispC变压器正、负序电抗值Xt：);disp(XT);dispC变压器零序电抗值XtO：*);disp(XT0);取Ub=220KV,SB=100MVA)%加分号避免出现过程项，下同取Ub=220KV,SB=100MVA)%加分号避免出现过程项，下同%线路各段电阻%输出R%线路各容纳之半%输出B%发电机暂态电抗有名值UGN=[16.5,18,13.8];fori=l:3K(i)=242/UGN(i);endfori=l:3Xdg(i)=j*l00/220A2*x(i)*K(i)*K(i);enddispC发电机暂态电抗Xdg=')；disp(Xdg);PL=[125,90,100];QL=[50,30,35];fori=l:3ZL(i)=i00/(PL(i)A2+QL(i)A2)*(PL(i)+j*QL(i));enddispC负荷节点阻抗Zl=');disp(ZL);%发电机暂态电抗%负荷节点阻抗计算结果为：线路电抗X=0.0860 0.0909 0.0810线路电阻!^=0.0101 0.0168 0.0160线路正序、负序阻抗Z1=Z2=Columns1through50.13470.03100.0101+0.0860iColumn60.0118+0.1008i线路容纳之半B/2=Columns1through50+0.0896iColumn60+0.105H线路零序阻抗Z0:Columns1through50.0101+0.2579iColumn60.0118+0.3025i0.0168+0.0909i0+0.0783i0.0168+0.2727i变压器正、负序电抗值Xt：0+0.0555i 0+0.0550i变压器零序电抗值XtO：0+0.0555i 0+0.0550i发电机暂态电抗Xdg=0+0.0270i 0+0.0447i负荷节点阻抗Zl=0.6897+0.2759i 1.0000+O.3333i0.07270.00850.0160+0.08lOi0+0.0759i0.0160+0.2430i0+0.0605i0+0.0605i0+0.1152i0.8909+0.3118i0.10080.01180.0310+0.1347i0+0.1420i0.0310+0.404li0.0085+0.0727i0+0.0758i0.0085+0.2182i支数参数形成的矩阵：EH=Columns1through51.00004.0000 0+0.0555i 0 1.00002.00007.0000 0+0.0550i 0 1.00003.00009.0000 0+0.0605i 0 1.00004.00005.0000 0.0101+0.0860i 0+0.0896i 1.00004.00006.0000 0.0168+0.0909i 0+0.0783i 1.00005.00007.0000 0.0160+0.08lOi 0+0.0759i 1.00006.00009.0000 0.0310+0.1347i 0+0.1420i 1.00007.00008.0000 0.0085+0.0727i 0+0.0758i 1.00008.00009.0000 0.0118+0.1008i 0+0.1051i 1.0000Column6000000000支数参数形成的矩阵：B2=Columns1through52.0000+1.2395i0 1.0500 1.0000 01.8000+1.1155i0 1.0000 1.0000 01.0000+0.6197i0 1.0000 1.0000 000 1.0000 0 00 1.2500+0.5000i 1.0000 0 00 0.9000+0.3000i 1.0000 0 00 01.0000 0 00 1.0000+O.35OOi 1.0000 0 00 01.0000 0 0Column61.00003.00003.00002.00002.00002.00002.00002.00002.0000

系统等值电路图的绘制根据以上计算结果，得到系统等值电路图如下:0.8909+J0.31180.8909+J0.3118E2J0.0447jaosso0.0310+J0X)7270.0118+J0.1008J0.06051J04051J0.1051371工JO0758 J0Q75^1-1-1J0,1278E2J0.0447jaosso0.0310+J0X)7270.0118+J0.1008J0.06051J04051J0.1051371工JO0758 J0Q75^1-1-1J0,12780.0160*J0.0810.II-.J01708

0.03310+JO.134^J0.1152J0.0759TIJ0.0759TIJ0J420II-0.6897+J0.2759—iiJ0.1060a0101+J0.0860J0.08960.6897+J0.2759—iiJ0.1060a0101+J0.0860J0.0896HRII—JO07190.0168+J0.0909ja0783-HI—1.0000+J0^333J0.0555J0.0270ElJ0.0555J0.0270El二、潮流计算.系统参数的设置设计中要求所有结点电压不得低于1.Op.U.,也不得高于1.05p.u.,若电压不符合该条件，可采取下面的方法进行调压：改变发电机的机端电压改变变压器的变比（即改变分接头）改变发电机的出力（4）在电压不符合要求的结点处增加无功补偿调压方式应属于逆调压。结点的分类：根据电力系统中各结点性质的不同，将结点分为三类：PQ结点、PV结点和平衡结点，在潮流计算中，大部分结点属于PQ结点，小部分结点属于PV结点，一般只设一个平衡结点。对于平衡结点，给定其电压的幅值和相位，整个系统的功率平衡由这一点承担。本设计中，选1号节点为平衡节点；2、3号节点为P、U节点；4、5、6、7、8、9号结点为P、Q节点。采用PQ分解法进行潮流的计算和分析。设计中，节点数:n=9,支路数:nl=9,平衡母线节点号:isb=l,误差精度:pr=O.OOOOl。由支路参数形成的矩阵：矩阵B1的每行是由卜列参数构成的：1、某支路的首段号p；2、末端号q,且p<q；3、支路的阻抗（R+Xi）；4、支路的对地导纳；5、支路的变比；6、折算到哪一侧的标志（支路首段p处于高压侧则输入“1”，否则为"0”）；Bl=[l40.05546i010;270.055i010;390.0605i010;450.01+0.0859510.1792i10;460.0168+0.091i0.1566i10;570.016+0.081i0.1518i10;690.031+0.1347i0.284i10;780.0085+0.0727i0.1516i10;890.0119+0.1008i0.2102i10]矩阵B2的每行是由下列参数构成的：1、节点所接发电机的功率SC；2、节点负荷功率S1;3、节点电压初始值；4、PU节点电压U的给定值：5、节点所接的无功补偿设备容量；6、节点分类标号（1—平衡节点；2—PQ节点；3—PU节点）本文采用的是自动筛选法，利用matlab选出合理的数据组成，然后验证合理性。首先采用只调节发电机，结果没有可选值。加上变压器后通过限制初次迭代值和节点电压值最终筛选出了结果。源程序如下clc;n=9;nl=9;isb=l;pr=O.0001;mm=0;Bl=[1.000004.00000+0.0555i01.00002.00007.00000+0.0550i01.000003.00009.00000+0.0605i01.000004.00005.00000.0101+0.0860i04-0.0896i1.000004.00006.00000.0168+0.0909i04-0.0783i1.000005.00007.00000.0160+0.0810i0+0.0759i1.000006.00009.00000.0310+0.1347i04-0.1420i1.000007.00008.00000.0085+0.0727i04-0.0758i1.000008.00009.00000.0118+0.1008i04-0.1051i1.00000]B2=[2.0000+1.2395i01.05001.000001.0000;1.8000+1.1155i01.00001.000003.0000;1.0000+0.6197i01.00001.000003.0000：001.0000002.0000;01.2500+0.5000i1.0000002.0000;00.9000+0.3000i1.0000002.0000;001.0000002.0000;01.0000+0.3500i1.0000002.0000;001.0000000000;];X=[l,0;2,0;3,0;4,0;5,0;6,0;7,0;8,0;9,0];foraa=0:5B2(2,4)=l+0.01*aa;forab=0:5B2(3,4)=l+0.01*ab;forba=0:5Bl(l,5)=0.95+0.025*ba;forbb=O:5Bl(2,5)=0.95+0.025*bb;forbc=0:5Bl(3,5)=0.95+0.025*bc;Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);0二zeros(1,n);Sl=zeros(nl)生成初始矩阵Y为n*n的全0矩阵efV0均为Ixn的全0矩阵fori=l:nifX(i,2)~=0;P=X(i,1);Y(p,p)=l./X(i,2);endendfori=l:nlifBl(i,6)=0p=Bl(i,1);q=Bl(i,2);elsep=Bl(i,2);q=Bl(i,1);end(p,q)=Y(p,q)-l./(Bl(i,3)*Bl(i,5));(q,p)=Y(p,q);(q,q)=Y(q,q)+l./(Bl(i,3)*Bl(i,5厂2)+Bl(i,4)./2;(p,p)=Y(p,p)+l./Bl(i,3)+Bl(i,4)./2;endG=real(Y);B=imag(Y);fori=l:ne(i)=real(B2(i,3));f(i)=imag(B2(i,3));V(i)=B2(i,4)；%判断X(i,2)不等于0%判断是否有接地阻抗支路%不含变压器的支路%含有变压器的支路%计算节点导纳矩阵%电导%电纳%电压fori=l:nS(i)=B2(i,l)-B2(i,2);B(i,i)=B(i,i)+B2(i,5);endP=real(S);Q=imag(S);ICT1=O;IT2=1;N0=2*n;N=NO+1;a=0;whileIT2'=0IT2=0;a=a+l;fori=l:nifi"=isbCG)=0;D⑴=0;forjl=l:nC(i)=C(i)+G(i,jl)*e(jl)-B(i,jl)*f(jl);D(i)=D(i)+G(i,jl)*f(jl)+B(i,jl)*e(jl);endPl=C(i)*e(i)+f(i)*D(i);Ql=f(i)*C(i)-D(i)*e(i);V2二e(i厂2+f(i)C2;ifB2(i,6)DP二P⑴-Pl;DQ=Q(i)-Ql;forjl=l:nifjP=isb&jl"=iXl=-G(i,jl)*e(i)-B(i,jl)*f(i)；X2=B(i,jl)*e(i)-G(i,jl)*f(i);X3=X2;X4=-X1;p=2*i-1;q=2*jl-l;J(p,q)=X3;J(p,N)=DQ;m=p+1;J(m,q)=Xl;J(m,N)=DP;q=q+l；J(p,q)=X4;J(m,q)=X2;elseifjl==i&jl"=isbXl=-C(i)-G(i,i)*e(i)-B(i,i)*f⑴；X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);X3=D(i)+B(i,i)*e(i)-G(i,i)*f(i);X4=-C(i)+G(i,i)*e(i)+B(i,i)*f(i);P=2*i-1;q=2*jl-l;J(p,q)=X3;J(p,N)=DQ;m=p+l;J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X4;J(m,q)=X2;endendelseDP=P(i)-Pl;DV=V(i)^2-V2;forjl=l:nifjT=isb&jr=iXl=-G(i,jl)*e(i)-B(i,jl)*f(i);X2=B(i,jl)*e(i)-G(i,jl)*f(i);X5=0;X6=0;P=2*i-1;q=2*jl-l;J(p,q)=X5;J(p,N)=DV;m=p+l;J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X6;J(m,q)=X2;elseifjl==i&jl~=isbXl=-C(i)-G(i,i)*e(i)-B(i,i)*f(i);X2=-D(i)+B(i,i)*e(i)-G(i,i)*f(i);X5=-2*e(i);X6=-2*f(i);p=2*i-1;q=2*j1-1;J(p,q)=X5;J(p,N)=DV;m=p+1;J(m,q)=Xl;J(m,N)=DP;q=q+l;J(p,q)=X6;J(m,q)=X2;endendendendendfork=3:N0kl=k+l;Nl=N;fork2=kl:NlJ(k,k2)=J(k,k2)./J(k,k);endJ(k,k)=l;ifk、=3k4=k-l;fork3=3:k4fork2=kl:N1J(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endifk==N0,break;endfork3=kl:N0fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endelsefork3=kl:N0fork2=kl:NlJ(k3,k2)=J(k3,k2)-J(k3,k)*J(k,k2);endJ(k3,k)=0;endendendfork=3:2:N0-lL=(k+l)./2;e(L)=e(L)-J(k,N);kl=k+l;f(L)=f(L)-J(kl,N);endfork=3:N0DET=abs(J(k,N));ifDET>=prIT2=IT2+1;endendICT2(a)=IT2;ICT1=ICT1+1;fork=l:ndy(k)=sqrt(e(k)*2+f(k)2);endfori=l:nDy(ICTl,i)=dy(i);endforvb=l:nifDy(l,vb)<=1.05&Dy(l,vb)>=1.00;zy(vb)=l;elsezy(vb)=0;endendend%用高斯消去法解〃w=-J*V〃fork=l:nV(k)=sqrt(e(k)"2+f(k)2);0(k)=atan(f(k)./e(k))*180./pi;endE=e+f*j;forva=l:nifV(va)<=1.05&V(va)>=1.000;yj(va)=l;elseyj(va)=0;endendifyjd)*yj(2)*yj(3)*yj(4)*yj(5)*yj(6)*yj(7)*yj(8)*yj(9)==1ifzy(l)*zy(2)*zy(3)*zy(4)*zy(5)*zy(6)*zy(7)*zy(8)*zy(9)==1disp('发电机2出口电压：')；disp(B2(2,4));disp('发电机3出口电压：')；disp(B2(3,4));disp('变压器1接头电压：')；disp(Bl(l,5));disp('变压器2接头电压：')；disp(Bl(2,5));disp('变压器3接头电压：')；disp(81(3,5));dispf各节点的电压大小V为(节点号从小到大排列)：’)；disp(V);disp('#################################################################')；endendendendendendenddisp('满足要求的数据组共有：’)；disp(mm)：dispC#################################################################')；运行后共筛选出了41组数据，经过改变电压范围最终确定了下面两组数据。第一组数据：节点数： 9节点支路数：9平衡节点： 1精度：1.0000e-004Bl=Columns1through51.00004.00000+O.O555i01.02502.00007.00000+0.0550i01.02503.00009.00000+0.0605i00.97504.00005.00000.0101+0.0860i0+0.089611.00004.00006.00000.0168+0.0909i0+0.078311.00005.00007.00000.0160+0.0810i0+0.0759i1.00006.00009.00000.0310+0.1347i0+0.1420i1.00007.00008.00000.0085+0.0727i0+0.0758i1.00008.00009.00000.0118+0.1008i0+0.1051i1.0000Column6000000000B2=Columns1through52.0000+1.2395i01.05001.000001.8000+1.1155i01.00001.030001.0000+0.6197i01.00001.03000001.00000001.2500+0.5000i1.00000000.9000+0.3000i1.000000001.00000001.0000+0.3500i1.000000001.000000Column61.00003.00003.00002.00002.000()2.00002.(XX)02.00002.0000迭代次数4没有达到精度要求的个数16 16 1 0各节点的实际电压标幺值E为(节点号从小到大排列)：Columns1through51.0500 1.0176+0.1595i 1.0241+O.HOOi1.0402•0.0213i 1.0065-

Columns6through91.0062-0.0383i 1.0299+0.0617i 1.0072+0.0163i 1.0116+0.051Oi各节点的电压大小V为（节点号从小到大排列）：Columns1through81.0500 1.0300 1.0300 1.0404 1.0069 1.0070 1.0317Column91.0129各节点的电压角O为（节点号从小到大排列）：Columns1through80 8.9093 6.1298 -1.1722 -1.5602 -2.1784 3.4298Column92.8886各节点的功率S为（节点号从小到大排列）：Columns1through50.3928+0.6657i 1.8000+0.5248i 1.0000-0.1229i 0.0000+O.OOOOiColumns6through9-0.9000-0.3000i-0.0000+O.OOOOi-1.0000-0.3500i 0各条支路的首段功率Si为（顺序同您输入Bl时•样）：0.3928+0.6657i1.8000+0.5248i1.0000-0.1229i0.1284+0.3424i1.00730.9292-1.2500-0.5000i0.2644+0.2933i1.00730.9292-1.2500-0.5000i-0.6384+0.0600i0.6569+0.2433i-0.3472-0.0628i各条支路的末段功率Sj为（顺序同您输入Bl时一样）:-0.3928-0.6356i-1.8000-0.3426i-1.0000+0.1808i-0.1268-0.4228i-0.2616-0.3600i1.1431+0.0992i0.6514-0.1484i-0.6528-0.2872i0.3486-0.0324i各条支路的功率损耗DS为（顺序同您输入Bl时一样）:0.0000+0.030li0+0.1822i-0.0000+0.0579i0.0016-0.0805i0.0028-0.0667i0.0199+0.0220i0.0130-0.0884i0.0041-0.0438i以下是各节点每次迭代后的电压值（如图所示）第二组数据：#################################################################节点数：9节点支路数：9平衡节点：1精度：Bl=Columns1through51.00004.00000+O.O555i01.02502.00007.00000+O.O55Oi01.02503.00009.00000+0.0605100.97504.00005.00000.0101+0.0860i0+0.0896i1.00004.00006.00000.0168+0.0909i0+0.0783i1.00005.00007.00000.0160+0.08lOi0+0.0759i1.00006.00009.00000.0310+0.1347i0+0.1420i1.00007.00008.00000.0085+0.0727i0+0.0758i1.00008.00009.00000.0118+0.100810+0.105Hl.(XX)OColumn6000000000Columns1through52.0000+1.2395i 01.05001.000001.8000+1.1155i01.00001.03000l.(XX)0+0.6197i01.00001.04000001.00000001.2500+0.5000i1.00000000.9000+0.3000i1.00000000l.OO(X)000l.OOOO+O.35OOi1.000000001.000000Column61.00003.00003.000()2.00002.00002.00002.00002.00002.0000迭代次数没有达到精度要求的个数16 16I0各节点的实际电压标幺值E为（节点号从小到大排列）：Columns1through51.0500 1.0177+0.1589i 1.0343+0.1086i 1.0418-0.0212i 1.0083-0.0275iColumns6through91.0101-0.0385i 1.0317+0.0613i 1.0112+0.0158i 1.0186+0.0499i各节点的电压大小V为（节点号从小到大排列）：Columns1through81.0500 1.0300 1.0400 1.0420 1.0086 1.0108 1.0335 1.0113Column91.0198各节点的电压角O为（节点号从小到大排列）：Columns1through80 8.8729 5.9939 -1.1681 -1.5623 -2.1817 3.4029 0.8959Column92.8057各节点的功率S为（节点号从小到大排列）：Columns1through50.3921+0.6364i 1.8000+0.492li1.0000-0.0747i-0.0000+O.OOOOi-1.2500-0.5000iColumns6through9-0.9000-0.3000i 0.0000+O.OOOOi-1.0000-O.35OOi0各条支路的首段功率Si为（顺序同您输入Bl时一样）：0.3921+0.6364i1.0000-0.0747i0.1298+0.3405i0.2623+0.2678i-1.1218-0.0786i-0.6402+0.0365i0.6584+0.2118i-0.3456-0.0928i各条支路的末段功率Sj为（顺序同您输入Bl时一样）:-0.3921-0.6083i-1.8000-0.3116i-1.0000+0.1310i-0.1282-0.4214i-0.2598-0.3365i1.1416+0.0998i0.6530-0.1273i-0.6544-0.2572i0.3470-0.0037i各条支路的功率损耗DS为（顺序同您输入Bl时一样）:0.0000+0.028li-0.0000+0.1805i0+0.0562i0.0016-0.0809i

0.0198+0.0212i0.0128-0.0908i0.0040-0.0454i0.0014-0.0965i以下是各节点每次迭代后的电压值（如图所示）Figure1,|n|Figure1FileEditViewInsertToolsDesktopWindowHelp ，口安。昌©0争@I变口国I・亘"电压迭代次数曲线迭代次数DS1=O.0428-0.0825iDS2=0.0422-0.0813i以上两组数据均满足要求，从减小线路损耗的角度看，选择了第二组数据。潮流分布图:bus-8Ol391140.3295!bus-8Ol391140.3295!-0.3900*Jft3270对运行结果的分析线路有功功率总是从电压相位超前的结点流向滞后的结点，无功功率的总是从电压幅值高的结点流向电压幅值低的结点。功率数值前的符号若为正表示和箭头的方向一致，若为负则表示和箭头的方向相反。在潮流分布图中，3号发电机发出的无功功率为负值，说明3号发电机从系统吸收无功。由于电网运行时负荷等的变化，引起系统中各节点电压的变化。由于发电机容量有限，而电压的调整又必须服从对电压质量的要求，调整的幅度不能太大，另一方面从电网运行经济性和安全等方面考虑，网络中的潮流往往需要控制。在实际的网络潮流控制中主要采用：串联电容（作用以容抗抵偿线路感抗）；串联电抗（作用在于限流）；附加串联加压器（作用在于产生一环流或强制循环功率，使强制循环功率与自然分布功率的叠加可达到理想值）等手段控制潮流。三、故障电流计算1.三相短路电流的计算利用结点阻抗矩阵和导纳矩阵都可以计算短路电流，其算法有所不同。利用结点阻抗阵时,只要形成了阻抗阵，计算网络中任意一点的对称短路电流和网络中电流、电压的分布非常方便，计算工作量小，但是，形成阻抗阵的工作量大，网络变化时的修改也比较麻烦，而且结点阻抗矩阵是满阵，需要计算机存储量较大。利用结点导纳矩阵计算短路电流，实质是先用它计算出与指定的短路点有关的结点阻抗矩阵的相应列的元素，然后利用公式/f=―2一,q=(/,(0)-zj=Ui~Uj(z,为连接节点ijz.+z * jij 4和j的支路阻抗)进行短路电流的计算。然而，导纳阵是对称、稀疏阵，极易形成，且网络结构变化时也易于修改。NF=1n=9nl=9B=Columns1through51.00004.00000+0.0555i01.02502.00007.00000+0.0550i01.02503.00009.00000+0.0605i00.97504.00005.00000.01014-0.0860i0+0.0896i1.00004.00006.00000.0168+0.0909i0+0.0783i1.00005.00007.00000.0160+0.0810i04-0.0759i1.00006.00009.00000.0310+0.1347i0+0.1420i1.00007.00008.00000.0085+0.0727i04-0.0758i1.00008.00009.00000.0118+0.1008i0+0.1051i1.0000Column600000000vo=1.05001.03001.04001.04201.00861.01081.03351.01131.0198D=7 07点短路时电流的标幺值If=0.5504-4.636li各节点的电压标幺值U为（节点号从小到大排）：1.05001.03001.04001.04200.3594-0.0296i1.01080+1.3878e-017i-0.0222+O.OOOOi-0.0137+0.0000i各支路短路电流的标幺值I为（顺序同您输入B时一样）:0-0.602li-0.0000-18.7273i-0.0000-17.4223i1.2588-7.7890i0.0613-0.3319i0.4920-4.3402i1.6624-7.2232i0.0352-0.3012i-0.0097+0.0832i2.单相短路电流的计算请输入短路点的数目：NF=1请输入短路类型（单相接地短路为1,两相短路为2,两相短路接地为3）:Lf=l正序、负序、零序节点数为：nl=9n2=9n0=9正序、负序、零序支路数为:nil=9nl2=9nlO=9短路节点号:f=7纵横故障标志：If=0接地阻抗:zf=0接地阻抗:zg=0正序支路参数:Bl=Columns1through51.00002.00004.00007.00000+0.0555i0+0.0550i001.02501.02503.00009.00000+0.0605i00.97504.00005.00000.0101+0.0860i0+0.0896i1.00004.00006.00000.0168+0.0909i0+0.0783i1.00005.00007.00000.0160+0.08lOi0+0.0759i1.00006.00009.00000.0310+0.1347i0+0.1420i1.00007.00008.00000.0085+0.0727i0+0.0758i1.00008.00009.00000.0118+0.1008i0+0.1051i1.0000Column6000000000负序支路参数B2=Columns1through51.00004.00000+0.0555i01.02502.00007.00000+0.0550i01.02503.00009.00000+0.0605i00.97504.00005.00000.0101+0.0860i0+0.0896i1.00004.00006.00000.0168+0.0909i0+0.0783i1.00005.00007.00000.0160+0.08lOi0+0.0759i1.00006.00009.00000.0310+0.1347i0+0.1420i1.00007.00008.00000.0085+0.0727i0+0.0758i1.0000

8.00009.00000.0118+0.1008i0+0.105H1.0000Column6000000000B0=Columns1through51.00004.0