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

1、电力系统分析课程设计复杂网络牛顿拉夫逊法潮流分析与计算的设计一、设计内容及要求复杂网络牛顿拉夫逊法潮流分析与计算的设计电力系统潮流计算是电力系统中一项最基本的计算,设计内容为复杂网络潮流计算的计算机算法-牛顿拉夫逊法。首先，根据给定的电力系统简图,通过手算完成计算机算法的两次迭代过程，从而加深对牛顿拉夫逊法的理解，有助于计算机编程的应用。其次，利用计算机编程对电力系统稳态运行的各参数进行解析和计算；编程完成复杂网络的节点导纳矩阵的形成；电力系统支路改变、节点增减的程序变化；编程完成各元件的功率损耗、各段网络的电压损耗、各点电压、功率大小和方向的计算。二、设计原始资料、给出一个六节点、环网、两电

2、源和多引出的电力系统；2、给出一个五节点、环网、两电源和多引出的电力系统；参数给定，可以选用直角坐标表示的牛拉公式计算,也可以选用极坐标表示的牛拉公式计算三、设计完成后提交的文件和图表1计算说明书部分设计报告和手算潮流的步骤及结果2。图纸部分：电气接线图及等值电路;潮流计算的计算机算法,即程序;运算结果等以图片的形式附在设计报告中.四、进程安排第一周:第一天上午：选题，查资料，制定设计方案;第一天下午第五天：复习潮流计算的计算机算法，完成给定网络的潮流计算;第二周:第一-四天：利用matlb编程完成潮流计算，并对照手算结果，分析误差第五天下午：答辩，交设计报告。五、主要参考资料电力系统分析(第

3、三版） 于永源主编，中国电力出版社,207年电力系统分析，何仰赞 温增银编著，华中科技大学出版社，002年版；电力系统分析，韩桢祥主编,浙江大学出版社,201年版;电力系统稳态分析,陈珩 编，水利电力出版社;复杂网络法潮流分析与计算的设计电力系统的潮流计算是电力系统分析课程基本计算的核心部分之一。它既有自身的独立意义,又有电力系统规划设计、运行和研究的理论基础,因此课程设计的重要性自不待言。一、 设计题目1.系统图的确定选择六节点、环网、两电源和多引出的电力系统，简化电力系统图如图所示,等值导纳图如图2所示。运用以直角坐标表示的牛顿拉夫逊计算如图所示系统中的潮流分布。计算精度要求各节点电压的误

4、差或修正量不大于。 图1 电力系统图图2 电力系统等值导纳图2。各节点的初值及阻抗参数该系统中,节点为平衡节点，保持=1。05+0为定值,节点为PV节点,其他四个节点都是P节点。给定的注入电压标幺值、线路阻抗标幺值、输出功率标幺值如下表注释。表1 各节点电压标幺值参数UUU101。0010。001。01。05表2 线路、变压器阻抗标幺值 线路L2LL4T1T22阻抗00j020。00。25008j0.300.1+。35j0。j0.1j0。25表3 节点输出功率节点功率2+j1。+j0。41。6+j0.8。7+。35注:各Q节点的电压取1是为了方便计算和最后验证程序的正确性.二、 N法的求解过程

5、1、 给定个节点电压初始值2、 将以上电压初始值代入下式（1）式，求出修正方程式常数项向量。3、 4、 将电压初始值代入下式（）式,求出修正方程式中系数矩阵（雅可比矩阵)的元素(为（n1）阶方阵)5、 解修正方程式，求出修正量方法：高斯消元法，即将雅可比矩阵化简为上三角阵或下三角阵求雅可比矩阵的逆矩阵 修正各节点电压6、将再代入（）式，求出7、校验是否收敛，其收敛条件为式中,为向量中大分量的绝对值。这个收敛条件比较直观，它可以直接显示出最终结果的功率误差、如果收敛，就进一步计算各段电力线路潮流和平衡节点功率，并打印出计算结果；如果不收敛，转回步进行下一次迭代计算，直到收敛为止。三、手算潮流计算

