从基因工程菌中释放提取白介素13的研究-北京石油化工学院.docVIP

基于MATLAB的电力系统潮流计算分析与仿真 北京石油化工学院2013届优秀毕业论文 PAGE 10 PAGE 9 基于MATLAB的电力系统潮流计算 分析与仿真 电气工程专业 黄 勐 哲 指导教师 张红洁 讲师 摘 要 本文首先对电力系统潮流计算进行了深入的研究，熟练掌握了牛顿－拉夫逊极坐标形式和直角坐标形式的两种算法的原理。接着，为电力系统网络中各元件建立了等效的数学模型：输电线路用π型等值电路等效，考虑了导线对地导纳；变压器用π型等值电路等效。进一步，为电力系统网络建立极坐标与直角坐标形式的牛顿－拉夫逊潮流计算数学模型。在此基础上，运用MATLAB语言编写了基于牛顿－拉夫逊方法的潮流计算程序，该程序可以实现原始数据输入、自动形成节点导纳矩阵、形成雅可比矩阵和求解满足精度要求的各节点电压和功率等功能。为验证程序的准确性，将该程序应用于某三节点环形电力网络潮流计算，并与手工算法求得的理论结果进行对比。最后，采用该程序对某复杂环形电力网络进行潮流计算。 关键词 潮流计算，牛顿－拉夫逊法，MATLAB，电力系统 1 前言 潮流计算是研究电力系统稳态运行情况的一种重要的分析计算,随着电力系统越来越复杂，使得计算机计算潮流成为主流，目前计算机计算潮流有三种算法：高斯赛德尔、牛顿-拉夫逊、PQ分解方法[1,2,3]。而牛顿-拉夫逊法是电力系统潮流计算的常用算法之一，它收敛性好，迭代次数少,精度高。且MATLAB使用方便，有着其他高级语言无法比拟的强大的矩阵处理功能[4]，所以本毕设研究的方向为基于MATLAB的电力系统牛顿-拉夫逊法潮流计算分析与仿真。 本毕设在实现潮流计算基本功能的同时，还扩展了一些新的功能。具体如下：根据从命令窗口输入i的值来选择算法，增加了程序算法可选择性；扩展了根据自己的需要规定相应的精度eps的功能；根据导纳矩阵形成原理自定义了函数Y=matrix(y,n,y0,Ltype,KT)来形成导纳矩阵，省去了以往手工算导纳矩阵花费的大量精力；增加了各支路间流动功率，各节点注入电流和各节点间线路功率损耗的计算功能；将已知线路串联的导纳y，线路并联的导纳y0，节点的功率S，线路类型Ltype（有变压器时置2，无变压器时置1）,变压器变比KT，存入相应矩阵中，并将矩阵保存为mat文件，程序运行前通过加载命令load调入。避免每次计算时重复输入线路参数。 2 潮流计算的数学模型 2.1 潮流计算的节点模型 根据电力系统中各节点性质的不同，很自然地把节点分成三类[1,2,3,5]： PQ节点 对这一类点，事先给定的是节点功率(P，Q)，待求的未知量是节点电压向量(U，)，所以叫PQ节点。 PV节点 这类节点给出的参数是该节点的有功功率P及电压幅值U，待求量为该节点的无功功率Q及电压向量的相角。 平衡节点 在潮流计算中，这类节点一般只设一个。对该节点，给定其电压值，并在计算中取该节点电压向量的方向作为参考轴，相当于给定该点电压向量的角度为零。也就是说，对平衡节点给定的运行参数是U和，因此待求量是该节点的P和Q，整个系统的功率平衡由这一节点承担。 以上三类节点4个运行参数P、Q、U、中，已知量都是两个，待求量也是两个，只是类型不同而已。 2.2 修正方程 采用牛顿法计算潮流时，需要对功率方程进行修改。下面将根据在不同坐标内的修改进行讨论： 当在极坐标系内时，计算的PQ节点功率与给定的PQ节点的注入功率的差值满足下面方程： （2－1） 式中：——i与j节点电压的相角差。 由方程（2－1）可以写出修正方程式如下: （2－2） 其中：H是阶方阵，其元素为；N是阶矩阵，其元素为；K是阶矩阵，其元素为；L是阶矩阵，其元素为。 对式 （2－1）求偏导数，可得雅克比矩阵元素的表达式如下： 非对角元素（） （2－3） 对角元素（） （2－4） 在直角坐标系内时，计算的PQ 节点功率给定的PQ 节点的注入功率的差值应满足以下方程： （2－5） 对于PV 节点，已知节点的注入有功功率及节点电压大小，记作，其节点的有功功率应满方程： （2－6） 对于平衡节点，因为其电压给定，故不需要迭代求解。 由方程（2－5）和（2－6）可以得到修正方程如下： （2－7） 式中 其中：雅克比矩阵的各元素可以对式（