Unit 19
Unit 19
lntroductio~i
Objectives
19.9 Summary
19.10 Key Words
19.1 1 Answers to SAQs
19.1 INTRODUCTION
The theory of stochastic process is now finding increasing applications to problems in
physical sciences which deal with process involving some random element in their
structure. The tendency to use the results of stochastic theory in building up of models in
applied science and engineering has increased considerably during the past decade. In the
present unit, we present an account of the general classification of stochastic processes. A
stochastic process, or sometimes random process, is the counterpart to a deterministic
process (or deterministic system) in probability theory. Instead of dealing with only one
possible 'reality' of how the process might evolve under time (as is the case, for example,
for solutions of an ordinary differential equation), in a stochastic or random process there
is some indeterminacy in its future evolution described by probability distributions. This
means that even if the initial portion (or starting point) is known, there are many
possibilities the process might go to, but some paths are more probable and others less.
Objectives
After studying this unit, you should be able to
b
explain 'stochastic process,
differentiate between various stochastic processes,
describe Markov Chain, and
define queuing model.
Thus we may regard a stochastic process either as a collection of random variables X (t)
which depend on the parameter t or as a collection of the realizations of the process X (t).
The value oEX (t) may be one dimensional, two dimensional or n-dimensional.
(When several quantities say XI (t), Xz (t) . . .X,, (t) are required for the complete
description of the state of the system at a fixed parameter point t).
For every fixed t say t = t,, X (w, t) is a random variable X(tl) = X (w, ti), the distribution Stochastic Models
of which, we denote by :
For every finite set of t-values say tl, . t, the corresponding random variables
X (1,) . . . X (t,,) have the joint n-dimensional distribution function :
.
Fn(x1,x2,..., ~ , ~ ; t...
~ t.n t) =~P r [ X (t,)_<x,... X(fn)5x,,]
One can express many essential properties of a stochastic process as properties of the
corresponding family of finite dimensional distributions.
SAQ 1
(a) Define stochastic process. What is its importance?
said to be evolutionary.
Gaussian or Normal Process
In these processes, the joint distribution function is multivariate normal. In the
Markov Process
future behavior of the process depends upon only on the present state but not on
the past. In formal terms a process is said to be Markovian, if
P r [ a , X ( t ) < b ( X , , = x , , X12 = x 2,..., X , , , ] = P r [ a < X ) ( t ) l b I X , , = x , , ]
Discrete Time
Continuous Value Process
Continuous Time
' Discrete Value Process
Discrete Time
Discrete Value Process
Example 19.3
Consider a simple experiment that there are k boxes and an infinitely large number
of identical balls which are thrown at random one by one into the boxes. The balls
thrown being equally likely to go into any one of the boxes. If X, represents the
I
number of filled boxes after :I throws, then {X,,, tz 2 1) constitutes a stochastic
II process.
Example 19.4
Suppose that X (t) is the random variable which represents the number of
increasing calls in an interval (0,t) of duration t units. Then the family
( X (t), t E 77( constitutes a stochastic process where T = [0, co).
SAQ 2
Describe various types of stochastic processes.
P, (t) = Pr {N (t) = n)
This probability is a function of the time I. Since the only possible values of n are
n = 0, 1, ...
Thus {p,,(t)) represents the probability distribution of the r. v. N ( t ) for every value of I.
The family of random variables {N (t), t 2 0) is a stochastic process. Here the time t is
continuous, the state space of N (t) is discrete and integral valued. Such a process is also
called a counting process. One of the most important integral valued processes is Poisson
process, which serves as a mathematical model for a wide range of empirical phenomena
with remarkable accuracy.
We now sho\v that undcr certain conditions, N (t) follows Poisson distribution with mean Stocl~asticModels
[ Theorem
Under the postulates I, II and 111, N ( t ) follows Poisson distribution with mean h t,
i.e. is given by the Poisson law
Corollary 1
E (N (t)) = h t
and var {N(t)) =h t
This is the manner in which customers form a queue (first come first served, last
come first served, random selection priority selection).
Notations and Assumptions
(a) Service is provided on a first come first served (FIFO) basis.
(b) Customers arrive completely at random but a certain average rate.
(c) The queuing system is in a steady state.
Theses three assumptions are valid in many realistic queuing systems and will
serve to illustrate the use of queuing theory.
Assumption (a) says that a customer who arrives first will be served first regardless
of whether he joins a queue or not. Assumption (b) states that an arrival is equally
likely to occur at any time and is independent of the time that has elapsed since the
last arrival. This is equivalent to saying that the number of arrivals per unit of time
is a random variable with a Poisson distribution.
That is, if
X = Number of arrivals per unit time
Then f (x) = Pr ( X = x ) = e LC x = 0, 1, 2, . .
x!
Where h is the average number of amvals per unit time? Other interesting result of
assumption (b) is that the time between consecutive arrivals T (also called the inter
arrival time) has an exponential distribution with the same parameter h.
Thus if T = time between consecutive arrivals
1
E ( t )= -
h
O p e r a t i u ~ ~Research
s In summary if the number of arrivals per unit of time has a Poisson diqlribution
Applications
with mean h , then the time between consecutive arrivals has an exponential
1
distribution with mean - . The queuing system in this situation is said to have a
h
Poisson input and customers are said to arrive according to a Poisson process.
Assuniption (c) means the queuing system has been operating long enough to be
independent of the initial state of the system and is independent of time. That Is the
system has reached a state of equilibrium with respect to time. The distribution of
the number of arrivals per unit of time and the distribution of the service time do
not change with tinie.
Suppose we let
Sn = Number of custoniers in the systeni,
PI,(t) = Probability of n customers in the system at tinie t.
A,, = Average arrival rate when n custoniers are in the systeni (both
waiting and being served), and
p,, = Average service rate when n custonlers are in the system.
The system being in the steady state does not imply that the arrival rate and service
rate are independent of the number of customers in the system. For finite queue.
Limited source and multiple servers' model A,, and p,, will be functions of the
number of customers in the system. In steady state we will use the notation
P,, = probability of n customers in the system at any point of t ~ m e ,
Queuing niodels with Poisson input-exponential service.
Here we assume
(a) Arrivals to tlle system occur con~pletelyat random.
(b) Arrivals form a single queue.
(c) First in first out discipline (FIFO).
(d) Departure from the system occur completely at random.
I The system can be in state Sf,at time t + A t in exactly four distinct ways :
(a) System is in state S,, at time t and no arrival or departure occur during
the interval A t.
(b) System is in state S,, - at time t and one arrival and no departure occur
during the interval A t.
(c) System is in state S, + at time t and no arrival and one departure occur
during the interval A t.
(d) One arrival and one departure or more that one arrival and/or
departure occur in the interval A t.
Since these four ways are al nii~tuallyexclusive, the probability of n custoniei-s in Stort~asticModels
the system at time t + A t with n > 0 may be expressed as the sum of three
independent compound probabilities (the probability of way 4 above is assumed to
be negligible)
el (t + A t ) = el (t) (1 - A l r At) (1 - p I r At) + Pn-, (t) ~1,-1 At (1 - pn-1 At)
P, (t + At) - 5
1 (1)
At
- e,-1(1) A,,-, - e, (1) (An + P,,) + el+l(1) PI?+,+ ta-mS in At
Taking limit as A t tends to zero.
We also assume here liniit P,,(t) = P,,as t tend to infinity.
Then the rate of change of P,,with respect to time is zero.
Thus, in steady state the differential equations become the difference equations
A,,-, e,-l + P I , + , e1+1 -A,, + P n ) P n =O n>O
- A , Po + t i l 4 =O 11=O
P2 =-A" A1 p,
PI P2
p? =
P1 P2 P3
PI P2 113 . . . PI1
Using the fact that
m
rr
C
=0
el=l
The parameter h, then is the average arrival rate per unit of time. Also let
-
T Time to service a customer.
Then g ( t ) = ~ l e - ~ t' t>O; p > 0
1
E (t) = -
Cl
We already know that if the number of amvals per unit of time has a Poisson
distribution with parameter h, then the inter arrival time the time between
consecutive arrivals) has an exponential distribution with the same parameter.
Likewise if the service time has an exponential distribution with parameter p
(completely random departure) then the number of custon~ersserviced per unit of
time has a Poisson distribution with parameter p.
Thus if X, = Number of customers serviced per unit of time.
S
Then F(X,)=~-~P~---! x S = 0 , 1 , 2,...
XY
.
E (Xs) = P
One is generally interested in modeling a given queuing system and then in
studying certain interesting characteristics to determine if measures should be
taken to modify the queuing system. Modification might take the form of adding
additional servers to the service facility, providing more space in the queuing area,
or buying a new piece of equipment to increase the service rate. Most decisions are
based on several useful quantities which can be obtained through the use of
equation for Po and P,,, namely
h
L = -= average number of customers in the system
P-A
h2
Lq = = average number of customers in the queue
P (Cl - A)
h
L, =- = average number of customers in nonempty queue
P-A
I
w=- = average time a custonlers spends in the system
P-1
h
TI = C l (11 - A) = average time a customers spends in the queue
1
w =- = average time a customers spends in the queue if he must wait
" Cl-h
'i ).'I/
-P\I
P(T>t)=e = probability the time in the system is greater than t
Example 19.5
A person repairing mobiles finds that the time spent on the mobile set has
exponential distribution which niean 20 minutes. If the mobiles are repaired in the
order in which they came in and their arrival is approximately Poisson with an
average rate of 15 for 8-hour day. What is the repaimlan's expected idle time each
day? How many jobs are ahead of the average set just brought in?
Solution
15 1
Arrival rate h = -- = - unitslmin
8 x 6 0 32
1
Service rate p = - unitslmin
20
h 5
Exponential No. of jobs in the system = L, = -- -
p-h 3
No. of hours for which the repairman remains busy in an
8h
8-hour day = - = 5 hours
P
:. Repairman's expected idle time each day = 8 - 5 = 3 hours.
SAQ 5
(a) What effect does doubling h and 11have on L, L,/, Wand W,in a single
channel queuing model.
(b) In a large maintenance department, fitters draw parts from the parts-stores
which is at present staffed by one storeman. The maintenance foreman is
concerned about the time spent by fitters getting parts and wants to know if
the employment of a stores labourer to assist the storeman would be
worthwhile on investigation it is found that
(i) a simple queue situation exists,
(ii) fitters cost Rs. 2.50 per hour,
(iii) the storeman costs Rs. 2 per hour and can deal, on the average, with
10 fitters per hour,
(iv) a labourer could be employed at Rs. 1.75 per hour and would increase
the service capacity of the stores to 12 per hour, and
(v) on the average 8 fitters visit the stores each hour.
A library wants to improve its service Eacilities in terms of the waiting time of its
borrowers. 'Tl~elibrary has to counters at present and borrowers arrive according to
Poisson distribution with arrival rate 1 every 6 minutes and service time follows
exponential distribution with a meall of I0 minutes. The library I ~ a srelaxed its
men~bersl~ip rules and a substantial increase in the nun~berof borrowers is
expected. Find the number of additional counters to be provided if the arrival rate
is expected to be twice the present value and the average waiting time of the
borrowers must be limited to half the present value.
Solution
1
h = - minutes = IOIhour
6
. . New arrival
1'= 20ihour
19.9 SUMMARY
Stochastic means being or having a random variable. A stochastic model is a tool for
estimating probability distribution of potential outcomes by allowing for random
variation in one or niore inputs over time. In this unit, you have studies about types of
stochastic niodels, specification of stochastic process and applications of stochastic
models.
-
19.10 KEY WORDS
Stochastic Model : A stochastic model is one that involves probability
or randomness.
Markovian Process : A Markovian process is a time-based process
where the probability of a state on the next trial
(time period) depends only on the state in the
current trial (time period) and not on how we got
to the present state.
Markovian Chain : A Markovian chain is a stochastic process with the
following properties :
(a) Discrete state space
(b) Markovian property
(c) One-step transition probability that remain
constant over tinie (termed stationary
transition probabilities).
Please refer the preceding text for all the Answers to SAQs.
FURTHER READINGS
Cooper, L. and Steinberg D., (1974), Methods arid Applications of Linear Prograrnnti~zg.
Saunders, Philadelphia, USA.
Mustafi, C. K., (1988), Operations Research Methods and Practice, Wiley Eastern
Limited, New Delhi.
Sasieni, M. A. Yaspan and L. Friedman, (1 9591, Operations Research Methods and
Problems, J . Wiley and Sons, New York, USA.
Hadly, G . , (1969), Linear Programming, Addison Wesley Reading Massachusetts, USA.
Mittal, K. V., (1976), Optimisation Methods and Operation Research and System
Analysis. Wiley Eastern Limited, New Delhi.