Announcements. Queueing Systems: Lecture 1. Lecture Outline. Topics in Queueing Theory

Size: px
Start display at page:

Download "Announcements. Queueing Systems: Lecture 1. Lecture Outline. Topics in Queueing Theory"

Transcription

1 Aoucemets Queueig Systems: Lecture Amedeo R. Odoi October 4, 2006 PS #3 out this afteroo Due: October 9 (graded by 0/23) Office hours Odoi: Mo. 2:30-4:30 - Wed. 2:30-4:30 o Oct. 8 (No office hrs 0/6) _ Or sed me a message Quiz #: October 25, ope book, i class Old quiz problems ad solutios: posted o 0/9 Topics i Queueig Theory Itroductio to Queues Little s Law Markovia Birth-ad-Death Queues The M/M/ ad Other Markovia Variatios The M/G/ Queue ad Extesios Priority Queues Some Useful Bouds ogestio Pricig Queueig Networks; State Represetatios Dyamic Behavior of Queues Lecture Outlie Itroductio to queueig s oceptual represetatio of queueig s odes for queueig models Termiology ad otatio Little s Law ad basic relatioships Birth-ad-death models The M/M/ queueig Referece: hapter 4, pp

2 Queues A Geeric Queueig System Queueig theory is the brach of operatios research cocered with waitig lies (delays/cogestio) A queueig cosists of a user source, a queue ad a service facility with oe or more idetical parallel servers A queueig etwork is a set of itercoected queueig s Fudametal parameters of a queueig : - Demad rate - apacity (service rate) - Demad iter-arrival times - Service times - Queue capacity ad disciplie (fiite vs. ifiite; FIFO/FFS, SIRO, LIFO, priorities) - Myriad details (feedback effects, balkig, jockeyig, etc.) Arrival poit at the Source of users/ customers Size of user source Arrivals process Queue Queue disciplie ad Queue capacity Service process Servers Number of servers Departure poit from the Queueig etwork cosistig of five queueig s I Queueig Queueig 2 Poit where users make a choice Queueig 4 Queueig 3 Poit where users merge + Queueig 5 Out Applicatios of Queueig Theory Some familiar queues: _ Airport check-i; aircraft i a holdig patter _ Automated Teller Machies (ATMs) _ Fast food restaurats _ Phoe ceter s lies _ Urba itersectio _ Toll booths _ Spatially distributed urba s ad services Level-of-service (LOS) stadards Ecoomic aalyses ivolvig trade-offs amog operatig costs, capital ivestmets ad LOS ogestio pricig

3 The Airside as a Queueig Network Queueig Models a Be Essetial i Aalysis of apital Ivestmets ost Total cost Optimal cost ost of buildig the capacity ost of losses due to waitig Optimal capacity Airport apacity Stregths ad Weakesses of Queueig Theory Queueig models ecessarily ivolve approximatios ad simplificatio of reality Results give a sese of order of magitude, chages relative to a baselie, promisig directios i which to move losed-form results essetially limited to steady state coditios ad derived primarily (but ot solely) for birth-ad-death s ad phase s Some useful bouds for more geeral s at steady state Numerical solutios icreasigly viable for dyamic s Huge umber of importat applicatios A ode for Queueig Models: A/B/m Distributio of service time Number of servers Queueig System / / S Queue ustomers S Service S facility S Distributio of iterarrival time Some stadard code letters for A ad B: _ M: Negative expoetial (M stads for memoryless) _ D: Determiistic _ E k :kth-order Erlag distributio _ G: Geeral distributio

4 Termiology ad Notatio Number i : umber of customers i queueig Number i queue or Queue legth : umber of customers waitig for service Total time i ad waitig time N(t) = umber of customers i queueig at time t P (t) = probability that N(t) is equal to at time t λ : mea arrival rate of ew customers whe N(t) = μ : mea (total) service rate whe N(t) = Termiology ad Notatio (2) Trasiet state: state of at t is iflueced by the state of the at t = 0 Steady state: state of the is idepedet of iitial state of the m: umber of servers (parallel service chaels) If λ ad the service rate per busy server are costats λ ad μ, respectively, the λ =λ, μ = mi (μ, mμ); i that case: _ Expected iter-arrival time = /λ _ Expected service time = /μ Some Expected Values of Iterest at Steady State Little s Law Give: _ λ = arrival rate _ μ = service rate per service chael Ukows: _ L = expected umber of users i queueig _ L q = expected umber of users i queue _ W = expected time i queueig per user (W = E(w)) _ W q = expected waitig time i queue per user (W q = E(w q )) 4 ukows We eed 4 equatios Number of users A(t): cumulative arrivals to the (t): cumulative service completios i the A(t) N(t) (t) t T T T N ( t) A T N ( t) 0 ( ) 0 L T = = = λt WT T T A ( T ) Time

