matpwer牛顿拉夫逊法计算机程序研究.doc

课程作业 牛顿-拉夫逊法计算机程序研究 课程名称：电力糸统分析 指导教师： 姓名： 姓名： 学号： 年级专业班级： 提交日期2014年1月12日 1概念 潮流计算是电力系统分析中的-?种最难本的计算,它的任务是对给定的运行条件确定系统的 运行状态，比如各母线上的电压幅值与相角、网络中的功率分布及功率损耗等。在简单电力 网络屮，一般可采取手工计算方法，如单端供电网络，给定酋端电压以及末端功率，从末端 向而推出功率损耗，再从首端向后推出电压损耗。然而，实际的电力系统十分复杂，少则儿 十个节点，多则上千节点，此时采取手工运算难木不可能实现运行要求，计算机代替手算的 方法由此产生。它服务于人系统，较之手算，速度快，结果精确，能够满足电力系统运行要 求。本文主要研究了基于matlab的牛顿-拉夫逊潮流计兑方法，结合书本例子，验证了该程 序的实川性。 2分析方法 网络方程式(如节点方程)是潮流计算的基础方程式。如果能够给出电压源(或电流源)， 直接求解网络方程就讨以求将网络内电流和电压的分布。但是在潮流计算中，4:网络的运行 状态确认以前，无论是电源的电势，还是节点的注入电流都是无法事先给定的。 对于一个三节点简单电力系统，其网络方程为： =rfai?a+r￡2^ + r^ = (i) 将节点电流川节点功率与电压表示后代入上式，这样n节点系统的潮流方程为 * * 巧 + 4 (i= 1，2,…，n) (2) 7=1 将上述方程的实部，虚部分开，对于每一个节点课的两个实数方程，似是变量仍有4个，P, Q，V，S。我们需要给定其中两个，这样方程就有解了。按给定量的不M，分为PQ W点(又 叫负份节点)，PV节点(乂叫电压控制节点)，还奋T衡节点(给定v, 3)。下面采川的方 法，是将节点电压表示为极叱标形式 V}=VZ6l=Vi(cos6i + ysin^.) (3) 将(3)带入(2)可得巧与化的表达式，即巾电压幅值，相角计算不平衡量的公式。 具体有K述两种方法： 直角坐标算法： 给定e(k)，f(k),算出AP(k)，AQ(k)，然p与预设偏差作比较，如果符合则进行全网潮流计算， 如果不符合，则用雅各比矩阵算出△e(k)，Af(k),修正得e(k+l)，f(k+l),重新计算△P(k+1)，A Q(k+1)并与预设值比较。 极少标算法： 给定V(k)，3 (k),算出AP(k)，AQ(k)，然后与预设偏差作比较，如果符合则进行全网潮流 计算，如果不符合，则用雅各比矩阵算出△〃(!)，△ 5⑻，修正得V(k+l),S(k+l),重新计算 AP（k+l）,AQ（k+l）并与预设值比较 3具体例子 以《电力系统分析》（何仰赞著）61页例11-5为例进行验证。采用极坐标法。 首先川busdata（卞点矩阵）矩阵描述1到4兮节点的数据： Busdata 从左到右依次为busj （节点号）type （节点类型）Pd （负荷需求柯功功率/MVA） Qd （符 介需求无功功率/MVar） Gs（并联电导的功率/MVA）Bs（并联电纳的功率/MVar）area （电力区 域号） Vm （ W点也压幅值）Va （节点屯压相角）baseKV （ W/、（处实阮也压/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]; 然P用gendata（电机矩阵）描述系统输入 Generator data bus （ W点号）Pg （发电机发出有功/MVA） Qg（发电机发岀无功/MVar） Qmax （S大发 出无功）Qmin （敁小发出无功）Vg （给定电压幅值标幺）mBase （电机容量/MVA） status（运行状态）Pmax （最人发出奋功/MVA） Pmin（最小发出奋功/MVA） Pcl（PQ性能|lh线 卜的实际输出功率）Pc2 （PQ性能曲线上的实际输出功率）Qclmin （pcf卜最小输出功率） Qc 1 max （pci下足人输fli功率）Qc2min （pc2下鐘小输出功率）Qc2max （pc2下姑人输fli 功率）ramp_agc （负荷跟踪/ AGC斜坡率） ramp_10 （10分钟储备斜坡率）ramp_30 （30 分钟储备斜坡率）ramp_q （2秒内功率曲线斜率）apf （区域参与因素） mpc.gen= [3 50 0 300 -300 1.1 100 1 250 10 00000000000 400