matpower在电力系统潮流计算中的应用,matpower

《matpower在电力系统潮流计算中的应用,matpower》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、MATPOWER在电力系统潮流计算中的应用MATPOWER是一个用MATLAB的M文件编写,用来解决电力系统潮流计算和优化潮流计算问题的软件包。MATPOWER特点是简单、易懂且程序代码公开，这为电力系统专业学生深入学习和理解掌握潮流计算中的难点（如节点导纳矩阵、算法及迭代过程等）提供了一个开放、便捷的平台。下面是MATPOWE4.0在电力系统潮流计算中的一个具体实例，仅供学习电力系统专业的读者参考。读者利用MATPOWER进行潮流计算时，结构体mpc各字段的形式可参照case2_5.m进行编写。图12机5节点系统对图1的2机5节点系统（具体参数

2、参见：现代电力系统分析-王锡凡），按“version2”格式编写成的case2_5.m的程序清单如下：functionmpc=case2_5%MATPOWERCaseFormat:Version2mpc.version='2';%%------PowerFlowData------%%%%systemMVAbasempc.baseMVA=100;%%Busdata%bus_itypePdQdGsBsareaVmVabaseKVzoneVmaxVminmpc.bus=[11160800011010011.10.94;2120010000110100

3、11.10.94;313701300011010011.10.94;42000011.05010011.10.94;53000011.05010011.10.94;];%%generatordata%busPgQgQmaxQminVgmbasestatusPmaxPminmpc.gen=[4500099990-99991.0510016000;50099990-99991.0510016000;];%%branchdata%fbustbusrxbrataArataBrataCratioanglestatusangminangmaxmpc.bran

4、ch=[210.040.250.5000001-360360;310.10.350000001-360360;320.080.30.5000001-360360;3500.0300001.0501-360360;2400.01500001.0501-360360;];Return当采用牛顿—拉夫逊法计算case2_5.m交流电网潮流时，在MATLAB的名令窗口输入以下命令即可>>runpf('case2_5')计算输出结果如下：MATPOWERVersion4.0,07-Feb-2011--ACPowerFlow(Newton)Newton'sm

5、ethodpowerflowconvergedin5iterations.Convergedin0.04seconds=====================================================================SystemSummary=====================================================================Howmany?Howmuch?P(MW)Q(MVAr)--------------------------------------

6、-------------------------------Buses5TotalGenCapacity1200.0-19998.0to199980.0Generators2On-lineCapacity1200.0-19998.0to199980.0CommittedGens2Generation757.9411.2Loads3Load730.0310.0Fixed3Fixed730.0310.0Dispatchable0Dispatchable0.0of0.00.0Shunts0Shunt(inj)0.00.0Branches5Losses

7、(I^2*Z)27.94204.78Transformers2BranchCharging(inj)——103.5Inter-ties0TotalInter-tieFlow0.00.0Areas1MinimumMaximum---------------------------------------------------------------------VoltageMagnitude0.862p.u.@bus11.078p.u.@bus2VoltageAngle-4.78deg@bus121.84deg@bus4PLosses(I^2*R

8、)-13.81MW@line3-2QLosses(I^2*X)-73.98MVAr@line2-1===================