5 Relatioships amog L, L q, W, W q Four ukows: L, W, L q, W q Need 4 equatios. We have the followig 3 equatios: _ L = λw (Little s law) _ L q = λw q _ W = W q + μ If we ca fid ay oe of the four expected values, we ca determie the three others The determiatio of L (or other) may be hard or easy depedig o the type of queueig at had L = P (P :probability that customers are i the ) =0 Birth-ad-Death Queueig Systems. m parallel, idetical servers. 2. Ifiite queue capacity (for ow). 3. Wheever users are i (i queue plus i service) arrivals are Poisso at rate of λ per uit of time. 4. Wheever users are i, service completios are Poisso at rate of μ per uit of time. 5. FFS disciplie (for ow). Time: t users The Fudametal Relatioship μ Δt λ Δt -(λ + μ )Δt Time: t+δt + users users - users P (t) = Prob [ users i at time t] P (t + Δt) = P + (t) μ + Δt + P (t) λ Δt + P (t) [ (μ + λ ) Δt] The differetial equatios that determie the state probabilities P (t + Δt) = P + (t) μ + Δt + P (t) λ Δt + P (t) [ (μ + λ ) Δt] After a simple maipulatio: dp (t) = (λ + μ ) P (t) + λ P (t) + μ + P + (t) () () applies whe =, 2, 3,.; whe = 0, we have: dp 0 (t) = λ0 P 0 (t) + μ P (t) (2) The of equatios () ad (2) is kow as the hapma-kolmogorov equatios for a birth-ad-death

6 The state balace equatios Birth-ad-Death System: State Trasitio Diagram We ow cosider the situatio i which the queueig λ 0 λ λ 2 λ m- λ m λ m+ has reached steady state, i.e., t is large eough to have P (t) = P, idepedet of t, or dp (t) = m m+ The, () ad (2) provide the state balace equatios: μ μ 2 μ μ m μ m+ μ m+2 3 λ 0 P 0 = μ P = 0 (3) We are iterested i the characteristics of the (λ + μ ) P = λ P + μ + P + =, 2, 3,.. (4) uder equilibrium coditios ( steady state ), i.e., whe The state balace equatios ca also be writte directly the state probabilities P (t) are idepedet of t for from the state trasitio diagram large values of t a write balace equatios ad obtai closed form expressios for P, L, W, L q, W q Solvig.. M/M/: Observig State Trasitio Diagram from Two Poits Solvig (3) ad (4), we have: From poit : λp 0 = μp (λ + μ)p = λp 0 + μp 2 (λ + μ)p = λp + μp + P = λ 0 P0 ; P 2 = λ P = λ λ 0 P0 etc. μ μ 2 μ 2 μ λ λ λ λ λ λ λ ad, i geeral, λ λ 2... λ λ 0 μ μ μ μ μ μ μ P = P 0 = K P 0 μ μ... μ 2 μ But, we also have: = P = P 0 ( + K ) =0 = Givig, P 0 = + K = From poit 2: oditio for steady state: K < = λp 0 = μp P μp2 λ λ λ λ λ λ λ μ λ = λ P μp = + μ μ μ μ μ μ

7 M/M/: Derivatio of P 0 ad P M/M/: Derivatio of L, W, W q, ad L q Step : P = λ P 0, P 2 = λ 2 P 0, L, P = λ μ μ μ Step 2: P =, P 0 λ = P = =0 =0 μ 0 P 0 λ =0 μ Step 3: ρ = λ, the λ = ρ = ρ = (Q ρ < ) μ =0 μ =0 ρ ρ Step 4: P 0 = ρ =0 = ρ ad P = ρ ( ρ) L = P = ρ ( ρ) = ( ρ) ρ = ( ρ)ρ ρ =0 =0 =0 = = ( ρ)ρ d ρ = ( ρ)ρ d dρ =0 dρ ρ λ = ( ρ)ρ ρ μ λ = = = ( ρ) 2 ( ρ) λ μ λ μ L λ W = = = λ μ λ λ μ λ λ W q = W = = μ μ λ μ μ(μ λ) λ λ 2 L q = λw q = λ = μ(μ λ) μ(μ λ) High Sesitivity of Delay at High Levels of Utilizatio Expected delay M/M/: A alterative, direct derivatio of L ad W For a M/M/, with FFS disciplie: apacity ρ = Demad W = ( +) P = E[ N + ] = E[N ] + = L + () =0 μ μ μ μ But from Little s theorem we also have: L = λ W (2) It follows from () ad (2) that, as before: λ L = ; W = μ λ μ λ Does the queueig disciplie matter?

8 Additioal importat M/M/ results The pdf for the total time i the, w, ca be computed for a M/M/ (ad FFS): f w (w) = ( ρ )μe ( ρ )μw = (μ λ )e (μ λ )w for w 0 Thus, as already show, W = /(μ -λ) = /[μ (-ρ)] The stadard deviatio of N, w, N q, w q are all proportioal to /(-ρ), just like their expected values (L, W, L q, W q, respectively) The expected legth of the busy period, E[B], is equal to /(μ -λ) M/M/: E[B], the expected legth of a busy period N t I B I B I B P 0 = B = busy period I = idle period E[legth Idle period ] E[legth Busy period ] + E[legth Idle period ] But, P 0 = ρ E[legth Idle period] = λ Therefore, E[B] = E[legth Busy period] = = μ ( ρ ) μ λ

B. Maddah ENMG 622 ENMG /27/07

B. Maddah ENMG 622 ENMG /27/07 B. Maddah ENMG 622 ENMG 5 3/27/7 Queueig Theory () What is a queueig system? A queueig system cosists of servers (resources) that provide service to customers (etities). A Customer requestig service will

More information

Queuing Theory. Basic properties, Markovian models, Networks of queues, General service time distributions, Finite source models, Multiserver queues

Queuing Theory. Basic properties, Markovian models, Networks of queues, General service time distributions, Finite source models, Multiserver queues Queuig Theory Basic properties, Markovia models, Networks of queues, Geeral service time distributios, Fiite source models, Multiserver queues Chapter 8 Kedall s Notatio for Queuig Systems A/B/X/Y/Z: A

