matlab电力系统潮流计算

1、华中科技大学信息工程学院课程设计报告书题目:电力系统潮流计算专业：电气工程及其自动化班级：学号：学生姓名：指导教师:2015年11月10日学号学生姓名专业（班级）电气工程及其自动化（）设计题目电力系统潮流计算设计技术参数1. 变压器的阻抗，励磁损耗。2. 线路的阻抗，无功损耗。3. 变压器及线路的功率损耗4. 变压器及线路的电压降落5. 各节点的功率和电压设计要求1. 计算过程包括理论计算和MATLA昉其沏大部分。2. 报告里面的理论分析（手工计算）过程必须要有详细推导过程，再把仿真过程放到相应的部分中去。工作量共24页；4000字工作计划第15周：对系统模型分析和资料收集，找到计算的思路，总

2、结计算方法；第16周：对整个模型进行详细计算，并完成word文档。参考资料1 何仰赞，温增银电力系统分析（第三版）M.华中科技大学，20022 刘同娟.MATLAB在电路分析中的应用.电气电子教学学报.20023 王守相，刘玉田电力系统潮流计算研究现状-山东电力技术1996年05期指导教师签字2015年11月12日信息工程学院课程设计成绩评定表学生姓名：学号：专业（班级）：班课程设计题目：电力系统潮流计算指导教师评语：成绩：指导教师：年月日电力系统稳态分析包括潮流计算和静态安全分析。本文主要运用的事潮流计算，潮流计算是电力网络设计与运行中最基本的运算，对电力网络的各种设计方案及各种运行方式进行

3、潮流计算，可以得到各种电网各节点的电压，并求得网络的潮流及网络中的各元件的电力损耗，进而求得电能损耗。本位就是运用潮流计算具体分析，并有MATLA所真。关键词：电力系统潮流计算MATLAB仿真AbstractElectricpowersystemsteadyflowcalculationandanalysisofthestaticsafetyanalysis.Thispaper,bymeansofthecalculation,flowcalculationisthetrendofthepowernetworkdesignandoperationofthemostbasicoperationso

4、felectricpowernetwork,variousdesignschemeandtheoperationwaystotidecomputation,cangetallkindsofeachnodeofthepowergridvoltageandseekthetrendofthenetworkandthenetworkofthecomponentsofthepowerloss,andgettingelectricpower.Thestandardistousethepowerflowcalculationandanalysis,thespecifichaveMATLABsimulatio

5、n.Keywords:Powersystem;Flowcalculation;MATLABsimulation1任务提出与方案论证2总体设计12.1潮流计算等值电路2.2建立电力系统系型2.3模型的调试与运行3详细设计3.1计算前提3.2手工计算4设计图及源程序4.1MATLA盼真4.2潮流计算源程序5总结错误！未定义书签。参考文献1任务提出与方案论证潮流计算是在给定电力系统网络结构、参数和决定系统运行状态的边界条件的情况下确定系统稳态运行状态的一种基本方法，是电力系统规划和运营中不可缺少的一个重要组成部分。可以说，它是电力系统分析中最基本、最重要的计算，是系统安全、经济分析和实时控制与调度的

6、基础。常规潮流计算的任务是根据给定的运行条件和网路结构确定整个系统的运行状态，如各母线上的电压（幅值及相角）、网络中的功率分布以及功率损耗等。潮流计算的结果是电力系统稳定计算和故障分析的基础。在电力系统运行方式和规划方案的研究中，都需要进行潮流计算以比较运行方式或规划供电方案的可行性、可靠性和经济性。同时，为了实时监控电力系统的运行状态，也需要进行大量而快速的潮流计算。因此，潮流计算是电力系统中应用最广泛、最基本和最重要的一种电气运算。在系统规划设计和安排系统的运行方式时，采用离线潮流计算；在电力系统运行状态的实时监控中，则采用在线潮流计算。是电力系统研究人员长期研究的一个课题。它既是对电力系

7、统规划设计和运行方式的合理性、可靠性及经济性进行定量分析的依据，乂是电力系统静态和暂态稳定计算的基础。潮流计算经历了一个由手工到应用数字电子计算机的发展过程，现在的潮流算法都以计算机的应用为前提用计算机进行潮流计算主要步骤在于编制计算机程序,这是一项非常复杂的工作。对系统进行潮流分析，本文利用MATLAB中的SimpowerSystemsX具箱设计电力系统，在simulink环境下，不仅可以仿真系统的动态过程，还可以对系统进行稳态潮流分析。2总体设计SimpowerSystems使用Simulink环境，可以将该系统中的发电机、变压器，线路等模型联结起来，形成电力系统仿真模拟图。在加人测量模块

