qnom
Exact analysis of open, multiple-class BCMP networks. The network can be made of queueing centers (FCFS, LCFS-PR or PS) or delay centers (IS). This function assumes a network with service centers and customer classes.
INPUTS
lambda(c)
If this function is invoked as qnom(lambda, S, V, …)
,
then lambda(c)
is the external arrival rate of class
customers (lambda(c) ≥ 0
). If this
function is invoked as qnom(lambda, S, P, …)
, then
lambda(c,k)
is the external arrival rate of class
customers at center (lambda(c,k)
≥ 0
).
S(c,k)
mean service time of class customers on the service center
(S(c,k)>0
). For FCFS nodes, mean service
times must be class-independent.
V(c,k)
visit ratio of class customers to service center
(V(c,k) ≥ 0
). If you pass this argument,
class switching is not allowed
P(r,i,s,j)
probability that a class job completing service at center
is routed to center as a class job.
If you pass argument P, class switching is allowed;
however, all servers must be fixed-rate or infinite-server nodes
(m(k) ≤ 1
for all ).
m(k)
number of servers at center . If m(k) < 1
,
enter is a delay center (IS); otherwise it is a regular
queueing center with m(k)
servers. Default is
m(k) = 1
for all .
OUTPUTS
U(c,k)
If is a queueing center, then U(c,k)
is the
class utilization of center . If is an IS
node, then U(c,k)
is the class defined as X(c,k)*S(c,k)
.
R(c,k)
class response time at center . The system
response time for class requests can be computed as
dot(R, V, 2)
.
Q(c,k)
average number of class requests at center . The
average number of class requests in the system Qc
can be computed as Qc = sum(Q, 2)
X(c,k)
class throughput at center .
NOTES
If the function call specifies the visit ratios V, class switching is not allowed. If the function call specifies the routing probability matrix P, then class switching is allowed; however, all nodes are restricted to be fixed rate servers or delay centers: multiple-server and general load-dependent centers are not supported. Note that the meaning of parameter lambda is different from one case to the other (see below).
REFERENCES
See also: qnopen, qnos, qnomvisits
P = zeros(2,2,2,2); lambda = zeros(2,2); S = zeros(2,2); P(1,1,2,1) = P(1,2,2,1) = 0.2; P(1,1,2,2) = P(2,2,2,2) = 0.8; S(1,1) = S(1,2) = 0.1; S(2,1) = S(2,2) = 0.05; rr = 1:100; Xk = zeros(2,length(rr)); for r=rr lambda(1,1) = lambda(1,2) = 1/r; [U R Q X] = qnom(lambda,S,P); Xk(:,r) = sum(X,1)'; endfor plot(rr,Xk(1,:),";Server 1;","linewidth",2, ... rr,Xk(2,:),";Server 2;","linewidth",2); legend("boxoff"); xlabel("Class 1 interarrival time"); ylabel("Throughput"); |