More information

First come, first served (FCFS) Batch

First come, first served (FCFS) Batch Queuig Theory Prelimiaries A flow of customers comig towards the service facility forms a queue o accout of lack of capacity to serve them all at a time. RK Jaa Some Examples: Persos waitig at doctor s

More information

CS/ECE 715 Spring 2004 Homework 5 (Due date: March 16)

CS/ECE 715 Spring 2004 Homework 5 (Due date: March 16) CS/ECE 75 Sprig 004 Homework 5 (Due date: March 6) Problem 0 (For fu). M/G/ Queue with Radom-Sized Batch Arrivals. Cosider the M/G/ system with the differece that customers are arrivig i batches accordig

More information

Simulation of Discrete Event Systems

Simulation of Discrete Event Systems Simulatio of Discrete Evet Systems Uit 9 Queueig Models Fall Witer 2014/2015 Uiv.-Prof. Dr.-Ig. Dipl.-Wirt.-Ig. Christopher M. Schlick Chair ad Istitute of Idustrial Egieerig ad Ergoomics RWTH Aache Uiversity

More information

A queueing system can be described as customers arriving for service, waiting for service if it is not immediate, and if having waited for service,

A queueing system can be described as customers arriving for service, waiting for service if it is not immediate, and if having waited for service, Queuig System A queueig system ca be described as customers arrivig for service, waitig for service if it is ot immediate, ad if havig waited for service, leavig the service after beig served. The basic

More information

Reliability and Queueing

Reliability and Queueing Copyright 999 Uiversity of Califoria Reliability ad Queueig by David G. Messerschmitt Supplemetary sectio for Uderstadig Networked Applicatios: A First Course, Morga Kaufma, 999. Copyright otice: Permissio

More information

TCOM 501: Networking Theory & Fundamentals. Lecture 3 January 29, 2003 Prof. Yannis A. Korilis

TCOM 501: Networking Theory & Fundamentals. Lecture 3 January 29, 2003 Prof. Yannis A. Korilis TCOM 5: Networkig Theory & Fudametals Lecture 3 Jauary 29, 23 Prof. Yais A. Korilis 3-2 Topics Markov Chais Discrete-Time Markov Chais Calculatig Statioary Distributio Global Balace Equatios Detailed Balace

More information

TUFTS UNIVERSITY DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING ES 152 ENGINEERING SYSTEMS Spring Lesson 14. Queue System Theory

