【基础讲解】基于matpower的电力系统潮流计算（Matlab代码实现）

 👨‍🎓个人主页：研学社的博客 

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述与安装

第一步：

第二步： 

 第三步：

第四步： 

 第五步：

第六步：

第七步：

📚2 入门讲解

主函数：

子函数：

🌈3 基于Matpower的电力系统潮流计算原理及仿真设计

🎉4 参考文献

---

​

💥1 概述与安装

matpower是基于 Matlab 语言的m文件组建包，作为电力系统潮流计算仿真软件。潮流计算仿真中，按照matpower中潮流计算程序的格式修改或重新编写一个新的程序“caseX.m”，在MATLAB软件中的命令窗口输入runpf（‘程序名’）就可以通过matpower已经编好的程序进行潮流计算，而通过调用函数runpf的参数计算潮流。
 

matpower安装包官网下载链接： Downloads – MATPOWER

第一步：

解压

​

我们可以看到matpower的安装包实际为一个Matlab程序包：

 ​

第二步： 

将matpower文件复制在C:\Program Files (x86)\MATLAB\R2016b\bin目录下。

 ​

 第三步：

找到菜单栏的【设置路径】

​

第四步： 

点击【添加并包含子文件】

​

 第五步：

选择所要添加的matpower文件夹，点击【选择文件夹】

​

原来的窗口中多了很多文件路径，点击【保存】，再点击【关闭】

第六步：

在MATLAB命令行中输入test_matpower

​

第七步：

完美，成功

！！！ 

​ 

📚2 入门讲解

以最简单的IEEE5节点为例说明：

主函数：

clc
clear
close all
runpf('IEEE5')

子函数：

function mpc = IEEE5
% MATPOWER Case Format : Version 2
mpc.version = '2';
%%-----  Power Flow Data  -----%%
%% system MVA base
mpc.baseMVA = 100;
 
%% bus data
%  bus_i   type    Pd  Qd  Gs  Bs  area    Vm  Va  baseKV  zone    Vmax    Vmin
mpc.bus = [
    1   1   160 80  0   0   1   1       0   100   1   1.1    0.94;
    2   1   200 100 0   0   1   1       0   100   1   1.1    0.94;
    3   1   370 130 0   0   1   1       0   100   1   1.1    0.94;
    4   2   0   0   0   0   1   1.050   0   100   1   1.1   0.94;
    5   3   0   0   0   0   1   1.050   0   100   1   1.1   0.94;
    ];
%% generator data
%   bus Pg  Qg  Qmax    Qmin    Vg  mBase   status  Pmax    Pmin
mpc.gen = [
    4   500 0   99990   -9999   1.050   100 1   600 0;
    5   0   0   99990   -9999   1.050   100 1   600 0;
    ];
 
%% branch data
%  fbus tbus  r   x   b   rateA   rateB   rateC   ratio   angle   status  angmin  angmax
mpc.branch = [
    2   1   0.04    0.25    0.5 0   0   0   0     0   1   -360    360;
    3   1   0.1     0.35    0   0   0   0   0      0   1   -360    360;
    3   2   0.08    0.3     0.5 0   0   0   0     0   1   -360    360;
    3   5   0     0.03    0   0   0   0   1.05    0   1   -360    360;
    2   4   0      0.015   0   0   0   0   1.05    0   1   -360    360;
    ];
return;

function mpc = IEEE5
% MATPOWER Case Format : Version 2
mpc.version = '2';
%%-----  Power Flow Data  -----%%
%% system MVA base
mpc.baseMVA = 100;

%% bus data
%  bus_i   type    Pd  Qd  Gs  Bs  area    Vm  Va  baseKV  zone    Vmax    Vmin
mpc.bus = [
    1   1   160 80  0   0   1   1       0   100   1   1.1    0.94;
    2   1   200 100 0   0   1   1       0   100   1   1.1    0.94;
    3   1   370 130 0   0   1   1       0   100   1   1.1    0.94;
    4   2   0   0   0   0   1   1.050   0   100   1   1.1   0.94;
    5   3   0   0   0   0   1   1.050   0   100   1   1.1   0.94;
    ];
%% generator data
%   bus Pg  Qg  Qmax    Qmin    Vg  mBase   status  Pmax    Pmin
mpc.gen = [
    4   500 0   99990   -9999   1.050   100 1   600 0;
    5   0   0   99990   -9999   1.050   100 1   600 0;
    ];

%% branch data
%  fbus tbus  r   x   b   rateA   rateB   rateC   ratio   angle   status  angmin  angmax
mpc.branch = [
    2   1   0.04    0.25    0.5 0   0   0   0     0   1   -360    360;
    3   1   0.1     0.35    0   0   0   0   0      0   1   -360    360;
    3   2   0.08    0.3     0.5 0   0   0   0     0   1   -360    360;
    3   5   0     0.03    0   0   0   0   1.05    0   1   -360    360;
    2   4   0      0.015   0   0   0   0   1.05    0   1   -360    360;
    ];
return;
 

🌈3 基于Matpower的电力系统潮流计算原理及仿真设计

这个博主总结得很棒，致谢之：

基于Matpower的电力系统潮流计算原理及仿真设计

​​​​​​​

🎉4 参考文献

部分理论来源于网络，如有侵权请联系删除。