case3bus_P6_6

CASE3BUS_P6_6  Case of 3 bus system.
   From Problem 6.6 in book 'Computational
   Methods for Electric Power Systems' by Mariesa Crow
   created by Rui Bo on 2007/11/12

0001 function [baseMVA, bus, gen, branch, areas, gencost] = case3bus_P6_6
0002 %CASE3BUS_P6_6  Case of 3 bus system.
0003 %   From Problem 6.6 in book 'Computational
0004 %   Methods for Electric Power Systems' by Mariesa Crow
0005 %   created by Rui Bo on 2007/11/12
0006 
0007 %   MATPOWER
0008 %   $Id: case3bus_P6_6.m,v 1.3 2010/04/19 14:40:12 ray Exp $
0009 
0010 %%-----  Power Flow Data  -----%%
0011 %% system MVA base
0012 baseMVA = 1000;
0013 
0014 %% bus data
0015 %    bus_i    type    Pd    Qd    Gs    Bs    area    Vm    Va    baseKV    zone    Vmax    Vmin
0016 bus = [
0017     1   3   350       100    0    0    1    1    0    230    1    1.00    1.00;
0018     2   2   400       250    0    0    1    1    0    230    1    1.02    1.02;
0019     3   2   250     100    0    0    1    1    0    230    1    1.02    1.02;
0020 ];
0021 
0022 %% generator data
0023 % Note:
0024 % 1)It's better of gen to be in number order, otherwise gen and genbid
0025 % should be sorted to make the lp solution output clearly(in number order as well)
0026 % 2)set Pmax to nonzero. set to 999 if no limit
0027 % 3)If change the order of gen, then must change the order in genbid
0028 % accordingly
0029 %    bus    Pg    Qg    Qmax    Qmin    Vg    mBase    status    Pmax    Pmin
0030 gen = [
0031     1    182.18   0    999    -999    1.00       100    1    600    0;
0032     2    272.77      0    999    -999    1.02       100    1    400    0;
0033     3    545.05   0    999    -999    1.02       100    1    100    0;
0034 ];
0035 %gen(:, 9) = 999; % inactive the Pmax constraints
0036 
0037 %% branch data
0038 %    fbus    tbus    r    x    b    rateA    rateB    rateC    ratio    angle    status
0039 branch = [
0040     1    2    0.01     0.1    0.050    999    100    100    0    0    1;
0041     1    3    0.05     0.1    0.025    999    100    100    0    0    1;
0042     2    3    0.05     0.1    0.025      999    100    100    0    0    1;
0043 ];
0044 
0045 %%-----  OPF Data  -----%%
0046 %% area data
0047 areas = [
0048     1    1;
0049 ];
0050 
0051 %% generator cost data
0052 %    2    startup    shutdown    n    c(n-1)    ...    c0
0053 gencost = [
0054     2    0    0    3    1.5     1    0;
0055     2    0    0    3    1       2    0;
0056     2    0    0    3    0.5      2.5    0;
0057 ];
0058 
0059 return;