6、用图1和上表的数据和等值网络计算节点导纳矩阵。 节点导纳矩阵 由图1可知，该系统以串联支路的导纳标幺值和对地并联导纳标幺值得等值电路如图2所示。可得图2的等值节点导纳矩阵。对角线上的元素为：非对角线上的元素为: Y=2。计算各节点功率的修正方程的初始值（不平衡量) 设各节点电压初始值为: 根据上述NR法的求解过程，将各节点电压的初始值代入（1）式和(）式,进行潮流计算,得节点功率和节点电压：雅可比矩阵:=修正各节点电压: 四、计算机算法潮流计算导纳矩阵=tme=1= = 100 1.712 1.02 0668 。66 1.1000= 0 。117 。24 0。355 -00391 .0358U

7、 =1。0 。171-0.14i 1。010.12 0。80.3595i1.16670。0391 1。1000+0。0358dU = 01 01174 0.01 0。182 0。0332 0.595 0。166 -00391 0 0。038PQ =-2。000 4。1863 .800 。00 -1000 0.-3。700 8。5726 5000 0tim = = 1.10 .108 0。957 0.807 317 .108F = 0 .1321 -0.605 0。3271 0。085 0054U=1。00 1.120-121i .9578-。65i 0。87-。7i1.1370。785i 1。

8、100.05iU 0。0 0147 00713 。09 0.1660 0324 。0350 -。9 0.000 0042 =0。356 -1。9 00288 0。1500.0666 0.51 。53 1。8333 0.37 0。0013time JJ =E = 1。10 1。12 .93 .7501 1.1271 1。099F = 0。351 。156 -0。321 -0.089 0.09U =1。100 1.1140。151i 0。9438。1560 0510。22 20。08i 1。99-0。169idU = -.0065 。000 .0 004 -007 。0050 -0006 0。10

9、5 -。008 0。0114Q 00 0。059 0。9 -。02 0。0154 。100 0。068 0。554 0.03 0。07te=4J=E = 1。0 。1136 0。9424 0.3 .1267 1。09F = 0。13550.55 -0。321 。090 0。081U 1000 1。136-。1355i 0。24。15 0。7430215 1。1270.001i .099-.0181d = 。006 -。0004 0.005 0。0005 0。05 0。006 -0.005 -0.0 。01 0。002PQ=000 -0。0 06 -0003 0.2 -0。011 0。011 0

10、003 。028 -0。000ime 5JJ =E= 。00 136 0。9423 0。43 1267 199 = 0 。 -.1554 0.325 。092 -0。0181U 。0 1。136-。5i。23-15i 0。4350.321i 116000i .99-0181idU 1。00 .10 0。006 0。055 000 0.1157 0。0 。00 。2 0。0 0.0207Q = 。0e-003 *0.0 -00010 -0.0059 。 -0。0439-0。1922 。 .002028 0.05ti=6J =E = 1.10 1。116 . 。435 。1267 1.0999 =

11、 。155 -0.55 0。321 -。0 0。08U =1。100 136-1355i 04230。1554i 0740.3215 1。1267。0902i 09990。181idU .0007 * .032 01 。80 0.086 0.367 0.03 .021 -00654 0。0030 。065PQ .0e-07 *00 -。00260.0163 00 -0.14 0090 0。49 0。0022 -0.07 。0043平衡节点功率S 4。7307 144i节点间线路上的流动功率矩阵 =输出电力线路功率总损耗at 个人心得经过两周的课程设计让我学到了很多以前没有学到的东西，首先让我明

12、白了电力系统的稳定运行，必须经过精密的设计和运算。在进行此次课程设计的过程中，加深了我对潮流计算的认识，尤其是对牛顿拉夫逊潮流计算的求解思路有了比较透彻的理解。同时整个求解过程中用到求节点导纳矩阵,求矩阵的逆阵等等，又让我们遇到了很多问题，我们经过查资料问同学、问老师,最后才把问题解决。使我对以前学过的知识有了一次很好的复习，同时也看到了我知识的欠缺，让我明白我应该学的东西还有很多，更是给我更多的激励让我以后会更加努力地学习.此次课程设计有让我学到了一本新的知识Matlb的使用，起初我们经过好大力气才把程序设计出来，可是设计的题目里的数据又不收敛,我们反复演算了好多次，最后在老师的指导下改了几

13、个数据程序终于收敛了。而且在此次课程设计中，我发现了自己的基础只是有很多的不足。这些基础的缺乏给我的设计造成了不小的障碍，今后我会努力学习来弥补这些不足。在这个过程中，我明白了，只要用心去做，认真去做，同学一起努力、团结，就会克服所有困难,取得最后的胜利。看到同学和自己设计的程序稳定运行,心中也有一种成功的喜悦。通过此次的实习让我发现手工的潮流计算和计算机潮流计算的优缺,对于小系统来说手工计算很简单，手工也很经济。可对于大系统来说手工计算是不可完成的任务,因此用计算机计算潮流是必不可少的，只要编好程序计算机就会给你一个自己想要的满意答案。而且方便、简单、准确.此次课程设计我能成功设计完成离不开

14、同学们的帮助、团队的合作以及老师细心的指导，在这里再次感谢老师的指导和同学的帮助。附录:y11=.03；y1=1/（1。05003j）；=0；y=0；y15=0；16=;y21=y1；y20.5j+1（6+0。25j)+1（.04+0.j)+1(（1。5)。0j）；y23=（0。06+0。0);y24=0；25=（0。04+0.2j)；y260；1=y3；y32y23；30。5j+1/（。06+.05j)+1/(0.8+。3)；4=1(0。03j）；y35=； y3=;y41y4；y4=y4；y4y4;4=0。j+1/（0.0+.3）+1(0.1.5j）；y45=1/(0.1+035)；y=;

15、y51y15；y5=y25；y533；y54=5；y5=0。5j+1（0。10。5）+/(0。+。2j）+1/（1.052）*0。15j)；y56=1(.00.015j)；y61；y2=0;y3=；4=；y6=y；y66=1/0。05j；%输入导纳参数YB=y 1 y1 14 y15 16; 21 y22 y3 y4 y2 y6; 32 y3y4 y y3； y41 y42 y4 4y45 y46；y 2 y3y4 y55y5;1 y62y63y4 y6 y6%形成导纳矩阵%*牛顿拉夫逊法潮流计算* E（2）=。0；（3）=1。0；E（4)=。0；E（）10;E(6）.1；()0;F（)=0;

16、F(）=0；F（）=0； F（6）=0； 设定电压初始值G=l（YB);B=imag(B)；%提取导纳实部,虚部分别为G，B()-1j；S(3）=-10.j；(）=1。6-0。8;S（5）=1.3;S（6)5；设定注入功率初始值Pe(）;=m(S）; %提取注入功率实部，虚部分别为，k=0;reiion=0。； 设定参数k为迭代次数初值为及精度大于0.0001N1=；%设定除平衡节点外节点数之和whlepiion 0。00001 E(1）=1。；()=0；输入平衡节点电压实部，虚部， r m:N1 for n1:1+1 t（)=(（）*（m,n)E（n)（,）（n)+F（m）（G(m，n)F(

17、）+（m,)E（）; Q（n)=（m）*（G（m，)E（n）B（，n）*（n）（)(G（,n)*F（n）+B（，n）*E（n)）；以设定值计算的注入功率 n dP(m）=P（m)-u（t）; (）=Q（m)-sum（Q）; 计算功率误差nf m2：N+ r n=1:N1+1 i（n）G（，n)(n）（m，)*F（n)； Bi（n)（m，)*F（）B（，n）*E（n）； nd N(m,m)=sum（Ai)+（m，m）*（m）+B（m，m）(）;计算d/de H（,）=su（B)-B（m，m）E（)+G（m，m）*F（m）;计算P/ (m，m) -sum(Bi）B（m,m)E（)G（m，m)F（）

