《电力系统潮流上机》课程设计报告院系：电气工程学院班级：电088班________学号：01学生姓名：刘东昇________指导教师：张新松设计周数：两周_________日期：2010年12月25日一、课程设计的目的与要求目的：培养学生的电力系统潮流计算机编程能力，掌握讣算机潮流讣算的相关知识要求：基本要求：1.编写潮流计算程序:2.在计算机上调试通过：3.运行程序并il•算岀正确结果：4.写出课程设计报告二、设计步骤：1.根据给泄的参数或工程具体要求（如图），收集和査阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。

勺节点1为平衡节点，节点2、3、4、5为PQ节点匚2.在给定的电力网络上画出等值电路图。

3.运用汁算机进行潮流计算。

4.编写设汁说明书。

三、设计原理1.牛顿•拉夫逊原理牛顿迭代法是取xO之后，在这个基础上，找到比xO更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。

牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f（x） = 0的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。

电力系统潮流il•算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡肖点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。

为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知肖点赋电压初值，一般为额泄电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。

牛顿一拉夫逊迭代法的一般步骤：（1）形成各节点导纳矩阵Y。

（2）设个节点电压的初始值U和相角初始值e还有迭代次数初值为0。

（3） II-算各个节点的功率不平衡量。

（4〉根据收敛条件判断是否满足，若不满足则向下进行。

（5）计算雅可比矩阵中的各元素。

（6）修正方程式个节点电压（7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。

（8）汁算平衡节点输出功率和各线路功率2.网络节点的优化1）静态地按最少出线支路数编号这种方法由称为静态优化法。

在编号以前。

首先统汁电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。

当由n个节点的出线支路相同时，则可以按任意次序对这n个节点进行编号。

这种编号方法的根据是导纳矩阵中，出线支路数最少的肖点所对应的行中非零元素也2 ）动态地按增加岀线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统讣出来的，在编号过程中认为固泄不变的，事实上，在廿点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化（增加，减少或保持不变）。

因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个盯点进行编号，就可以预期得到更好的效果，动态按最少岀线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中％巧点出线支路数目的变动情况。

3.MATLAB编程应用Matlab是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。

由于使用Matlab编程运算与人进行科学汁算的思路和表达方式完全一致，所以不像学习髙级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一髙效的科研助手。