8、，并对各元件的参数进行设置后，用measurement和sink中的仪器可以观察各元件的电压、电流、功率的大小。2.1潮流计算等值电路2.2建立电力系统模型在Simulink中按照电力系统原型选择元件进行建模。所建立的模型和建立的方法在详细设计中详述。在电力系统模型的建立工程中主要涉及到的是：元器件的选择及其参数的设置；发电机选型；变压器选择；线路的选择；负荷模型的选择；母线选择。2.3模型的调试与运行建立系统模型，并设置好参数以后，就可以在Simulink环境下进行仿真运行。运行的具体结果和分析也在详细设计中详述。3详细设计3.1计算前提首先是发电机的参数计算，先对5个发电厂简化为5台发电机

9、来计算。发电机G1:发电机G2:发电机G3:发电机G4:发电机G5:其次是变电站的参数计算，我们还是对7个变电站简化为7台变压器来计算。变压器T1:变压器T2:（双并联）变压器T3:（四并联）变压器T4:（双并联）变压器T5:变压器T6:（两个三绕组变压器并联）变压器T7:（双并联）再次是传输线参数计算，5条传输线的具体计算如下。根据教材查得r0=0.21，.i/kmxo=04/kmb0=2.810S/km线路L1:线路L2:线路L3:（双回路）线路L4:线路L5:（双回路）3.2手工计算FLR1:P2102.S1=(RT1jXT1)=(3.450j74.406)=(0.0285j0.6562

10、)MVAVn110Sa=10MWS1.jQL1=(10.0442j0.1142)MVA_2_2_2_2P2Q210.044220.11422§1=-(%jX)=而(8.4j16)=(0.070j0.1334)MVA§2=(禺2jXT2)=(1.346j31.7625)=(0.4032j9.5156)MVAVn110Sb=SG1-20-§2=60j45-20-0.4032-j9.5156=(39.5968j35.4844)MVAFLR2Sc=Sb-Sa-25-jQL1-SL1=(4.4826j35.9144)MVA_2_2_2.2-PQ252156S3=(Rt3jX

11、T3)=2(0.092j5.042)=(0.6679j36.6024)MVAVn1103:Sc=(4.4931j34.1048)MVAP2Q2.4.4931234.1048SL2=-TT2(RL2,jXL2)-L2L2.Vn110_，一一一一一Sd=宅2Sc-120-§3-：S03-jQL2-SL22(27.3j52)=(2.67j5.0854)MVA=(132.9792j149.229)MVAFLR2_222cPQ627Vn2S4=(RmjXT4)=(1.725j39.703)=(0.1091j2.5101)MVA1102Sd=(133.5955j149.9956)MVA_2_2一

12、_2_一.2-PQ133.5955149.99561102.：SL3=(RL3jXL3)=2(7.35j14)=(24.51j46.682)MVAVnSe=SG3Sd一30-25-.：&4一£。4一jQL3£L3=(89.945j130.0151)MVAFLR4：P2Q2502312(Rt5jXT5)(0.368j20.168)=(0.1052j5.7687)MVA_，Se=(92.7481j133.9937)MVAP2Q292.74872133.99372S4=.,2(Rl4jXL4)=2(12.6j24)=(27.654j52.674)MVA-,：ST5Vn2_

13、21102Vn2_-'.-.-_一_Sf=Sg4Se-80-$5-S05-jQL4-Sl4=(34.9449j107.3469)MVAFLR5:计算每一个FLR的功率分布和电压分布计算如下:FLR1:FLR2:功率分布:*QZT3Q(Vb-Vn)、/Sl2一*Sd*VnZL2ZT3ZL2ZT3=(4.8812-j13.8097)MVA*"SdVnZL2ZT3ZL2ZT3=(108.687-j122.62)MVA(0.092-j5.042)(132.9792j149.229)-1418.6727.392-j57.042(27.3-j52)(132.9792j149.229)-1

14、418.6727.392-j57.042电压分布:Sc1=SL2：SL2=(4.8812-j13.8097)(2.67j5.0854)=(7.5512-j8.7243)MVAFLR3:、，7.551227.3-8.724352VL2=-2.424kV102.1030V=Vb_VL2=102.103-(-2.424)=104.527kV功率分布:*"*Zt*4&"Vd)VnZl3Zt4Zl3Zt4=(59.444一j16.846)MVA*SZl3s(VbM)V54一*De*VnZL3，ZT4ZL3，Zt4=(31.811j60.1256)MVA(1.725-j39.7

15、03)(89.945j130.0151)-1037.9279.075-j53.73(7.35-j14)(89.945j130.0151)1037.9279.075-j53.73电压分布:鬲=Sl3Sl3=(59.444-j19.846)(24.51j46.682)=(83.954j26.836)MVA、，83.9547.3526.83614.VL3=9.404kVFLR4:105.5643M=Vd-VL3=96.16kV功率分布:*z_S*fZL4ZT5(Vs3-Vd)_(0.368-j20.168)(34.9449j107.3469)-1037.927-*Vn=ZL4ZT512.968-j4

