1、.电力系统分析潮流计算程序设计报告题目：13 节点配电网潮流计算学 院 电气工程学院 专业班级 学生姓名 学 号 班内序号 指导教师 房大中 提交日期 2015 年 05 月 04 日 .目 录一、程序设计目的 2二、程序设计要求 4三、13 节点配网潮流计算 .43.1 主要流程 .43.1.1 第一步的前推公式如下（1-1）-（1-5）： .43.1.2 第二步的回代公式如下（1-6） （1-9）： .43.2 配网前推后代潮流计算的原理 .53.3 配网前推后代潮流计算迭代过程 .53.3 计算原理 .6四、计算框图流程 7五、确定前推回代支路次序 7六、前推回代计算输入文件 8主程序：
2、 8输入文件清单： 9计算结果： 10数据分析： 10七、配电网潮流计算的要点 11八、自我总结 11九、参考文献 12附录一 MATLAB 的简介 .12.一、程序设计目的开式网络潮流计算：配电网的结构特点呈辐射状，在正常运行时是开环的；配电网的潮流计算采用的方法是前推回代法，本程序利用前推回代法的基本原理、收敛性。(1)在电网规划阶段，通过潮流计算，合理规划电源容量及接入点，合理规划网架，选择无功补偿方案，满足规划水平年的大、小方式下潮流交换控制、调峰、调相、调压的要求。(2)在编制年运行方式时，在预计负荷增长及新设备投运基础上，选择典型方式进行潮流计算，发现电网中薄弱环节，供调度员日常调
3、度控制参考，并对规划、基建部门提出改进网架结构，加快基建进度的建议。(3)正常检修及特殊运行方式下的潮流计算，用于日运行方式的编制，指导发电厂开机方式，有功、无功调整方案及负荷调整方案，满足线路、变压器热稳定要求。及电压质量要求。 (4)预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。 1234567891012图 1 13 节点配电网结构图.表 1 系统支路参数支路 R（ ） X（） B/2(TK) （S ）12 3.367 3.685 0.023 2.356 2.541 0.034 1.145 1.28 0.045 4.524 5.04 0.026 0.856 1.
4、14 0.067 2.745 2.965 0.028 3.743 4.251 0.089 2.237 2.756 0.0310 4.14 4.696 0.0311 1.328 1.763 0.01112 2.436 2.866 0.0413 3.521 3.966 0.0表 2 系统负荷参数节点编号 节点类型 节点初始电（kV）Pi（MVA）Qi（MVA）1 根节点 10.4 0 02 中间节点 10.0 0.0342 0.03013 中间节点 10.0 0.0693 0.06424 中间节点 10.0 0.0845 0.07635 叶节点 10.0 0.0295 0.02616 中间节点 1
5、0.0 0.0474 0.04097 叶节点 10.0 0.1176 0.09578 中间节点 10.0 0.0946 0.08579 叶节点 10.0 0.0916 0.085910 叶节点 10.0 0.0271 0.0229.11 中间节点 10.0 0.0696 0.064312 叶节点 10.0 0.0676 0.057913 叶节点 10.0 0.0298 0.0242二、程序设计要求1. 看懂前推回代法计算程序；2. 报告叙述计算原理及计算流程；3. 绘制计算流程框图；4. 确定前推回代支路次序（广度优先，或深度优先）并编写前推回代计算输入文件，然后进行潮流计算；5. 整理专利计
6、算结果；6. 总结配电网潮流计算的要点；三、设计内容1.根据电力系统网络推导电力网络数学模型，写出节点导纳矩阵； 2.赋予各节点电压变量（直角坐标系形式）初值后，求解不平衡量；3.形成雅可比矩阵； 4.求解修正量后，重新修改初值，从2开始重新循环计算； 5.求解的电压变量达到所要求的精度时，再计算各支路功率分布、功率损耗和节点 6.上机编程调试；连调； 7.计算分析给定系统潮流分析并与手工计算结果作比较分析。 8.准备计算机演示答辩，书写该课程设计说明书（必须计算机打印）四、13 节点配网潮流计算4.1 牛拉法的原理及其基本方程牛顿迭代法（Newtons method）又称为牛顿-拉夫逊（拉弗
7、森）方法（Newton-Raphson method），它是牛顿在 17 世纪提出的一种在实数域和复数域上近似求解方程的方法。设 r 是 f(x) = 0 的根，选取 x0 作为 r 初始近似值，过点（x0,f(x0)）做曲线 y = f(x)的切线 L，L 的方程为 y = f(x0) f(x0)(x-x0)，求出 L 与 x 轴交点的横坐标 x1 = x0-f(x0)/f(x0)，称 x1 为 r 的一次近似值。过点（ x1,f(x1)）做曲线 y = f(x)的.切线，并求该切线与 x 轴的横坐标 x2 = x1-f(x1)/f(x1)，称 x2 为 r 的二次近似值。重复以上过程，得
8、r 的近似值序列，其中 x(n+1)=x(n)f(x(n)/f(x(n)，称为 r的 n+1 次近似值，上式称为牛顿迭代公式。解非线性方程 f(x)=0 的牛顿法是把非线性方程线性化的一种近似方法。把 f(x)在 x0 点附近展开成泰勒级数 f(x) = f(x0)+(xx0)f(x0)+(xx0)2*f(x0)/2! + 取其线性部分，作为非线性方程 f(x) = 0 的近似方程，即泰勒展开的前两项，则有 f(x0)+f(x0)(xx0)=f(x)=0 设 f(x0)0 则其解为 x1=x0f(x0)/f(x0) 这样，得到牛顿法的一个迭代序列：x(n+1)=x(n)f(x(n)/f(x(n
9、)。4.2 PQ 分解法的原理及其基本方程PQ 分解法的基本思想是根据电力系统实际运行特点：通常网络上的电抗远大于电阻值，则系统母线电压幅值的微小变化 对母线有功功率的改变 影响很小。同样，UP母线电压相角的少许改变 ，也不会引起母线有功功率的明显改变 ，因此，节 Q点功率方程在用极坐标形式表示时，它的修正方程式可简化为（2-0/PHQLU1）这就是把 2（n-1）阶的线性方程组变成了两个 n-1 阶的线性方程组，将 P 和 Q 分开来进行迭代计算，因而大大地减少了计算工作量。但是 H、L 在迭代过程中仍然在不断的变化而且又都是不对称矩阵。对牛顿法的进一步简化，即把式（2-1）中的系数矩阵简化
10、为在迭代过程中不变的对称矩阵。在一般情况下，线路两端的电压相角 是不大的（不超过 1020），因此，ij可以认为 （2-cos1inijjjijGB2）此外，与系统各节点无功功率相应的导纳 远小于该节点自导纳的虚部，即LDi.2iLDi iiQBBU因而 （2-ii3）考虑到以上关系，式（2-1）的系数矩阵中的个元素可以表示为( i , j=1,2, n-1)ijijiHUB( i , j=1,2, m)ijijiL而系数矩阵 H 和 L 则可以分别写成：1121,12 2,1,11,21,1.nnnnnUBBBUU ,2122,12 21,21, 1. .nnn nBBUU(2-4)11DU
11、B2121 212.mmmmUBULUB 11212 2212. .mmmmBBUU(2-5)22“DB.将式（2-4）和（2-5）代人式 中，得到(cosin)ijijiji iiYGBUe11DPUB2“Q用 和 分别左乘以上两式，便得1D12D(2-6)11UPBU(2-7)2“DQ这就是简化了的修正方程式，它们也可以展开写成(2-8) 1211,2,1,11,2221,121 nnnnnn UBBUPP (2-9)mmmm UBBUQ 21,2,1, ,221,121在这两个修正方程式中系数矩阵元素就是系统导纳矩阵的虚部，因而系数矩阵是对称矩阵，且在迭代过程中保持不变。这就大大减少了计
12、算工作量。用极坐标表示的节点功率增量为(2-10)11(cosin)0iniisijiijjjniisijiijjjPUGBQ式（2-8）、（2-9）和（2-10）构成了 PQ 分解法迭代过程的基本方程式。.4.3 配网前推后代潮流计算的原理前推回代法在配电网潮流计算中简单实用，所有的数据都是以矢量形式存储，因此节省了大量的计算机内存，对于任何种类的配电网只要有合理的 R/X 值，此方法均可保证收敛。算法的稳定性也是评价配电网潮流算法的重要指标。一般情况下，算法的收敛阶数越高，算法的稳定性越差，前推回代法的收敛阶数为一阶，因此它也具有较好的稳定性。比较而言，前推回代法充分利用了网络呈辐射状的结
13、构特点，数据处理简单，计算效率高，具有较好的收敛性，被公认是求解辐射状配电网潮流问题的最佳算法之一。4.4 配网前推后代潮流计算迭代过程配网潮流前推回代潮流算法第 n+1 步的迭代公式如下：节点 i 的前推计算公式为由图 1 可知节点 i+1 的注入有功功率和无功功率分别为：NB NB-1P(i+1)= PL(j)+ LP(j)Q(i+1)= QL(j)+ LQ(j)其中：i=1 ， 2， 。 。 。NB-1，NB 为节点数，PL(i)为第 i 节点所带负荷有功功率；QL(i)为第 i 节点所带负荷无功功率；LP(i) 为第 i 条线段上的有功功率损耗；LQ(i)为第i 条线段上的无功功率损耗
14、。收敛条件以前后两次迭代的电压偏差作为迭代收敛条件，(8)式表明节点电压幅值最大偏差小于设定阈值，即认为迭代收敛，则迭代过程结束。max|Vk(i)-Vk-1(i)|0.0001PQ2=PQ; %每一次迭代各节点的注入有功和无功相同PL=0.0;for i=1:NBkf=FT(i,1); %前推始节点号kt=FT(i,2); %前推终节点号x=(PQ2(kf,2)2+PQ2(kf,3)2)/(V(kf)2); %计算沿线电流平方 APQ1(i,1)= RX(i,1)*x+PQ2(kf,2); %计算支路首端有功/MW RX(i,1)RPQ1(i,2)= RX(i,2)*x+PQ2(kf,3);
15、 %计算沿支路的无功损耗/Mvar RX(i,2)XPQ2(kt,2)= PQ2(kt,2)+PQ1(i,1); %用 PQ1 去修正支路末端节点的有功 P 单位 MWPQ2(kt,3)= PQ2(kt,3)+PQ1(i,2); %用 PQ1 去修正支路末端节点的有功 Q 单位 MvarPL=PL+RX(i,1)*x; %累积计算整个网络的功率损耗endangle(1)=0.0;for i=NB:-1:1kf=FT(i,2); %回代始节点号 kt=FT(i,1); %回代终节点号dv1=(PQ1(i,1)*RX(i,1)+PQ1(i,2)*RX(i,2)/V(kf); %计算支路电压损耗的纵
16、分量 dv1dv2=(PQ1(i,1)*RX(i,2)-PQ1(i,2)*RX(i,1)/V(kf); %计算支路电压损耗的横分量 dv2V2(kt)=sqrt(V(kf)-dv1)2+(dv2)2); %计算支路末端电压/kVangle(kt)= atand(dv2/(V(kf)-dv1)+angle(kf); %计算支路电压相角（单位度）endmaxd=abs(V2(2)-V(2);V2(1)=V(1);for i=3:1:NNif abs(V2(i)-V(i)maxd;maxd=abs(V2(i)-V(i); end endmaxdk=k+1PQ1 %潮流分布 即支路首端潮流 MVAV=
17、V2 %节点电压模计算结果 kV.angle %节点电压角度计算结果单位度PL %网损单位 MWendclear输入文件清单：function PQ,FT,RX=case115()PQ=%节点电压 有功 无功10.4 0 010.0 0.0342 0.030110.0 0.0693 0.064210.0 0.0845 0.076310.0 0.0295 0.026110.0 0.0474 0.040910.0 0.1176 0.095710.0 0.0946 0.085710.0 0.0916 0.085910.0 0.0271 0.022910.0 0.0696 0.064310.0 0.0
18、676 0.057910.0 0.0298 0.0242;FT=%首端 末端 5 413 412 114 310 3 11 37 69 83 26 28 22 1;RX=% R X 4.524 5.043.521 3.966.2.436 2.8661.145 1.284.14 4.6961.328 1.7632.745 2.9652.237 2.7562.356 2.5410.856 1.143.743 4.2513.367 3.685;计算结果：k =8PQ1 = 0.0296 0.02620.0299 0.02430.0678 0.05810.1444 0.12730.0272 0.023
19、00.1379 0.12310.1183 0.09640.0920 0.08640.3852 0.34450.1661 0.13790.1891 0.17500.8116 0.7278V =10.4000 9.8795 9.6991 9.6652 9.6377 9.8492 9.7872 9.7326 9.6870 9.6764 9.6578 9.6235 9.6444angle = 0 0.3011 0.4011 0.4250 0.4438 0.3431 0.3942 0.3899 0.4266 0.4209 0.4498 0.4823 0.4452PL =0.0488数据分析：经过 8 次
20、迭代，各节点的电压误差均在 0.0001kv 以内，计算到此结束。.七、配电网潮流计算的要点1. 需要注意系统中所有母线的电压是否在允许的范围以内。2. 系统中各种元件(线路、变压器等)是否会出现过负荷。3. 在进行电压和功率分布计算以前，先要对网络等值电路作些简化处理。4. 确立前推回代的计算次序。5. 在开始编写程序之前，要了解 matlab 仿真软件，包括其程序的编写，程序的运行，我们都要清楚的知道6. 为了提高计算精度，重复计算，给出一个容许误差作为计算的收敛的判据。八、自我总结此次课程设计首先让我明白了要使电力系统运行的稳定，必须经过精密的设计和计算。在进行课题设计的过程中，加深了我
21、对潮流计算的认识，尤其是对牛顿拉夫逊潮流计算的求解思路有了比较透彻的理解。在此次程序设计中，我在 MATLAB 编程，13 节点配电网潮流计算，WORD 文档的编辑方面均有所提高，之前学过的潮流计算仅限于书本，通过完成这个作业使我对这块知识了解得更加全面，但也暴露出了一些问题：理论知识还是不足，基本知识不充足，做起这个大作业难度还是很大的，对 MATLAB 的性能和特点还不能有一个全面的把握，相信通过以后的学习能弥补这些不足，从而达到一个新的层次。在过程设计中我发现了一些要点，要求得支路的首端功率，从与叶节点联接的支路开始，该支路的末端功率即等于叶节点功率，利用这个功率和对应的节点电压计算支路
22、功率损耗，才能求首段功率。还有在迭代计算开始之前，先要确定支路的计算顺序问题，即确定好是深度优先还是广度优先。要求得各支路终结点的电压的话利用上述所得的支路首端功率和本步骤刚算出的本支路始节点的电压，从电源点开始逐条支路进行计算。而且在此次课程设计中，我发现了自己的基础知识有很多的不足。这些基础的缺乏给我的设计计划造成了不小的障碍。在这个过程中，我明白了，只要用心去做，认直去做，持之以恒，就会有新的发现，有意外的收获。虽然很多错误被克服了，但是系统中难免还有很多不足之处，有些地方感觉做的.不是很正确，希望在以后的学习中能够在同学们的帮助和老师的指导下做得更好。九、参考文献1何仰赞 温增银. 电
23、力系统分析. 华中科技大学出版社，20022蔡旭晖、刘卫国、蔡立燕. MATLAB 基础与应用教程D. 北京：人民邮电出版社 2009 附录一 MATLAB 的简介MATLAB 是由美国 mathworks 公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如 C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MA
24、TLAB 和 Mathematica、Maple 并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB 可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、.信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 MATLAB 的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用 MATLAB 来解算问题要比用 C，FORTRAN 等语言完成相同的事情简捷得多，并且 mathwork 也吸收了像 Maple 等软件的优点,使 MATLAB 成为一个强大的数学软件。在新的版本中
25、也加入了对 C，FORTRAN ，C+ ，JAVA 的支持。可以直接调用,用户也可以将自己编写的实用程序导入到 MATLAB 函数库中方便自己以后调用，此外许多的 MATLAB 爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。MATLAB 产品族可以用来进行以下各种工作： 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真 数字图像处理 技术 数字信号处理 技术 通讯系统设计与仿真 财务与金融工程 MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用 MATLAB 函数集）扩展了 MATLAB 环境，以解决这些应用领域内特定类型问题。
