matpower牛顿拉夫逊法计算机程序研究分析.docVIP

课程作业 牛顿-拉夫逊法计算机程序研究 课程名称：电力系统分析 指导教师： 姓 名： ___ 学号：____________ 年级专业班级：_____________ 提交日期 2014年1月12日 1概念 潮流计算是电力系统分析中的一种最基本的计算，它的任务是对给定的运行条件确定系统的运行状态，比如各母线上的电压幅值与相角、网络中的功率分布及功率损耗等。在简单电力网络中，一般可采取手工计算方法，如单端供电网络，给定首端电压以及末端功率，从末端向前推出功率损耗，再从首端向后推出电压损耗。然而，实际的电力系统十分复杂，少则几十个节点，多则上千节点，此时采取手工运算基本不可能实现运行要求，计算机代替手算的方法由此产生。它服务于大系统，较之手算，速度快，结果精确，能够满足电力系统运行要求。本文主要研究了基于matlab的牛顿-拉夫逊潮流计算方法，结合书本例子，验证了该程序的实用性。 2分析方法 网络方程式（如节点方程）是潮流计算的基础方程式。如果能够给出电压源（或电流源），直接求解网络方程就可以求得网络内电流和电压的分布。但是在潮流计算中，在网络的运行状态确认以前，无论是电源的电势，还是节点的注入电流都是无法事先给定的。 对于一个三节点简单电力系统，其网络方程为: （1） 将节点电流用节点功率与电压表示后代入上式，这样n节点系统的潮流方程为 = (i=1,2,…,n) （2） 将上述方程的实部，虚部分开，对于每一个节点课的两个实数方程，但是变量仍有4个，P，Q,V, 。我们需要给定其中两个，这样方程就有解了。按给定量的不同，分为PQ节点（又叫负荷节点），PV节点（又叫电压控制节点），还有平衡节点（给定V，）。下面采用的方法，是将节点电压表示为极坐标形式 （3） 将（3）带入（2）可得与的表达式，即由电压幅值，相角计算不平衡量的公式。 具体有下述两种方法： 1）直角坐标算法： 给定e(k),f(k)，算出△P(k),△Q(k),然后与预设偏差作比较，如果符合则进行全网潮流计算，如果不符合，则用雅各比矩阵算出△e(k),△f(k),修正得e(k+1),f(k+1)，重新计算△P(k+1),△Q(k+1）并与预设值比较。 2）极坐标算法： 给定V（k）,δ（k），算出△P(k),△Q(k),然后与预设偏差作比较，如果符合则进行全网潮流计算，如果不符合，则用雅各比矩阵算出△V(k),△δ(k),修正得V(k+1),δ(k+1)，重新计算△P(k+1),△Q(k+1）并与预设值比较。 3具体例子 以《电力系统分析》（何仰赞著）61页例11-5为例进行验证。采用极坐标法。 等效电力网络图 首先用busdata(节点矩阵)矩阵描述1到4号节点的数据： Busdata 从左到右依次为 bus_i（节点号）type（节点类型） Pd（负荷需求有功功率/MVA）Qd（符合需求无功功率/MVar）Gs(并联电导的功率/MVA) Bs(并联电纳的功率/MVar) area（电力区域号） Vm（节点电压幅值）Va（节点电压相角）baseKV（节点处实际电压/kV） zone(损耗区域) Vmax（最大电压幅值） Vmin（最小电压幅值） mpc.bus=[1 1 30 18 0 0 1 1 0 110 1 1.1 0.9 2 1 55 13 0 0 1 1 0 110 1 1.1 0.9 3 2 0 0 0 0 1 1.1 0 110 1 1.1 0.9 4 3 0 0 0 0 1 1.05 0 110 1 1.1 0.9]; 然后用gendata(电机矩阵)描述系统输入 Generator data bus（节点号） Pg（发电机发出有功/MVA） Qg(发电机发出无功/MVar) Qmax（最大发出无功） Qmin（最小发出无功） Vg（给定电压幅值标幺） mBase（电机容量/MVA） status(运行状态) Pmax（最大发出有功/MVA）Pmin(最小发出有功/MVA) Pc1(PQ性能曲线下的实际输出功率) Pc2（PQ性能曲线上的实际输出功率） Qc1min（pc1下最小输出功率） Qc1max（pc1下最大输出功率） Qc2min（pc2下最小输出功率） Qc2max（pc2下最大输出功率） ram