case6bus

CASE6BUS  6-bus system
   From in problem 3.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] = case6bus
0002 %CASE6BUS  6-bus system
0003 %   From in problem 3.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: case6bus.m,v 1.3 2010/04/19 14:40:12 ray Exp $
0009 
0010 %%-----  Power Flow Data  -----%%
0011 %% system MVA base
0012 baseMVA = 100;
0013 
0014 %% bus data
0015 %    bus_i    type    Pd    Qd    Gs    Bs    area    Vm    Va    baseKV    zone    Vmax    Vmin
0016 bus = [
0017     1    3    0.25    0.1    0    0    1    1    0    230    1    1.1    0.9;
0018     2    2    0.15    0.05    0    0    1    1    0    230    1    1.1    0.9;
0019     3    1    0.275    0.11    0    0    1    1    0    230    1    1.1    0.9;
0020     4    1    0    0    0    0    1    1    0    230    1    1.1    0.9;
0021     5    1    0.15    0.09    0    0    1    1    0    230    1    1.1    0.9;
0022     6    1    0.25    0.15    0    0    1    1    0    230    1    1.1    0.9;
0023 ];
0024 
0025 bus(:, 3) = bus(:, 3)*baseMVA;
0026 bus(:, 4) = bus(:, 4)*baseMVA;
0027 
0028 %% generator data
0029 % Note:
0030 % 1)It's better of gen to be in number order, otherwise gen and genbid
0031 % should be sorted to make the lp solution output clearly(in number order as well)
0032 % 2)set Pmax to nonzero. set to 999 if no limit
0033 % 3)If change the order of gen, then must change the order in genbid
0034 % accordingly
0035 %    bus    Pg    Qg    Qmax    Qmin    Vg    mBase    status    Pmax    Pmin
0036 gen = [
0037     1    0    0    100    -100    1.05    100    1    100    0;
0038     2    0.5    0    100    -100    1.05    100    1    100    0;
0039 ];
0040 
0041 gen(:, 2) = gen(:, 2)*baseMVA;
0042 
0043 %% branch data
0044 %    fbus    tbus    r    x    b    rateA    rateB    rateC    ratio    angle    status
0045 branch = [
0046     1    4    0.020    0.185    0.009    999    100    100    0    0    1;
0047     1    6    0.031    0.259    0.010    999    100    100    0    0    1;
0048     2    3    0.006    0.025    0    999    100    100    0    0    1;
0049     2    5    0.071    0.320    0.015    999    100    100    0    0    1;
0050     4    6    0.024    0.204    0.010    999    100    100    0    0    1;
0051     3    4    0.075    0.067    0    999    100    100    0    0    1;
0052     5    6    0.025    0.150    0.017    999    100    100    0    0    1;
0053 ];
0054 
0055 %branch(:, 3) = branch(:, 3)*1.75;
0056 
0057 return;