-
Notifications
You must be signed in to change notification settings - Fork 149
/
t_runmarket.m
119 lines (102 loc) · 4.2 KB
/
t_runmarket.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
function t_runmarket(quiet)
% t_runmarket - Tests for runmkt, smartmkt and auction.
% MATPOWER
% Copyright (c) 2005-2024, Power Systems Engineering Research Center (PSERC)
% by Ray Zimmerman, PSERC Cornell
%
% This file is part of MATPOWER.
% Covered by the 3-clause BSD License (see LICENSE file for details).
% See https://matpower.org for more info.
if nargin < 1
quiet = 0;
end
n_tests = 20;
t_begin(n_tests, quiet);
[PQ, PV, REF, NONE, BUS_I, BUS_TYPE, PD, QD, GS, BS, BUS_AREA, VM, ...
VA, BASE_KV, ZONE, VMAX, VMIN, LAM_P, LAM_Q, MU_VMAX, MU_VMIN] = idx_bus;
[GEN_BUS, PG, QG, QMAX, QMIN, VG, MBASE, GEN_STATUS, PMAX, PMIN, ...
MU_PMAX, MU_PMIN, MU_QMAX, MU_QMIN, PC1, PC2, QC1MIN, QC1MAX, ...
QC2MIN, QC2MAX, RAMP_AGC, RAMP_10, RAMP_30, RAMP_Q, APF] = idx_gen;
if ~have_feature('smartmarket')
t_skip(n_tests, 'smartmarket code not available');
else
mpc = loadcase('t_auction_case');
mpopt = mpoption('opf.ac.solver', 'MIPS', 'out.lim.all', 1, 'out.branch', 0, 'out.sys_sum', 0, 'out.all', 0, 'verbose', 1);
% mpopt = mpoption('out.gen', 1, 'out.branch', 0, 'out.sys_sum', 0);
offers.P.qty = [
12 24 24;
12 24 24;
12 24 24;
12 24 24;
12 24 24;
12 24 24;
];
offers.P.prc = [
20 50 60;
20 40 70;
20 42 80;
20 44 90;
20 46 75;
20 48 60;
];
bids.P.qty = [
10 10 10;
10 10 10;
10 10 10;
];
bids.P.prc = [
100 70 60;
% 100 64.3 20;
% 100 30.64545 0;
100 50 20;
100 60 50;
];
offers.Q.qty = [ 60; 60; 60; 60; 60; 60; 0; 0; 0 ];
offers.Q.prc = [ 0; 0; 0; 0; 0; 3; 0; 0; 0 ];
bids.Q.qty = [ 15; 15; 15; 15; 15; 15; 15; 12; 7.5 ];
% bids.Q.prc = [ 0; 0; 0; 0; 0; 0; 0; 83.9056; 0 ];
bids.Q.prc = [ 0; 0; 0; 0; 0; 0; 0; 20; 0 ];
t = 'marginal Q offer, marginal PQ bid, auction_type = 5';
mkt = struct( 'auction_type', 5, ...
't' , [], ...
'u0', [], ...
'lim', [] );
[r, co, cb, f, dispatch, success, et] = runmarket(mpc, offers, bids, mkt, mpopt);
co5 = co;
cb5 = cb;
% [ co.P.qty co.P.prc ]
% [ cb.P.qty cb.P.prc ]
% [ co.Q.qty co.Q.prc ]
% [ cb.Q.qty cb.Q.prc ]
i2e = r.bus(:, BUS_I);
e2i = sparse(max(i2e), 1);
e2i(i2e) = (1:size(r.bus, 1))';
G = find( ~isload(r.gen) ); %% real generators
L = find( isload(r.gen) ); %% dispatchable loads
Gbus = e2i(r.gen(G,GEN_BUS));
Lbus = e2i(r.gen(L,GEN_BUS));
t_is( co.P.qty, ones(6, 1) * [12 24 0], 2, [t ' : gen P quantities'] );
t_is( co.P.prc(1,:), 50.1578*ones(1,3), 3, [t ' : gen 1 P prices'] );
t_is( cb.P.qty, [10 10 10; 10 0.196 0; 10 10 0], 2, [t ' : load P quantities'] );
t_is( cb.P.prc(2,:), 56.9853*ones(1,3), 4, [t ' : load 2 P price'] );
t_is( co.P.prc(:,1), r.bus(Gbus, LAM_P), 8, [t ' : gen P prices'] );
t_is( cb.P.prc(:,1), r.bus(Lbus, LAM_P), 8, [t ' : load P prices'] );
t_is( co.Q.qty, [4.2722; 11.3723; 14.1472; 22.8939; 36.7886; 12.3375; 0; 0; 0], 2, [t ' : Q offer quantities'] );
t_is( co.Q.prc, [0;0;0;0;0;3; 0.4861; 2.5367; 1.3763], 4, [t ' : Q offer prices'] );
t_is( cb.Q.qty, [0;0;0;0;0;0; 15; 4.0785; 5], 2, [t ' : Q bid quantities'] );
t_is( cb.Q.prc, [0;0;0;0;0;3; 0.4861; 2.5367; 1.3763], 4, [t ' : Q bid prices'] );
t_is( co.Q.prc, r.bus([Gbus; Lbus], LAM_Q), 8, [t ' : Q offer prices'] );
t_is( cb.Q.prc, co.Q.prc, 8, [t ' : Q bid prices'] );
t = 'marginal Q offer, marginal PQ bid, auction_type = 0';
mkt.auction_type = 0;
[r, co, cb, f, dispatch, success, et] = runmarket(mpc, offers, bids, mkt, mpopt);
t_is( co.P.qty, co5.P.qty, 8, [t ' : gen P quantities'] );
t_is( cb.P.qty, cb5.P.qty, 8, [t ' : load P quantities'] );
t_is( co.P.prc, offers.P.prc, 8, [t ' : gen P prices'] );
t_is( cb.P.prc, bids.P.prc, 8, [t ' : load P prices'] );
t_is( co.Q.qty, co5.Q.qty, 8, [t ' : gen Q quantities'] );
t_is( cb.Q.qty, cb5.Q.qty, 8, [t ' : load Q quantities'] );
t_is( co.Q.prc, offers.Q.prc, 8, [t ' : gen Q prices'] );
t_is( cb.Q.prc, bids.Q.prc, 8, [t ' : load Q prices'] );
end
t_end;