四、设计内容1.设计流程图2.程序clear;clc%重新编号，把原题中的节点1,2,3,4,5重新依次编号为54,2,3,4.其中1-4号为PQ节点，5号为平衡节点y 二0;%输入原始数据，求节点导纳矩阵y (i,2)=i/+； y(13)=1/+; y ⑴ 4)=i/+；y ⑴ 5)=1/+;y(2/3)=l/+；y(2/5)=l/+;y(3z4)=l/+;y(4,5)=0;for i=l:5for j=i:5y(j』=y(i,j);endendY=0;%求互导纳for i=l:5for j=l:5if r=jY(i,j)=-y(i,j)；endendend%求自导纳for i=l:5Y(i/i)=sum(y(i/:));endY%Y为导纳矩阵G=real(Y);B=imag(Y);%原始节点功率S(l)=+;S(2)=賦初值U 二ones(l,5);U(5)=;e=zeros(：U5);ox=on€s(8j);fx=ones(8j);count=0 %计算迭代次数while max(fx)>le-5for i=l:4for j=l:4H(i,j)=0;N(i』=0;M(i』=0;L(i,j)=0;oP(i)=0;oQ(i)=0；endendfor i=l:4for j=l:5oP(i)=oP(i)-U(i)#U(j)*(G(ij)*cos(e(i)-e(j))+B(iJ)#sin(e(i)-e(j))); oQ(i)=oQ(i)-U(i)*UO)*(G(i,j)*sin(e(i)-e(j))-B(i z j)*cos(e(i)-e(j)));endoP(i)=oP(i)+P(i); oQ(i)=oQ(i)+Q(i);endfx=[oP/oQ]';%求雅克比矩阵%当中时候求H,NML如下:for i=l:4for j=l:4 ifr=jH(i,j)=-U(i)*U0)*{G(i,j)*sin(e(i)-e(j))-B(i/j)*cos(e(i)-e(j)))；N(i/j)=-U(i)*U(j)*(G(ij)*cos{e(i)-e0))+B(iJ)*sin(e(i)-eO))); L(i,j)=H(i,j);M(iJ)=-N(iJ);endendendH,NML%当冃时HNM丄如下:for i=l:4for j=l:5if r=jH{i,i)=H(i z i)+U(i)*U0)*{G{i/j)*sin(e(i)-e(j))-B(i z j)*cos (e(i)-e(j)));N(i/i)=N(i,i)-U(i)*U0)*(G(i/j)*cos(e(i)-e(j))+B(i/j)*sin(e(i)-e(j))); M(i/i)=M(i,i)-U(i)*UO)*{G(i,j)*cos(e(i)-e(j))+B(i/j)*sin(e(i)-e(j))); L(i,i)=L(i/i)-U(i)*UO)*(G(i,j)*sin(e(i)-e(j))-B(i z j)*cos(e(i)-e(j))); endendN(i/i)=N(i z i)-2*(U(i))^2*G(i/i);L(i,i)=L(i/i)+2*(U(i))A2^B(i/i)；endJ=(H,N;M,L] %J为雅克比矩阵ox=-((inv(J))*fx);for i=l:4oe(i)=ox(i); oU(i)=ox(i+4)*U(i);endfor i=l:4e(i)=e(i)+oe(i); U(i)=U(i)+oU(i);endcount=co un t+1;endOX,U QCO unt%求节点注入的净功率i=5;for j=l:5P(i)=U(i)*U(j)*(G(i,j)*cos(e(i)-e(j))+B(i,j)*sin(e(i)-e(j)))+P(i); Q(i)=U(i)*UO)*(G(i,j)*sin(e(i)-e(j))-B(iJ)*cos(e(i)-e(j)))+Q(i); endS(5)=P(5)+Q(5)*sqrt(-l);S%求节点注入电流l=Y*U'3.运行结果Y值：(p New to MATLAB? Vifatch this Video, see Demos, or read Gettina Started.10. 8333 -32.50001 -1.6667 + 5-OOOOi -1. 6667 + S.OOOOi -2.5000 + 7.5000i -5.0000 +15.0000i-1. 6667 + 5.0000i 12.9167 -38.7500i -10. 0000 +30. 00001 0 ・ L.2500 + 3.75001 -1. 6667 + 5.00001 -10.0000 +30.0000i 12. 9167 -38.75001 -1.2500 + 3. 7500i 0-2. 5000 + 7.5000i 0 -1.2500 + 3. 7500i 3.7500 -11.2500i 0-5.0000 +15.00001 -1.2500 + 3.7500i 0 0 6.2500 -18. 7500i 迭代过程:count =0 5. 0000 5. OOOO7. 50005. 0000 0 30. 0000 05. 0000 30.0000 0 3. 75007. 5000 0 3. 7500 00 1. 6667 1. 6667 2・ 50001.6667 0 10. 0000 01. 6667 10.0000 0 1. 25002. 5000 O 1.2500 00 -1. 6667 -1. 6667 -2. 50001.6667 0-10. 000001. 6667 -10.0000 0 一 1・25002. 5000 0 -1. 2500 00 5. 00005. 0000 7. 50005. 000030. 0000 05.0000 30.00000 3.75007.50003. 7500 0-33.4000 5. 0000 5. 0000 7. 5000 -10. 5333 1. 66671. 66672. 50005. 0000 -38. 9750 30. 0000 0 1. 6667 -12. 8417 10. 0000 0 5・ 0000 30. 0000 -38. 7500 3・ 7500 1. 6667 10. 0000-12. 9167 1. 2500 7. 5000 03. 7500 -11. 2500 2. 50000 1. 2500 -3. 7500 11. 1333-1.6667-1.6667 -2.5000 -31.60005.00005.0000 7. 5000 -1.6667 12.9917 -10.0000 05.0000 -38.5250 30. 0000 0-1.6667 -10.0000 12.9167 -1.2500 5.0000 30.0000 -38.7500 3. 7500 -2.5000-1.25003. 75007.50003.7500-11. 2500M =-1・ 5572-L ・ 5235-2. 1511 一1・ 97010 一 1O ・ 1 13*7 0一1・ 9985 -IO ・ 48030 -1. 2203一3・ 2290 -L ・ 3384 00 5. 35995. 36228. 03965. 2222 0 30. 95205・ 2038 30. 8298 0 3. 85767. 7223 0 3. 8183 0-35. 0648 5. 3599 5. 3622 8. 0396-12. 0310 1. 5572 1. 5235 2. 15115. 2222 -40. 0794 30. 9520 0 1. 9701 -12. S627 10. 1137 0 5. 2038 30. 8298 -39. 8913 3. 85761. 9985 10. 4803 -12. 86821. 2203 7. 7223 03. 8183 -11.5406 3. 1029 01. 3384. -3. 1996 1 1. 5373-1. 5572 -1. 5235 -2. 1511-35. 6400 5. 3599 5. 36228. 0396-1. 9701 L3. 7720 -10. 1 1375. 2222 -39. 8247 dO ・ 9520 0-1. 9986 -L0. 4S03 13. 6991 -1. 2203 5. 203S 30. 8298 -39. S105 3. S576-3. 1029 0 -1. 33844. 44127. 72230 3. 8183 -11. 38185 35995・ 2222 5・ 2038 7・ 722336 32985・ 3622 30・ 95203・ & 1838.03963. 8576O1・ 9701 1・ 9985 3・ 10291. 5572 4803L ・ 5235 10・ 1 137 0L ・ 33842. 1511O1. 22035. 29045・ 2J912 7・ 924(3峪.1585 O30・ 5426石.1397 30・ 42673・ 80L8T ・ 62 17 O 3. 7C44 O0 1・ 54361・ 5 1 122・ 13681・ 9393 O9・ 93761・ Q6DS 1 O. 3300O1 ・ 2OE5O3・ 0463O1・ 3X710 一1・ 5436一1・ 5 112一2・1368 一 1・ 0303O —Q. QS76O一 1・ 9658 -10・ 3366一1・ 2080一3・ 0453O 一1・ 3 1715. 29045. 2912 7. 92465. 158530. 54265. 1397了0. 42673. 80187. 62173. 7644 0-1・ 5436 —1.・ 511 1 -2-136S 一 1・ 93910 -9・ 9867 0 一1. 9655 -10.3343 0 一1. 2049一3・ 04.4S 0 -L ・ 3169 0行.28995・ 2907 92375・ 1S80 030. 5395O5・ 1392 30・ 42360 3. 80137・ 62100 3・ 7640 0-34.7163 5. 2904 5. 2912 7. 9246 5. 1585 -39.5849 30. 5426 05. 1397 30. 4267 -39. 36813. 8018 7・ 62173・ 76<14-LI ・ 3S611 1・ 4384 -1.5436-1. 5112 -2. 1368 -1. 9393 13.5950 -9. 9876 0-1. 9658 -10.335513. 5062-1. 2050-3. 0453-1. 31714. 3623H =5・ 2899 5・ 29077. 92375・ 1580 036 53955・ 1392 30・ 42363.80137・ 62103. 7640N =0 1・ 5436 L ・ 5111 2.13681・ 9391 0 9・ 9867 0 1. 9555 1O.3343 0 1-2049 3・ oaa8 0 L ・ 3169 0-11.8394 1. 5436 1. 51122. 13681.9393 -12. 6950 9・ 9876 0 1.9658 10.3355 -12. 70591. 20503・ 0<1531. 3171 7・ 1619 -36.1173 5. 29045. 2912 7.9246 氐 1585 -39. 2853 30.5426 0 5. 1397 30. 4267-39. 2683 3. 80187.62173. 7644 -1 1. 1866-34.7136 5. 2899 5.2907 7. 9237 -H.837S 1.5436 1. 5111 2.13685.1580 -39.5812 30.5395 0 1. 9391 -12.6937 9. 9867 05.1392 30. 4236 -39.3642 3.8013 1. 9655 10.3343 -12. 7047 1.20497.6210 0 3.7640 -11.3850 3. 0448 0 1.3169 -3.161711.4379 -1. 5436 -1.5111 -2.136S -35. 1136 5.2899 5・ 2907 7.9237-1.9391 13. 5937 -9.9867 0 5. 1580 ・ 39.2812 30. 5395 0-1.9655 -10. 3343 13.5047 -1.2049 5. 1392 30.4236-39. 2642 3.8013-3.04480 -1.3169 4.3617 7. 62100 3. 7640 -11.1850电压值:ox =1.0 e-OO8 *-O.0004O.08350.1153O-2L57-O.2645-O.3080-O-3215-0.39281.0365 1・ 0088 1・ 0073L・ 00 16 1・ 060。