电力网潮流电压计算例题与matlab程序分解.docxVIP

 PAGE \* MERGEFORMAT 9 电力网潮流电压计算例题与MATLAB程序 编写 佘名寰 本文介绍了两个电力网潮流计算例题。一例为5个节点5条支路；另一例有6个节点7条支路，有PQ负荷节点也有PV发电机节点，变压器支路标么变比不为1。本文给出了完整的计算程序和计算结果。程序包括计算网络节点导纳矩阵和潮流电压两部分，程序中引用的M函数见参考文献1。例题选自研究生教材，比较典型，计算结果与教材核对基本一致。本文可供电力专业师生和基层技术人员潮流计算时参考。 【例1.3】如图 1-4所示，五个节点三条母线简单系统，系统参数均用标么值（pu.）表示。节点4、5通过变压器与节点2、3相连，变比为1.05:1，节点1、2、3为PQ负荷节点，节点4为PV发电机节点，节点5为参考节点，节点4、5电压均为1.05，图中箭头表示各支路的正方向，各个节点功率如表2.1所示,计算各节点电压、发电机无功功率和系统潮流分布： 4 Z23 Z13 Z21 YC2 YC3 YC1 Z42 Z53 1:1.05 1.05:1 4 2 3 1 5 S3 S2 S1 Z21=0.04+J0.25 Z23=0.08+J0.30 Z13=0.1+J0.35 Z42=J0.015 Z53=J0.03 YC1=J0.25 YC2=J0.50 YC3=J0.25 S1=1.6+J0.8 S2=2+J1 S3=3.7+J1.3 P4=5.0 V4=1.05 V5=1.05 P4=5.0 V4=1.05v V5=1.05 图1-4 [例1.3]网络接线图 表 1-1 母线有功负载无功负载发电机有功发电机无功1PD1=1.6QD1=0.82PD2=2.0QD2=1.03PD3=3.7QD3=1.34 5PG4=5.0 PG5=？ QG4=未定 QG5=？ 【例1.3】 程序说明： 本例题与参考文献2的【例2.3.1】不同之处在于有两个变比不为1的变压器支路，并且引入了PV发电机节点。由于网络结构与例图1-1完全相同，网络节点导纳矩阵已在文献1之1.3节求出不再赘述。潮流计算程序参考【例2.3.1】 节点电压符号矩阵：一共五个节点，独立节点数N1=4，PV节点数N2=1 节点电压幅值符号矩阵为 u=sym([u1,u2,u3,u4,u5]); 节点电压相角符号矩阵为 delt=sym([d1,d2,d3,d4,d5]); 节点功率给定值为：三个PQ节点，一个PV 节点 p=[-1.6,-2,-3.7,5];q=[-0.8,-1.0,-1.3]; 【例1.3】 源程序 NU114N.m %******* NU114N.m , example 2-3-2, fig 2-2 ******** % %The following Program for load flow calculation is based on MATLAB2007 clear % bus 1,2,3 is PQ bus,bus 4 is PV bus,bus 5 is slack bus global Np Nb bt p m k % Np is number of node point,Nb is number of braches, % p the transformer turns ratio,with off-nominal tap-setting,fig 2-4 ,PSCLF p6, % bt is two-dimensional array ,line one is point 1,line two is point 2, % m is number of node point, k is number of transformers Np=5;Nb=5; bt=[2,3;4,5];p=[1/1.05,1/1.05];m=5;k=2; % nstart--the start point of branches ,nend -- the end point, % mm -- network incidence matrix nstart=[2,2,1,4,5]; nend=[1,3,3,2,3]; mm=ffm(nstart,nend); % zb1,the series impedances of transmission line % yb1,the series admittances of transmission line zb1r=[0.04