Gideon Weiss University of Haifa. Joint work with students: Anat Kopzon Yoni Nazarathy. Stanford University, MSE, February, 2009
|
|
- Sophia Rosalyn Gilbert
- 5 years ago
- Views:
Transcription
1 Optimal Finite Horizon Control of Manufacturing Systems: Fluid Solution by SCLP (separated continuous LP) and Fluid Tracking using IVQs (infinite virtual queues) Stanford University, MSE, February, 29 Gideon Weiss University of Haifa Joint work with students: Anat Kopzon Yoni Nazarathy Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29
2 Control of computer chip semi-conductor wafer fab Cost: 3x1 9 $ Return: 3 years Cycle time: 6 weeks WIP: 6, wafers, 18x1 6 $ Challenge: Control the queues at ~5 work steps No steady state, Finite horizon control, Rolling horizon update Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 1
3 Finite horizon control of multi class queueing networks MCQN Control MCQN (Q(t),T(t)), over <t<t objective Q k (t) = Q k ()! S k (T k (t)) + " k 'k (S k ' (T k ' (t))) k #K $ k ' min "! k T k (T ) + c k Q k (t)dt k # T Fluid Problem T V * = min ("! u(s) + c! q(s) )ds s.t. # t # Ru(s)ds + q(t) = q(), Au(t) $ 1, u(t),q(t) %, t &[,T ] Optimal Fluid Solution Tracking of Fluid: - Model Deviations from Fluid as MCQN with IVQ - Use Max Pressure to keep Deviations Stable Asymptotically Optimal Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 2
4 Continuous linear programming Bellman (1953) Economic models, Dantzig, Perold, Anstreicher. T " max c!(s)u(s)ds H (t) u(t) + Anderson (1975) job shop scheduling, " G(s,t)u(s)ds # a(t) CLP Anderson, Philpott & Nash papers and Book (1987),. t u(t) $, < t < T T # max c(t! s )" u(s)ds t # G u(s)ds + x(t) = a(t) H u(t) x, u %, < t < T $ b(t) SCLP Pullan (1993) Dual problem, Strong duality, convergent (not finite) algorithm Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 3
5 SCLP for MCQN T ( ) max $ ("! + (T # t) c! )u(t) + d! x(t) dt s.t. t $ G u(s)ds + Fx(t) % & + at H u(t) % b u(t), x(t) ' < t < T Variables: K+J+I+L States Controls G KxJ F KxL H IxJ slacks x k : k=1,,k, k=k +1,,K + L, slacks u j : j=1,,j, j= J +1,,J + I Structure: Time horizon T Control cost! Inventory cost c d System G I F a input rate! initial state H I b control capacity Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 4
6 Solving SCLP Discretization vs Simplex Discrete Time: CLP and SCLP become LP (1) LP is large (2) approximation is doubtful (3) problem structure is lost Most past approaches (including Pullan) used discrete approximation SCLP provides excellent models, but because it needed discretization, problems were modeled as dictrete time multi-period LPs We solve SCLP in continuous time, exactly, in a finite number of steps. Our algorithm performs simplex steps, in function space. We gain much insight to the solution, e.g. sensitivity analysis in time and space Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 5
7 A simplex approach Solution features: partition of time horizon =t <t 1 < <t N =T piecewise constant controls u(t) continuous piecewise linear x(t) Symmetric dual Extreme points are sequences of bases Edges via validity regions SCLP Pivots Parametric sequence of steps T ( ) max $ ("! + (T # t) c! )u(t) + d! x(t) dt s.t. t $ G u(s)ds + Fx(t) % & + at H u(t) u(t), x(t) ' % b < t < T Simplex in the space of bounded measurable (density) controls Essential Assumption! a # $ is in general " b% &! G I F$ position to # " H I & % Problem Non-Degenerate:! c # $ is in general " d% &! G' (I H ' $ position to # " F ' (I & % Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 6
8 Symmetric Duality & Complementary Slackness max $ ("! + (T # t) c! )u(t) + d! x(t) dt s.t. T t ( ) $ G u(s)ds + Fx(t) % & + at H u(t) % b u(t), x(t) ' < t < T T " x k (t)p k (T! t)dt = " q j (T! t)u j (t)dt = primal states primal controls x k (t) u j (t) T T (( ) p(t) + b! q(t) ) dt min $ "! + (T # t) a! t s.t. $ G! p(s)ds + H! q(t) % & + ct F! p(t) % d dual states dual controls q j (T! t) p k (T! t) p(t),q(t) % < t < T Dual runs in reversed time We produce primal dual feasible complementary slack solutions Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 7
9 Complementary slack optimal solutions x k!x k 1!x k 2 p k 3 = t t 1 t 2 t N = T t 3 T! t = T T! t 1 T! t 2 T! t N = T! t 3 u j 1!q j 2!q j 3 q j Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 8
10 Boundary and Rates LP/LP* Stage 1 Solve for T = (horizon dt) Determine x(), q() Boundary LP max d! x s.t. Fx " #, x $ T! c d G I F H I min b! q N a b! s.t. H! q N " #, q N " T max $ (("! + (T # t) c! )u(t) + d! x(t) ) dt t s.t. $ G u(s)ds + Fx(t) % & + at H u(t) % b u(t), x(t) ' < t < T x(), q() remain unchanged for all time horizons Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 9
11 Boundary and Rates LP/LP* Boundary LP max d! x s.t. Fx " #, x $ min b! q N s.t. H! q N " #, q N " Determine x(), q() Stage 2 Determine values of u(t), p(t-t) slopes of x(t), q(t-t) for t n-1 <t<t n T! c d G I F H I a b! T max $ (("! + (T # t) c! )u(t) + d! x(t) ) dt t s.t. $ G u(s)ds + Fx(t) % & + at H u(t) % b u(t), x(t) ' < t < T Rates LP max c! u n + d!!x n s.t. Gu n + F!x n " a Hu n " b min a! p n + b!!q n s.t. G! p n + H!!q n " c F! p n " d Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 1
12 Boundary and Rates LP/LP* Boundary LP max d! x s.t. Fx " #, x $ min b! q N s.t. H! q N " #, q N " Determine x(), q() Rates LP max c! u n + d!!x n s.t. Gu n + F!x n " a Hu n " b min a! p n + b!!q n s.t. G! p n + H!!q n " c F! p n " d Determine values of u(t), p(t-t) slopes of x(t), q(t-t) for t n-1 <t<t n Sign constraints for the interval t n-1 <t<t n n n If x k (t n!1 ) > then!x k unrestricted, pk = n n If q j (T! t n ) > then!q j unrestricted, u j = n Else :!x k, u n j,!q n j, n pk " Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 11
13 Solution Structure theorem SCLP solution is determined by partition of time horizon =t <t 1 <...<t N =T and a sequence of bases B 1,,B N B n :!x n,u n,!q n, p n B n! B n+1 adjacent " n = t n # t n#1 B n! B n+1 : &!x k leaves the basis " = x k (t n ) = x k + # m ( $ m!x k ( m=1 ' N ( N u j leaves the basis " = q j (T % t n ) = q j + # m $ m!q ( j ) m=n+1 n Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 12
14 The Structure Theorem Equations and inequalities for interval lengths and slacks " 1!1 % $ '"( % $ A ' $ ' = # $ B!I&' #) & " T% $ ' $ g' # $ h& ' Theorem If a base-sequence B 1,, B N is: Admissible (u,p ), adjacent, consistent with the boundary (x,q N ) and has interval lengths and slacks!," >, then the solution u, x, p,q is optimal Conversely: almost all solutions are of this form These Base-sequences are the extreme points Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 13
15 Validity Regions Theorem An optimal base-sequence B 1,, B N is optimal for a convex polyhedral cone of boundary values x,q N,T. We call this the Validity Region of the base sequence Definition: Two base sequences whose validity regions touch are Neighboring base sequences. This defines an Edge Moving along an edge is an SCLP Pivot. Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 14
16 Parametric Algorithm x t, qt N,Tt! (R)! (1) x s, qs N,Ts Algorithm pivots at! " <! (1) <!<! (R) < 1"! Moving from N x s,qs,ts to N x t,qt,tt Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 15
17 SCLP pivot T < T (r) T = T (r) A necessary condition for pivoting is that the problem is completely nondegenerate:! a$ "# b %&! c$ "# d% & is in general position to is in general position to! G I F$ "# H I % &! G ' (I H ' $ "# F ' (I %& If that condition does not hold we have no way to find where to `repair solution after T (r] T > T (r) Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 16
18 Simple SCLP example max 6 " t ((8! t)u 1 (t)! x 2 (t)) dt s.t. " u 1 (s)ds + x 2 (t) # 3 + t u 1 (t) # 2 u, x $ Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 17
19 Simple SCLP example 1 T " t ( ) dt max (!4 + 2(T! t) )u 1 (t)! x 2 (t) s.t. " u 1 (s)ds + x 2 (t) # 3 +t u 1 (t) # 2 u,x $ Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 18
20 Simple SCLP example 2 T " t ( ) dt max (!4 + 2(T! t) )u 1 (t)! x 2 (t) s.t. " u 1 (s)ds + x 2 (t) # 3 +t u 1 (t) # 2 u,x $ T " ( ) dt min (3+ T! t)p 1 (t) + 2q 2 (t) t s.t. " p 1 (s)ds + q 2 (t) #!4 + 2t p 1 (t) #!1 p,q # Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 19
21 Simple SCLP example 3 T " t ( ) dt max (!4 + 2(T! t) )u 1 (t)! x 2 (t) s.t. " u 1 (s)ds + x 1 (t) + x 2 (t) = 3 +t u 1 (t) + u 2 (t) = 2 u,x # T " ( ) dt min (3+ T! t)p 1 (t) + 2q 2 (t) t s.t. " p 1 (s)ds! q 1 (t) + q 2 (t) =!4 + 2t p 1 (t)! p 2 (t) =!1 p,q # Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 2
22 Simple SCLP example 4 T " t ( ) dt max (!4 + 2(T! t) )u 1 (t)! x 2 (t) s.t. " u 1 (s)ds + x 1 (t) + x 2 (t) = 3 +t u 1 (t) + u 2 (t) = 2 u,x # T " min (3+ T! t)p 1 (t) + 2q 2 (t) ( ) dt t s.t. " p 1 (s)ds! q 1 (t) + q 2 (t) =!4 + 2t p 1 (t)! p 2 (t) =!1 p,q # Boundary LP T= max! x 2 s.t. x 1 + x2 = 3 x 1 min 2q 2 N s.t.! q 1 N + q2 N =!4 x " Solution x 1 = 3 T! t = T t N = T q N " N Solution q 1 = 4 q 1 N Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 21
23 Simple SCLP example 7 Boundary LP max! x 2 min 2q 2 N s.t. x 1 + x2 = 3 x " Solution x 1 = 3 x 1 x 1 s.t.! q 1 N + q2 N =!4 q N " Solution q 1 N = 4 Rates LP <T<2 max 2u 1!!x 2 s.t. u 1 +!x 1 +!x 2 = 1 u 1 + u 2 = 2 min p 1 + 2!q 2 s.t. p 1!!q 1 +!q 2 = 2 p 1! p 2 =!1!x 1 U, u 1 Z,!x 2,u 2 P Solution u 2 = 2,!x 1 = 1 q 1 q 1 N!q 1 U, p 1 Z, p 2,!q 2 P Solution p 2 = 1,!q 1 =!2 Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 22
24 Simple SCLP example 7 Boundary LP max! x 2 min 2q 2 N s.t. x 1 + x2 = 3 x " Solution x 1 = 3 x 1 x 1 s.t.! q 1 N + q2 N =!4 q N " Solution q 1 N = 4 Rates LP <T<2 max 2u 1!!x 2 s.t. u 1 +!x 1 +!x 2 = 1 u 1 + u 2 = min p 1 + 2!q 2 s.t. p 1!!q 1 +!q 2 = 2 p 1! p 2 =!1!x 1 U, u 1 Z,!x 2,u 2 P Solution u 2 = 2,!x 1 = 1 q 1 q 1 N!q 1 U, p 1 Z, p 2,!q 2 P Solution p 2 = 1,!q 1 =!2 Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 23
25 Simple SCLP example 7 Boundary LP max! x 2 min 2q 2 N s.t. x 1 + x2 = 3 s.t.! q 1 N + q2 N =!4 x " Solution x 1 = 3 x 1 x 1 q N " Solution q 1 N = 4 Rates LP <T<2 max 2u 1!!x 2 s.t. u 1 +!x 1 +!x 2 = 1 u 1 + u 2 = min p 1 + 2!q 2 s.t. p 1!!q 1 +!q 2 = 2 p 1! p 2 =!1!x 1 U, u 1 Z,!x 2,u 2 P Solution u 2 = 2,!x 1 = 1 q 1!q 1 U, p 1 Z, p 2,!q 2 P Solution p 2 = 1,!q 1 =!2 q 1 N Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 24
26 Simple SCLP example 8 Boundary LP max! x 2 min 2q 2 N s.t. x 1 + x2 = 3 s.t.! q 1 N + q2 N =!4 x " Solution x 1 = 3 x 1 x 1 q N " Solution q 1 N = 4 Rates LP <T<2 max 2u 1!!x 2 s.t. u 1 +!x 1 +!x 2 = 1 u 1 + u 2 = min p 1 + 2!q 2 s.t. p 1!!q 1 +!q 2 = 2 p 1! p 2 =!1!x 1 U, u 1 Z,!x 2,u 2 P Solution u 2 = 2,!x 1 = 1 q 1!q 1 U, p 1 Z, p 2,!q 2 P Solution p 2 = 1,!q 1 =!2!x 1 U,!x 2,u 1,u 2 P q 1 N p 1 Z,!q 1, p 2,!q 2 P Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 25
27 Simple SCLP example 8 Boundary LP max! x 2 min 2q 2 N s.t. x 1 + x2 = 3 s.t.! q 1 N + q2 N =!4 x " Solution x 1 = 3 x 1 x 1 q N " N Solution q 1 = 4 Rates LP 2<T max 2u 1!!x 2 s.t. u 1 +!x 1 +!x 2 = 1 u 1 + u 2 = 2 q 2 min p 1 + 2!q 2 s.t. p 1!!q 1 +!q 2 = 2 p 1! p 2 =!1!x 1 U, u 1 Z,!x 2,u 2 P!q 1 U, p 1 Z, p 2,!q 2 P Solution u 2 = 2,!x 1 = 1!x 1 U,!x 2,u 1,u 2 P Solution u 1 = 2,!x 1 =!1 q 1 q 1 N Solution p 2 = 1,!q 1 =!2 p 1 Z,!q 1, p 2,!q 2 P Solution p 2 = 1,!q 2 = 2 Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 26
28 Simple SCLP example 9 Boundary LP max! x 2 s.t. x 1 + x2 = 3 x " Solution x 1 = 3 x 1 min 2q 2 N N s.t.! q 1 + N q2 =!4 q N " N Solution q 1 = 4 Rates LP 2<T<5 max 2u 1!!x 2 s.t. u 1 +!x 1 +!x 2 = 1 u 1 + u 2 = 2 5 x 1 x 1 5 min p 1 + 2!q 2 s.t. p 1!!q 1 +!q 2 = 2 p 1! p 2 =!1!x 1 U, u 1 Z,!x 2,u 2 P!q 1 U, p 1 Z, p 2,!q 2 P Solution u 2 = 2,!x 1 = 1!x 1 U,!x 2,u 1,u 2 P Solution u 1 = 2,!x 1 =!1!x 1,!x 2,u 1,u 2 P q 2 q 1 q 1 N Solution p 2 = 1,!q 1 =!2 p 1 Z,!q 1, p 2,!q 2 P Solution p 2 = 1,!q 2 = 2 p 1,!q 1, p 2,!q 2 P Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 27
29 Simple SCLP example 1 Boundary LP max! x 2 s.t. x 1 + x2 = 3 x " Solution x 1 = 3 x 1 min 2q 2 N N s.t.! q 1 + N q2 =!4 q N " N Solution q 1 = 4 Rates LP 5<T< max 2u 1!!x 2 s.t. u 1 +!x 1 +!x 2 = 1 u 1 + u 2 = 2 6 x 1 x 1 6 min p 1 + 2!q 2 s.t. p 1!!q 1 +!q 2 = 2 p 1! p 2 =!1!x 1 U, u 1 Z,!x 2,u 2 P!q 1 U, p 1 Z, p 2,!q 2 P Solution u 2 = 2,!x 1 = 1 q 1 Solution p 2 = 1,!q 1 =!2!x 1 U,!x 2,u 1,u 2 P q 2 N p 1 Z,!q 1, p 2,!q 2 P q 1 Solution u 1 = 2,!x 1 =!1 Solution p 2 = 1,!q 2 = 2!x 1,!x 2,u 1,u 2 P Solution u 1 = 1, u 2 = 1 p 1,!q 1, p 2,!q 2 P Solution p 1 = 2, p 2 = 3 Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 28
30 Simple SCLP example 11 Boundary LP max! x 2 s.t. x 1 + x2 = 3 x " Solution x 1 = 3 Rates LP 2<T<5 max 2u 1!!x 2 s.t. u 1 +!x 1 +!x 2 = 1 u 1 + u 2 = 2!x 1 U, u 1 Z,!x 2,u 2 P Solution u 2 = 2,!x 1 = 1!x 1 U,!x 2,u 1,u 2 P Solution u 1 = 2,!x 1 =!1!x 1,!x 2,u 1,u 2 P Solution u 1 = 1, u 2 = 1 x 1 6 x 1 p 2 u 1 q 2 p 2 p 1 u 1 u 2 min p 1 + 2!q 2 s.t. p 1!!q 1 +!q 2 = 2 p 1! p 2 =!1!q 1 U, p 1 Z, p 2,!q 2 P Solution p 2 = 1,!q 1 =!2 p 1 Z,!q 1, p 2,!q 2 P Solution p 2 = 1,!q 2 = 2 p 1,!q 1, p 2,!q 2 P Solution p 1 = 2, p 2 = 3 Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, p 2 u 2 q 1 x 1 6 q 1 N min 2q 2 N s.t.! q 1 N + q2 N =!4 q N " Solution q 1 N = 4
31 SCLP algorithm demo x q Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 3
32 SCLP algorithm demo x q Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 31
33 SCLP algorithm demo x q Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 32
34 SCLP algorithm demo x q Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 33
35 SCLP algorithm demo x q Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 34
36 SCLP algorithm demo x q Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 35
37 SCLP algorithm demo x q Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 36
38 SCLP algorithm demo Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 37
39 Finite horizon control of multi class queueing networks MCQN Control MCQN (Q(t),T(t)), over <t<t objective Q k (t) = Q k ()! S k (T k (t)) + " k 'k (S k ' (T k ' (t))) k #K $ k ' min "! k T k (T ) + c k Q k (t)dt k # T Fluid buffers can have fluid with positive flow In real system: Standard Qs k!k Optimal Fluid Solution Or they have fluid in whole interval: In real system: IVQs k!k " IVQs have nominal outflow rates, from fluid solution:! k = u k : k "K # Tracking of Fluid: - Model Deviations from Fluid as MCQN with IVQ - Use Max Pressure to keep Deviations Stable Asymptotically Optimal Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 38
40 Multi Class Queueing Networks MCQN (Harrison 88, Dai 94, Q k (t) = Q k () + A k (t)! S k (T k (t)) + % k '#$ " k 'k (S k ' (T k ' (t))) Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 39
41 MCQN with Infinite Virtual Queues (Massey 84, W 6) Two types of classes Standard Queues Infinite Virtual Queues! 6 K = {1,2,3,5} K! = {4,6} No exogenous arrivals Nominal Input Rate 5 4! k!k : Q k (t) = Q k () " S k (T k (t)) + % # k 'k (S k ' (T k ' (t))) & k '!$ k!k ' : Q k (t) = Q k () + ( k t " S k (T k (t)) Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 4
42 Infinite Virtual Queues - Queueing vs Manufacturing Single server queue Machine with controlled input! µ min(!,µ)! µ µ A tandem of queues! µ 1 µ 2 µ 3! Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 41
43 A new paradigm for balanced heavy traffic Kumar-Seidman Rybko-Stolyar network (KSRS) 1 2 " 1! 1 µ 1 Heavy traffic when! 1!,! 2! µ 2! 2 " 2 Balanced heavy traffic! 1 " # 1,! 2 " # 2 Push-Pull network - full utilization, ρ = 1, no congestion Both machines work non-stop, and no flow accumulates implies:!! 1 µ 1 " Q 1 µ 2 Q 2! 2 " 1! 1 = " 1 # 1 = (1$ " 2 )µ 1 = # 1µ 1 (µ 2 $ # 2 ) µ 1 µ 2 $ # 1 # 2! 2 = " 2 # 2 = (1$ " 1 )µ 2 = # 2µ 2 (µ 1 $ # 1 ) µ 1 µ 2 $ # 1 # 2 Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 42
44 Finite horizon control of multi class queueing networks MCQN Control MCQN (Q(t),T(t)), over <t<t objective Q k (t) = Q k ()! S k (T k (t)) + " k 'k (S k ' (T k ' (t))) k #K $ k ' min "! k T k (T ) + c k Q k (t)dt k # T Fluid buffers can have fluid with positive flow In real system: Standard Qs Or they have fluid in whole interval: In real sysatem: IVQs k!k " IVQs have nominal outflow rates, from fluid solution:! k = u k : k "K # k!k Optimal Fluid Solution Tracking of Fluid: - Model Deviations from Fluid as MCQN with IVQ - Use Max Pressure to keep Deviations Stable Asymptotically Optimal Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 43
45 Example Network 1 2 µ 1 Stacked Queue level representation: µ 3 µ 2 Stacked Queue Levels Finished Jobs Trajectory of a single job Resource 2 is the bottleneck 1 > µ 2 µ 1 µ 3 Q2 Q3 How to control queues 1,3 Q1 Schedule resource 1, for optimal draining time T Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 44
46 Fluid formulation 1 2 min T! q 1 (t) + q 2 (t) + q 3 (t)dt s.t : q 1 (t) = q 1 () " µ 1 u 1 (s)ds t! t q 2 (t) = q 2 () +! µ 1 u 1 (s)ds "! µ 2 u 2 (s)ds t q 3 (t) = q 3 () +! µ 2 u 2 (s)ds "! µ 3 u 3 (s)ds t t µ 1 µ 3 µ 2 u 1 (t) + u 3 (t) # 1 u 2 (t) # 1 u,q $ t %[,T ] This is a Separated Continuous Linear Program (SCLP) Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 45
47 Fluid policies: LBFS, min Makespan, Optimal LBFS V=376 min Makespan V=36 Optimal V= µ 1 µ 3 µ 2 Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 46
48 Time intervals of Optimal Fluid Solution! = [, 4)! = [4,8)! = [8, 24)! = [24, 4) Allocations u 1 = u 1 = u 1 =.25 u 1 =.25 u 2 = u 2 = 1 u 2 = 1 u 2 = 1 u 3 = 1 u 3 = 1 u 3 =.75 u 3 =.25 Define: n " = { k qk ( t) =, t #! n} n # = { k q ( t) $, t %! } " k n Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 47
49 Maximum pressure policy (Tassiulas, Stolyar, Dai, Lin, Ata) Max pressure achieves stability for any system with offered load < 1! d dt d dt Q(t) =! " Ru(t)! 2 d Q k (t) = dt Q(t)iQ(t) = 2Q(t)i(" # Ru(t)) k rate stability if offered load = 1 Consider MCQN with fluid dynamics described by where R is the input output matrix, u(t) is the machine allocation, and is the input rate. The machine allocations (controls ) are subject to resource constraints Max pressure attempts to maximize the gradient of the sum of squares of queue lengths At any time t choose allocation u(t) such that max Q(t)iRu(t) s.t. Au(t)! 1, u(t) ", u(t) is available In balanced heavy traffic it may optimize the diffusion approximation Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 48
50 Maximum Pressure (Tassiulas, Stolyar, Dai-Lin) for MCQN with IVQs Theorem (Dai and Lin 25): MCQN under Max Pressure, when ρ 1 is rate stable: Q lim k (t) t!" t = Adaptation to MCQN with IVQs Theorem: Same holds for MCQN with IVQ, k!k : Q k (t) = Q k () " S k (T k (t)) + % # k 'k (S k ' (T k ' (t))) & k '!$ k!k ' : Q k (t) = Q k () + ( k t " S k (T k (t)) At any time t choose allocation u(t) such that max Q(t)iRu(t) s.t. Au(t)! 1, u(t) ", u(t) is available u(t) = T!(t)! k come from the fluid solution Q k (t) Measure deviations from fluid solution Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 49
51 Input-Output matrixes when Infinite Virtual Queues exist Time intervals 1 and 2 Time interval 3 Time interval R! µ " # $ = R = µ 2 # µ $ % 3 & R 3! µ 1 " # $ = % µ 1 µ 2 # µ $ & 3 ' R 4 = " µ 1 % $ ' $!µ 1 µ 2 ' # $!µ 2 µ 3 & ' Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 5
52 Implementing Maximum Pressure Pressure for buffer k based on k and k+1 current downstream k + 1! K k + 1" K! k! K µ Q µ k ( Qk! Qk + 1) k k k " K! µ k (! k t " D k " Q k +1 ) µ k (! k t " D k ) Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 51
53 Using Max Pressure to track fluid solution by MCQN-IVQ Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 52
54 Asymptotic Optimality Theorem: Let Q(t) be the queue length process of a finite horizon MCQN. Let Q N (t) be N scalings of Q(t), with Q N ()=NQ(), and with N-fold speed Of processing. Let q(t) be the optimal fluid solution and let V f be its Objective value. (i) Let V N denote the objective values of Q N (t) for any general policy. Then: lim inf N!" 1 N V N # V f a.s. (ii) and Under max pressure tracking of the fluid solution lim N!" 1 N Q(t)N = q(t) a.s. uniformly on < t < T lim N!" 1 N V N = V f a.s. Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 53
55 Summary We considered control of MCQN over a finite time horizon We formulated a fluid approximation and formulated SCLP We described a simplex type algorithm to solve SCLP Solution of SCLP provides optimal fluid solution We Introduced MCQN with IVQs IVQs add another level of control to MCQN IVQs provide a new paradigm for balanced heavy traffic. We fitted a MCQN with IVQ model to each interval of the fluid solution We adapted Max Pressure policy to track fluid solution This provides asymptotically optimal policy Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 54
56 References (1) A Simplex based algorithm to solve SCLP W, Math Prog A, 115: , 28. (2) Near optimal control of queueing networks over a finite time horizon Yoni Nazarathy and Gideon Weiss, Annals of OR online (3) A Push-Pull network with infinite supply of work Anat Kopzon, Yoni Nazarathy and Gideon Weiss, Under review, QUESTA (4) Positive Harris recurrent and diffusion scale analysis of a Push-Pull queueing network with infinite supply of work Yoni Nazarathy and Gideon Weiss, To be presented at Value Tools 28, Performance Evaluation, to appear Challenges - Can we track fluid solution so as to be optimal on the diffusion scale? - Can we extend the policy for Push-Pull to general networks with infinite supply at all heavy traffic nodes? Gideon Weiss, University of Haifa, Finite Horizon Manufacturing - SCLP & IVQs, 29 55
A Semiconductor Wafer
M O T I V A T I O N Semi Conductor Wafer Fabs A Semiconductor Wafer Clean Oxidation PhotoLithography Photoresist Strip Ion Implantation or metal deosition Fabrication of a single oxide layer Etching MS&E324,
More informationControl of computer chip semi-conductor wafer fabs
Conrol of compuer chip semi-conducor wafer fabs Cos: x9 $ Reurn: years Cycle ime: 6 weeks WIP: 6, wafers, 8x6 $ Challenge: Conrol he queues a ~5 work seps Conrol of Manufacuring Sysems:! Sochasic racking
More informationPositive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network. Haifa Statistics Seminar May 5, 2008
Positive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network Yoni Nazarathy Gideon Weiss Haifa Statistics Seminar May 5, 2008 1 Outline 1 Preview of Results 2 Introduction Queueing
More informationThis lecture is expanded from:
This lecture is expanded from: HIGH VOLUME JOB SHOP SCHEDULING AND MULTICLASS QUEUING NETWORKS WITH INFINITE VIRTUAL BUFFERS INFORMS, MIAMI Nov 2, 2001 Gideon Weiss Haifa University (visiting MS&E, Stanford)
More informationControl of multi-class queueing networks with infinite virtual queues
Control of multi-class queueing networks with infinite virtual queues Erjen Lefeber (TU/e) Workshop on Optimization, Scheduling and Queues Honoring Gideon Weiss on his Retirement June 8, 2012 Where innovation
More information5 Lecture 5: Fluid Models
5 Lecture 5: Fluid Models Stability of fluid and stochastic processing networks Stability analysis of some fluid models Optimization of fluid networks. Separated continuous linear programming 5.1 Stability
More informationMaximum pressure policies for stochastic processing networks
Maximum pressure policies for stochastic processing networks Jim Dai Joint work with Wuqin Lin at Northwestern Univ. The 2011 Lunteren Conference Jim Dai (Georgia Tech) MPPs January 18, 2011 1 / 55 Outline
More informationPositive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network
Positive Harris Recurrence and Diffusion Scale Analysis of a Push Pull Queueing Network Yoni Nazarathy a,1, Gideon Weiss a,1 a Department of Statistics, The University of Haifa, Mount Carmel 31905, Israel.
More informationLinear Model Predictive Control for Queueing Networks in Manufacturing and Road Traffic
Linear Model Predictive Control for ueueing Networks in Manufacturing and Road Traffic Yoni Nazarathy Swinburne University of Technology, Melbourne. Joint work with: Erjen Lefeber (manufacturing), Hai
More informationStability and Asymptotic Optimality of h-maxweight Policies
Stability and Asymptotic Optimality of h-maxweight Policies - A. Rybko, 2006 Sean Meyn Department of Electrical and Computer Engineering University of Illinois & the Coordinated Science Laboratory NSF
More informationA Push Pull Network with Infinite Supply of Work
DOI 10.1007/s11134-009-9121-1 A Push Pull Network with Infinite Supply of Work Anat Kopzon Yoni Nazarathy Gideon Weiss Received: 31 January 2007 / Revised: 28 July 2008 Springer Science+Business Media,
More informationDynamic Matching Models
Dynamic Matching Models Ana Bušić Inria Paris - Rocquencourt CS Department of École normale supérieure joint work with Varun Gupta, Jean Mairesse and Sean Meyn 3rd Workshop on Cognition and Control January
More informationScheduling Multiclass Queueing Networks via Fluid Models
Scheduling Multiclass Queueing Networks via Fluid Models John Hasenbein OR/IE, UT-Austin Ron Billings OR/IE, UT-Austin Leon Lasdon MSIS, UT-Austin Gideon Weiss Statistics, Univ of Haifa www.me.utexas.edu/~has
More informationA fluid approach to job shop scheduling: theory, software and experimentation
A fluid approach to job shop scheduling: theory, software and experimentation Yoni Nazarathy, TU/e Eindhoven, The Department of Mechanical Engineering and EURANDOM, Gideon Weiss, The University of Haifa,
More informationt dt t SCLP Bellman (1953) CLP (Dantzig, Tyndall, Grinold, Perold, Anstreicher 60's-80's) Anderson (1978) SCLP
Coninuous Linear Programming. Separaed Coninuous Linear Programming Bellman (1953) max c () u() d H () u () + Gsusds (,) () a () u (), < < CLP (Danzig, yndall, Grinold, Perold, Ansreicher 6's-8's) Anderson
More informationMIT Manufacturing Systems Analysis Lectures 19 21
MIT 2.852 Manufacturing Systems Analysis Lectures 19 21 Scheduling: Real-Time Control of Manufacturing Systems Stanley B. Gershwin Spring, 2007 Copyright c 2007 Stanley B. Gershwin. Definitions Events
More informationOperations Research Letters. Instability of FIFO in a simple queueing system with arbitrarily low loads
Operations Research Letters 37 (2009) 312 316 Contents lists available at ScienceDirect Operations Research Letters journal homepage: www.elsevier.com/locate/orl Instability of FIFO in a simple queueing
More informationOn the Resource/Performance Tradeoff in Large Scale Queueing Systems
On the Resource/Performance Tradeoff in Large Scale Queueing Systems David Gamarnik MIT Joint work with Patrick Eschenfeldt, John Tsitsiklis and Martin Zubeldia (MIT) High level comments High level comments
More informationOptimeringslära för F (SF1811) / Optimization (SF1841)
Optimeringslära för F (SF1811) / Optimization (SF1841) 1. Information about the course 2. Examples of optimization problems 3. Introduction to linear programming Introduction - Per Enqvist 1 Linear programming
More informationLinear Programming. H. R. Alvarez A., Ph. D. 1
Linear Programming H. R. Alvarez A., Ph. D. 1 Introduction It is a mathematical technique that allows the selection of the best course of action defining a program of feasible actions. The objective of
More informationYinyu Ye, MS&E, Stanford MS&E310 Lecture Note #06. The Simplex Method
The Simplex Method Yinyu Ye Department of Management Science and Engineering Stanford University Stanford, CA 94305, U.S.A. http://www.stanford.edu/ yyye (LY, Chapters 2.3-2.5, 3.1-3.4) 1 Geometry of Linear
More informationSupply chains models
Supply chains models Maria Ivanova Technische Universität München, Department of Mathematics, Haupt Seminar 2016 Garching, Germany July 3, 2016 Overview Introduction What is Supply chain? Traffic flow
More informationConvexity Properties of Loss and Overflow Functions
Convexity Properties of Loss and Overflow Functions Krishnan Kumaran?, Michel Mandjes y, and Alexander Stolyar? email: kumaran@lucent.com, michel@cwi.nl, stolyar@lucent.com? Bell Labs/Lucent Technologies,
More informationLinear Programming: Chapter 5 Duality
Linear Programming: Chapter 5 Duality Robert J. Vanderbei September 30, 2010 Slides last edited on October 5, 2010 Operations Research and Financial Engineering Princeton University Princeton, NJ 08544
More information1 Markov decision processes
2.997 Decision-Making in Large-Scale Systems February 4 MI, Spring 2004 Handout #1 Lecture Note 1 1 Markov decision processes In this class we will study discrete-time stochastic systems. We can describe
More informationWritten Examination
Division of Scientific Computing Department of Information Technology Uppsala University Optimization Written Examination 202-2-20 Time: 4:00-9:00 Allowed Tools: Pocket Calculator, one A4 paper with notes
More information(P ) Minimize 4x 1 + 6x 2 + 5x 3 s.t. 2x 1 3x 3 3 3x 2 2x 3 6
The exam is three hours long and consists of 4 exercises. The exam is graded on a scale 0-25 points, and the points assigned to each question are indicated in parenthesis within the text. Problem 1 Consider
More informationMODELING WEBCHAT SERVICE CENTER WITH MANY LPS SERVERS
MODELING WEBCHAT SERVICE CENTER WITH MANY LPS SERVERS Jiheng Zhang Oct 26, 211 Model and Motivation Server Pool with multiple LPS servers LPS Server K Arrival Buffer. Model and Motivation Server Pool with
More informationA Note on the Event Horizon for a Processor Sharing Queue
A Note on the Event Horizon for a Processor Sharing Queue Robert C. Hampshire Heinz School of Public Policy and Management Carnegie Mellon University hamp@andrew.cmu.edu William A. Massey Department of
More informationDynamic Power Allocation and Routing for Time Varying Wireless Networks
Dynamic Power Allocation and Routing for Time Varying Wireless Networks X 14 (t) X 12 (t) 1 3 4 k a P ak () t P a tot X 21 (t) 2 N X 2N (t) X N4 (t) µ ab () rate µ ab µ ab (p, S 3 ) µ ab µ ac () µ ab (p,
More informationOPEN MULTICLASS HL QUEUEING NETWORKS: PROGRESS AND SURPRISES OF THE PAST 15 YEARS. w 1. v 2. v 3. Ruth J. Williams University of California, San Diego
OPEN MULTICLASS HL QUEUEING NETWORKS: PROGRESS AND SURPRISES OF THE PAST 15 YEARS v 2 w3 v1 w 2 w 1 v 3 Ruth J. Williams University of California, San Diego 1 PERSPECTIVE MQN SPN Sufficient conditions
More informationMATH2070 Optimisation
MATH2070 Optimisation Linear Programming Semester 2, 2012 Lecturer: I.W. Guo Lecture slides courtesy of J.R. Wishart Review The standard Linear Programming (LP) Problem Graphical method of solving LP problem
More informationThe Strong Duality Theorem 1
1/39 The Strong Duality Theorem 1 Adrian Vetta 1 This presentation is based upon the book Linear Programming by Vasek Chvatal 2/39 Part I Weak Duality 3/39 Primal and Dual Recall we have a primal linear
More informationTutorial: Optimal Control of Queueing Networks
Department of Mathematics Tutorial: Optimal Control of Queueing Networks Mike Veatch Presented at INFORMS Austin November 7, 2010 1 Overview Network models MDP formulations: features, efficient formulations
More informationMaximum Pressure Policies in Stochastic Processing Networks
OPERATIONS RESEARCH Vol. 53, No. 2, March April 2005, pp. 197 218 issn 0030-364X eissn 1526-5463 05 5302 0197 informs doi 10.1287/opre.1040.0170 2005 INFORMS Maximum Pressure Policies in Stochastic Processing
More informationSTABILITY OF MULTICLASS QUEUEING NETWORKS UNDER LONGEST-QUEUE AND LONGEST-DOMINATING-QUEUE SCHEDULING
Applied Probability Trust (7 May 2015) STABILITY OF MULTICLASS QUEUEING NETWORKS UNDER LONGEST-QUEUE AND LONGEST-DOMINATING-QUEUE SCHEDULING RAMTIN PEDARSANI and JEAN WALRAND, University of California,
More informationMATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem
MATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem Pulemotov, September 12, 2012 Unit Outline Goal 1: Outline linear
More informationManagement of demand-driven production systems
Management of demand-driven production systems Mike Chen, Richard Dubrawski, and Sean Meyn November 4, 22 Abstract Control-synthesis techniques are developed for demand-driven production systems. The resulting
More informationBRAVO for QED Queues
1 BRAVO for QED Queues Yoni Nazarathy, The University of Queensland Joint work with Daryl J. Daley, The University of Melbourne, Johan van Leeuwaarden, EURANDOM, Eindhoven University of Technology. Applied
More informationMidterm Review. Yinyu Ye Department of Management Science and Engineering Stanford University Stanford, CA 94305, U.S.A.
Midterm Review Yinyu Ye Department of Management Science and Engineering Stanford University Stanford, CA 94305, U.S.A. http://www.stanford.edu/ yyye (LY, Chapter 1-4, Appendices) 1 Separating hyperplane
More informationSection Notes 9. Midterm 2 Review. Applied Math / Engineering Sciences 121. Week of December 3, 2018
Section Notes 9 Midterm 2 Review Applied Math / Engineering Sciences 121 Week of December 3, 2018 The following list of topics is an overview of the material that was covered in the lectures and sections
More informationHyperbolic Models for Large Supply Chains. Christian Ringhofer (Arizona State University) Hyperbolic Models for Large Supply Chains p.
Hyperbolic Models for Large Supply Chains Christian Ringhofer (Arizona State University) Hyperbolic Models for Large Supply Chains p. /4 Introduction Topic: Overview of conservation law (traffic - like)
More informationStability and Rare Events in Stochastic Models Sergey Foss Heriot-Watt University, Edinburgh and Institute of Mathematics, Novosibirsk
Stability and Rare Events in Stochastic Models Sergey Foss Heriot-Watt University, Edinburgh and Institute of Mathematics, Novosibirsk ANSAPW University of Queensland 8-11 July, 2013 1 Outline (I) Fluid
More informationPart 1. The Review of Linear Programming
In the name of God Part 1. The Review of Linear Programming 1.5. Spring 2010 Instructor: Dr. Masoud Yaghini Outline Introduction Formulation of the Dual Problem Primal-Dual Relationship Economic Interpretation
More informationEnergy Optimal Control for Time Varying Wireless Networks. Michael J. Neely University of Southern California
Energy Optimal Control for Time Varying Wireless Networks Michael J. Neely University of Southern California http://www-rcf.usc.edu/~mjneely Part 1: A single wireless downlink (L links) L 2 1 S={Totally
More informationECE7850 Lecture 9. Model Predictive Control: Computational Aspects
ECE785 ECE785 Lecture 9 Model Predictive Control: Computational Aspects Model Predictive Control for Constrained Linear Systems Online Solution to Linear MPC Multiparametric Programming Explicit Linear
More informationDeterministic Dynamic Programming
Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0
More informationIntroduction to Mathematical Programming IE406. Lecture 10. Dr. Ted Ralphs
Introduction to Mathematical Programming IE406 Lecture 10 Dr. Ted Ralphs IE406 Lecture 10 1 Reading for This Lecture Bertsimas 4.1-4.3 IE406 Lecture 10 2 Duality Theory: Motivation Consider the following
More informationChapter 1 Linear Programming. Paragraph 5 Duality
Chapter 1 Linear Programming Paragraph 5 Duality What we did so far We developed the 2-Phase Simplex Algorithm: Hop (reasonably) from basic solution (bs) to bs until you find a basic feasible solution
More information4. Duality and Sensitivity
4. Duality and Sensitivity For every instance of an LP, there is an associated LP known as the dual problem. The original problem is known as the primal problem. There are two de nitions of the dual pair
More informationModeling of Material Flow Problems
Modeling of Material Flow Problems Simone Göttlich Department of Mathematics University of Mannheim Workshop on Math for the Digital Factory, WIAS Berlin May 7-9, 2014 Prof. Dr. Simone Göttlich Material
More information[4] T. I. Seidman, \\First Come First Serve" is Unstable!," tech. rep., University of Maryland Baltimore County, 1993.
[2] C. J. Chase and P. J. Ramadge, \On real-time scheduling policies for exible manufacturing systems," IEEE Trans. Automat. Control, vol. AC-37, pp. 491{496, April 1992. [3] S. H. Lu and P. R. Kumar,
More information56:171 Operations Research Midterm Exam - October 26, 1989 Instructor: D.L. Bricker
56:171 Operations Research Midterm Exam - October 26, 1989 Instructor: D.L. Bricker Answer all of Part One and two (of the four) problems of Part Two Problem: 1 2 3 4 5 6 7 8 TOTAL Possible: 16 12 20 10
More informationOn the Flow-level Dynamics of a Packet-switched Network
On the Flow-level Dynamics of a Packet-switched Network Ciamac Moallemi Graduate School of Business Columbia University ciamac@gsb.columbia.edu Devavrat Shah LIDS, EECS Massachusetts Institute of Technology
More informationScheduling: Queues & Computation
Scheduling: Queues Computation achieving baseline performance efficiently Devavrat Shah LIDS, MIT Outline Two models switched network and bandwidth sharing Scheduling: desirable performance queue-size
More informationThe Dual Simplex Algorithm
p. 1 The Dual Simplex Algorithm Primal optimal (dual feasible) and primal feasible (dual optimal) bases The dual simplex tableau, dual optimality and the dual pivot rules Classical applications of linear
More informationDEPARTMENT OF STATISTICS AND OPERATIONS RESEARCH OPERATIONS RESEARCH DETERMINISTIC QUALIFYING EXAMINATION. Part I: Short Questions
DEPARTMENT OF STATISTICS AND OPERATIONS RESEARCH OPERATIONS RESEARCH DETERMINISTIC QUALIFYING EXAMINATION Part I: Short Questions August 12, 2008 9:00 am - 12 pm General Instructions This examination is
More informationDynamic Scheduling of Multiclass Queueing Networks. Caiwei Li
Dynamic Scheduling of Multiclass Queueing Networks A Thesis Presented to The Academic Faculty by Caiwei Li In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy in Industrial
More information56:171 Operations Research Fall 1998
56:171 Operations Research Fall 1998 Quiz Solutions D.L.Bricker Dept of Mechanical & Industrial Engineering University of Iowa 56:171 Operations Research Quiz
More informationSubject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)
Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must
More informationUNIT-4 Chapter6 Linear Programming
UNIT-4 Chapter6 Linear Programming Linear Programming 6.1 Introduction Operations Research is a scientific approach to problem solving for executive management. It came into existence in England during
More informationLecture 10: Linear programming. duality. and. The dual of the LP in standard form. maximize w = b T y (D) subject to A T y c, minimize z = c T x (P)
Lecture 10: Linear programming duality Michael Patriksson 19 February 2004 0-0 The dual of the LP in standard form minimize z = c T x (P) subject to Ax = b, x 0 n, and maximize w = b T y (D) subject to
More informationCOT 6936: Topics in Algorithms! Giri Narasimhan. ECS 254A / EC 2443; Phone: x3748
COT 6936: Topics in Algorithms! Giri Narasimhan ECS 254A / EC 2443; Phone: x3748 giri@cs.fiu.edu https://moodle.cis.fiu.edu/v2.1/course/view.php?id=612 Gaussian Elimination! Solving a system of simultaneous
More informationAllocating Resources, in the Future
Allocating Resources, in the Future Sid Banerjee School of ORIE May 3, 2018 Simons Workshop on Mathematical and Computational Challenges in Real-Time Decision Making online resource allocation: basic model......
More information4.6 Linear Programming duality
4.6 Linear Programming duality To any minimization (maximization) LP we can associate a closely related maximization (minimization) LP Different spaces and objective functions but in general same optimal
More informationTales of Time Scales. Ward Whitt AT&T Labs Research Florham Park, NJ
Tales of Time Scales Ward Whitt AT&T Labs Research Florham Park, NJ New Book Stochastic-Process Limits An Introduction to Stochastic-Process Limits and Their Application to Queues Springer 2001 I won t
More informationA Shadow Simplex Method for Infinite Linear Programs
A Shadow Simplex Method for Infinite Linear Programs Archis Ghate The University of Washington Seattle, WA 98195 Dushyant Sharma The University of Michigan Ann Arbor, MI 48109 May 25, 2009 Robert L. Smith
More informationSolving Dual Problems
Lecture 20 Solving Dual Problems We consider a constrained problem where, in addition to the constraint set X, there are also inequality and linear equality constraints. Specifically the minimization problem
More information56:270 Final Exam - May
@ @ 56:270 Linear Programming @ @ Final Exam - May 4, 1989 @ @ @ @ @ @ @ @ @ @ @ @ @ @ Select any 7 of the 9 problems below: (1.) ANALYSIS OF MPSX OUTPUT: Please refer to the attached materials on the
More information(includes both Phases I & II)
Minimize z=3x 5x 4x 7x 5x 4x subject to 2x x2 x4 3x6 0 x 3x3 x4 3x5 2x6 2 4x2 2x3 3x4 x5 5 and x 0 j, 6 2 3 4 5 6 j ecause of the lack of a slack variable in each constraint, we must use Phase I to find
More informationmin 4x 1 5x 2 + 3x 3 s.t. x 1 + 2x 2 + x 3 = 10 x 1 x 2 6 x 1 + 3x 2 + x 3 14
The exam is three hours long and consists of 4 exercises. The exam is graded on a scale 0-25 points, and the points assigned to each question are indicated in parenthesis within the text. If necessary,
More informationLecture: Algorithms for LP, SOCP and SDP
1/53 Lecture: Algorithms for LP, SOCP and SDP Zaiwen Wen Beijing International Center For Mathematical Research Peking University http://bicmr.pku.edu.cn/~wenzw/bigdata2018.html wenzw@pku.edu.cn Acknowledgement:
More information5. Duality. Lagrangian
5. Duality Convex Optimization Boyd & Vandenberghe Lagrange dual problem weak and strong duality geometric interpretation optimality conditions perturbation and sensitivity analysis examples generalized
More informationMATH 373 Section A1. Final Exam. Dr. J. Bowman 17 December :00 17:00
MATH 373 Section A1 Final Exam Dr. J. Bowman 17 December 2018 14:00 17:00 Name (Last, First): Student ID: Email: @ualberta.ca Scrap paper is supplied. No notes or books are permitted. All electronic equipment,
More informationChapter 1: Linear Programming
Chapter 1: Linear Programming Math 368 c Copyright 2013 R Clark Robinson May 22, 2013 Chapter 1: Linear Programming 1 Max and Min For f : D R n R, f (D) = {f (x) : x D } is set of attainable values of
More informationResource Pooling for Optimal Evacuation of a Large Building
Proceedings of the 47th IEEE Conference on Decision and Control Cancun, Mexico, Dec. 9-11, 28 Resource Pooling for Optimal Evacuation of a Large Building Kun Deng, Wei Chen, Prashant G. Mehta, and Sean
More informationAdvanced Computer Networks Lecture 3. Models of Queuing
Advanced Computer Networks Lecture 3. Models of Queuing Husheng Li Min Kao Department of Electrical Engineering and Computer Science University of Tennessee, Knoxville Spring, 2016 1/13 Terminology of
More informationLecture 9: Deterministic Fluid Models and Many-Server Heavy-Traffic Limits. IEOR 4615: Service Engineering Professor Whitt February 19, 2015
Lecture 9: Deterministic Fluid Models and Many-Server Heavy-Traffic Limits IEOR 4615: Service Engineering Professor Whitt February 19, 2015 Outline Deterministic Fluid Models Directly From Data: Cumulative
More informationA Strongly Polynomial Simplex Method for Totally Unimodular LP
A Strongly Polynomial Simplex Method for Totally Unimodular LP Shinji Mizuno July 19, 2014 Abstract Kitahara and Mizuno get new bounds for the number of distinct solutions generated by the simplex method
More informationInfinite-dimensional nonlinear predictive controller design for open-channel hydraulic systems
Infinite-dimensional nonlinear predictive controller design for open-channel hydraulic systems D. Georges, Control Systems Dept - Gipsa-lab, Grenoble INP Workshop on Irrigation Channels and Related Problems,
More informationOn the equivalence of economic lot scheduling and switched production systems
On the equivalence of economic lot scheduling and switched production systems K.G.M. Jacobs, Dieter Armbruster,, Erjen Lefeber and J.E. Rooda January, 2 Abstract Scheduling production of several product
More informationLecture 2: The Simplex method. 1. Repetition of the geometrical simplex method. 2. Linear programming problems on standard form.
Lecture 2: The Simplex method. Repetition of the geometrical simplex method. 2. Linear programming problems on standard form. 3. The Simplex algorithm. 4. How to find an initial basic solution. Lecture
More informationDuality Theory, Optimality Conditions
5.1 Duality Theory, Optimality Conditions Katta G. Murty, IOE 510, LP, U. Of Michigan, Ann Arbor We only consider single objective LPs here. Concept of duality not defined for multiobjective LPs. Every
More informationΩ R n is called the constraint set or feasible set. x 1
1 Chapter 5 Linear Programming (LP) General constrained optimization problem: minimize subject to f(x) x Ω Ω R n is called the constraint set or feasible set. any point x Ω is called a feasible point We
More informationConvex Optimization of Graph Laplacian Eigenvalues
Convex Optimization of Graph Laplacian Eigenvalues Stephen Boyd Stanford University (Joint work with Persi Diaconis, Arpita Ghosh, Seung-Jean Kim, Sanjay Lall, Pablo Parrilo, Amin Saberi, Jun Sun, Lin
More informationAlgorithms and Theory of Computation. Lecture 13: Linear Programming (2)
Algorithms and Theory of Computation Lecture 13: Linear Programming (2) Xiaohui Bei MAS 714 September 25, 2018 Nanyang Technological University MAS 714 September 25, 2018 1 / 15 LP Duality Primal problem
More informationSummary of the simplex method
MVE165/MMG630, The simplex method; degeneracy; unbounded solutions; infeasibility; starting solutions; duality; interpretation Ann-Brith Strömberg 2012 03 16 Summary of the simplex method Optimality condition:
More informationLinear conic optimization for nonlinear optimal control
Linear conic optimization for nonlinear optimal control Didier Henrion 1,2,3, Edouard Pauwels 1,2 Draft of July 15, 2014 Abstract Infinite-dimensional linear conic formulations are described for nonlinear
More informationAsymptotically Optimal Inventory Control For Assemble-to-Order Systems
Asymptotically Optimal Inventory Control For Assemble-to-Order Systems Marty Reiman Columbia Univerisity joint work with Mustafa Dogru, Haohua Wan, and Qiong Wang May 16, 2018 Outline The Assemble-to-Order
More informationLinear and Combinatorial Optimization
Linear and Combinatorial Optimization The dual of an LP-problem. Connections between primal and dual. Duality theorems and complementary slack. Philipp Birken (Ctr. for the Math. Sc.) Lecture 3: Duality
More informationQUALIFYING EXAM IN SYSTEMS ENGINEERING
QUALIFYING EXAM IN SYSTEMS ENGINEERING Written Exam: MAY 23, 2017, 9:00AM to 1:00PM, EMB 105 Oral Exam: May 25 or 26, 2017 Time/Location TBA (~1 hour per student) CLOSED BOOK, NO CHEAT SHEETS BASIC SCIENTIFIC
More informationMinimum cost transportation problem
Minimum cost transportation problem Complements of Operations Research Giovanni Righini Università degli Studi di Milano Definitions The minimum cost transportation problem is a special case of the minimum
More informationMS-E2140. Lecture 1. (course book chapters )
Linear Programming MS-E2140 Motivations and background Lecture 1 (course book chapters 1.1-1.4) Linear programming problems and examples Problem manipulations and standard form problems Graphical representation
More informationAn introductory example
CS1 Lecture 9 An introductory example Suppose that a company that produces three products wishes to decide the level of production of each so as to maximize profits. Let x 1 be the amount of Product 1
More informationControl of Fork-Join Networks in Heavy-Traffic
in Heavy-Traffic Asaf Zviran Based on MSc work under the guidance of Rami Atar (Technion) and Avishai Mandelbaum (Technion) Industrial Engineering and Management Technion June 2010 Introduction Network
More informationA Parametric Simplex Algorithm for Linear Vector Optimization Problems
A Parametric Simplex Algorithm for Linear Vector Optimization Problems Birgit Rudloff Firdevs Ulus Robert Vanderbei July 9, 2015 Abstract In this paper, a parametric simplex algorithm for solving linear
More information1. Algebraic and geometric treatments Consider an LP problem in the standard form. x 0. Solutions to the system of linear equations
The Simplex Method Most textbooks in mathematical optimization, especially linear programming, deal with the simplex method. In this note we study the simplex method. It requires basically elementary linear
More informationCO759: Algorithmic Game Theory Spring 2015
CO759: Algorithmic Game Theory Spring 2015 Instructor: Chaitanya Swamy Assignment 1 Due: By Jun 25, 2015 You may use anything proved in class directly. I will maintain a FAQ about the assignment on the
More informationStochastic process. X, a series of random variables indexed by t
Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,
More informationAsymptotic Coupling of an SPDE, with Applications to Many-Server Queues
Asymptotic Coupling of an SPDE, with Applications to Many-Server Queues Mohammadreza Aghajani joint work with Kavita Ramanan Brown University March 2014 Mohammadreza Aghajanijoint work Asymptotic with
More information