Compute Balanced System Bounds for closed, multiclass networks with K service centers and C customer classes. Only single-server nodes are supported.
INPUTS
N(c)
number of class c requests in the system (vector of length C).
D(c, k)
class c service demand at center k (C \times K
matrix, D(c,k) ≥ 0
).
S(c, k)
mean service time of class c
requests at center k (C \times K matrix, S(c,k) ≥ 0
).
V(c,k)
average number of visits of class c
requests to center k (C \times K matrix, V(c,k) ≥ 0
).
OUTPUTS
Xl(c)
Xu(c)
Lower and upper class c throughput bounds (vector of length C).
Rl(c)
Ru(c)
Lower and upper class c response time bounds (vector of length C).
See also: qncsbsb.
The following code
S = [10 7 5 4; ... 5 2 4 6]; NN=20; Xl = Xu = Rl = Ru = Xmva = Rmva = zeros(NN,2); for n=1:NN N=[n,10]; [a b c d] = qncmbsb(N,S); Xl(n,:) = a; Xu(n,:) = b; Rl(n,:) = c; Ru(n,:) = d; [U R Q X] = qncmmva(N,S,ones(size(S))); Xmva(n,:) = X(:,1)'; Rmva(n,:) = sum(R,2)'; endfor subplot(2,2,1); plot(1:NN,Xl(:,1), 1:NN,Xu(:,1), 1:NN,Xmva(:,1),";MVA;", "linewidth", 2); ylim([0, 0.2]); title("Class 1 throughput"); legend("boxoff"); subplot(2,2,2); plot(1:NN,Xl(:,2), 1:NN,Xu(:,2), 1:NN,Xmva(:,2),";MVA;", "linewidth", 2); ylim([0, 0.2]); title("Class 2 throughput"); legend("boxoff"); subplot(2,2,3); plot(1:NN,Rl(:,1), 1:NN,Ru(:,1), 1:NN,Rmva(:,1),";MVA;", "linewidth", 2); ylim([0, 700]); title("Class 1 response time"); legend("location", "northwest"); legend("boxoff"); subplot(2,2,4); plot(1:NN,Rl(:,2), 1:NN,Ru(:,2), 1:NN,Rmva(:,2),";MVA;", "linewidth", 2); ylim([0, 700]); title("Class 2 response time"); legend("location", "northwest"); legend("boxoff");
Produces the following figure
Figure 1 |
---|
Package: queueing