TUFTS UNIVERSITY DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING ES 152 ENGINEERING SYSTEMS Spring Lesson 14. Queue System Theory TUFTS UNIVERSITY DEARTMENT OF CIVI AND ENVIRONMENTA ENGINEERING ES 52 ENGINEERING SYSTEMS Sprig 2 esso 4 Queue System Theory There exists a cosiderable body of theoretical aalysis of ueues. (Chapter 7

More information

On Poisson Bulk Arrival Queue: M / M /2/ N with. Balking, Reneging and Heterogeneous servers

On Poisson Bulk Arrival Queue: M / M /2/ N with. Balking, Reneging and Heterogeneous servers Applied Mathematical Scieces, Vol., 008, o. 4, 69-75 O oisso Bulk Arrival Queue: M / M // with Balkig, Reegig ad Heterogeeous servers M. S. El-aoumy Statistics Departmet, Faculty of Commerce, Dkhlia, Egypt

More information

Generalized Semi- Markov Processes (GSMP)

Generalized Semi- Markov Processes (GSMP) Geeralized Semi- Markov Processes (GSMP) Summary Some Defiitios Markov ad Semi-Markov Processes The Poisso Process Properties of the Poisso Process Iterarrival times Memoryless property ad the residual

More information

SINGLE-CHANNEL QUEUING PROBLEMS APPROACH

SINGLE-CHANNEL QUEUING PROBLEMS APPROACH SINGLE-CHANNEL QUEUING ROBLEMS AROACH Abdurrzzag TAMTAM, Doctoral Degree rogramme () Dept. of Telecommuicatios, FEEC, BUT E-mail: xtamta@stud.feec.vutbr.cz Supervised by: Dr. Karol Molár ABSTRACT The paper

More information

Queueing Systems: Lecture 3. Amedeo R. Odoni October 18, Announcements

Queueing Systems: Lecture 3. Amedeo R. Odoni October 18, Announcements Queueing Systems: Lecture 3 Amedeo R. Odoni October 18, 006 Announcements PS #3 due tomorrow by 3 PM Office hours Odoni: Wed, 10/18, :30-4:30; next week: Tue, 10/4 Quiz #1: October 5, open book, in class;

More information

Control chart for number of customers in the system of M [X] / M / 1 Queueing system

Control chart for number of customers in the system of M [X] / M / 1 Queueing system Iteratioal Joural of Iovative Research i Sciece, Egieerig ad Techology (A ISO 3297: 07 Certified Orgaiatio) Cotrol chart for umber of customers i the system of M [X] / M / Queueig system T.Poogodi, Dr.

More information

4.6 M/M/s/s Loss Queueing Model

4.6 M/M/s/s Loss Queueing Model 4.6 M/M// Lo Queueig Model Characteritic 1. Iterarrival time i expoetial with rate Arrival proce i Poio Proce with rate 2. Iterarrival time i expoetial with rate µ Number of ervice i Poio Proce with rate

More information

Birth-Death Processes. Outline. EEC 686/785 Modeling & Performance Evaluation of Computer Systems. Relationship Among Stochastic Processes.

Birth-Death Processes. Outline. EEC 686/785 Modeling & Performance Evaluation of Computer Systems. Relationship Among Stochastic Processes. EEC 686/785 Modelig & Perforace Evaluatio of Couter Systes Lecture Webig Zhao Deartet of Electrical ad Couter Egieerig Clevelad State Uiversity webig@ieee.org based o Dr. Raj jai s lecture otes Relatioshi

More information

Chapter 13: Tests of Hypothesis Section 13.1 Introduction

Chapter 13: Tests of Hypothesis Section 13.1 Introduction Chapter 13: Tests of Hypothesis Sectio 13.1 Itroductio RECAP: Chapter 1 discussed the Likelihood Ratio Method as a geeral approach to fid good test procedures. Testig for the Normal Mea Example, discussed

More information

Topic 5 [434 marks] (i) Find the range of values of n for which. (ii) Write down the value of x dx in terms of n, when it does exist.

Topic 5 [434 marks] (i) Find the range of values of n for which. (ii) Write down the value of x dx in terms of n, when it does exist. Topic 5 [44 marks] 1a (i) Fid the rage of values of for which eists 1 Write dow the value of i terms of 1, whe it does eist Fid the solutio to the differetial equatio 1b give that y = 1 whe = π (cos si

More information

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering CEE 5 Autum 005 Ucertaity Cocepts for Geotechical Egieerig Basic Termiology Set A set is a collectio of (mutually exclusive) objects or evets. The sample space is the (collectively exhaustive) collectio

More information

Introduction to probability Stochastic Process Queuing systems. TELE4642: Week2

Introduction to probability Stochastic Process Queuing systems. TELE4642: Week2 Itroductio to probability Stochastic Process Queuig systems TELE4642: Week2 Overview Refresher: Probability theory Termiology, defiitio Coditioal probability, idepedece Radom variables ad distributios

More information

Weight Moving Average = n

Weight Moving Average = n 1 Forecastig Simple Movig Average- F = 1 N (D + D 1 + D 2 + D 3 + ) Movig Weight Average- Weight Movig Average = W id i i=1 W i Sigle (Simple) Expoetial Smoothig- F t = F t 1 + α(d t 1 F t 1 ) or F t =

More information

Castiel, Supernatural, Season 6, Episode 18

Castiel, Supernatural, Season 6, Episode 18 13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio

More information

Queueing Theory II. Summary. M/M/1 Output process Networks of Queue Method of Stages. General Distributions

Queueing Theory II. Summary. M/M/1 Output process Networks of Queue Method of Stages. General Distributions Queueig Theory II Suary M/M/1 Output process Networks of Queue Method of Stages Erlag Distributio Hyperexpoetial Distributio Geeral Distributios Ebedded Markov Chais 1 M/M/1 Output Process Burke s Theore:

More information

B. Maddah ENMG 622 ENMG /20/09

B. Maddah ENMG 622 ENMG /20/09 B. Maddah ENMG 6 ENMG 5 5//9 Queueig Theory () Distributio of waitig time i M/M/ Let T q be the waitig time i queue of a ustomer. The it a be show that, ( ) t { q > } =. T t e Let T be the total time of

More information

Sigma notation. 2.1 Introduction

Sigma notation. 2.1 Introduction Sigma otatio. Itroductio We use sigma otatio to idicate the summatio process whe we have several (or ifiitely may) terms to add up. You may have see sigma otatio i earlier courses. It is used to idicate

More information

Queueing theory and Replacement model

Queueing theory and Replacement model Queueig theory ad Replacemet model. Trucks at a sigle platform weigh-bridge arrive accordig to Poisso probability distributio. The time required to weigh the truck follows a expoetial probability distributio.

More information

Queueing Theory (Part 3)

Queueing Theory (Part 3) Queueig Theory art 3 M/M/ Queueig Sytem with Variatio M/M/, M/M///K, M/M//// Queueig Theory- M/M/ Queueig Sytem We defie λ mea arrival rate mea ervice rate umber of erver > ρ λ / utilizatio ratio We require

More information

1 Introduction to reducing variance in Monte Carlo simulations

1 Introduction to reducing variance in Monte Carlo simulations Copyright c 010 by Karl Sigma 1 Itroductio to reducig variace i Mote Carlo simulatios 11 Review of cofidece itervals for estimatig a mea I statistics, we estimate a ukow mea µ = E(X) of a distributio by

More information

http://www.xelca.l/articles/ufo_ladigsbaa_houte.aspx imulatio Output aalysis 3/4/06 This lecture Output: A simulatio determies the value of some performace measures, e.g. productio per hour, average queue

More information

ECE 308 Discrete-Time Signals and Systems

ECE 308 Discrete-Time Signals and Systems ECE 38-5 ECE 38 Discrete-Time Sigals ad Systems Z. Aliyazicioglu Electrical ad Computer Egieerig Departmet Cal Poly Pomoa ECE 38-5 1 Additio, Multiplicatio, ad Scalig of Sequeces Amplitude Scalig: (A Costat

More information

Analysis of a queueing model with service threshold

Analysis of a queueing model with service threshold Rochester Istitute of Techology RIT Scholar Works Theses Thesis/Dissertatio Collectios 2-7-2 alysis of a queueig model with service threshold Kathry Graf Follow this ad additioal works at: http://scholarworks.rit.edu/theses

More information

CHAPTER 10 INFINITE SEQUENCES AND SERIES

CHAPTER 10 INFINITE SEQUENCES AND SERIES CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece

More information

Optimization Methods MIT 2.098/6.255/ Final exam

Optimization Methods MIT 2.098/6.255/ Final exam Optimizatio Methods MIT 2.098/6.255/15.093 Fial exam Date Give: December 19th, 2006 P1. [30 pts] Classify the followig statemets as true or false. All aswers must be well-justified, either through a short

More information

6.3 Testing Series With Positive Terms

6.3 Testing Series With Positive Terms 6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial

More information

Chapter 3. Strong convergence. 3.1 Definition of almost sure convergence

Chapter 3. Strong convergence. 3.1 Definition of almost sure convergence Chapter 3 Strog covergece As poited out i the Chapter 2, there are multiple ways to defie the otio of covergece of a sequece of radom variables. That chapter defied covergece i probability, covergece i

More information

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015 ECE 8527: Itroductio to Machie Learig ad Patter Recogitio Midterm # 1 Vaishali Ami Fall, 2015 tue39624@temple.edu Problem No. 1: Cosider a two-class discrete distributio problem: ω 1 :{[0,0], [2,0], [2,2],

More information

Approximations and more PMFs and PDFs

Approximations and more PMFs and PDFs Approximatios ad more PMFs ad PDFs Saad Meimeh 1 Approximatio of biomial with Poisso Cosider the biomial distributio ( b(k,,p = p k (1 p k, k λ: k Assume that is large, ad p is small, but p λ at the limit.

More information

INTRODUCTORY MATHEMATICAL ANALYSIS

INTRODUCTORY MATHEMATICAL ANALYSIS INTRODUCTORY MATHEMATICAL ANALYSIS For Busiess, Ecoomics, ad the Life ad Social Scieces Chapter 4 Itegratio 0 Pearso Educatio, Ic. Chapter 4: Itegratio Chapter Objectives To defie the differetial. To defie

More information

This is an introductory course in Analysis of Variance and Design of Experiments.

This is an introductory course in Analysis of Variance and Design of Experiments. 1 Notes for M 384E, Wedesday, Jauary 21, 2009 (Please ote: I will ot pass out hard-copy class otes i future classes. If there are writte class otes, they will be posted o the web by the ight before class

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 19 11/17/2008 LAWS OF LARGE NUMBERS II THE STRONG LAW OF LARGE NUMBERS

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 19 11/17/2008 LAWS OF LARGE NUMBERS II THE STRONG LAW OF LARGE NUMBERS MASSACHUSTTS INSTITUT OF TCHNOLOGY 6.436J/5.085J Fall 2008 Lecture 9 /7/2008 LAWS OF LARG NUMBRS II Cotets. The strog law of large umbers 2. The Cheroff boud TH STRONG LAW OF LARG NUMBRS While the weak

More information

FLUID LIMIT FOR CUMULATIVE IDLE TIME IN MULTIPHASE QUEUES. Akademijos 4, LT-08663, Vilnius, LITHUANIA 1,2 Vilnius University

FLUID LIMIT FOR CUMULATIVE IDLE TIME IN MULTIPHASE QUEUES. Akademijos 4, LT-08663, Vilnius, LITHUANIA 1,2 Vilnius University Iteratioal Joural of Pure ad Applied Mathematics Volume 95 No. 2 2014, 123-129 ISSN: 1311-8080 (prited versio); ISSN: 1314-3395 (o-lie versio) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v95i2.1

More information

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows

More information

Chapter 11 Output Analysis for a Single Model. Banks, Carson, Nelson & Nicol Discrete-Event System Simulation

Chapter 11 Output Analysis for a Single Model. Banks, Carson, Nelson & Nicol Discrete-Event System Simulation Chapter Output Aalysis for a Sigle Model Baks, Carso, Nelso & Nicol Discrete-Evet System Simulatio Error Estimatio If {,, } are ot statistically idepedet, the S / is a biased estimator of the true variace.

More information

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

Chapter 8: Estimating with Confidence

Chapter 8: Estimating with Confidence Chapter 8: Estimatig with Cofidece Sectio 8.2 The Practice of Statistics, 4 th editio For AP* STARNES, YATES, MOORE Chapter 8 Estimatig with Cofidece 8.1 Cofidece Itervals: The Basics 8.2 8.3 Estimatig

More information

If we want to add up the area of four rectangles, we could find the area of each rectangle and then write this sum symbolically as:

If we want to add up the area of four rectangles, we could find the area of each rectangle and then write this sum symbolically as: Sigma Notatio: If we wat to add up the area of four rectagles, we could fid the area of each rectagle ad the write this sum symbolically as: Sum A A A A Liewise, the sum of the areas of te triagles could

More information

MATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4

MATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4 MATH 30: Probability ad Statistics 9. Estimatio ad Testig of Parameters Estimatio ad Testig of Parameters We have bee dealig situatios i which we have full kowledge of the distributio of a radom variable.

More information

Modeling and Performance Analysis with Discrete-Event Simulation

Modeling and Performance Analysis with Discrete-Event Simulation Simulatio Modelig ad Performace Aalysis with Discrete-Evet Simulatio Chapter 5 Statistical Models i Simulatio Cotets Basic Probability Theory Cocepts Useful Statistical Models Discrete Distributios Cotiuous

More information

Kinetics of Complex Reactions

Kinetics of Complex Reactions Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet

More information

( θ. sup θ Θ f X (x θ) = L. sup Pr (Λ (X) < c) = α. x : Λ (x) = sup θ H 0. sup θ Θ f X (x θ) = ) < c. NH : θ 1 = θ 2 against AH : θ 1 θ 2

( θ. sup θ Θ f X (x θ) = L. sup Pr (Λ (X) < c) = α. x : Λ (x) = sup θ H 0. sup θ Θ f X (x θ) = ) < c. NH : θ 1 = θ 2 against AH : θ 1 θ 2 82 CHAPTER 4. MAXIMUM IKEIHOOD ESTIMATION Defiitio: et X be a radom sample with joit p.m/d.f. f X x θ. The geeralised likelihood ratio test g.l.r.t. of the NH : θ H 0 agaist the alterative AH : θ H 1,

More information

Integer Programming (IP)

Integer Programming (IP) Iteger Programmig (IP) The geeral liear mathematical programmig problem where Mied IP Problem - MIP ma c T + h Z T y A + G y + y b R p + vector of positive iteger variables y vector of positive real variables

More information

Ma 530 Infinite Series I

Ma 530 Infinite Series I Ma 50 Ifiite Series I Please ote that i additio to the material below this lecture icorporated material from the Visual Calculus web site. The material o sequeces is at Visual Sequeces. (To use this li

More information

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

Increasing timing capacity using packet coloring

Increasing timing capacity using packet coloring 003 Coferece o Iformatio Scieces ad Systems, The Johs Hopkis Uiversity, March 4, 003 Icreasig timig capacity usig packet colorig Xi Liu ad R Srikat[] Coordiated Sciece Laboratory Uiversity of Illiois e-mail:

More information

Summary: CORRELATION & LINEAR REGRESSION. GC. Students are advised to refer to lecture notes for the GC operations to obtain scatter diagram.

Summary: CORRELATION & LINEAR REGRESSION. GC. Students are advised to refer to lecture notes for the GC operations to obtain scatter diagram. Key Cocepts: 1) Sketchig of scatter diagram The scatter diagram of bivariate (i.e. cotaiig two variables) data ca be easily obtaied usig GC. Studets are advised to refer to lecture otes for the GC operatios

More information

Chapter 4. Fourier Series

Chapter 4. Fourier Series Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,

More information

lim za n n = z lim a n n.

lim za n n = z lim a n n. Lecture 6 Sequeces ad Series Defiitio 1 By a sequece i a set A, we mea a mappig f : N A. It is customary to deote a sequece f by {s } where, s := f(). A sequece {z } of (complex) umbers is said to be coverget

More information

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense, 3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [

More information

ENGI Series Page 6-01

ENGI Series Page 6-01 ENGI 3425 6 Series Page 6-01 6. Series Cotets: 6.01 Sequeces; geeral term, limits, covergece 6.02 Series; summatio otatio, covergece, divergece test 6.03 Stadard Series; telescopig series, geometric series,

More information

Most text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t

Most text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said

More information

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.

More information

Ma 530 Introduction to Power Series

Ma 530 Introduction to Power Series Ma 530 Itroductio to Power Series Please ote that there is material o power series at Visual Calculus. Some of this material was used as part of the presetatio of the topics that follow. What is a Power

More information

Chapter 10: Power Series

Chapter 10: Power Series Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because

More information

Mark Lundstrom Spring SOLUTIONS: ECE 305 Homework: Week 5. Mark Lundstrom Purdue University

Mark Lundstrom Spring SOLUTIONS: ECE 305 Homework: Week 5. Mark Lundstrom Purdue University Mark udstrom Sprig 2015 SOUTIONS: ECE 305 Homework: Week 5 Mark udstrom Purdue Uiversity The followig problems cocer the Miority Carrier Diffusio Equatio (MCDE) for electros: Δ t = D Δ + G For all the

More information

The statistical pattern of the arrival can be indicated through the probability distribution of the number of the arrivals in an interval.

The statistical pattern of the arrival can be indicated through the probability distribution of the number of the arrivals in an interval. Itroductio Queuig are the most freuetly ecoutered roblems i everyday life. For examle, ueue at a cafeteria, library, bak, etc. Commo to all of these cases are the arrivals of objects reuirig service ad

More information

NOTES ON DISTRIBUTIONS

NOTES ON DISTRIBUTIONS NOTES ON DISTRIBUTIONS MICHAEL N KATEHAKIS Radom Variables Radom variables represet outcomes from radom pheomea They are specified by two objects The rage R of possible values ad the frequecy fx with which

More information

Convergence of a queueing system in heavy traffic with general patience-time distributions

Convergence of a queueing system in heavy traffic with general patience-time distributions Available olie at www.sciecedirect.com Stochastic Processes ad their Applicatios 121 (211) 257 2552 www.elsevier.com/locate/spa Covergece of a queueig system i heavy traffic with geeral patiece-time distributios

More information

There is no straightforward approach for choosing the warmup period l.

There is no straightforward approach for choosing the warmup period l. B. Maddah INDE 504 Discrete-Evet Simulatio Output Aalysis () Statistical Aalysis for Steady-State Parameters I a otermiatig simulatio, the iterest is i estimatig the log ru steady state measures of performace.

More information

The Poisson Distribution

The Poisson Distribution MATH 382 The Poisso Distributio Dr. Neal, WKU Oe of the importat distributios i probabilistic modelig is the Poisso Process X t that couts the umber of occurreces over a period of t uits of time. This

More information

PC5215 Numerical Recipes with Applications - Review Problems

PC5215 Numerical Recipes with Applications - Review Problems PC55 Numerical Recipes with Applicatios - Review Problems Give the IEEE 754 sigle precisio bit patter (biary or he format) of the followig umbers: 0 0 05 00 0 00 Note that it has 8 bits for the epoet,

More information

MCT242: Electronic Instrumentation Lecture 2: Instrumentation Definitions

MCT242: Electronic Instrumentation Lecture 2: Instrumentation Definitions Faculty of Egieerig MCT242: Electroic Istrumetatio Lecture 2: Istrumetatio Defiitios Overview Measuremet Error Accuracy Precisio ad Mea Resolutio Mea Variace ad Stadard deviatio Fiesse Sesitivity Rage

More information

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018 CSE 353 Discrete Computatioal Structures Sprig 08 Sequeces, Mathematical Iductio, ad Recursio (Chapter 5, Epp) Note: some course slides adopted from publisher-provided material Overview May mathematical

More information

Course 4: Preparation for Calculus Unit 1: Families of Functions

Course 4: Preparation for Calculus Unit 1: Families of Functions Course 4: Preparatio for Calculus Uit 1: Families of Fuctios Review ad exted properties of basic fuctio families ad their uses i mathematical modelig Develop strategies for fidig rules of fuctios whose

More information

Sequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

Sequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet

More information

M/M/1 queueing system with server start up, unreliable server and balking

M/M/1 queueing system with server start up, unreliable server and balking Research Joural of Mathematical ad Statistical Scieces ISS 232-647 Vol. 6(2), -9, February (28) Res. J. Mathematical ad Statistical Sci. Aalysis of a -policy M/M/ queueig system with server start up, ureliable

More information

c. Explain the basic Newsvendor model. Why is it useful for SC models? e. What additional research do you believe will be helpful in this area?

c. Explain the basic Newsvendor model. Why is it useful for SC models? e. What additional research do you believe will be helpful in this area? 1. Research Methodology a. What is meat by the supply chai (SC) coordiatio problem ad does it apply to all types of SC s? Does the Bullwhip effect relate to all types of SC s? Also does it relate to SC

More information

ECO 312 Fall 2013 Chris Sims LIKELIHOOD, POSTERIORS, DIAGNOSING NON-NORMALITY

ECO 312 Fall 2013 Chris Sims LIKELIHOOD, POSTERIORS, DIAGNOSING NON-NORMALITY ECO 312 Fall 2013 Chris Sims LIKELIHOOD, POSTERIORS, DIAGNOSING NON-NORMALITY (1) A distributio that allows asymmetry differet probabilities for egative ad positive outliers is the asymmetric double expoetial,

More information

A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece,, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet as

More information

Math 312 Lecture Notes One Dimensional Maps

Math 312 Lecture Notes One Dimensional Maps Math 312 Lecture Notes Oe Dimesioal Maps Warre Weckesser Departmet of Mathematics Colgate Uiversity 21-23 February 25 A Example We begi with the simplest model of populatio growth. Suppose, for example,

More information

Expectation and Variance of a random variable

Expectation and Variance of a random variable Chapter 11 Expectatio ad Variace of a radom variable The aim of this lecture is to defie ad itroduce mathematical Expectatio ad variace of a fuctio of discrete & cotiuous radom variables ad the distributio

More information

0, otherwise. EX = E(X 1 + X n ) = EX j = np and. Var(X j ) = np(1 p). Var(X) = Var(X X n ) =

0, otherwise. EX = E(X 1 + X n ) = EX j = np and. Var(X j ) = np(1 p). Var(X) = Var(X X n ) = PROBABILITY MODELS 35 10. Discrete probability distributios I this sectio, we discuss several well-ow discrete probability distributios ad study some of their properties. Some of these distributios, lie

More information

Machine Learning Brett Bernstein

Machine Learning Brett Bernstein Machie Learig Brett Berstei Week 2 Lecture: Cocept Check Exercises Starred problems are optioal. Excess Risk Decompositio 1. Let X = Y = {1, 2,..., 10}, A = {1,..., 10, 11} ad suppose the data distributio

More information

: Transforms and Partial Differential Equations

: Transforms and Partial Differential Equations Trasforms ad Partial Differetial Equatios 018 SUBJECT NAME : Trasforms ad Partial Differetial Equatios SUBJECT CODE : MA 6351 MATERIAL NAME : Part A questios REGULATION : R013 WEBSITE : wwwharigaeshcom

More information

4. Basic probability theory

4. Basic probability theory Cotets Basic cocepts Discrete radom variables Discrete distributios (br distributios) Cotiuous radom variables Cotiuous distributios (time distributios) Other radom variables Lect04.ppt S-38.45 - Itroductio

More information

ADVANCED SOFTWARE ENGINEERING

ADVANCED SOFTWARE ENGINEERING ADVANCED SOFTWARE ENGINEERING COMP 3705 Exercise Usage-based Testig ad Reliability Versio 1.0-040406 Departmet of Computer Ssciece Sada Narayaappa, Aeliese Adrews Versio 1.1-050405 Departmet of Commuicatio

More information

Statistical Inference (Chapter 10) Statistical inference = learn about a population based on the information provided by a sample.

Statistical Inference (Chapter 10) Statistical inference = learn about a population based on the information provided by a sample. Statistical Iferece (Chapter 10) Statistical iferece = lear about a populatio based o the iformatio provided by a sample. Populatio: The set of all values of a radom variable X of iterest. Characterized

More information

Session 5. (1) Principal component analysis and Karhunen-Loève transformation

Session 5. (1) Principal component analysis and Karhunen-Loève transformation 200 Autum semester Patter Iformatio Processig Topic 2 Image compressio by orthogoal trasformatio Sessio 5 () Pricipal compoet aalysis ad Karhue-Loève trasformatio Topic 2 of this course explais the image

More information

Properties and Hypothesis Testing

Properties and Hypothesis Testing Chapter 3 Properties ad Hypothesis Testig 3.1 Types of data The regressio techiques developed i previous chapters ca be applied to three differet kids of data. 1. Cross-sectioal data. 2. Time series data.

More information

A Feedback Retrial Queuing System with Starting Failures and Single Vacation

A Feedback Retrial Queuing System with Starting Failures and Single Vacation Tamkag Joural of Sciece ad Egieerig, Vol., No. 3, pp. 8392 (27) 83 A Feedback Retrial Queuig System with Startig Failures ad Sigle Vacatio G. S. Mokaddis, S. A. Metwally ad B. M. Zaki Ai Shams Uiversity,

More information

AMS 216 Stochastic Differential Equations Lecture 02 Copyright by Hongyun Wang, UCSC ( ( )) 2 = E X 2 ( ( )) 2

AMS 216 Stochastic Differential Equations Lecture 02 Copyright by Hongyun Wang, UCSC ( ( )) 2 = E X 2 ( ( )) 2 AMS 216 Stochastic Differetial Equatios Lecture 02 Copyright by Hogyu Wag, UCSC Review of probability theory (Cotiued) Variace: var X We obtai: = E X E( X ) 2 = E( X 2 ) 2E ( X )E X var( X ) = E X 2 Stadard

More information

ECONOMIC OPERATION OF POWER SYSTEMS

ECONOMIC OPERATION OF POWER SYSTEMS ECOOMC OEATO OF OWE SYSTEMS TOUCTO Oe of the earliest applicatios of o-lie cetralized cotrol was to provide a cetral facility, to operate ecoomically, several geeratig plats supplyig the loads of the system.

More information

EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY

EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA, 016 MODULE : Statistical Iferece Time allowed: Three hours Cadidates should aswer FIVE questios. All questios carry equal marks. The umber

More information

SECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES

SECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES SECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES Read Sectio 1.5 (pages 5 9) Overview I Sectio 1.5 we lear to work with summatio otatio ad formulas. We will also itroduce a brief overview of sequeces,

More information

PBW 654 Applied Statistics - I Urban Operations Research

PBW 654 Applied Statistics - I Urban Operations Research PBW 654 Applied Statistics - I Urban Operations Research Lecture 2.I Queuing Systems An Introduction Operations Research Models Deterministic Models Linear Programming Integer Programming Network Optimization

More information

If a subset E of R contains no open interval, is it of zero measure? For instance, is the set of irrationals in [0, 1] is of measure zero?

If a subset E of R contains no open interval, is it of zero measure? For instance, is the set of irrationals in [0, 1] is of measure zero? 2 Lebesgue Measure I Chapter 1 we defied the cocept of a set of measure zero, ad we have observed that every coutable set is of measure zero. Here are some atural questios: If a subset E of R cotais a

More information

Discrete Mathematics and Probability Theory Spring 2016 Rao and Walrand Note 19

Discrete Mathematics and Probability Theory Spring 2016 Rao and Walrand Note 19 CS 70 Discrete Mathematics ad Probability Theory Sprig 2016 Rao ad Walrad Note 19 Some Importat Distributios Recall our basic probabilistic experimet of tossig a biased coi times. This is a very simple

More information

4.1 Sigma Notation and Riemann Sums

4.1 Sigma Notation and Riemann Sums 0 the itegral. Sigma Notatio ad Riema Sums Oe strategy for calculatig the area of a regio is to cut the regio ito simple shapes, calculate the area of each simple shape, ad the add these smaller areas

More information

NANYANG TECHNOLOGICAL UNIVERSITY SYLLABUS FOR ENTRANCE EXAMINATION FOR INTERNATIONAL STUDENTS AO-LEVEL MATHEMATICS

NANYANG TECHNOLOGICAL UNIVERSITY SYLLABUS FOR ENTRANCE EXAMINATION FOR INTERNATIONAL STUDENTS AO-LEVEL MATHEMATICS NANYANG TECHNOLOGICAL UNIVERSITY SYLLABUS FOR ENTRANCE EXAMINATION FOR INTERNATIONAL STUDENTS AO-LEVEL MATHEMATICS STRUCTURE OF EXAMINATION PAPER. There will be oe 2-hour paper cosistig of 4 questios.

More information

2C09 Design for seismic and climate changes

2C09 Design for seismic and climate changes 2C09 Desig for seismic ad climate chages Lecture 02: Dyamic respose of sigle-degree-of-freedom systems I Daiel Grecea, Politehica Uiversity of Timisoara 10/03/2014 Europea Erasmus Mudus Master Course Sustaiable

More information

STOCHASTIC NETWORKS EXAMPLE SHEET 2 SOLUTIONS

STOCHASTIC NETWORKS EXAMPLE SHEET 2 SOLUTIONS STOCHASTIC NETWORKS EXAMPLE SHEET 2 SOLUTIONS ELENA YUDOVINA Exercise. Show that the traffic equatios for a ope migratio process have a uique solutio, ad that this solutio is positive. [Hit: From the irreducibility

More information