👨🎓个人主页:研学社的博客
💥💥💞💞欢迎来到本博客❤️❤️💥💥
🏆博主优势:🌞🌞🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。
⛳️座右铭:行百里者,半于九十。
📋📋📋本文目录如下:🎁🎁🎁
目录
💥1 知识回顾
【基础讲解】基于matpower的电力系统潮流计算(Matlab代码实现)
📚2 case文件格式
我们打开一个case文件,例如‘case5’:
在这里,我们定义了一个mpc结构体,名称为‘case5’,往下滑动页面可以看到这个结构体中的所有字段:
version:数据文件格式的版本,共有两种——version1和version2。MATPOWER3.0及以前的版本采用version1。现在一般都采用version2。
🌈3 Matlab代码实现
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;