OR Merged
OR Merged
LINEAR PROGRAMMING -1
Formulation &
Graphical Solution
12/1/2020 1
LINEAR PROGRAMMING I
Examples of modelling business
problems for linear programming
problems (LPP)
Basic Requirements for applying LPP
Basic assumptions of LPP
Graphical solution of LPP for two
variables
Technical issues involved in linear
programming
Applications of LPP
12/1/2020 2
LEARNING OBJECTIVES
After this session, you should be able :
◦ To formulate linear programming models for
resource allocation problems
◦ To understand the requirements of a linear
programming problem
◦ To solve a linear programming maximization
problem involving two variables graphically
◦ To solve a linear programming minimization
problem involving two variables graphically
◦ To appreciate the technical issues involved in
linear programming
LINEAR PROGRAMMING
A problem solving approach developed to
help managers make decisions.
TS1-10
REQUIREMENTS FOR APPLYING LPP
1. An objective which should be clearly
identifiable and measurable in quantitative terms
such as maximization of sales, maximization of
profits, minimization of costs etc.
TS1-19
LP Model Formulation: Example- 1
Product mix problem
Beaver Creek Pottery Company
How many bowls and mugs should be produced to
maximize profits given labor and materials
constraints?
Product resource requirements and unit profit:
Resource Requirements
Bowl 1 4 40
Mug 2 3 50
Product mix problem
SOLUTION (Product mix problem)
Resource 40 hrs of labor per day
Availability: 120 lbs of clay
Decision x1 = number of bowls to produce per day
Variables: x2 = number of mugs to produce per day
Objective Maximize Z = $40x1 + $50x2
Function: Where Z = profit per day
Resource 1x1 + 2x2 40 hours of labor
Constraints: 4x1 + 3x2 120 pounds of clay
Non-Negativity x1 0; x2 0
Constraints:
SOLUTION (Product mix problem)
Complete Linear Programming Model:
x1, x2 0
EXAMPLE II (Media Selection)
An advertising company is planning a media
campaign for a client, willing to spend Rs. 20
lakhs to promote a new fuel economy model
of a pressure cooker.
Objective Function
Objective is to maximize total number of
housewives (z) to be exposed to the
advertisements
Maximize z = 1,00,000x1 + 5,00,000x2 + 10,00,000x3
SOLUTION (Media Selection)
Constraints
Maximum total expense on advertising not to
exceed Rs.20,00,000
» 20,000x1+40,000x2+1,00,000x3≤ 20,00,000
» x1 + 2x2 + 5x3 ≤ 100
At least 50,00,000 housewives to be exposed
to TV advertising
» 10,00,000x3 ≥ 50,00,000
» x3 ≥ 5
subject to constraints
x1 ≥ 0, x2 ≥ 0, x3 ≥ 0 (Non-negativity constraints)
Example III (Diet Mix)
A kitchen manager at Calcutta Medical Collage Hospital
has to decide food mix for patients.
Dietary instructions are that each patient must get at
least: i) One gm of protein ii) one gm of fat iii) three gm
of carbohydrates
Additional instructions are that in no case
carbohydrates should exceed 6 gm per patient.
Availability of protein, fat and carbohydrates in gms per
kg of chicken, rice and bread is given as follows:
Protein Fat Carbohydrates Price/Kg
Chicken 10 2 0 160
Rice 2 1 15 30
Bread 2 0 10 20
Formulate a suitable mathematical model for the above
diet mix assuming 100 patients on that day.
SOLUTION (Diet Mix)
Decision Variables
Problem is to find the quantity of chicken,
rice and bread to be used for diet of 100
patients.
x1 : quantity (kg) of chicken to be used,
x2 : quantity (kg) of rice to be used and
x3 : quantity (kg) of bread to be used
Objective Function
Objective is to minimize the cost of diet mix
Minimize z = 160x1 + 30x2 + 20x3
SOLUTION (Diet Mix)
Constraints
Food to contain at least 100gm (= 1gm per
patient Х 100 patients) of protein
» 10x1 + 2x2 + 2x3 ≥ 100
subject to constraints
10x1+2x2+2x3 ≥ 100 (Minimum Protein
requirements)
2x1+ x2 ≥ 100 (Minimum Fat requirements)
15x2+10x3 ≥ 300 (Minimum Carbohydrate
requirements)
15x2+10x3 ≤ 600 (Maximum limit of
Carbohydrates)
x1 ≥ 0, x2 ≥ 0, x3 ≥ 0 (Non-negativity constraints)
Example IV (By-product Generation)
In a chemical Industry, two products A & B are
made involving two operations.
subject to constraints
3x1+4x2 ≤18 (Capacity constraint of Operation I)
4x1+5x2 ≤ 21 (Capacity constraint of Operation II)
x3 ≤ 5 (Sale of by-product C)
3x2– x3– x4=0 (Generation of by-product C)
Objective Function
◦ Objective is to minimize manpower
requirements
◦ i.e. no. of waiters for all six periods
together must be minimum possible i.e.
Minimize z = x1 + x2 + x3 + x4 + x5 + x6
SOLUTION (Staffing Problem)
Constraints
◦ At any 4 hrs period, no. of waiters added
& continuing from previous shift should
be at least equal to their requirement
» x6 + x1 ≥ 2
x1 + x2 ≥ 3
x2 + x3 ≥ 5
x3 + x4 ≥ 4
x4 + x5 ≥ 7
x5 + x6 ≥ 3
X1, ………. , X6 ≥ 0 and have integer values
Financial Model: Rate of Return
An investor has $100,000 to invest in the stock
market and is considering investments in six stocks
with the following current projected rates of return:
Tech Stocks: Cisco(8%), Microsoft(6%) & Intel(5%)
Bank Stocks: Bank of America(7%), First Bank (4%)
Money Market: ING(2%)
Constraints
At least $20,000 is to be invested in ING
At least 25% of funds in tech stocks must be in Intel
At least as much in bank stocks as tech stocks
No more than $25,000 in investments with projected
returns of 5% or less
◦ Formulate a suitable LP model
CASE – Regional Planning
Southern confederation of Kibbutzim is a group of
3 Kibbutzim (communal farming communities) in
Israel.
LINEAR PROGRAMMING -1
Graphical Solution
12/4/2020 1
LINEAR PROGRAMMING I
Special
types of linear programming
problems
Applications of LPP
12/4/2020 2
LEARNING OBJECTIVES
After this session, you should be able :
◦ To solve a linear programming maximization
problem involving two variables graphically
Bowl 1 4 40
Mug 2 3 50
x1, x2 0
Graphical Solution of Maximization
Model: Coordinate Axes
X2: mugs
Maximize
Z = $40x1 + $50x2
subject to:
1x1 + 2x2 40
4x2 + 3x2 120
x1, x2 0
subject to:
1x1 + 2x2 40
4x2 + 3x2 120
x1, x2 0
Maximize
Z = $40x1 + $50x2
subject to:
1x1 + 2x2 40
4x1+ 3x2 120
x1, x2 0
Maximize
Z = $40x1 + $50x2
subject to:
1x1 + 2x2 40
4x1 + 3x2 120
x1, x2 0
Maximize
Z = $40x1 + $50x2
subject to:
1x1 + 2x2 40
4x2 + 3x2 120
x1, x2 0
Maximize
Z = $40x1 + $50x2
subject to:
1x1 + 2x2 40
4x2 + 3x2 120
x1, x2 0
Maximize
Z = $40x1 + $50x2
subject to:
1x1 + 2x2 40
4x2 + 3x2 120
x1, x2 0
subject to:
1x1 + 2x2 40
4x2 + 3x2 120
x1, x2 0
Maximize
Z = $40x1 + $50x2
subject to:
1x1 + 2x2 40
4x2 + 3x2 120
x1, x2 0
Maximize
Z = $40x1 + $50x2
subject to:
1x1 + 2x2 40
4x2 + 3x2 120
x1, x2 0
Maximize
Z = $40x1 + $50x2
subject to:
1x1 + 2x2 40
4x2 + 3x2 120
x1, x2 0
subject to:
1x1 + 2x2 40
4x2 + 3x2 120
x1, x2 0
Set up equations No
of objective function Is the objective
function maximized?
Identify nearest vertex
Set up equations Yes among vertices in feasible
of constraints
region in relation to
Identify farthest vertex among
objective function line
vertices in feasible region in
Write complete relation to objective function line
LPP formulation
Fertilizing
farmer’s field
Chemical Contribution
Nitrogen Phosphate
Brand
(lb/bag) (lb/bag)
Super-gro 2 4
Crop-quick 4 3
LP Model Formulation – Minimization
Decision Variables:
x1 = bags of Super-gro
x2 = bags of Crop-quick
Model Constraints:
2x1 + 4x2 16 (nitrogen constraint)
4x1 + 3x2 24 (phosphate constraint)
x1, x2 0 (non-negativity constraint)
Constraint Graph – Minimization
Minimize
Z = $6x1 + $3x2
subject to:
2x1 + 4x2 16
4x1 + 3x2 24
x1, x2 0
Minimize
Z = $6x1 + $3x2
subject to:
2x1 + 4x2 16
4x2 + 3x2 24
x1, x2 0
Minimize
Z = $6x1 + $3x2
subject to:
2x1 + 4x2 16
4x2 + 3x2 24
x1, x2 0
◦ Subject to constraints
6x + 4y ≤ 3600
2x + 4y ≤ 2000
3.6x + 4y ≤ 3600
0 ≤ x ≤ 500 0 ≤ y ≤ 400
x2 SOLUTION TO EXAMPLE
Feasible solution region
1000
OABCDE shows every straight
900
(0, 900) line in the graph except (3), i.e.
every constraint in this problem
4
800 except constraint 3.6x+4y≤3600
defines a certain portion of
700 boundary of feasible solution
region and hence it is a
600
redundant constraint.
(0, 500)
500
5
(0, 400) B
400
A
C
300
200
D
100
(600,0) (1000,0)
E
(500,0)
0 100 200 300
400 500 600 700 800 900 1000 2
1 3
SOLUTION TO EXAMPLE
Co-ordinates of extreme points of feasible region are
O(0,0), A(0,400), B(200,400), C(400,300), D(500,150),
E(500,0)
Step 4: Z values corresponding to extreme points are
Extreme point (x,y) Z=25x+20y
O (0,0) 0
A (0,400) 8000
B (200,400) 13000
C (400,300) 16000
D (500,150) 15500
E (500,0) 12500
•Max value of Z occurs at extreme point C(400,300) »
company should produce 400 & 300 nos. of chairs of type
A and type B respectively.
REDUNDANT CONSTRAINTS
•Max value of Z occurs at extreme point C(400,300) »
company should produce 400 & 300 nos. of chairs of
type A and type B respectively.
Maximize
Z=$40x1 + 30x2
subject to:
1x1 + 2x2 40
4x2 + 3x2 120
x1, x2 0
x1 = number of bowls
x2 = number of mugs Multiple Optimal Solutions
CONDITIONS FOR EXISTENCE OF MULTIPLE
OPTIMAL SOLUTION
Unbounded Problem
UNBOUNDED SOULTION
A LPP is unbounded if value of solution is made
infinitely large in a maximization LPP or
infinitely small in minimization LPP
without violating any of the constraints.
12/9/2020
LINEAR PROGRAMMING II
2
Subject to
gi(x1, x2 , …, xn)
{ }
=
bi , ,, i =1,2,…,m
xj ≥ 0, j = 1,…,n
,…,m
MODEL COMPONENTS
{}
• •gi (x1, x2 ,…, xn)
=
bi
: structural or functional or technological
constraints
• xj ≥ 0 : non-negativity constraints
• f(x1, x2 ,…, xn) :objective function
x1
.
A feasible solution x = .
xn )
satisfies all the
constraints (both structural and non negativity)
Objective function ranks feasible solutions
LP – A MATHEMATICAL PROGRAM
A linear program is a special case of a
mathematical program where f and g1 ,…, gm are
linear functions
Linear Program:
Maximize/Minimize z = c1x1 + c2x2 + • • • + cnxn
{}
Subject to ≥
ai1x1 + ai2x2 + • • • + ainxn ≤
bi m ,
,i = 1,…,m =
xj ≥ 0, j = 1,…,n
LP MODEL COMPONENTS
goal is to find x
LP model : Canonical form
xi >= 0 , (i =1,2,…,n)
STANDARD FORM OF LPP
,,i = 1,…,m
x1, x2 0
Slack Variables
Max
Z=40x1 + 50x2 +0. s1 +
0.s2
subject to:
1x1 + 2x2 + s1 = 40
4x2 + 3x2 + s2 = 120
x1, x2, s1, s2 0
Where:
x1 = number of bowls
x2 = number of mugs
s1, s2 are slack
variables
Solution Points A, B, and C with Slack
LP Model – FERTILIZER EXAMPLE
Decision Variables:
x1 = bags of Super-gro
x2 = bags of Crop-quick
Model Constraints:
2x1 + 4x2 16 (nitrogen constraint)
4x1 + 3x2 24 (phosphate constraint)
x1, x2 0 (non-negativity constraint)
Surplus Variables
Model Constraints:
2x1 + 4x2 16 (nitrogen constraint)
4x1 + 3x2 24 (phosphate constraint)
x1, x2 0 (non-negativity constraint)
Minimize
Z = $6x1 + $3x2 + 0s1 +
0s2
subject to:
2x1 + 4x2 – s1 = 16
4x2 + 3x2 – s2 = 24
x1, x2, s1, s2 0
Graph of Fertilizer Example
LP PROBLEM 1
LPP model :
Max Z = 40x1 + 35x2 (Profit)
subject to constraints
2x1 + 3x2 ≤ 60 (Raw material constraint)
4x1 + 3x2 ≤ 96 (Labour hour constraint)
x1 ≥ 0, x2 ≥ 0 (Non-negativity restrictions)
STEP 1: STANDARDIZATION OF PROBLEM
Standard form of LPP (introducing slack
variables x3 & x4):
Max Z = 40x1 + 35x2 + 0.x3 + 0.x4
subject to constraints
2x1 + 3x2 + x3 + 0.x4 = 60
4x1 + 3x2 + 0.x3 + x4 = 96
xj ≥ 0, j = 1, 2, 3, 4
Initial basic feasible solution:
x1 = 0, x2 = 0, x3 = 60, x4 = 96
Value of objective function:
Z=40.0 + 35.0 + 0.60 + 0.96= 0
IBFS is not to produce any of the products A and
B with a zero profit.
JUSTIFICATION & SIGNIFICANCE
OF SLACK VARIABLES
2x1 + 3x2 = 0 » x3 = 60
& 2x1 + 3x2 = 60 » x3 = 0
Assumption:
Unused resources have no cost and do not affect
profit.
Coefficients assigned is zero.
SESSION 5
1
12/14/2020
LP PROBLEM 1
LPP model :
Max Z = 40x1 + 35x2 (Profit)
subject to constraints
2x1 + 3x2 ≤ 60 (Raw material constraint)
4x1 + 3x2 ≤ 96 (Labour hour constraint)
x1 ≥ 0, x2 ≥ 0 (Non-negativity restrictions)
STEP 1: STANDARDIZATION OF PROBLEM
Standard form of LPP (introducing slack
variables x3 & x4):
Max Z = 40x1 + 35x2 + 0.x3 + 0.x4
subject to constraints
2x1 + 3x2 + x3 + 0.x4 = 60
4x1 + 3x2 + 0.x3 + x4 = 96
xj ≥ 0, j = 1, 2, 3, 4
Initial basic feasible solution:
x1 = 0, x2 = 0, x3 = 60, x4 = 96
Value of objective function:
Z=40.0 + 35.0 + 0.60 + 0.96= 0
IBFS is not to produce any of the products A and
B with a zero profit.
JUSTIFICATION & SIGNIFICANCE
OF SLACK VARIABLES
2x1 + 3x2 = 0 » x3 = 60
& 2x1 + 3x2 = 60 » x3 = 0
Assumption:
Unused resources have no cost and do not affect
profit.
Coefficients assigned is zero.
STEP 2: SIMPLEX TABLE 1
(NON-OPTIMAL SOLUTION)
Step 2:
Simplex Table 1 (Non-optimal solution)
Cj 40 35 0 0 Replacement
Basis Ratio
CB XB Xj X1 X2 X3 X XBi/air, air > 0
4
C3=0 X3=60 2 3 1 0 60/2 = 30
Zj 40 30 0 10
Δj =Cj-ZJ 0 5 0 -10
STEP 3: SIMPLEX TABLE 2
(NON-OPTIMAL SOLUTION)
To obtain improved iteration (Simplex Table 2),
following operations are used:
Max Z = 35 X 8 + 40 X 18 = 1000
JUSTIFICATION AND SIGNIFICANCES
XB values:
In Table 1, containing initial Simplex Table 1 (Non-optimal solution)x
solution, Product A (x1)
yielded a profit of Rs. 40/unit Cj 40 35 0 0 Repla
Basis ceme
while product B(x2) would nt
contribute Rs. 35/unit & Ratio
hence x1 is added to the new Xj X1 X2 X3 X4 XBi/air,
solution. CB XB air > 0
With focus on product A, the
quantity to be added could C3=0 X3=60 2 3 1 0 60/2=
be determined as: 30
restriction
C3=0 X3=12 0 3/2 1 -1/2 12/(3/2)
• Production of 24 units of =8
product A will exhaust all the
labour hours, while this C1=40 X1=24 1 3/4 0 1/4 24/(3/4)
= 32
output would require 24x2 =
48 kg of raw material leaving Zj 40 30 0 10
60 – 48 = 12 kg unused.
Δj =Cj-ZJ 0 5 0 -10
• XB values X3 = 12 and X1 =
24 in Simplex table 2 are
indicative of this
SUBSTITUTION RATES (Column headed X4)
simplex table 2
Since one unit of product A requires 4 hours of
labour, the solution in table 2 indicates complete
use of available 96 labour hours.
For x3, net change is zero, while for x4, adding a unit
to the solution would involve a net loss of Rs. 10.
LINEAR PROGRAMMING II
12/17/2020
LINEAR PROGRAMMING II
2
12/17/2020
LEARNING OBJECTIVES
After this session, you should be able :
LP formulation:
Minimize Z = 40x1 + 24x2 (Total Cost)
subject to
subject to constraints
20x1 + 50x2 - x3 + 0.x4 + A1 + 0.A2 = 4800
80x1 + 50x2 + 0.x3 - x4 + 0.A1 + A2 = 7200
xj ≥ 0, j = 1, 2, 3, A1 , A2 ≥ 0
Δj =Cj-ZJ 40 - 24 - M M 0 0
100M 100M
SOLUTION TO MINIMIZATION CASE
Step 3:
Simplex Table 2 (Non-optimal solution)
Cj 40 24 0 0 M Replacem
Basis ent Ratio
CB XB Xj X1 X2 X3 X4 A2 XBi/air,
air > 0
C2=24 X2 = 96 2/5 1 -1/50 0 0 96/(2/5)
= 240
C6=M A2=2400 60 0 +1 -1 1 2400/60
= 40
Zj 60M+ 24 M- -M M
48/5 12/25
Δj =Cj-ZJ 152/5 0 12/25 - M 0
– 60M M
SOLUTION TO MINIMIZATION CASE
Step 4:
Simplex Table 3 (Non-optimal solution)
Cj 40 24 0 0 Replacem
Basis ent Ratio
CB X Xj X1 X2 X3 X4 XBi/air,
B air > 0
C2=24 X2 = 80 0 1 -2/75 1/150
C1=40 X1=40 1 0 1/60 -1/60 40/(1/60)
=2400
Zj 40 24 2/75 -38/75
Δj =Cj-ZJ 0 0 -2/75 38/75
SOLUTION TO MINIMIZATION CASE
:optimal table
Cj 40 24 0 0 Replacemen
Basis t Ratio
CB XB Xj X1 X2 X3 X4 XBi/air, air > 0
subject to constraints
x1 + x2 ≤ 200
x2 ≤ 125
3x1 + 6x2 ≤ 900
x1 ≥ 0, x2 ≥ 0 (Non-negativity restrictions)
MULTIPLE OPTIMAL SOLUTIONS
OPTIMAL SIMPLEX TABLE
Cj 8 16 0 0 0
Basis
Xj X1 X2 X3 X4 X5
CB XB
C4=0 X3 = 25 0 0 1 1 -1/3
C2=16 X2 = 125 0 1 0 1 0
Zj 8 16 0 0 8/3
Δj =Cj-ZJ 0 0 0 0 -8/3
MULTIPLE OPTIMAL SOLUTIONS
Optimal solution is X1=50, X2=125, Zmax= 2400
Xj X1 X2 X3 X4 X5
CB XB
C5=0 X4 = 25 0 0 1 1 -1/3
Zj 8 16 0 0 8/3
Δj =Cj-ZJ 0 0 0 0 -8/3
MULTIPLE OPTIMAL SOLUTIONS
subject to constraints
2x1 + x2 ≤ 40
4x1 - x2 ≤ 20
x1 ≥30
x1 ≥ 0, x2 ≥ 0 (Non-negativity restrictions)
INFEASIBLE SOLUTIONS
(OPTIMAL SIMPLEX TABLE)
Cj 20 30 0 0 0 -M
Basis
CB XB Xj X1 X2 X3 X4 X5 A1
Maximize Z = 2x1 + x2
subject to
x1 - x2 ≤ 10
2x1 ≤ 40
x1≥0, x2≥0, (Non-negativity constraints)
UNBOUNDED SOLUTIONS
(NON OPTIMAL SIMPLEX TABLE-1)
Cj 2 1 0 0
Basis
Replacement
Xj X1 X2 X3 X4 Ratio
CB XB XBi/air, air > 0
C3 =0 X3= 10 1* -1 1 0 10
C4=0 X4= 40 2 0 0 1 20
Zj 0 0 0 0
Δj =Cj-ZJ 2 1 0 0
UNBOUNDED SOLUTIONS
Cj 2 1 0 0
Basis Replacement
CB XB X j X1 X2 X3 X4 Ratio
XBi/air, air > 0
C1 =2 X1= 10 1 -1 1 0
C4=0 X4= 20 0 2* -2 1 10
Zj 2 -2 2 0
Δj =Cj-ZJ 0 3 -2 0
UNBOUNDED SOLUTIONS
(NON OPTIMAL SIMPLEX TABLE-3)
Cj 2 1 0 0
Basis Replacement Ratio
CB XB Xj X1 X2 X3 X4 X Bi/ air, air > 0
subject to constraints
6x1 + 3x2 ≤ 18
3x1 + x2 ≤ 8
4x1 + 5x2 ≤ 30
x1 ≥ 0, x2 ≥ 0 (Non-negativity restrictions)
DEGENERACY
IBFS: x3 = 18, x4 = 8, x5 = 30
Non-basic variables: x1 = 0, x2 = 0
DEGENERACY
Step 2:
Simplex Table 1 (Non-optimal solution)
Cj 28 30 0 0 0 Replacement
Basis Ratio
CB XB j X X 1 X 2 X3 X 4 X XBi/air, air > 0
5
C3=0 X3 =18 6 3 1 0 0 18/3 = 6
C5=0 X5 = 30 4 5 0 0 1 30/5 = 6
Zj 0 0 0 0 0
Δj =Cj-ZJ -28 -30 0 0 0
DEGENERACY
Step 3:
Simplex Table 2 (Optimal solution)
Cj 28 30 0 0 0 Replacement
Basis Ratio
XBi/air, air > 0
CB XB Xj X1 X2 X3 X4 X5
C2=30 X2 = 6 2 1 1/3 0 0
C4=0 X4 = 2 1 0 -1/3 1 0
C5=0 X5 = 0 -6 0 -5/3 0 1
Zj 60 30 10 0 0
Δj =Cj-ZJ -32 0 -10 0 0
DEGENERACY
LINEAR PROGRAMMING
III
1 12/28/2020
LINEAR PROGRAMMING III
2
Formulation of dual
Duality and simplex method
subject to constraints
ai1x1+ ai2x2+ ………. +ainxn = bi , i = 1,2,…….,m
j = 1,2,…….,n
xj ≥ 0
DUAL PROBLEM (MINIMIZATION)
6
Minimize Z* = b1w 1 + b2w 2 + ………… + bmw m
subject to constraints
a1jw1+ a2jw2+ ………. +amjwm ≥ cj,
wi unrestricted, i = 1,2,…….,m
j = 1,2,…….,n
subject to constraints
ai1x1+ ai2x2+ ………. +ainxn = bi ,
i = 1,2,…….,m
j = 1,2,…….,n
xj ≥ 0
DUAL PROBLEM (MAXIMIZATION)
8
Maximize Z* = b1w 1 + b2w 2 + ………… + bmw m
subject to constraints
a1jw1+ a2jw2+ ………. +amjwm ≤ cj ,
wi unrestricted i = 1,2,…….,m
j = 1,2,…….,n
subject to constraints
3x1 – x2 + 2x3 ≥ 7
2x1 – 4x2 ≥ 12
– 4x1 +3x2 + 8x3 = 10
x1, x2 ≥ 0 , x3 is unrestricted.
SOLUTION TO PROBLEM 2
16
Step 1: Writing x3 = xt – xs where xt ≥ 0 , xs ≥ 0 we get
Minimize Z = x1 – 3x2 – 2xt + 2xs
Subject to constraints
3x1 - x2 + 2xt - 2xs ≥ 7
2x1 – 4x2 + 0.xt + 0.xs ≥ 12
– 4x1 + 3x2 + 8xt – 8xs = 10
xj ≥ 0, j = 1, 2, t, s
Subject to constraints
3w1 + 2w2 - 4w3 ≤ 1
w1 + 4w2 – 3w3 ≥ 3
- 2w1 – 8w3 = 2
w1 ≥ 0, w2 ≥ 0 & w3 unrestricted
DUALITY AND SIMPLEX METHOD
19
Since any LPP can be solved by using simplex
method, the method is applicable to both the
primal as well as to its dual.
Standard form :
Max Z = 25x1 + 20x2 + 0.x3 + 0.x4
Optimal Solution:
Producing 18 units of A & 8 units of B per week
would yield maximum profit of Rs.1000.
28
OPTIMAL SOLUTION)
CB XB Xj X1 X2 X3 X4
C2=35 X2=8 0 1 2/3 -1/3
Zj 40 35 10/3 25/3
Δj =Cj-ZJ 0 0 -10/3 -25/3
ECONOMIC INTERPRETATION OF DUAL
Let the firm be approached by another firm who
would like to rent facilities of firm for one week.
2y1 + 4y2 ≥ 40
3y1 + 3y2 ≥ 35
Also, rentals can not be negative. 31
ECONOMIC INTERPRETATION OF DUAL
COMPLETE FORM OF PROBLEM:
Minimize G = 60y1 + 96y2
subject to 2y1 + 4y2 ≥ 40
3y1 + 3y2 ≥ 35, y1, y2 ≥ 0
Same as dual to the given problem.
Rental rates y1 and y2 obtained from solutions of dual:
y1 = 10/3 y2 = 25/3
Can also be obtained from Δj row of optimal primal
simplex table.
Objective function value of primal and dual would be
identical.
i.e. Minimum total rental acceptable by farm is equal
to maximum profit that it can earn by producing
output itself using given resources. 32
ECONOMIC INTERPRETATION OF DUAL
Shadow/ Imputed Price: Individual rental rates y1 &
y2 (dual variables corresponding to slack variables)
Indicates unit worth of a resource- change in optimal
objective value per unit change in availability of the
resource .
Prices assigned to services of two resources, raw
materials and labour hours.
Imputed from profit obtained from utilizing their
services.
Bear no relationship with original costs of the two.
So long as the total revenue from all the activities is
less than the worth of resources , corresponding
primal & dual solutions are not optimal.
Optimality (maximum revenue) is reached only when
the resources have been exploited completely which
can happen only when the input (worth of
resources) equals the output (revenue).
33
ECONOMIC INTERPRETATION OF DUAL
Each unit of product A contributes Rs40 to the profit.
Imputed prices of material and labour Rs10/3 per kg and
Rs25/3 per hour respectively.
P1 P2 P3 P4
R1 4 4 3 7 90 Kgs
R2 6 3 5 4 120 Kgs
L1 5 2 3 3 60 Hrs
L2 6 5 1 2 100 Hrs
Determine the optimal product mix for the firm and
carryout the sensitivity analysis for changes in the
objective funtion coefficients and RHS values of the
constraint.
FORMULATION OF LPP
x1, x2, x3 & x4: No. of units to be produced.
The LP model:
x1 x2 x3 x4 x5 x6 x7 x8
cj - zj -115/2 0 0 -85/3 -35/6 0 -125/6 0
Coeff 2 0 1 -1/3 -1/3 0 2/3 0
in
Row
x3
Ratio -115/4 - 0 85 35/2 - -125/4) -
(Least (Least
-ve) +ve
Allowable increase (Least +ve ratio) = 35/2 = 17.50
Allowable decrease (Least -ve ratio) = 115/4 = 28.75
SENSITIVITY ANALYSIS-Changes in
coefficient in the objective function
Changes in coefficient in the objective function for
Non basic variables are x1 & x4
x5 x6 x7 x8 xB x5 x6 x7 x8
2 120 infinity 25 95 to ∞
4 100 infinity 15 85 to ∞
INTRODUCING A NEW PRODUCT
•Suppose that the firm is contemplating introduction of a
new product P5 with a likely profit margin of Rs. 60.
•Data shows that product P5 shall consume
3 kg of R1 and 8 kg R2 (raw materials) &
4 hr of L1 and 6 hr of L2 (labour hours)
Should the firm decide for the introduction of the product?
Value of resources (RM & LH) consumed by Product P5
3 kg of R1 @ Rs. 35/6 Rs.3 Х 35/6 = 35/2 = 17.50
8 kg of R2 @ Rs. 0 Rs.8 Х 0 = 0
4 hr of L1 @ of 125/6 Rs.4 Х 125/6 = 250/3 = 83.33
6 hr of L2 @ of 0 Rs.6 Х 0 = 0
TOTAL 100.83
Since profit per unit of P5 (Rs60) is less than value of
resources (Rs100.83) consumed by it, the product should
not be included in the production list of the firm.
Decision Theory and
Game Theory
Maximin
The maximin strategy is a conservative one; it consists of identifying
the worst (minimum) payoff for each alternative and then selecting the
alternative that has the best (maximum) of the worst payoffs. In effect,
the decision maker is setting a floor for the potential payoff; the actual
payoff cannot be less than this amount.
Table 5 Maximax Solution for Real Estate Problem
Maximax
The maximax approach is the opposite of the previous one: The
best payoff for each alternative is identified, and the alternative with
the maximum of these is the designated decision.
Table 6 Payoff Table with Similar Maximum Payoffs
Minimax Regret
An approach that takes all payoffs into account. To use this approach,
it is necessary to develop an opportunity loss table that reflects the
difference between each payoff and the best possible payoff in a
column (i.e., given a state of nature). Hence, opportunity loss
amounts are found by identifying the best payoff in a column and then
subtracting each of the other values in the column from that payoff.
Table 7 Opportunity Loss Table for Real Estate Problem
Table 8 Identifying the Minimax Regret Alternative
Table 9 Minimax Regret Can Lead in a Poor Decision
The Hurwicz (Realism) Criterion (Weighted
Average or Realism Criterion)
• The approach offers the decision maker a compromise
between the maximax and the maximin criteria.
– Requires the decision maker to specify a degree of optimism, in
the form of a coefficient of optimism α, with possible values of
α ranging from 0 to 1.00.
– The closer the selected value of α is to 1.00, the more optimistic
the decision maker is, and the closer the value of α is to 0, the
more pessimistic the decision maker is.
Table 10 Equal Likelihood Criterion
Table 11 Summary of Methods for Decision Making under Complete
Uncertainty
A Retail store desires to determine the optimal daily
order size for a perishable item. The store buys the
perishable item at the rate of Rs. 80per Kg. and sells at
the rate Rs. 100 per Kg. if the order size is more than
the demand, the excess quantity can be sold at Rs. 70
per Kg in a secondary market; otherwise the opportunity
cost for the store is Rs. 15 per Kg for the unsatisfied
portion of the demand. Based on the past experience, it
is found that demand varies from 50 Kg to 250Kg, in
steps of 50. The possible value of the order size are
from 75Kg to 300Kg in steps of 75Kg. Determine the
optimal order size which will maximize the daily profit of
the store.
Game Theory
– By number of players
• Two players - Chess
• Multiplayer – Poker
– By total return
• Zero Sum - the amount won and amount lost by all
competitors are equal (Poker among friends)
• Nonzero Sum -the amount won and the amount lost by all
competitors are not equal (Poker In A Casino)
– By sequence of moves
• Sequential - each player gets a play in a given sequence.
• Simultaneous - all players play simultaneously.
Dominance property :
Reduce the size of payoff matrix by eliminating
redundant rows or columns
• Data
– The weekly percentage gain in market share for
IGA, as a function of advertising emphasis.
Sentry's Emphasis
Meat Produce Grocery Bakery
IGA's Meat 2 2 -8 6
Emphasis Produce -2 0 6 -4
Grocery 2 -7 1 -3
• Decision variables
– X1 = the probability IGA’s advertising focus is on
meat.
– X2 = the probability IGA’s advertising focus is on
produce.
– X 3 = the probability IGA’s advertising focus is on
groceries.
• Constraints
– IGA’s market share increase for any given
advertising focus selected by Sentry, must be at
least V.
• The model
IGA’s expected change
Max V in market share.
S.T.
Meat 2X1 – 2X2 + 2X3 V
Sentry’s Produce 2X1 – 7 X3 V
advertising Groceries -8X1 – 6X2 + X3 V
emphasis Bakery 6X1 – 4X2 – 3X3 V
Probability X1 + X2 + X3 = 1
Sentry’s Linear Programming Model
• Decision variables
– Y1 = the probability Sentry’s advertising focus is on meat.
– Y2 = the probability Sentry’s advertising focus is on
produce.
– Y 3 = the probability Sentry’s advertising focus is on
groceries.
– Y4 = the probability Sentry’s advertising focus is on
bakery.
• Constraints
– Sentry’s market share decrease for any given advertising
focus selected by IGA, must not exceed V.
• The Model
Min V
S.T.
2Y1 + 2Y2 – 8Y3 + 6Y4 V
-2Y1 + 6Y3 – 4Y4 V
2Y1 – 7Y2 + Y3 – 3Y4 V
Y1 + Y2 + Y3 + Y4 = 1
• For IGA
– X1 = 0.3889; X2 = 0.5; X3 = 0.1111
• For Sentry
– Y1 = .3333; Y2 = 0; Y3 = .3333; Y4 = .3333
▪ Introduction
▪ Setting Up a Transportation Problem
▪ Developing an Initial Solution: Northwest
Corner Rule
▪ Stepping-Stone Method: Finding a Least-
Cost Solution
▪ Modified Distribution (MODI) Method (also
known as U-V method)
▪ Vogel’s Approximation Method (VAM):
Another Way to Find an Initial Solution
▪ Unbalanced Transportation Problems
Session Outline
Factories Warehouses
(Sources) (Destinations)
Initial Solution
NorthWest Corner Rule
Least Cost Method
Vogel’s Approximate Method
Optimum Solution
Stepping Stone Method
MODI/UV/RK Method
Special-purpose
algorithms
Although standard LP methods can be
used to solve transportation problems,
special-purpose algorithms have been
developed that are more efficient.
They still involve finding an initial
solution and developing improved
solutions until an optimal solution is
reached.
They are fairly simple in terms of
computation.
Common Techniques
Figure 1
Setting Up a Transportation
Problem
Geographical locations of Executive Furniture’s
factories and warehouses
Boston
Cleveland
Factory
Des Moines
Evanston Warehouse
Albuquerque
Fort Lauderdale
Figure 2
Setting Up a Transportation
Problem
Production costs are the same at the three
factories so the only relevant costs are
shipping from each source to each destination.
Costs are constant no matter the quantity
shipped.
The transportation problem can be described
as how to select the shipping routes to be used
and the number of desks to be shipped on
each route so as to minimize total
transportation cost.
Restrictions regarding factory capacities and
warehouse requirements must be observed.
Setting Up a Transportation
Problem
The first step is setting up the transportation table.
Its purpose is to summarize all the relevant data
and keep track of algorithm computations.
EVANSVILLE $8 $4 $3
FORT LAUDERDALE $9 $7 $5
Setting Up a Transportation
Problem
Transportation table for Executive FurnitureDes Moines
capacity
constraint
DES MOINES $5 $4 $3
100
FACTORY
EVANSVILLE $8 $4 $3
300
FACTORY
FORT LAUDERDALE $9 $7 $5
300
FACTORY
WAREHOUSE
300 200 200 700
REQUIREMENTS
Cell representing a
Table 2 Total supply source-to-destination
Cost of shipping 1 unit from Cleveland (Evansville to Cleveland)
and demand
Fort Lauderdale factory to warehouse shipping assignment
Boston warehouse demand that could be made
Setting Up a Transportation
Problem
DES MOINES $5 $4 $3
100 100
(D)
EVANSVILLE $8 $4 $3
300
(E)
FORT LAUDERDALE $9 $7 $5
300
(F)
WAREHOUSE
300 200 200 700
REQUIREMENTS
Developing an Initial Solution:
Northwest Corner Rule
2. Assign 200 units from Evansville to
Albuquerque. This meets Albuquerque’s
demand. Evansville has 100 units remaining
so we move to the right to the next column of
the second row.
TO ALBUQUERQUE BOSTON CLEVELAND FACTORY
FROM (A) (B) (C) CAPACITY
DES MOINES $5 $4 $3
100 100
(D)
EVANSVILLE $8 $4 $3
200 300
(E)
FORT LAUDERDALE $9 $7 $5
300
(F)
WAREHOUSE
300 200 200 700
REQUIREMENTS
Developing an Initial Solution:
Northwest Corner Rule
3. Assign 100 units from Evansville to Boston.
The Evansville supply has now been
exhausted but Boston is still 100 units short.
We move down vertically to the next row in
the Boston column.
TO ALBUQUERQUE BOSTON CLEVELAND FACTORY
FROM (A) (B) (C) CAPACITY
DES MOINES $5 $4 $3
100 100
(D)
EVANSVILLE $8 $4 $3
200 100 300
(E)
FORT LAUDERDALE $9 $7 $5
300
(F)
WAREHOUSE
300 200 200 700
REQUIREMENTS
Developing an Initial Solution:
Northwest Corner Rule
DES MOINES $5 $4 $3
100 100
(D)
EVANSVILLE $8 $4 $3
200 100 300
(E)
FORT LAUDERDALE $9 $7 $5
100 300
(F)
WAREHOUSE
300 200 200 700
REQUIREMENTS
Developing an Initial Solution:
Northwest Corner Rule
5. Assign 200 units from Fort Lauderdale to
Cleveland. This exhausts Fort Lauderdale’s
supply and Cleveland’s demand. The initial
shipment schedule is now complete.
TO ALBUQUERQUE BOSTON CLEVELAND FACTORY
FROM (A) (B) (C) CAPACITY
DES MOINES $5 $4 $3
100 100
(D)
EVANSVILLE $8 $4 $3
200 100 300
(E)
FORT LAUDERDALE $9 $7 $5
100 200 300
(F)
WAREHOUSE
300 200 200 700
REQUIREMENTS
Developing an Initial Solution:
Northwest Corner Rule
We can easily compute the cost of this
shipping assignment
ROUTE
UNITS PER UNIT TOTAL
FROM TO SHIPPED x COST ($) = COST ($)
D A 100 5 500
E A 200 8 1,600
E B 100 4 400
F B 100 7 700
F C 200 5 1,000
4,200
TO FACTORY
A B C
FROM CAPACITY
$5 $4 $3
D 100 100
$8 $4 $3
E 100 200 300
$9 $7 $5
F 100 200 300
WAREHOUSE
300 200 200 700
REQUIREMENTS
ROUTE
DESKS PER UNIT TOTAL
FROM TO SHIPPED x COST ($) = COST ($)
D A 100 5 500
E B 200 4 800
E C 100 3 300
F A 200 9 1,800
F C 100 5 500
3,900
Obtaining an Improved
Solution
◼ Third and optimal solution
TO FACTORY
A B C
FROM CAPACITY
$5 $4 $3
D 100 100
$8 $4 $3
E 200 100 300
$9 $7 $5
F 200 100 300
WAREHOUSE
300 200 200 700
REQUIREMENTS
Summary of Steps in
Transportation Algorithm
(Minimization)
1. Set up a balanced transportation table.
2. Develop initial solution using either the northwest
corner method or Vogel’s approximation method.
3. Calculate an improvement index for each empty
cell using either the stepping-stone method or the
MODI method (also known as U-V method). If
improvement indices are all nonnegative, stop as
the optimal solution has been found. If any index
is negative, continue to step 4.
4. Select the cell with the improvement index
indicating the greatest decrease in cost. Fill this
cell using the stepping-stone path and go to step
3.
Optimal Solution
II. MODI Method (or U-V Method; also known as R-K
Method)
MODI Method (or U-V Method)
$8 $4 $3
R2 E 200 100 300
$9 $7 $5
R3 F 100 200 300
WAREHOUSE
300 200 200 700
REQUIREMENTS
Solving the Executive Furniture
Corporation Problem with MODI
The first step is to set up an equation for each
occupied square
By setting R1 = 0 we can easily solve for K1, R2,
K2, R3, and K3
(1) R1 + K1 = 5 0 + K1 = 5 K1 = 5
(2) R2 + K1 = 8 R2 + 5 = 8 R2 = 3
(3) R2 + K2 = 4 3 + K2 = 4 K2 = 1
(4) R3 + K2 = 7 R3 + 1 = 7 R3 = 6
(5) R3 + K3 = 5 6 + K3 = 5 K3 = –1
Solving the Executive Furniture
Corporation Problem with MODI
The next step is to compute the improvement
index for each unused cell using the formula
Improvement index (Iij) = Cij – Ri – Kj
We have
Des Moines- IDB = C12 – R1 – K2 = 4 – 0 – 1
Boston index = +$3
Des Moines- IDC = C13 – R1 – K3 = 3 – 0 – (–1)
Cleveland index = +$4
Evansville- IEC = C23 – R2 – K3 = 3 – 3 – (–1)
Cleveland index = +$1
Fort Lauderdale- IFA = C31 – R3 – K1 = 9 – 6 – 5
Albuquerque index = –$2
Solving the Executive Furniture
Corporation Problem with MODI
Table 1
Transportation Tableau
Table 2
Model
Optimal Solution of MG Auto
in Network Framework
<Please refer to the class notes
for further details>
Unbalanced Transportation
Problems
$8 $4 $3 0
E 50 200 50 300
$9 $7 $5 0
F 150 150 300
WAREHOUSE
300 200 200 150 850
REQUIREMENTS
▪ Introduction
▪ Network
▪ Setting up an LP Transshipment Problem
Introduction
5 200
600 1 2 2
3 6
3
3 6 150
6
4
3 4
7 350
1 6
4
5
400 2
8 300
Atlanta 3 1 Kansa 2 6 3 6
s
Louisv 4 4 6 5
ille
Model Formulation
Rearranging we get
-x13-x23+x35+x36+x37+x38=0
Similarly, the constraint corresponding to node 4 is
-x14-x24+x45+x46+x47+x48=0
To develop the constraints associated with the
destination nodes, we recognize that for each node the
amount shipped to the destination must equal to the
demand. For example, to satisfy the demand for 200
units at node 5 (the Detroit retail outlet), we write
X35+x45=200
Model Formulation
Minimize 2x13+3x14+3x23+1x24+2x35+6x36+3x37+6x38+4x45+
4x46+6x47+5x48
Subject to
x13+x14 <=600
x23+x24 <=400
-x13 -x23 +x35+x36+x37+x38 =0
-x14 -x24 +x45+x46+x47+x48 =0
x35 +x45 =200
x36 +x46 =150
x37 x47 =350
x38 +x48 =300
xij>=0 for all i and j
General Transshipment
Problem
Min c xi
all arcs
ij ij
s.t.
x
arcs out
ij − x
arcs in
ij si Origin nodes i
x
arcs out
ij − x
arcs in
ij = 0 Transshipment nodes
x
arcs in
ij − x
arcs out
ij = dj Destination nodes j
xij 0 i and j
where
5 200
600 1 2 2
3 6
3
3 6 150
6
4
3 4
7 350
1 6
4
5 1
400 2
4
8 300
Can study systems with preliminary data and critical data elements
can be identified
Simulation generates data by using existing data decreasing need
to collect additional data
Methods are often easier to apply than other methods
Analytical techniques are great, but usually require set of
simplifying assumptions that cannot always be justified
Simplifying assumptions are often not required
Importance of
Simulation…(cont’d.)
Continuous system
A system where state variables change continuously
over time
Many continuous systems are modeled as discrete and
vice versa for tractability
Deterministic and Stochastic
System
Deterministic system
A system where all state variables change only deterministically.
Stochastic system
A system where some state variables evolves probabilistically over
time.
Time Dependent and Time
Independent Simulations
Capstone is a rich,
complex business
simulation designed to
teach strategy,
competitive analysis,
finance, cross-functional
alignment, and the
selection of tactics to
build a successful and
focused company.
Business Simulation
Games
Finance Economics
Financial risk Economic indicators:
management
capital-labor ratio,
Option pricing exchange rates
Marketing Environmental sciences
Network systems Pollution transport and
Communication systems: prediction
Data flow under different
protocols Corrosion analysis
Energy Systems
Biological Systems
Epidemiology Groundwater and surface
Genetics water modeling
Meteorology Chemistry
Need of Models to Understand
System
Input Output
Advantages of Simulation
Define a problem
Introduce the variables associated with the problem
Construct a numerical model
Set up possible courses of action for testing
Run the experiment
Consider the results
Decide what courses of action to be taken
Process of Simulation
Define Problem
Introduce Important
Variables
Construct Simulation
Model
Specify Values of
Variables to Be Tested
Conduct the
Simulation
Examine the
Results
Validation
Check model results against historical or expectations of
system performance.
Review accuracy of predictions and plan experiments
accordingly.
Verify logic and data.
Building Blocks of Simulation
Models
A definition of the state of the system (the
number of customers in a queueing system)
Identification of possible states of the system
that can occur.
Identification of possible events (e.g., arrivals
and service completions in a queueing system)
that would change the state of the system.
Simulation clock (passage of time)
A method for randomly generating the events
A formula for identifying state transitions that
are generated by the various kinds of events.
Requirements of a
Simulation Model
Create random numbers
Create random variables
Record data for output
Perform statistical analysis on recorded data
Arrange and carry out computations using the output in
specified formats
Detect and report logical inconsistencies and other error
and provide correction
Two Statistical Theories
Inverse method
Acceptance-rejection method
Inverse Method
1
b a for a xb
f ( x)
0 for all other values
Area = 1
a b
Uniform Distribution
Probability
P( x X x ) x x
2 1
1 2
ba
45 42 1
P( 42 X 45)
47 41 2 Area
= 0.5
42 45
Uniform Distribution
Mean and Standard Deviation
Mean Mean
a+b 41 + 47 88
= = 44
2 2 2
Area = 1
Uniform Distribution
Mean and Standard Deviation
Mean Mean
a+b 39 + 27
= = 33
2 2
Uniform Distribution of
Assembly of Plastic Modules
35 30 5
P(30 X 35) 0.4167
39 27 12
35 30
0.4167
30 39 27
f (x) 42 45
27 35 39 x
Uniform Distribution of
Assembly of Plastic Modules
30 27 3
P( X 30) 0.2500
39 27 12
30
f ( x)
27 39 x
Uniform Distribution
Continuous
Family of distributions
Skewed to the right
X varies from 0 to infinity
Apex is always at X = 0
Steadily decreases as X gets larger
Probability function
X
f (X) e for X 0, 0
Different Exponential
Distributions
Exponential Distribution:
Computation of Probability
1.2
X 0
1.0
P X X 0 e
(12
P X 2| 12
. e
. )(2)
0.8
.0907
0.6
0.4
0.2
0.0
0 1 2 3 4 5
Generating Exponential
Random Variable
Inventory demand
Lead time for inventory
Times between machine breakdowns
Times between arrivals
Service times
Times to complete project activities
Number of employees absent
Simulation of a Queuing
Problem
Simulation of a Queuing Problem
Modeling waiting lines is an important application
of simulation
The assumptions of queuing models are quite
restrictive
Sometimes simulation is the only approach that fits
In a Port, fully loaded container ships arrive at
night for unloading
The number of ships per night varies from 0 - 5
The number of ships vary from day to day
Port of New Orleans
Table 10
Port of New Orleans
Table 11
Port of New Orleans
Simulation languages
Simula
Simscript
GPSS
Simio
Software packages
Crystal ball
@ Risk
GoldSim
If any error or inconsistencies are
found in this presentation, pl.
report it to [email protected]
Thank You
Queuing Theory Models
Dr. Prabir K. Das
Professor (QT-BRM/Analytics)
Indian Institute of Foreign Trade
Learning Objectives
Introduction
Waiting Line Costs
Characteristics of a Queuing System
Poisson process and exponential
distribution
Single-Channel Queuing Model with
Poisson Arrivals and Exponential Service
Times (M/M/1)
Multichannel Queuing Model with
Poisson Arrivals and Exponential Service
Times (M/M/m)
Introduction
* Service Level
Optimal
Figure 1 Service
Level
Elements of a Queuing
Model
Arrivals (from a source) (Customers are also
known as calling population)
Interarrival time: The arrival of customers is
represented by interarrival time between
successive customers
Service facilities
Service: is described by service time per customer
The actual waiting line
Characteristics of a Queuing
System
e x
P( X x) for x 0, 1, 2, 3, 4,...
x!
where
e x
P( x)
x!
e 2 20 0.1353(1)
P (0 ) 0.1353 14%
0! 1
e 2 21 e 2 2 0.1353(2)
P (1) 0.2706 27%
1! 1 1
e 2 22 e 2 4 0.1353( 4)
P (2 ) 0.2706 27%
2! 2(1) 2
Computation of
Probability
The probability of arrival (0, 1, 2, 3, …) can be easily computed using the
function “POISSON” in MS Excel.
Computation of
Probability
Characteristics of a Queuing
System
Two examples of the Poisson distribution for
arrival rates
0.30 – 0.25 –
0.25 – 0.20 –
0.20 –
0.15 –
Probability
Probability
0.15 –
0.10 –
0.10 –
0.05 – 0.05 –
0.00 –| | | | | | | | | | 0.00 –| | | | | | | | | |
0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9
X X
= 2 Distribution = 4 Distribution
Figure 2
Characteristics of a Queuing
System
Behavior of arrivals
Most queuing models assume customers are patient and
will wait in the queue until they are served and do not
switch lines
Customer may switch from one Q to another in the hope
of reducing waiting time. It is known as jockey
Balking refers to customers who refuse to join the queue
Reneging customers enter the queue but become
impatient and leave without receiving their service
That these behaviors exist is a strong argument for the
use of queuing theory to managing waiting lines
Characteristics of a Queuing
System
Queue
Service Departures
Arrivals Facility after Service
Queue
Type 1 Type 2 Departures
Arrivals Service Service after Service
Facility Facility
Figure 3 (a)
Characteristics of a Queuing
System
Four basic queuing system configurations
Service
Facility Departures
Queue
1
Service
Arrivals Facility after
2
Service
Facility Service
3
Figure 3 (b)
Characteristics of a Queuing
System
Four basic queuing system configurations
Type 1 Type 2
Queue Service Service
Facility Facility
1 1 Departures
Arrivals after
Service
Type 1 Type 2
Service Service
Facility Facility
2 2
Figure 3 (c)
Characteristics of a Queuing
System
Average
–
Service Time of
–
20 Minutes
Average Service Time of 1 Hour
Figure 4 |– | | | | | |
0 30 60 90 120 150 180 Service Time (Minutes)
The Role of the Exponential
Distribution
The operating characteristics of queuing
systems are determined largely by the
probability distribution of interarrival times
The probability distribution of service times
For real queuing systems these distributions
can take any form
To formulate a queuing theory model as a
representation of the real system, it is
necessary to specify the assumed form of
each of these distributions
The Role of the Exponential
Distribution
To be useful, the assumed form should be
sufficiently realistic that the model provides
reasonable predictions while, at the same time,
being sufficiently simple that the model is
mathematically tractable
Based on these considerations, the most
important probability distribution in queuing
theory is the exponential distribution
Let a random variable X represents interarrival
or service times
End of these times are considered as events
The Role of the Exponential
Distribution
A random variable X is said to have an
exponential distribution with
parameter α if its probability density
function is
f (x) = αe–αx for x ≥ 0
0 for x<0
and α> 0
a/b/c:d/e/f
a=Arrival b=Service time c=Number of service
distribution distribution channels open
L W
and
Lq Wq
W = Wq + 1/α
Single-Channel Model, Poisson
Arrivals, Exponential Service Times
(M/M/1)
1. The average time a customer spends in the
system, W W 1
L
3. The average number of customers in the
queue, Lq
2
Lq
( )
Single-Channel Model, Poisson
Arrivals, Exponential Service Times
(M/M/1)
4. The average time a customer spends waiting in
the queue, Wq
Wq
( )
Single-Channel Model, Poisson
Arrivals, Exponential Service Times
(M/M/1)
6. The percent idle time, P0, the probability no
one is in the system
P0 1
Pw An arriving customer has to wait
7. The probability that the number of customers
in the system is greater than k, Pn>k
k 1
Pn k
Case: Arnold’s Muffler Shop
Lq
2
22
4
( ) 3(3 2) 3(1)
1.33 cars waiting in line
on the average
Wq
2
hour
( ) 3
40 minutes average
waiting time per car
in the Q
2
0.67
3
percentage of time
mechanic is busy
Case: Arnold’s Muffler
Shop
P0 1
2
1 0.33
3
probability that there are 0
cars in the system
k Pn>k = (2/3)k+1
where
m = number of channels
Cs = service cost of each channel
Case : Arnold’s Muffler Shop
Waiting cost when the cost is based on time in
the system
Total waiting cost (Total time spent waiting by all arrivals) x
= (Cost of waiting)
Arnold has identified the mechanics wage $7 per hour as the service
cost
L
2
2 1 car in the system on the
42 2
average
W
1
1 1/2 hour that an average car spends in
42 the system
Case : Arnold’s Muffler Shop
Wq
1
hour 15 minutes average waiting time
( ) 4 per car
2
0.5 percentage of time mechanic is busy
4
P0 1
2
1 0.5 probability that there are 0 cars in
4 the system
Case : Arnold’s Muffler Shop
k Pn>k = (2/4)k+1
0 0.500
1 0.250
2 0.125
3 0.062
4 0.031
5 0.016
6 0.008
7 0.004
Case : Arnold’s Muffler Shop
The customer waiting cost is the same $10 per
hour
Total daily waiting
cost = (8 hours per day) WqCw
= (8)(2)(1/4)($10) = $40.00
1
P0 for m
nm1 1 n
1 mm
n0 n! m! m
Multichannel Model, Poisson
Arrivals, Exponential Service
Times (M/M/m)
2. The average number of customer in the system
( / )m
L P
(m 1)!(m ) 2 0
3. The average time a unit spends in the waiting line or being served, in
the system
( / ) m 1 L
W P0
(m 1)!(m ) 2
Multichannel Model, Poisson
Arrivals, Exponential Service
Times (M/M/m)
4. The average number of customers or units in
line waiting for service
Lq L
1 Lq
Wq W
m
Arnold’s Muffler Shop
Revisited
Arnold wants to investigate opening a second
garage bay
He would hire a second worker who works at the
same rate as his first worker
The customer arrival rate remains the same
1
P0 for m
nm1 1 n 1
m
m
n0 n!
m! m
P0 0.5
probability of 0 cars in the system
Arnold’s Muffler Shop
Revisited
Average number of cars in the system
( / )m
L P 0.75
(m 1)!(m ) 2 0
L 3
W hours 22 1 minutes
8 2
Arnold’s Muffler Shop
Revisited
Average number of cars in the queue
3 2 1
Lq L 0.083
4 3 12
1 Lq 0.083
Wq W 0.0415 hour 2 1 minutes
2 2
Arnold’s Muffler Shop
Revisited
Effect of service level on Arnold’s operating
characteristics
LEVEL OF SERVICE
ONE TWO ONE FAST
OPERATING MECHANIC MECHANICS MECHANIC
CHARACTERISTIC α= 3 α= 3 FOR BOTH α=4
Probability that the system
0.33 0.50 0.50
is empty (P0)
Average number of cars in
2 cars 0.75 cars 1 car
the system (L)
Average time spent in the
60 minutes 22.5 minutes 30 minutes
system (W)
Average number of cars in
1.33 cars 0.083 car 0.50 car
the queue (Lq)
Average time spent in the
40 minutes 2.5 minutes 15 minutes
queue (Wq)
Table 2
Arnold’s Muffler Shop
Revisited
Adding the second service bay reduces the waiting
time in line but will increase the service cost as a
second mechanic needs to be hired
Total daily waiting cost = (8 hours per day) WqCw
= (8)(2)(0.0415)($10) = $6.64
Table 3
Port of New Orleans
Table 4
Port of New Orleans