16、4.168=(20.843j19.689)MVA*QZl4u_(VG3-Vd)/_(12.6-j24)(34.9449j107.3469)-1037.927为4*Sf*VN=ZL4ZT5ZL4ZT5=(1.398j44.389)MVA12.968-j44.168电压分布:FLR5:这里我们先将f点和发电机G5当做电源，经过ZT61和ZT63构成两端供电网络以g点作为运算负荷进行计算。京=(0.386-的(20.2656"砸制(22.0938-37)35=(3.900.j25.1175)MVA0.772-j10.718St4(0.386-j6.584)(20.2656j70.9293)

17、-(22.0938-37)350.772-j10.718=(16.5061j91.7905)MVA电压分布:4设计图及源程序4.1MATLAB"相关的原始数据输入格式如下：1、B1是支路参数矩阵，第一列和第二列是节点编号。节点编号由小到大编写。2、对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点编号，将变压器的串联阻抗置于低压侧处理，第三列为支路的串列阻抗参数，第四列为支路的对地导纳参数，第五烈为含变压器支路的变压器的变比，第六列为变压器是否是否含有变压器的参数，其中1”为含有变压器，0”为不含有变压器。3、B2为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列

18、为节点负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中1”为平衡节点，2”为PQ节点，3”为PV节点参数。4、X为节点号和对地参数矩阵。其中第一列为节点编号，第二列为节点对地参数。4.2潮流计算源程序closeallclearalln=input('请输入节点数:n=');n1=input('请输入支路数：n1=');isb=input('请输入平衡节点号：isb=');pr=input('请输入误差精度:pr=');B1=input('请输入支路参数:B1=');B2=input('请输入节

19、点参数:B2=');X=input('节点号和对地参数:X=');Y=zeros(n);Times=1;%置迭代次数为初始值%建节点导纳矩阵fori=1:n1ifB1(i,6)=0%不含变压器的支路p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/B1(i,3);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4);Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4);else%含有变压器的支路p=B1(i,1);q=B1(i,2);Y(p,q)=Y(p,q)-1/(B1(i,3)*B

20、1(i,5);Y(q,p)=Y(p,q);Y(p,p)=Y(p,p)+1/B1(i,3);Y(q,q)=Y(q,q)+1/(B1(i,5)A2*B1(i,3);endYOrgS=zeros(2*n-2,1);DetaS=zeros(2*n-2,1);%将OrgS、DetaS初始化唯，建OrgS,用于存储初始功率参数h=0;j=0;fori=1:n%对PQ节点的处理ifi=isb&B2(i,6)=2h=h+1;forj=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*ima

21、g(B2(j,3)+imag(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3);OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)-real(B2()*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3);endendendfori=1:n%对PV节点的处理，注意这时不可再将h初始化为0ifi=isb&B2(i,6)=3h=h+1;forj=1:nO

22、rgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)+imag(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3);OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)-real(B2()*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j

23、,3);endendendOrgS嗾J建PVU用于存储PV节点的初始电压PVU=zeros(n-h-1,1);t=0;fori=1:nifB2(i,6)=3t=t+1;PVU(t,1)=B2(i,3);endendPVU%建DetaS,用于存储有功功率、无功功率和电压幅值的不平衡量h=0;fori=1:n%对PQ节点的处理ifi=isb&B2(i,6)=2h=h+1;DetaS(2*h,1)=imag(B2(i,2)-OrgS(2*h,1);endendt=0;fori=1:n%对PV节点的处理，注意这时不可再将h初始化为0ifi=isb&B2(i,6)=3h=h+1;t=t+

24、1;DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1);DetaS(2*h,1)=real(PVU(t,1)A2+imag(PVU(t,1)A2-real(B2(i,3)A2-imag(B2(i,3)A2;endendDetaS%建I,用于存储节点电流参数i=zeros(n-1,1);h=0;fori=1:nifi=isbh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1)/conj(B2(i,3);endendI%建Jacbi(雅可比矩阵)Jacbi=zeros(2*n-2);h=0;k=0;fori=1:n%对PQ

25、节点的处理ifB2(i,6)=2h=h+1;forj=1:nifj=isbk=k+1;ifi=j%对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3)+imag(I(h,1);Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+real(I(h,1);Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-

26、2*imag(I(h,1);else%非对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3);Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3);Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);endifk=(n-1)%将用于内循环的指针置于初始值，以确保雅可比矩阵换行k=0;endendendendk=0;fori=1:n%

27、对PV节点的处理ifB2(i,6)=3h=h+1;forj=1:nifj=isbk=k+1;ifi=j%对角元素的处理Jacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3)+imag(I(h,1);Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+real(I(h,1);Jacbi(2*h,2*k-1)=2*imag(B2(i,3);Jacbi(2*h,2*k)=2*real(B2(i,3);else%非对角元素的处理Jacbi

28、(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3);Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3);Jacbi(2*h,2*k-1)=0;Jacbi(2*h,2*k)=0;endifk=(n-1)%将用于内循环的指针置于初始值，以确保雅可比矩阵换行k=0;endendendendendJacbi%求解修正方程，获取节点电压的不平衡量DetaU=zeros(2*n-2,1);DetaU=inv(Jacbi)*DetaS;DetaU%修

29、正节点电压j=0;fori=1:n%对PQ节点处理ifB2(i,6)=2j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endendfori=1:n%对PV节点的处理ifB2(i,6)=3j=j+1;B2(i,3)=B2(i,3)+DetaU(2*j,1)+DetaU(2*j-1,1)*sqrt(-1);endB2初始循环*whileabs(max(DetaU)>prOrgS=zeros(2*n-2,1);%!初始功率参数在迭代过程中是不累加的，所以在这里必须将其初始化为零矩阵h=0;j=0;fori=1:nifi=i

30、sb&B2(i,6)=2h=h+1;forj=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)+imag(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3);OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)-real(B2()*(real(Y(i,j)*ima

31、g(B2(j,3)+imag(Y(i,j)*real(B2(j,3);endendendfori=1:nifi=isb&B2(i,6)=3h=h+1;forj=1:nOrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3)*(real(Y(i,j)*real(B2(j,3)-imag(Y(i,j)*imag(B2(j,3)+imag(B2(i,3)*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3);OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3)*(real(Y(i,j)*real(B2

32、(j,3)-imag(Y(i,j)*imag(B2(j,3)-real(B2()*(real(Y(i,j)*imag(B2(j,3)+imag(Y(i,j)*real(B2(j,3);endendendOrgS%建DetaSh=0;fori=1:nifi=isb&B2(i,6)=2h=h+1;DetaS(2*h-1,1)=real(B2(i,2)-OrgS(2*h-1,1);DetaS(2*h,1)=imag(B2(i,2)-OrgS(2*h,1);endendt=0;fori=1:nifi=isb&B2(i,6)=3h=h+1;t=t+1;DetaS(2*h,1)=real(

33、PVU(t,1)A2+imag(PVU(t,1)A2-real(B2(i,3)A2-imag(B2(i,3)A2;endendDetaS唯，建Ii=zeros(n-1,1);h=0;fori=1:nifi=isbh=h+1;I(h,1)=(OrgS(2*h-1,1)-OrgS(2*h,1)*sqrt(-1)/conj(B2(i,3);endendI%建JacbiJacbi=zeros(2*n-2);h=0;k=0;fori=1:nifB2(i,6)=2h=h+1;forj=1:nifj=isbk=k+1;ifi=jJacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(

34、i,3)+real(Y(i,j)*imag(B2(i,3)+imag(I(h,1);Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+real(I(h,1);Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k)+2*real(I(h,1);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1)-2*imag(I(h,1);elseJacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3);Jacbi(2*

35、h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3);Jacbi(2*h,2*k-1)=-Jacbi(2*h-1,2*k);Jacbi(2*h,2*k)=Jacbi(2*h-1,2*k-1);endifk=(n-1)k=0;endendendendendk=0;fori=1:nifB2(i,6)=3h=h+1;forj=1:nifj=isbk=k+1;ifi=jJacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3)+imag(I(h,1);Jacbi(2*h-1,2*k)=real(Y(i,j)*real(B2(i,3)+imag(Y(i,j)*imag(B2(i,3)+real(I(h,1);Jacbi(2*h,2*k-1)=2*imag(B2(i,3);Jacbi(2*h,2*k)=2*real(B2(i,3);elseJacbi(2*h-1,2*k-1)=-imag(Y(i,j)*real(B2(i,3)+real(Y(i,j)*imag(B2(i,3);Jacbi(2