《2022年电力系统分析潮流计算课程序设计及其MATLAB程序设计 .pdf》由会员分享，可在线阅读，更多相关《2022年电力系统分析潮流计算课程序设计及其MATLAB程序设计 .pdf（17页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、电力系统分析潮流计算程序设计报告题目： 13 节点配电网潮流计算学院电气工程学院专业班级学生姓名学号班内序号指导教师房大中提交日期 2015年 05 月 04 日名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 17 页 - - - - - - - - - 1 目 录一、程序设计目的. 2二、程序设计要求. 4三、 13 节点配网潮流计算. 43.1 主要流程 . 错误！未定义书签。3.1.1第一步的前推公式如下（1-1）- （1-5 ）： . 错误！未定义书签。3.1.2
2、第二步的回代公式如下（1-6）（ 1-9）： . 错误！未定义书签。3.2 配网前推后代潮流计算的原理 . 73.3 配网前推后代潮流计算迭代过程 . 83.3 计算原理 . 9四、计算框图流程. 10五、确定前推回代支路次序. 错误！未定义书签。六、前推回代计算输入文件. 11主程序： . 11输入文件清单：. 12计算结果： . 13数据分析： . 13七、配电网潮流计算的要点. 14八、自我总结. 14九、参考文献. 15附录一 MATLAB 的简介 . 15名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 -
3、 - - - - - - 第 2 页，共 17 页 - - - - - - - - - 华南理工大学广州学院电气工程学院课程设计报告2 一、程序设计目的开式网络潮流计算：配电网的结构特点呈辐射状，在正常运行时是开环的；配电网的潮流计算采用的方法是前推回代法，本程序利用前推回代法的基本原理、收敛性。(1) 在电网规划阶段，通过潮流计算，合理规划电源容量及接入点，合理规划网架，选择无功补偿方案，满足规划水平年的大、小方式下潮流交换控制、调峰、调相、调压的要求。 (2)在编制年运行方式时，在预计负荷增长及新设备投运基础上，选择典型方式进行潮流计算，发现电网中薄弱环节，供调度员日常调度控制参考，并对规
4、划、基建部门提出改进网架结构，加快基建进度的建议。 (3)正常检修及特殊运行方式下的潮流计算，用于日运行方式的编制，指导发电厂开机方式，有功、无功调整方案及负荷调整方案，满足线路、变压器热稳定要求。及电压质量要求。(4) 预想事故、设备退出运行对静态安全的影响分析及作出预想的运行方式调整方案。12345678910111213图 1 13 节点配电网结构图表 1 系统支路参数名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 17 页 - - - - - - - - - 华南
5、理工大学广州学院电气工程学院课程设计报告3 支路R（）X（ ）B/2(TK) （S）12 3.367 3.685 0.0 23 2.356 2.541 0.0 34 1.145 1.28 0.0 45 4.524 5.04 0.0 26 0.856 1.14 0.0 67 2.745 2.965 0.0 28 3.743 4.251 0.0 89 2.237 2.756 0.0 310 4.14 4.696 0.0 311 1.328 1.763 0.0 1112 2.436 2.866 0.0 413 3.521 3.966 0.0 表 2 系统负荷参数节点编号节点类型节点初始电（ kV）
6、Pi（MVA）Qi（MVA）1 根节点10.4 0 0 2 中间节点10.0 0.0342 0.0301 3 中间节点10.0 0.0693 0.0642 4 中间节点10.0 0.0845 0.0763 5 叶节点10.0 0.0295 0.0261 6 中间节点10.0 0.0474 0.0409 7 叶节点10.0 0.1176 0.0957 8 中间节点10.0 0.0946 0.0857 9 叶节点10.0 0.0916 0.0859 10 叶节点10.0 0.0271 0.0229 11 中间节点10.0 0.0696 0.0643 12 叶节点10.0 0.0676 0.0579
7、 13 叶节点10.0 0.0298 0.0242 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 17 页 - - - - - - - - - 华南理工大学广州学院电气工程学院课程设计报告4 二、程序设计要求1.看懂前推回代法计算程序；2.报告叙述计算原理及计算流程；3.绘制计算流程框图；4.确定前推回代支路次序（广度优先，或深度优先）并编写前推回代计算输入文件，然后进行潮流计算；5.整理专利计算结果；6.总结配电网潮流计算的要点；三、设计内容 1.根据电力系统网络推导
8、电力网络数学模型，写出节点导纳矩阵； 2.赋予各节点电压变量（直角坐标系形式）初值后，求解不平衡量； 3.形成雅可比矩阵； 4.求解修正量后，重新修改初值，从2开始重新循环计算； 5.求解的电压变量达到所要求的精度时，再计算各支路功率分布、功率损耗和节点 6.上机编程调试；连调； 7.计算分析给定系统潮流分析并与手工计算结果作比较分析。 8.准备计算机演示答辩，书写该课程设计说明书（必须计算机打印）四、13 节点配网潮流计算4.1 牛拉法的原理及其基本方程牛顿迭代法（Newtons method） 又称为牛顿 - 拉夫逊 （拉弗森）方法 （Newton-Raphson method），它是牛顿
9、在 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 的二次近似值。重复以上过程，得 r 的近似值序列，其中x(n+1)=x(n) f(x(n)/
10、f(x(n)，称为 r 的 n+1次近似值，上式称为牛顿迭代公式。解非线性方程f(x)=0的牛顿法是把非线性方程线性化的一种近似方法。 把 f(x) 在 x0 点附 近展开成泰勒级数 f(x) = f(x0)+(x名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 17 页 - - - - - - - - - 华南理工大学广州学院电气工程学院课程设计报告5 x0)f(x0)+(xx0)2*f(x0)/2! + 取其线性部分，作为非线性方程f(x) = 0的近似方程，即泰勒展开
11、的前两项，则有f(x0)+f(x0)(xx0)=f(x)=0 设 f(x0)0 则其解 为 x1=x0 f(x0)/f(x0) 这 样， 得到 牛顿 法 的 一个 迭 代 序列： x(n+1)=x(n)f(x(n)/f(x(n)。4.2 PQ 分解法的原理及其基本方程PQ分解法的基本思想是根据电力系统实际运行特点：通常网络上的电抗远大于电阻值，则系统母线电压幅值的微小变化U对母线有功功率的改变P影响很小。同样，母线电压相角的少许改变，也不会引起母线有功功率的明显改变Q，因此，节点功率方程在用极坐标形式表示时，它的修正方程式可简化为00/PHQLUU（2-1）这就是把 2（n-1）阶的线性方程组
12、变成了两个n-1 阶的线性方程组， 将 P和 Q分开来进行迭代计算，因而大大地减少了计算工作量。但是H、L 在迭代过程中仍然在不断的变化而且又都是不对称矩阵。对牛顿法的进一步简化，即把式（2-1）中的系数矩阵简化为在迭代过程中不变的对称矩阵。在一般情况下，线路两端的电压相角ij是不大的（不超过1020），因此，可以认为cos1sinijijijijGB（2-2）此外，与系统各节点无功功率相应的导纳LDiB远小于该节点自导纳的虚部，即2iLDiiiiQBBU因而2iiiiQU B（2-3）考虑到以上关系，式（ 2-1）的系数矩阵中的个元素可以表示为ijijijHU U B ( i , j=1,2
13、, n-1) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 17 页 - - - - - - - - - 华南理工大学广州学院电气工程学院课程设计报告6 ijijijLU U B ( i , j=1,2, m) 而系数矩阵 H和 L 则可以分别写成：1111112211,112211222222,1111,1111,2211,11.nnnnnnnnnnnnU B UU B UU BUU B UU B UU BUHUBUUBUUBU11121, 11121222,1221
14、,11,21, 111.nnnnnnnnBBBUUBBBUUBBBUU11DDUBU (2-4) 111111221122112222221122.mmmmmmmmmmmmU B UU B UU B UU B UU B UU BULUB UU BUUBU1111211221222212.mmmmmmmmUBBBUUBBBUUBBBU22DDUBU (2-5) 将式（ 2-4）和（ 2-5）代人式(cossin )iijijijjiiiiiYGjBUU eUj中，得到11DDPUBU2DQUBU用11DU和12DU分别左乘以上两式，便得111DDUPBU (2-6) 名师资料总结 - - -精品
15、资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 17 页 - - - - - - - - - 华南理工大学广州学院电气工程学院课程设计报告7 12DUQBU (2-7) 这就是简化了的修正方程式，它们也可以展开写成1122111, 12, 11 , 11,222211,11211112211nnnnnnnnnnUUUBBBBBBBBBUPUPUP (2-8) mmmmmmmmmUUUBBBBBBBBBUQUQUQ21,2,1 ,22221,112112211 (2-9) 在这两个修正方程式中系数矩阵
16、元素就是系统导纳矩阵的虚部，因而系数矩阵是对称矩阵，且在迭代过程中保持不变。这就大大减少了计算工作量。用极坐标表示的节点功率增量为11(cossin)0(cossin)0niisijijijijjniisijijijijjPPUUGBQQUUGB (2-10) 式（2-8）、（ 2-9）和（ 2-10）构成了 PQ分解法迭代过程的基本方程式。4.3 配网前推后代潮流计算的原理前推回代法在配电网潮流计算中简单实用，所有的数据都是以矢量形式存储，因此节省了大量的计算机内存，对于任何种类的配电网只要有合理的 R/X 值，此方法均可名师资料总结 - - -精品资料欢迎下载 - - - - - - -
17、- - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 17 页 - - - - - - - - - 华南理工大学广州学院电气工程学院课程设计报告8 保证收敛。算法的稳定性也是评价配电网潮流算法的重要指标。一般情况下，算法的收敛阶数越高，算法的稳定性越差，前推回代法的收敛阶数为一阶，因此它也具有较好的稳定性。比较而言，前推回代法充分利用了网络呈辐射状的结构特点，数据处理简单，计算效率高，具有较好的收敛性， 被公认是求解辐射状配电网潮流问题的最佳算法之一。4.4 配网前推后代潮流计算迭代过程配网潮流前推回代潮流算法第n+1 步的迭代公式如下：节点i
18、的前推计算公式为由图 1 可知节点 i+1 的注入有功功率和无功功率分别为：NB NB-1 P(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 条线段上的无功功率损耗。收敛条件以前后两次迭代的电压偏差作为迭代收敛条件，(8)式表明节点电压幅值最大偏差小于设定阈值，即认为迭代收敛，则迭代过程结束。max|Vk(i)-Vk- 1(i)|.(8)在开始计算的时候先把
19、各节点上的LP(i) 、LQ(i) 分别设为零，利用式(1)计算出P(i+1)、Q(i+1)，分别代入式 (6)计算各节点电压幅值，作为初始化的节点电压值.然后把此电压值代入公式 (2)来修正线路的损耗LP(i)、LQ(i)，用经过修正的线路损耗求新的节点注入功率，再根据节点注入功率求新的节点电压，根据前后两次迭代的电压偏差是否小于设定阈值判断是否收敛， 如果满足， 则输出结果， 如果不满足 ,必须重复上述过程直到满足收敛条件为止。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9
20、页，共 17 页 - - - - - - - - - 华南理工大学广州学院电气工程学院课程设计报告9 4.5 计算原理第一步，从与叶节点联系的支路开始，该支路的末端功率即等于叶节点功率，利用这个功率和对应的节点电压计算支路功率损耗，求得支路的首端功率。当以某节点为始节点的各支路都计算完毕后，便想象将这些支路都拆去，是的该节点成为新的叶节点，其节点功率等于原有的负荷功率与该节点为始节点的各支路首端功率之和。于是计算便可延续下去，直到全部支路计算完毕。第二步，利用第一步所得的支路首端功率和本步骤刚算出的本支路始节点的电压（对电源点为已知电压） ，从电源点开始逐条支路进行计算，求得各支路终节点的电压
21、。在迭代计算开始之前，要先处理好支路的计算顺序问题。介绍两种确定之路计算顺序的方法。第一种方法是，按与叶节点联接的之路顺序，并将已排序的之路拆除，在此过程中将不断出现新的叶节点，而与其联接的支路又加入排序行列。这样就可以全部排列好从叶节点向电源点计算功率损耗的支路顺序。其逆序就是进行电压计算的支路顺序。第二种是逐条追加支路的方法。首先从根节点（电源点）开始接出第一条支路，引出一个新节点，以后每次追加的支路都必须从已出现的节点接出遵循这个原则逐条追加支路，直到全部支路追加完毕。所得到的支路追加顺序即是进行电压计算的支路顺序，其逆序便是功率损耗计算的支路顺序。名师资料总结 - - -精品资料欢迎下
22、载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 17 页 - - - - - - - - - 华南理工大学广州学院电气工程学院课程设计报告10 五、计算框图流程六、启动输入原始数据形成导纳矩阵 Y给定电压初值 e0，f0置K0利用式 (4-46)(4-47) 计算deltP(k),deltQ(k),deltU(k)IdeltP(k),deltQ(k)I0.0001 PQ2=PQ; %每一次迭代各节点的注入有功和无功相同PL=0.0; for i=1:NB kf=FT(i,1); %前推始节点号kt=FT(i
23、,2); %前推终节点号x=(PQ2(kf,2)2+PQ2(kf,3)2)/(V(kf)2); %计算沿线电流平方A PQ1(i,1)= RX(i,1)*x+PQ2(kf,2); %计算支路首端有功/MW RX(i,1)R PQ1(i,2)= RX(i,2)*x+PQ2(kf,3); %计算沿支路的无功损耗/Mvar RX(i,2)X PQ2(kt,2)= PQ2(kt,2)+PQ1(i,1); %用 PQ1 去修正支路末端节点的有功P单位 MW PQ2(kt,3)= PQ2(kt,3)+PQ1(i,2); %用 PQ1 去修正支路末端节点的有功Q 单位 Mvar PL=PL+RX(i,1)*
24、x; %累积计算整个网络的功率损耗end angle(1)=0.0; for i=NB:-1:1 kf=FT(i,2); %回代始节点号kt=FT(i,1); %回代终节点号dv1=(PQ1(i,1)*RX(i,1)+PQ1(i,2)*RX(i,2)/V(kf); %计算支路电压损耗的纵分量dv1 dv2=(PQ1(i,1)*RX(i,2)-PQ1(i,2)*RX(i,1)/V(kf); %计算支路电压损耗的横分量dv2 V2(kt)=sqrt(V(kf)-dv1)2+(dv2)2); %计算支路末端电压/kV angle(kt)= atand(dv2/(V(kf)-dv1)+angle(kf
25、); %计算支路电压相角（单位度）end maxd=abs(V2(2)-V(2); V2(1)=V(1); for i=3:1:NN if abs(V2(i)-V(i)maxd; maxd=abs(V2(i)-V(i); end end maxd k=k+1 PQ1 %潮流分布即支路首端潮流MV A V=V2 %节点电压模计算结果kV 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 17 页 - - - - - - - - - 华南理工大学广州学院电气工程学院课程设计报
26、告12 angle %节点电压角度计算结果单位度PL %网损单位MW end clear 输入文件清单：function PQ,FT,RX=case115() PQ= %节点电压有功无功10.4 0 0 10.0 0.0342 0.0301 10.0 0.0693 0.0642 10.0 0.0845 0.0763 10.0 0.0295 0.0261 10.0 0.0474 0.0409 10.0 0.1176 0.0957 10.0 0.0946 0.0857 10.0 0.0916 0.0859 10.0 0.0271 0.0229 10.0 0.0696 0.0643 10.0 0.0
27、676 0.0579 10.0 0.0298 0.0242 ; FT= %首端末端5 4 13 4 12 11 4 3 10 3 11 3 7 6 9 8 3 2 6 2 8 2 2 1 ; RX= % R X 4.524 5.04 3.521 3.966 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 17 页 - - - - - - - - - 华南理工大学广州学院电气工程学院课程设计报告13 2.436 2.866 1.145 1.28 4.14 4.696 1.
28、328 1.763 2.745 2.965 2.237 2.756 2.356 2.541 0.856 1.14 3.743 4.251 3.367 3.685 ; 计算结果：k =8 PQ1 = 0.0296 0.0262 0.0299 0.0243 0.0678 0.0581 0.1444 0.1273 0.0272 0.0230 0.1379 0.1231 0.1183 0.0964 0.0920 0.0864 0.3852 0.3445 0.1661 0.1379 0.1891 0.1750 0.8116 0.7278 V =10.4000 9.8795 9.6991 9.6652 9
29、.6377 9.8492 9.7872 9.7326 9.6870 9.6764 9.6578 9.6235 9.6444 angle = 0 0.3011 0.4011 0.4250 0.4438 0.3431 0.3942 0.3899 0.4266 0.4209 0.4498 0.4823 0.4452 PL =0.0488数据分析：经过 8 次迭代，各节点的电压误差均在0.0001kv 以内，计算到此结束。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 17 页
30、 - - - - - - - - - 华南理工大学广州学院电气工程学院课程设计报告14 七、配电网潮流计算的要点1.需要注意系统中所有母线的电压是否在允许的范围以内。2.系统中各种元件 ( 线路、变压器等 )是否会出现过负荷。3.在进行电压和功率分布计算以前，先要对网络等值电路作些简化处理。4.确立前推回代的计算次序。5.在开始编写程序之前，要了解matlab 仿真软件，包括其程序的编写，程序的运行，我们都要清楚的知道6.为了提高计算精度，重复计算，给出一个容许误差作为计算的收敛的判据。八、自我总结此次课程设计首先让我明白了要使电力系统运行的稳定，必须经过精密的设计和计算。在进行课题设计的过程
31、中，加深了我对潮流计算的认识，尤其是对牛顿拉夫逊潮流计算的求解思路有了比较透彻的理解。在此次程序设计中，我在MATLAB 编程，13 节点配电网潮流计算， WORD 文档的编辑方面均有所提高，之前学过的潮流计算仅限于书本，通过完成这个作业使我对这块知识了解得更加全面，但也暴露出了一些问题：理论知识还是不足，基本知识不充足，做起这个大作业难度还是很大的，对MATLAB 的性能和特点还不能有一个全面的把握，相信通过以后的学习能弥补这些不足，从而达到一个新的层次。在过程设计中我发现了一些要点，要求得支路的首端功率，从与叶节点联接的支路开始，该支路的末端功率即等于叶节点功率，利用这个功率和对应的节点电
32、压计算支路功率损耗，才能求首段功率。 还有在迭代计算开始之前， 先要确定支路的计算顺序问题，即确定好是深度优先还是广度优先。要求得各支路终结点的电压的话利用上述所得的支路首端功率和本步骤刚算出的本支路始节点的电压，从电源点开始逐条支路进行计算。而且在此次课程设计中，我发现了自己的基础知识有很多的不足。这些基础的缺乏给我的设计计划造成了不小的障碍。在这个过程中，我明白了，只要用心去做，认直去做，持之以恒，就会有新的发现，有意外的收获。虽然很多错误被克服了，但是系统中难免还有很多不足之处，有些地方感觉做的不是很正确，希望在以后的学习中能够在同学们的帮助和老师的指导下做得更好。名师资料总结 - -
33、-精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 17 页 - - - - - - - - - 华南理工大学广州学院电气工程学院课程设计报告15 九、参考文献1 何仰赞温增银 . 电力系统分析 . 华中科技大学出版社，2002 2 蔡旭晖、刘卫国、蔡立燕. MATLAB基础与应用教程D. 北京：人民邮电出版社 2009 附录一 MATLAB的简介MATLAB是由美国 mathworks 公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化
34、以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB 和 Mathematica、Maple 并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。名师资料总结
35、 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页，共 17 页 - - - - - - - - - 华南理工大学广州学院电气工程学院课程设计报告16 MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB 来解算问题要比用C，FORTRAN 等语言完成相同的事情简捷得多，并且 mathwork 也吸收了像 Maple 等软件的优点 ,使 MATLAB 成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C+ ，JAVA 的支持。可
36、以直接调用,用户也可以将自己编写的实用程序导入到MATLAB 函数库中方便自己以后调用，此外许多的 MATLAB 爱好者都编写了一些经典的程序， 用户可以直接进行下载就可以用。MATLAB 产品族可以用来进行以下各种工作： 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与仿真 数字图像处理技术 数字信号处理技术 通讯系统设计与仿真 财务与金融工程MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用 MATLAB 函数集） 扩展了 MATLAB 环境， 以解决这些应用领域内特定类型问题。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 17 页，共 17 页 - - - - - - - - -