现代电力系统潮流计算作业.doc

现代电力系统潮流计算作业 作业要求A组0.1 调用matpower中的runpf函数，分析输入文件中各矩阵定义； 0.2 调用某一个算例，输出潮流结果，并分析。 0.3 完成0.1和0.2的基础上，分析matpower中牛顿法和快速解耦法给出流程图写出newtonpf和fdpf函数每行程序定义 0.4 完成0 1 分析输入文件中各矩阵的定义 1.1 MATPOWER的安装 MATPOWER工具箱的安装步骤如下： 1） 下载matpower压缩包。官方下载网址：/matpower/， 目前最新版本为 2） 解压压缩包，得到文件夹matpower5.1，并将文件夹移动到MATLAB所在路径的toolbox文件夹下 C:\Program Files\MATLAB\R2016a\toolbox。 3） 添加地址到MATLAB路径。打开MATLAB，点击“文件”(“设置路径(“添加并包含子文件夹…找到power5.1所在的位置点击确定(“关闭 4） 测试matpower工具是否安装成功。在MATLAB命令行窗口输入test_matpower”，出现一系列的测试均显示ok”，最后显示“All tests successful (3256 passed, 682 skipped of 3938)”，则表示安装成功 1.2 矩阵的定义 打开文档“caseformat.m”，或者在MATLAB命令行窗口中输入“help caseformat”，可以得到关于输入矩阵的数据定义当然也可以参考docs文件夹下的manual文档其中进行了详细说明 Data），发电机参数矩阵（Generator Data），支路参数矩阵（Branch Data）。为了进行最优潮流的相关计算输入矩阵还包含发电机费用矩阵generator cost data）。以下对三种基本的输入参数矩阵数据格式进行详细说明。 表 列 名称 含义 单位 1 bus number 母线编号，从1开始，正整数 2 bus type 节点类型，1为PQ节点注入负荷的有功功率注入负荷的功功率 ar 5 Gs 与母线并联的电导 与母线并联的电纳Var 7 area number 电网断面号，一般设置为1 8 Vm 母线电压的幅值.u. 9 Va 母线电压的相 度 10 baseKV 基准电压 zone 损耗分区号，一般设置为1 12 Vmax 工作时母线的最高电压幅值 工作时母线的最电压幅值 表 列 名称 含义 单位 1 bus 发电机所在母线的编号 2 Pg 接入发电机的有功功率 MW 3 Qg 接入发电机的无功功率 MVar 4 Qmax 发电机的最大输出无功功率 MVar 5 Qmin 发电机的最小输出无功功率 MVar 6 Vg 发电机的工作电压幅值 p.u. 7 mBase 发电机的功率基准值，默认为baseMVA MVA 8 status 发电机是否投入运行，大于0表示投入运行，小于或等于0表示退出运行 9 Pmax 发电机的最大输出有功功率 MW 10 Pmin 发电机的最小输出有功功率 MW 表 列 名称 含义 单位 1 fbus 起始母线编号 2 tbus 终止母线编号 3 r 支路电阻 p.u. 4 x 支路电抗 p.u.支路充电电纳 p.u.rateA 支路长期运行允许的功率 MVA 7 rateB 支路短期运行允许的功率 MVA 8 rateC 支路紧急运行允许的功率 MVA 9 ratio 支路上变压器的变比如果支路元件变压器则该值为支路变压器的转角如果支路元件变压器则该值为 度 11 status 支路的初始工作状态 12 angmin 支路最小相角差 支路最相角差 1.3 case9数据分析 根据以上分析ase9，进行分析ase9.m文件包含两个变量和四个矩阵其中baseMVA=100表示功率的基准值为 表9的母线参数矩阵 2 3 4 5 6 7 8 9 10 11 12 13 bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin 1 3 0 0 0 0 1 1 0 345 1 1.1 0.9 2 2 0 0 0 0 1 1 0 345 1 1.1 0.9 3 2 0 0 0 0 1 1 0 345 1 1.1 0.9 4 1 0 0 0 0 1 1 0 345 1 1.1 0.9 5 1 90 30 0 0 1 1 0 345 1 1.1 0.9 6 1 0 0 0 0 1 1 0 345 1 1.1 0.9 7 1 100 35 0 0 1 1 0 345 1 1.1 0.9 8 1 0 0