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
|
## Copyright (C) 2012, 2016 Moreno Marzolla
##
## This file is part of the queueing toolbox.
##
## The queueing toolbox is free software: you can redistribute it and/or
## modify it under the terms of the GNU General Public License as
## published by the Free Software Foundation, either version 3 of the
## License, or (at your option) any later version.
##
## The queueing toolbox is distributed in the hope that it will be
## useful, but WITHOUT ANY WARRANTY; without even the implied warranty
## of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
## General Public License for more details.
##
## You should have received a copy of the GNU General Public License
## along with the queueing toolbox. If not, see <http://www.gnu.org/licenses/>
1; # not a function file
# figure (1, "visible", "off"); # do not show plot window
page_screen_output(0); # avoid output pagination
N = 100; # total population size
beta = linspace(0.1,0.9,18); # population mix for class 1
D = [12 14 23 20 80 31; ...
2 20 14 90 30 33 ];
V = ones(size(D));
X1 = X1 = XX = zeros(size(beta));
R1 = R2 = RR = zeros(size(beta));
for i=1:length(beta)
pop = [fix(beta(i)*N) N-fix(beta(i)*N)];
[U R Q X] = qncmmva( pop, D, V );
X1(i) = X(1,1) / V(1,1);
X2(i) = X(2,1) / V(2,1);
XX(i) = X1(i) + X2(i);
R1(i) = dot(R(1,:), V(1,:));
R2(i) = dot(R(2,:), V(2,:));
RR(i) = N / XX(i);
endfor
## Plot throughput and response times
## set(gcf,"paperorientation","landscape");
## papersize=[4 3] * 2; margin=[0 0];
## set(gcf,"papersize",papersize);
## set(gcf,"paperposition", [margin papersize-margin*2]);
subplot(2,1,1);
plot(beta, X1, "--;Class 1;", "linewidth", 2, ...
beta, X2, ":;Class 2;", "linewidth", 2, ...
beta, XX, "-;System;", "linewidth", 2 );
ylabel("Throughput");
legend("location", "south", "orientation", "horizontal");
legend("boxoff");
title("Throughput and response time vs population mix");
subplot(2,1,2);
plot(beta, R1, "--;Class 1;", "linewidth", 2, ...
beta, R2, ":;Class 2;", "linewidth", 2, ...
beta, RR, "-;System;", "linewidth", 2 );
ax = axis();
ax(3) = 0;
axis(ax);
legend("location", "south", "orientation", "horizontal"); legend("boxoff");
xlabel("Class 1 population mix");
ylabel("Response Time");
print("web.eps", "-deps2", "-mono", "-tight");
clf;
## Plot system power
## set(gcf,"paperorientation","landscape");
## papersize=[4 2.5] * 1.2; margin=[0 0];
## set(gcf,"papersize",papersize);
## set(gcf,"paperposition", [margin papersize-margin*2]);
plot(beta, X1./R1, "--;Class 1;", "linewidth", 2, ...
beta, X2./R2, ":;Class 2;", "linewidth", 2, ...
beta, XX./RR, "-;System;", "linewidth", 2);
legend("location","south", "orientation", "horizontal"); legend("boxoff");
xlabel("Class 1 population mix");
ylabel("Power");
title("Power as a function of the population mix");
print("power.eps", "-deps2", "-mono", "-tight");
|