18、;计算dQ/de J(，m)= sum(A)-G(m，m）*()B(m,m）F(m）;计算Q/df endor m=2：N+1 JJ（2m3，*m）=N(,）; J(2*m3,2）=H（m，m）; JJ(2*m2，2*3)=（m，m）； J（*-,2*m）J（m,m); end 给雅可比矩阵的对角元素赋值for =：1 or n=2：N1+ =n ele N(,)(，n）*E(m)+B（，）F（m）；计算Pd (,n）=B（m,）*E（m）+G(,n)（m）;计算dPf L(，)B(，n）*E()+G(m,）F（）；%计算Q J(m,n）=-B(，n)（m)(m，n）*E（）;计算dQ/df*

19、给雅可比矩阵的非对角元素赋值* JJ（2m,2*n3)=N（m，）； （2m-3，2*n-)=(，n）； J（2m-2，2n-）L（m，n); JJ（2*-2,*n2）J(m,n）; nd en edf：8J(10，i)=0；修正雅可比矩阵第十行,到元素为dJJ(1，9)2(); %修正雅可比矩阵第十行9元素为2e(0，10)=2F（）;修正雅可比矩阵第十行0元素为2 or =1:1 PQ（m1）=d(m+1); PQ（2）=dQ（m+1);d写Q矩阵，奇数行为P，偶数行为QPQ（）=1(E(6)2+F（6)2);%修正Q矩阵第十个元素为u2 dUin（J)PQ；求dU为的逆矩阵乘以PQ pe

20、cin=x(abs（Q）； %取精度为PQ的绝对值的最大值 orn1:N1 E（n+1）=E(n+1)d（2n）； F（n1）=F（n+1）+d（2n）；修正各节点电压 end fr n=：+1 （)（)（)）*j；%把电压写成e+jf的形式 ed=k+1;迭代次数加 ti=k，J, E ，F， U , dU, 输出k，JJ， ，F，,，sif k9；bekndd*功率及功率损耗的计算*for=1:N11 （m）(，m)（）;计算平衡节点电流ndS=U（1)*u（nj（)）; 计算平衡节点功率y0=。5/（1。03j）;y20=05(（1.052）0。03j）+。5j；30=0。5j；y40=0。;y5=。5-0。0/(（102)0。015);6=0。05/(1。5*0.015j)；=10 y20 y30 y40 y0 y6； %输入各节点对地导纳矩阵r :6fon=: if n当m=n时 S（，）0;流动节点功率为0 es%否则 S(m,）=U(m）*（onj（）*（)）（oj(）)- cn（U(n）)on（YB（，)）； 节点间线路上的流动功率矩阵 ndendnSSaaS=S+S输出电力线路功率总损耗datS及平衡节点功率s及节点间线路上的流动功率矩阵s