电力系统分析课程设计复杂网络牛顿—拉夫逊法潮流分析与计算的设计.docx

一、设计内容及要求复杂网络牛顿一拉夫逊法潮流分析与计算的设计电力系统潮流计算是电力系统中一项最基本的计算,设计内容为复杂网络潮流计算的计算机算法一牛顿-拉夫逊法。首先，根据给定的电力系统简图,通过手算完成计算机算法的两次迭代过程，从而加深对牛顿一拉夫逊法的理解，有助于计算机编程的应用。其次，利用计算机编程对电力系统稳态运行的各参数进行解析和计算；编程完成复杂网络的节点导纳矩阵的形成；电力系统支路改变、节点增减的程序变化；编程完成各元件的功率损耗、各段网络的电压损耗、各点电压、功率大小和方向的计算。二、设计原始资料1、给出一个六节点、环网、两电源和多引出的电力系统；2、给出一个五节点、环网、两电源和多引出的电力系统；参数给定，可以选用直角坐标表示的牛拉公式计算,也可以选用极坐标表示的牛拉公式计算.三、设计完成后提交的文件和图表1.计算说明书部分设计报告和手算潮流的步骤及结果2o图纸部分：电气接线图及等值电路；潮流计算的计算机算法,即程序;运算结果等以图片的形式附在设计报告中.四、进程安排第一周：第一天上午：选题，查资料，制定设计方案；第一天下午-第五天：复习潮流计算的计算机算法，完成给定网络的潮流计算;"f三l弟一J第一一-四天：利用matlab编程完成潮流计算，并对照手算结果，分析误差第五天下午：答辩，交设计报告。五、主要参考资料电力系统分析（第三版）于永源主编，中国电力出版社,2007年电力系统分析，何仰赞温增银编著，华中科技大学出版社，2002年版；电力系统分析，韩桢祥主编,浙江大学出版社,2001年版；电力系统稳态分析，陈玲编，水利电力出版社；复杂网络NR法潮流分析与计算的设计电力系统的潮流计算是电力系统分析课程基本计算的核心部分之一。它既有自身的独立意义,又有电力系统规划设计、运行和研究的理论基础,因此课程设计的重要性自不待言。一、设计题目1.系统图的确定选择六节点、环网、两电源和多引出的电力系统，简化电力系统图如图1所示,等值导纳图如图2所示。运用以直角坐标表示的牛顿一拉夫逊计算如图1所示系统中的潮流分布。计算精度要求各节点电压的误差或修正量不大于*=10-5o3.7÷j13图1电力系统图该系统中,节点为平衡节点，保持U三lo05+j0为定值,节点为PV节点其他四个节点都是PQ节点。给定的注入电压标幺值、线路阻抗标幺值、输出功率标幺值如下表注释。uIU2U3U4U5U61.05Io001.001O00Io00Io05各节点电压标幺值参数表1表2线路、变压器阻抗标幺值线路L2L3L4L5TlT2Y/2阻抗0.06+j.25Oo04+jOo250.08+j0.300.1+j0o35j0o03j0.015j0o25表3节点输出功率节点25功率2÷jl1o8+j0o4Io6+j0.83o7+j1035注:各PQ节点的电压取1是为了方便计算和最后验证程序的正确性二、N-R法的求解过程1、给定个节点电压初始值。(。)、/(°)2、将以上电压初始值代入下式(1)式，求出修正方程式常数项向量P(O)>Q(O)、U2(0)oP=P-西e(Ge-Bf+Be)lA,*J"jijiiii)UJJ3×2=e-G-/)-Xg+Be)lLjgji/jiijjijj、J=1U2=U2-(e2+f2),'ii4、将电压初始值代入下式(2)式,求出修正方程式中系数矩阵(雅可比矩阵)的元素(为2(n-1)阶方阵).当=i时，对角元素是 P = P-斗(Ge -Bf )+(g +Be )i i Liiii iij i ij i ij je.i=.1Gwj-Bijfj-Gaei-BiifiP)1.=乙SfBeBe-Gff.可JJuJ«"“Qj五(Gf+Be)+Be-Gfe.1Vjij)"i«*Q()A、=乙Se-Bfj+Ge+BfofUjUJ«»iiIij=T吗de工0笠ZLi时，矩阵中非对角元素是 p Q (r n _； =-L=AG e +B f /方 efij i ij id=-dIFe'de如2-rSf.J=B e -G f5、解修正方程式，求出修正量 函)、/(0)方法：高斯消元法，即将雅可比矩阵化简为上三角阵或下三角阵 求雅可比矩阵的逆矩阵修正各节点电压"= e(°)_Ae(°)/=/(0)- (0)6、将纳、/再代入(1 )式，求出P、。、¢/2(1)7、校验是否收敛，其收敛条件为(x)M=IA),Qm”式 AAA中JPw,Q(T为向量Rk)、Q（八）中大分量的绝对值。这个收敛条件比较直观，它可以直接显示出最终结果的功率误差A8、如果收敛，就进一步计算各段电力线路潮流和平衡节点功率，并打印出计算结果；如果不收敛，转回3步进行下一次迭代计算，直到收敛为止。三、手算潮流计算用图1和上表的数据和等值网络计算节点导纳矩阵yB10节点导纳矩阵yB由图1可知，该系统以串联支路的导纳标幺值和对地并联导纳标幺值得等值电路如图2所示。可得图2的等值节点导纳矩阵。对角线上的元素为：I K-Iy=_r÷Y=-733.33II TlK17T1I-Ay=_r+Y+y+y+y+y=1.53-j37.4222KTlK2TI20302325y=Jz+7'+y+y=1.74-y6.393330402334y=y+y.+y+y=1.58-5.254440503445y=-Y+U+y,+y+y+y=1.37-166.5155KTIK2T220502545y=Ijf+/=-j66.6766KTlKT2非对角线上的元素为:U(0)=1.00+j,u(O)=LO0+jU(O)=I.00+川U(O)=L00+jU(O)=I.05+j23456=/31.75y=y=0;y=y=0;y=y=0;y=y=0*z13J31，14*z41J1551*z1661%=%f91+378y=y=0;y=y=0,z24,422662yy0.62÷/3.9025524=%=-°83+Rlly=y=0;y=y=035533663y=y=-0.75+/2.644554J九二匕=04664y=y=/63.495665-j33331.75;31.751.53-137.420-0.91+73.78000-0.62+73.900100-0.91+73.781.74-y6.39-0.83+73.1100Y二00-0.83+3.111.58-75.25-0.75+;2.6400-0.62+73.900-0.75+/2.641.37-766.51/63.491。000;63.49;66.6712o计算各节点功率的修正方程的初始值（不平衡量）设各节点电压初始值为：u=.05÷yoI根据上述N-R法的求解过程，将各节点电压的初始值代入（1）式和（2）式,进行潮流计算,得节点功率和节点电压：P=-2;Q2=2.61;P=-i.8j=o.;33P=-1.6j(?=-0.30;4 4P=-3.7j=-5.40;5 5P=5.0;t2=0.00;66雅可比矩阵：-1.53-41.020.913.78000.623.9000-33.811.533.78-0.91003.90-0.64000.913.78-1.74-6.890.833.1100003.78-0.91-5.891.743.11-0.830000000.833.11-1.58-5.750.752.6400003.11-0.83-4.751.582.64-0.75000.623.90000.752.641.38-73.2106373.90-0.62002.64-0.75-59.811.3863.490000000066.670-63.4900000000-2.10修正各节点电压:(0) =-0.0952f(0) = 0.1232() = 0.0613f(o) = 0.6043e(o) = 0.1314f(0)= 0.6384() = - 0.0965f(0)= 0.1595e(o) = 06f(0) = 0.886四、计算机算法潮流计算U"-0.123e( D = 0 .9393f(1)= - 0 .6403e<D =0 .8694f(i)=-0.6384e<D = 1.0965f(i) =-0.1595e( 1)= 16f (I) =-0.886导纳矩阵YR=O0-33.3333i0+31.746000000+31.460i14.8252-39.55161-14.2012÷5.91720-0.624J÷3.9J0200-14.2012+5.9172115.0311-8.5292i-0.8299+3.112010000-0.8299+3.1120i1.5846-5.25351-0.7547÷2.6415100-0.6240+3.900210-0.754742.64161L3787-665103i0+63.4921i00000+63.492li0-66.666Htime=11. 1 7121.02 91 O 1 66 714.825244.7380-14.2012-5.91T200-0.6240-3.90020034.3653-14.8252-5.917214.201200-3.90020.624000-14.2012-5.917215.031i9.0292-0.8299-3.11200000-5.917214.20128.0292-15.0311-3.11200.8299000000-0.8299-3.11201.58465.7635-0.7647-2.64150000-3.11200.82994.7535-1.5846-2.64150.754700-0.6240-3.900200-0.7547-2.64151.378776.38290-63.4921-3.90020.624000-2.64150.764756,6377-1.3787-63.492100000000-69.8413063.4921000000002.200001.10000o117 4-O o 1 8 24-Oo 35 9 5-O . 0391O.0358Io 1 OOOo 171 2 -0.1 17 4iIo 02 910.18 24 1Oo 9 6 680.3595i1.1667Oo 03911Io 1000+0o 03581dU=0.1712-0.1174-Oo1824-O003320.3595Oo1667-0.03910Oo0358PQ=0 00 T 6 000 - 0. 3 0 0 0-2o00004o18631.8000Qo-3o70008o57265.0000time=2J20.615744.6939-17.3269-5.262700-1.1888-4.49450044.4695-23.3979-5.262717.326900-4.49451.188800-15.6934-3.498815.31516.4208-1.4217-3.05110000-3.498815,69345.6499-18.7327-3.05111.4217000000-1.9210-2.71042.00035.3261-1.6792-2.28250000-2.1041.92103.6928-4.8408-2.28251.679200-0.8805-4.526100-0.9838-3.05240.893277.1997-2.4809-74.0782-4.62610.880500-3.05240.983877.8918-7.6216-74.07822.48090000002.2711-69.84132.480974.0782000000002.20000.07151.10001.10081.1208Oo95780.80071.1317-0. 1 605 -Oo 3271 -Oo 0 7 8500.1321-0.0054U=Io 1 0 00 1.1208-0 . 13 21i0.9578-0o1605iOo8007-Oo3271i1.1317-000785iIo1008-0.0054id=-Oo05040.0147-0.0713Oo0219-0.16600.0324-0。0350-Oo03940.0008PQ二-Oo3576-Io03190.0288-Oo31500.06660.5211Oo1533-Io83330.3787-Oo0013time=320.179343.4131-16.6982-4.755800-1.2146-4.28880041.3278-23.5025-4.755816.698200-4.28881.214600-14.5516-3.388414.01226.4469-1.2943-2.8450000-3.388414.55165.0673-17.5190-2.84751.2943000000-1.6824-2.22051.59545.1261-1.4683-1.86830000-2.22051.68242.2509-4.3789-1.86831.46830Q-1.0123-4.365000-1.0615-2.93033.542176.4887-4.9828-71.8569-4.36501.012300-2.93031.061573.8403-10.0180-71.85694.9828000000-0.3458-69.88944.982871.8569000000002.2015-0.0109E：二Io1000lo11420.94380.75011.1271Io0999F：二0-Oo1351Oo1560-Oo3221-0.08890.0169Io10001.1142Oo1351iOo9438Oo1560i0.750I-Oo3221i1.1271-Oo0889ilo0999-0oO169idU二-O.0065-Oo0030-0.0140O.0045-O.0507Oo0050-O.0046-Oo0105-Oo0008-Oo0114PQ=O.OOOO-Oo0509-Oo0099-Oo02O7-Oo0154-Oo1040Oo0685-Oo0554-0.0935-OoOO17time=4JJ电力系统分析课程设计复杂网络牛顿一拉夫逊法潮流分析与计算的设计20.200543.1659-16.6230-4.674600-1.2222-4.26140040.9686-23.5242-4.674616.623000-4.26141.222200-14.3264-3.370013.72946.4249-1.2686-2.80780000-3.370014.32644.9870-17.3047-2.80781.2686000000-1.6248-2.06701.46355.0904-1.4168-1.73830000-2.06701.62481.770!-4.2976-1.73831.416800-1.0502-4.340400-1.0856-2.91024.295876.2461-5.6463-71.5628-4.34041.050200-2.91021.085673.4380-10.6417-71.56285.6463000000-1.0715-69.83685.646371.5628000000002.1999-0.0338E=IoIo099910001O1136Oo94240.74361.1267F=10Oo0181Oo13550,1555-Oo'32150O090U=1.100OIo1136-0o1355iOo9424Oo1555i0o7436-0.3215iIo12670.09Oli1.0999-0.0181i-0o0006-0o0004-0.0015Oo0005-Oo0065O00006-0.0005-0.0012-Oo0001-Oo0012PQ二1 o 1 0 001 . 1 1 36 Oo 9423Oo 7 43 51 . 12670.0002-Oo00040.0006-0.，0003o0001-0.0024-Oo0111OoOOll0.0003-Oo0028-0time=5JJ=20.205643.1382-16.6166-4.665900-1.2233-4.25880040.9379-23.5298-4.665916.616600-4.25881.223300-14.3024-3.368513.69926.4209-1.2658-2.80360000-3.368514.30244.9810-17.2818-2.80361.2658000000-1.6175-2.04721.44605.0870-1.4103-1.72150000-2.04721.61751.7069-4.2882-1.72151.410300-1.0546-4.337900-1.0884-2.90814.377276.2180-5.7232-71.5341-4.33791.054600-2.90811.088473.4028-10.7200-71.53415.7232000000-1.1491-69.83195.723271.53410E=00000002.1997-0.0362F=0-Oo1355-0.1554-0.3215-Oo0902-Oo0181U=1o1000Io1136-0o1355iOo9423-0.1554iOo74350.3215i1.12670.0902i1.0999-0.0181idU=lo0e-003大-0.0104-Oo0063-Oo02550.0090-0.1157O00106-O00077-Oo0207-Oo00100.0207PQ=1oOe-003*0. 0 0 0 4-0 .0010-0.0059 0o 0 0 0 3 -Oo0439 -Oo 1922Oo0.00 0 20 28 90.0 0 1 5timeJJ20.205743.1377-16.6165-4.665700-1.2233-4.25880040.9374-23,5299-4.665716.616500-4.25881.223300-14.3020-3.368513.69876.4209-1.2658-2.80360000-3.368514.30204.9809-17.2814-2.80361.2658000000-1.6174-2.04691.44576.0869-1.4102-1.72120000-2.04691.61741.7058-4.2880-1.72121.410200-1.0547-4.337900-1.0885-2.90804.378676.2175-5.7245-71.5336-4.33791.054700-2.90801.088573.4022-10.7214-71.53365.7245000000-1.1504-69.83185.724571.5336000E=1.100011.0999F=00.O1136Oo135500.9-0.0423155400Oo7435-Oo3215-2.1997一0.03621O1267Oo0902U=Io10001136-0.1355i0.9423-Oo1554i0.74350.3215iIo1267-0。0902i1.0999-Oo0181idU=1.Oe-007*0.03280.0198-Oo08090.02860.36740.03350.0241-0.0654-Oo0030-Oo0652PQ=1.Oe-OO7*0.0007-0o0026-0.01630.0023-0.14210.6090Oo0149-Oo0022-0.0872-Oo0043平衡节点功率S1=4。7307+1.4446i节点间线路上的流动功率矩阵S二0-4.7307+2.3966i0+1.9206i0+1.9206104-l.920i0+1.9206i4.7307-3.5963iO-2.9395-2.999110-2.531810.2089-2.531810-2.6318i0-0,4561i2.5172-0.1647iO-0.7172-0.B0350-0.450Ii0-0.4561i0-0,3280iO-0.3280i0.6615-0.1895100.9385÷0.0064i0-0.3280iO-4,5011i-0.2103-4.5098i0-4.60111-1.0897-5.363705.0000-6.83091O+3,8413i0+3.8413i0+3.841310+3.84131-5.0000+5.79590刊U山七尔所午心仪枇dataS=00+5.9929i0+2.3767i0+2,2487i0+6,4217i0-1.920610-59929i0-0.4224-2.8344i0-2,2038i-0.D014+i.978ii0-63731i0-2.3767i-0.4224+2.83440-0.0557-0.6141i0+4,04500-4.29730-2.2487i0+2.2038-0.055740.61410-0.1512+5.36910-4.16930-6.4217i-0.0014-1.97810-4.0450-0.1512-5.369Ii00-12.62680+1,9206i0+6.373110+4.2973i0+4,1693i0+12.6268i0个人心得经过两周的课程设计让我学到了很多以前没有学到的东西，首先让我明白了电力系统的稳定运行，必须经过精密的设计和运算。在进行此次课程设计的过程中,加深了我对潮流计算的认识，尤其是对牛顿拉夫逊潮流计算的求解思路有了比较透彻的理解。同时整个求解过程中用到求节点导纳矩阵求矩阵的逆阵等等，又让我们遇到了很多问题，我们经过查资料问同学、问老师,最后才把问题解决。使我对以前学过的知识有了一次很好的复习，同时也看到了我知识的欠缺，让我明白我应该学的东西还有很多，更是给我更多的激励让我以后会更加努力地学习.此次课程设计有让我学到了一本新的知识Matlab的使用，起初我们经过好大力气才把程序设计出来，可是设计的题目里的数据又不收敛,我们反复演算了好多次，最后在老师的指导下改了几个数据程序终于收敛了。而且在此次课程设计中，我发现了自己的基础只是有很多的不足。这些基础的缺乏给我的设计造成了不小的障碍，今后我会努力学习来弥补这些不足。在这个过程中，我明白了，只要用心去做，认真去做，同学一起努力、团结，就会克服所有困难,取得最后的胜利。看到同学和自己设计的程序稳定运行,心中也有一种成功的喜悦。通过此次的实习让我发现手工的潮流计算和计算机潮流计算的优缺,对于小系统来说手工计算很简单，手工也很经济。可对于大系统来说手工计算是不可完成的任务,因此用计算机计算潮流是必不可少的，只要编好程序计算机就会给你一个自己想要的满意答案。而且方便、简单、准确.此次课程设计我能成功设计完成离不开同学们的帮助、团队的合作以及老师细心的指导，在这里再次感谢老师的指导和同学的帮助。附录：yll=1/0.03j；yl2=1/(Io05*0.03j)；y13=0;y14=0；yl5=0;y16=0;y21=yl2;y22=0.5j+l/(0.06+0o025j)+l/(0.04+0,25j)+l(lo052)*Oo03j)；y23=-1/(0o06+0o025j);y24=0；y25=-1/(0o04+0.25j)；y2=0；y3l=y13；y32=y23；y33=0o5j+l(0o06+0.025j)+l(0.08+0o3j)；y34=-1/(0o08+0.3j)；y35=0;y36=0;y41=y14；y42=y24；y43=y34;y44=0o5j+l(0.08+0,3j)+l(0.1÷0.35j)；y45=-1/(0.1+0.35j)；y46=0;y51=yl5；y52=y25；y53=y35；y54=y45；y55=0o5j+l(0ol+0o35j)+1/(0o04+0o25j)+1/(1.052)*0o015j)；y56=-1/(1.05*0.015j)；y61=0；y62=0;y63=0;y64=0;y65=y56;y66三l0o015j；%输入导纳参数YB=y11y12yl3y14yl5y16;y21y22y23y24y25y26;y31y32y33y34y35y36；y41y42y43y44y45y46；y51y52y53y54y55y56;y61y62y63y64y65y66%形成导纳矩阵%*牛顿一拉夫逊法潮流计算夫夫夫夫*E(2)=1o00；E(3)=Io00；E(4)=1。00；E(5)=1.00;E(6)=1.1；F(2)=0jF(3)=0;F(4)=0；F(5)=0；F(6)=0；%设定电压初始值G=real(YB)B=imag(YB)；%提取导纳实部,虚部分别为G1BS(2)=-2-lj;S(3)=-1.8-0.4j;S(4)=Io6-0。8j;S(5)=-3.7-1.3j;S=5;%设定注入功率初始值P=rea1(三)mag(三);%提取注入功率实部，虚部分别为P,Qk=0;preciSion=Oo1；%设定参数k为迭代次数初值为0及精度大于0.00001Nl=5；驰设定除平衡节点外节点数之和wh1IepreCiSion>Oo00001E(I)=I01；F(I)=0;%输入平衡节点电压实部，虚部e.fform=2:N1÷1forn=l:Nl+lPt(n)=(E(m)*(G(m,n)*E()一B(m,n)*F(n)+F(m)*(G(m1n)*F(n)+B(m1n)*E(n);Qt()=(F(m)*(G(m,n)*E()B(r,n)*F(n)-E(m)*(G(m,n)*F(n)÷B(m,n)*E(n)；%以设定值计算的注入功率enddP(m)=P(m)-Sm(Pt);dQ(m)=Q(m)-sum(Qt)；%计算功率误差。end。form=2:N1+1fOrn=l:Nl+l。Ai(n)=G(mtn)*E(n)B(m,)*F()；Bi()=G(m1n)*F(n)+B(m,)*E()；endN(m,m)=sum(Ai)+G(m1m)*E(m)+B(mtm)*F(m);%计算dP/deH(m1m)=SUm(Bi)-B(m,m)*E(m)+G(m,m)*F(m);%计算dP/df1.(m,m)=-Sum(Bi)一B(m,m)*E(m)÷G(m,m)*F(r);%计算dQ/deJ(m,m)=sum(A1)-G(m,m)*E(m)B(m,m)*F(m);%计算dQ/dfendform=2:N1+1JJ(2*m3,2*-3)=N(m1m);JJ(2*m-3,2*m-2)=H(m,m);JJ(2*m-2,2*m一3)=L(m,m)；JJ(2*m-212*r2)=J(m,m);end%给雅可比矩阵的对角元素赋值form=2:Nl+1forn=2:Nl+1-jfm=nelSeN(mtn)=G(m,n)*E(m)+B(m,n)*F(m);%计算dP/deH(m,n)=B(m1)*E(m)+G(m,)*F(m);%计算dPdf1.(m,n)=B(m,)*E(m)+G(m,)*F(m);%计算dQdeJ(m,n)=-B(m,)*F(m)G(m,)*E(m);%计算dQ/df%*给雅可比矩阵的非对角元素赋值*JJ(2*m-3,2*n-3)=N(mtn);JJ(2*m-3,2*-2)=H(mr);JJ(2*m-2,2*n-3)=L(m,);JJ(2*m-212*n-2)=J(m,n);endendendfor1=1:8JJ(10,i)=0；%修正雅可比矩阵第十行1,到8元素为OendJJ(10,9)=2*E(6);驰修正雅可比矩阵第十行9元素为2eJJ(10,10)=2*F(6);%修正雅可比矩阵第十行10元素为2fform=l:NlPQ(2*m1)=dP(m+l);PQ(2*m)=dQ(m+l);end%写PQ矩阵，奇数行为dP,偶数行为dQPQ(1。)=1.r2-(E(6)2+F八