Construction of Control Chart for Random Queue Length for (M / M / c): ( / FCFS) Queueing Model Using Skewness

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

B. Maddah ENMG 622 ENMG /20/09

Optimal design of N-Policy batch arrival queueing system with server's single vacation, setup time, second optional service and break down

First come, first served (FCFS) Batch

ON STEREOGRAPHIC CIRCULAR WEIBULL DISTRIBUTION

Modeling and Performance Analysis with Discrete-Event Simulation

THE MEASUREMENT OF THE SPEED OF THE LIGHT

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

Optimal Management of the Spare Parts Stock at Their Regular Distribution

Statistical Fundamentals and Control Charts

Chapter 8 Hypothesis Testing

Nonparametric Goodness-of-Fit Tests for Discrete, Grouped or Censored Data 1

B. Maddah ENMG 622 ENMG /27/07

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,

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

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

Summation Method for Some Special Series Exactly

Simulation of Discrete Event Systems

SINGLE-CHANNEL QUEUING PROBLEMS APPROACH

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

Certain inclusion properties of subclass of starlike and convex functions of positive order involving Hohlov operator

R. van Zyl 1, A.J. van der Merwe 2. Quintiles International, University of the Free State

Basic Probability/Statistical Theory I

Queueing theory and Replacement model

What is a Hypothesis? Hypothesis is a statement about a population parameter developed for the purpose of testing.

International Journal of Advanced Research in Computer Science and Software Engineering

Confidence Interval for Standard Deviation of Normal Distribution with Known Coefficients of Variation

Pass-Fail Testing: Statistical Requirements and Interpretations

CONSIDERATION OF SOME PERFORMANCE OF CONTAINERS FLOW AT YARD

The standard deviation of the mean

(# x) 2 n. (" x) 2 = 30 2 = 900. = sum. " x 2 = =174. " x. Chapter 12. Quick math overview. #(x " x ) 2 = # x 2 "

CONTROL CHARTS FOR THE LOGNORMAL DISTRIBUTION

Chapter 22. Comparing Two Proportions. Copyright 2010 Pearson Education, Inc.

Topic 1 2: Sequences and Series. A sequence is an ordered list of numbers, e.g. 1, 2, 4, 8, 16, or

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

Sample Size Requirements For Stratified Random Sampling Of Pond Water With Cost Considerations Using A Bayesian Methodology

IE 230 Probability & Statistics in Engineering I. Closed book and notes. No calculators. 120 minutes.

Probability & Statistics Chapter 8

An Introduction to Randomized Algorithms

BIOS 4110: Introduction to Biostatistics. Breheny. Lab #9

Using the IML Procedure to Examine the Efficacy of a New Control Charting Technique

Approximate Confidence Interval for the Reciprocal of a Normal Mean with a Known Coefficient of Variation

Record Values from T-X Family of. Pareto-Exponential Distribution with. Properties and Simulations

Chapter 22. Comparing Two Proportions. Copyright 2010, 2007, 2004 Pearson Education, Inc.

(Dependent or paired samples) Step (1): State the null and alternate hypotheses: Case1: One-tailed test (Right)

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

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

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

Estimation of Gumbel Parameters under Ranked Set Sampling


Estimation of Population Mean Using Co-Efficient of Variation and Median of an Auxiliary Variable

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

Increasing timing capacity using packet coloring

Economics 250 Assignment 1 Suggested Answers. 1. We have the following data set on the lengths (in minutes) of a sample of long-distance phone calls

MATH/STAT 352: Lecture 15

Chapter 6 Sampling Distributions

NCSS Statistical Software. Tolerance Intervals

STAT 515 fa 2016 Lec Sampling distribution of the mean, part 2 (central limit theorem)

Direction: This test is worth 150 points. You are required to complete this test within 55 minutes.

Mixed Acceptance Sampling Plans for Multiple Products Indexed by Cost of Inspection

Modeling and Performance Analysis with Discrete-Event Simulation

Statistics 511 Additional Materials

ENGI 4421 Confidence Intervals (Two Samples) Page 12-01

EconoQuantum ISSN: Universidad de Guadalajara México

NOTES ON DISTRIBUTIONS

Review Questions, Chapters 8, 9. f(y) = 0, elsewhere. F (y) = f Y(1) = n ( e y/θ) n 1 1 θ e y/θ = n θ e yn

MATHEMATICS Paper 2 22 nd September 20. Answer Papers List of Formulae (MF15)

The Poisson Distribution

E. Teimoury *, I.G. Khondabi & M. Fathi

COMP26120: Introducing Complexity Analysis (2018/19) Lucas Cordeiro

A statistical method to determine sample size to estimate characteristic value of soil parameters

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering

THE POTENTIALS METHOD FOR A CLOSED QUEUEING SYSTEM WITH HYSTERETIC STRATEGY OF THE SERVICE TIME CHANGE

Optimization design in Wind Turbine Blade Based on Specific Wind Characteristics

Bayesian Control Charts for the Two-parameter Exponential Distribution

Average Number of Real Zeros of Random Fractional Polynomial-II

Test of Statistics - Prof. M. Romanazzi

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

STA Learning Objectives. Population Proportions. Module 10 Comparing Two Proportions. Upon completing this module, you should be able to:

Math Third Midterm Exam November 17, 2010


Central Limit Theorem the Meaning and the Usage

STAT 350 Handout 19 Sampling Distribution, Central Limit Theorem (6.6)

PH 425 Quantum Measurement and Spin Winter SPINS Lab 1

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.

4. Basic probability theory

Lecture 5. Materials Covered: Chapter 6 Suggested Exercises: 6.7, 6.9, 6.17, 6.20, 6.21, 6.41, 6.49, 6.52, 6.53, 6.62, 6.63.

Soo King Lim Figure 1: Figure 2: Figure 3: Figure 4: Figure 5: Figure 6: Figure 7:

Statistical Intervals for a Single Sample

Estimation for Complete Data

Parameter, Statistic and Random Samples

If, for instance, we were required to test whether the population mean μ could be equal to a certain value μ

Closed book and notes. No calculators. 60 minutes, but essentially unlimited time.

MATH CALCULUS II Objectives and Notes for Test 4

Modified Ratio Estimators Using Known Median and Co-Efficent of Kurtosis

Reliability and Queueing

Lecture 24 Floods and flood frequency

SOLUTIONS. 1. a) Let X and Y be random variables with. σ Y = 9. σ X 30 σ XY σ X 30 ρ σ X σ Y σ X + 24 σ XY

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

Confidence interval for the two-parameter exponentiated Gumbel distribution based on record values

Transcription:

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember ISSN 5-5 Costrutio of Cotrol Chart for Radom Queue Legth for (M / M / ): ( / FCFS) Queueig Model Usig Skewess Dr.(Mrs.) A.R. Sudamai Ramaswamy, Mrs. B.Veila *Departmet of Mathematis, Aviashiligam Istitute for HomeSiee ad Higher Eduatio for Wome, Coimbatore-4. **Departmet of Mathematis, Sri Ramakrisha Egieerig College, Coimbatore-. Dr. (Mrs.) A.R. Sudamai Ramaswamy, Ph.D, arsudamairamaswamy@hotmail.om Correspodee Author: Mrs. B.Veila, M.Phil, veila.bs@gmail.om ABSTRACT I this paper, first-ome, first served, multiple hael Poisso/epoetial queueig system M / M / with ifiite waitig apaity is osidered. For this model, we itrodue a ew proedure for ostrutio of i) Shewhart s otrol hart C, ii) Cotrol hart C for radom queue legth N usig the method o skewess suggested by Shore() ad the otrol harts C ad C are ompared. Ide Terms: Radom queue legth, Average ru legth, False alarm rate, Average queue legth.. INTRODUCTION Waitig i lie for servie is oe of the most upleasat eperiees of life o this world. Barrer (957) says, I ertai queueig proesses a potetial ustomer is osidered lost if the system is busy at the time servie is demaded. If ot served durig this time, the ustomer leaves the system ad is osidered lost. I queueig system the ustomer satisfatio a be ireased by ostrutig otrol harts for N ad providig otrol limits for N so as to make effetive utilizatio of time. If otrol limits are displayed so that ustomer a have prior idea about otrol limits. With this fator, a attempt is made to fid otrol limits for radom queue legth N for (M/M/) : ( /FCFS) queueig model. The pioeerig work i this diretio was made by Haim Shore i. Two otrol harts C ad C are ostruted for radom queue legth N for (M/M/) : ( /FCFS)

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember ISSN 5-5 queueig model by Khaparde ad Dhabe () ad also C ad C 4 are ostruted for radom queue legth N for (M/M/) : ( /FCFS) queueig model by Khaparde ad Dhabe () where Cotrol hart C : is the Shewhart otrol hart, Cotrol hart C : is the otrol hart usig method of weighted variae, Cotrol hart C : is the otrol hart for radom queue legth N based o skewess ad Cotrol hart C 4 : is the otrol hart usig Nelso s power trasformatio. I this paper ostrutio of otrol hart C for radom queue legth N for (M/M/) : ( /FCFS) queueig model is doe usig method based o skewess. Performae of otrol hart C is ompared with otrol hart C.. (M / M / ) : ( / FCFS) QUEUEING MODEL I this model arrival follows Poisso distributio with parameter ad servie time follows idepedetly ad idetially distributed epoetial distributio with parameter. There are idetial servers preset i the system. Figure I represet the state-trasitio diagram for this model. + () FIGURE I : STATE-TRANSITION RATE DIAGRAM FOR M / M / MODEL Let radom variable N deote the umber of ustomers i the system both waitig ad i servie ad P deote the steady state distributio of radom variable N give by P = P,! P,! () where P =,!! is the mea arrival rate ad is the mea servie rate. Let radom variable W S deote the waitig time spet i the system by the ustomer. This iludes both the waitig time ad servie time. The p.d.f. of radom variable W S is give by

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember ISSN 5-5 ( ) WS f(w S ) = ( ) e, W S > () Thus W S follows a epoetial distributio with mea E(W S ) = = W S f(ws ) d WS ( )W S ( ) e d WS W S E(W S ) = ( ) () E(W S ) = ( ) (4) Usig (.) ad (.4), the variae is give by V(W S ) = ( ) (5) The distributio futio of W S is F() = P(W S ) = f(w ) d S W S = ( )W S ( ) e d W S = e ( ), >. F() = e ( ), >. The probability distributio of N is a geometri distributio with parameter ( ). The steady state distributio of the radom variable N depeds o two parameters ad oly through their ratio. Usig (), the average umber of ustomers ad variae of queue legth i the system is give by E(N) = = P + P = P +!! P

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember 4 ISSN 5-5 = P )! ( )! ( Here P =!! E(N) = =!!!! (6) E(N ) = = ) P )P (( = )P ( + )P ( + P + P = ) (! P + ) (! P +! P +! P E(N ) = =!!!! (7) Usig (6) ad (7), the variae of queue legth i the system is give by V(N) =!!!!

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember 5 ISSN 5-5!! )! ( )! ( (8) Usig (8), the stadard deviatio is alulated as = /!! )! ( )! (!!!! (9) E(N ) = = )P ) ( ) )( (( = P! ) )( ( +! ) )( ( +! ) ( +! ) ( +!! =!!!! ()

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember 6 ISSN 5-5 Usig (), the third etral momet is give by =!!!! -!!!!! ) ( )! ( +!! =!!!! )! ( )! (!!!!!! () Let S k (N) deote the traditioal measure of skewess for radom queue legth amely, S k (N) =. where is the third etral momet.

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember 7 ISSN 5-5 Usig (9) ad (), the skewess for radom queue legth is give by S k (N) = / )! ( )! (!!!! )! ( )! (!! )! ( )! (!!!!!! (). SHEWHART S CONTROL CHART C Whe it is possible to speify stadard values for the proess mea ad stadard deviatio, we may use these stadards to establish the otrol hart. The the parameters of the hart are UCL = + A CL = LCL = A where the quatity = A, say, is a ostat that depeds o. Usig epressios (6) ad (9), the otrol limits for (M / M / ) : ( / FCFS) model are alulated as follows.

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember 8 ISSN 5-5 UCL =!! )! ( )! (!!!! A )! ( )! ( / CL =!! )! ( )! ( LCL =!! )! ( )! (!!!! A )! ( )! ( / TABLE I MEAN, STANDARD DEVIATION, UPPER CONTROL LIMIT FOR THE RANDOM VARIABLE N FOR DIFFERENT VALUES OF USING SHEWHART S CONTROL CHART C WITH A = AND C =

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember 9 ISSN 5-5 S.No. CL UCL u ARL..5.48.9.448.97.54...44..64.55.8..5.579.6.7979.86.6 4. 4..6.4..96.5 5. 5.5.675.5.75.444.5 6. 6..79.6.9.598.96 7. 7.5.7569.7.4569.5947.68 8. 8.4.7857.8.4857.778.9 9. 9.45.85.9.75.855.. CONSTRUCTION OF CONTROL CHART C FOR RANDOM VARIABLE N Before ostrutig the hart, it is essetial to otie that the probability futio of a geometri distributio is a mootoously dereasig futio of ad sie otrol limits for queueig system are to be obtaied, there is o meaig i defiig lower otrol limit (LCL) for this distributio. Also upper otrol limit (UCL) obviously makes sese sie it will eable moitorig improbable high value of N.. CONTROL CHART C I the queueig model (M / M / ) : ( / FCFS), the radom variable N deote the umber of ustomers i the system both waitig ad i servie. The steady state distributio of radom variable N is give by epressio (). I order to obtai otrol limits usig Shore s method (), the first three momets of the radom variable N give by epressios (6), (8) ad () respetively, are eeded. The epressio () shows that skewess of the distributio of N depeds o,, ad. Takig ito aout the mea, stadard deviatio ad skewess of the uderlyig geometri distributio of N the epressio for the ew otrol limits for the umber of ustomers i (M / M / ) : ( / FCFS) system usig Shore s method () are as follows : Cotrol Limits for Cotrol Chart C whe S k <.5

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember ISSN 5-5 UCL = E(N) + V(N) +.4 V(N) S k (N). CL = E(N) LCL = E(N) V(N) +.4 V(N) S k (N) +. Cotrol Limits for Cotrol Chart C whe S k >.5 UCL = E(N) +.64 V(N) +.946 V(N) S k (N). CL = E(N) LCL = E(N).64 V(N) + (.4) (.946) V(N) S k (N) +. Cotrol Limits for Cotrol Chart C for (M / M / ) : ( / FCFS) Model whe S k <.5 UCL = )! ( )! (!! +!!!! )! ( )! (!!

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember ISSN 5-5 +.4!!!! )! ( )! (!!!!!!!! )! ( )! (!! + )! ( )! ( CL = )! ( )! (

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember ISSN 5-5!! LCL = )! ( )! (!!!!!! )! ( )! (!! +.4!!!! )! ( )! (!!

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember ISSN 5-5!!!!!! )! ( )! (!! + )! ( )! ( Cotrol Limits for Cotrol Chart C for (M / M / ) : ( / FCFS) Model whe S k >.5 UCL = )! ( )! (!! +.64!!!!

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember 4 ISSN 5-5 )! ( )! (!! + (.946)!!!! )! ( )! (!!!!!!!! )! ( )! (!!

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember 5 ISSN 5-5 + )! ( )! ( CL = )! ( )! (!! LCL = )! ( )! (!!.64!!!! )! ( )! (!! + (.4) (.946)!!!!

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember 6 ISSN 5-5 )! ( )! (!!!!!!!! )! ( )! (!! + )! ( )! ( Table. shows that the miimum value of skewess is early. for high values of. TABLE II SKEWNESS, CONTROL LIMITS FOR THE RANDOM VARIABLE N FOR DIFFERENT VALUES OF USING CONTROL CHART C WITH L = ad C = S.No. Skewess CL UCL u ARL..5..9.45.9 9.8...657..87.99.4

Upper Cotrol Limit (UCL) Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember 7 ISSN 5-5..5.6.6..6 8.4 4. 4..66.4.5.56 7.86 5. 5.5.6.5.47.9. 6. 6..654.6.9.8 5.49 7 7.5.5597.7.58.85.5 8 8.4.65.8.54.456. 9 9.45.88.9.58.6896.45 4.5 Cotrol hart C Cotrol hart C.5.5.5....4.5 Traffi itesity (ρ) FIGURE II :COMPARISON OF UPPER CONTROL LIMIT (UCL) OF CONTROL CHARTS C AND C WITH L = AND C = FOR DIFFERENT VALUES OF It is otied that the upper otrol limit of otrol hart C outperforms the otrol hart C for all values of. Upper otrol limit values i otrol hart C are very very large tha that of hart C. Therefore for skewed populatio the use of otrol hart C is advisable. This a be observed from Figure II

Iteratioal Joural of Sietifi ad Researh Publiatios, Volume, Issue, Deember 8 ISSN 5-5 FIGURE III: COMPARISON OF AVERAGE RUN LENGTH OF CONTROL CHARTS C AND C WITH L = AND C = FOR DIFFERENT VALUES OF CONCLUSON The urves of ARL values are displayed i Figure III for harts C ad C. The horizotal ais of spas.5 to.45, Figure III demostrates that the ARL produed by hart C is superior to the hart C. REFERENCES. Arold, A.O. (978). Probability, Statistis ad Queueig Theory with Computer Siee Appliatios, Seod Editio, Aademi Press Limited, New York.. Barrer, D.Y. (957). Queueig with Impatiet Customers ad Idifferet Clerks, Operatios Researh, 5, 644-649.. Khaparde, M.V. ad Dhabe, S.D. (). Cotrol Charts for Radom Queue Legth N for (M / M / ) : ( / FCFS) Model, Iteratioal Joural of Agriultural ad Statistial Siees,, 9-4. 4. Khaparde, M.V. ad Dhabe, S.D. (). Cotrol Charts for Radom Queue Legth N for (M / M / ) : ( / FCFS) Queueig Model Usig Skewess ad Power Trasformatio, Bulleti of Pure ad Applied Siees, E, 7-8. 5. Motgomery, D.C. (5). Itrodutio to Statistial Quality Cotrol, Fifth Editio, Joh Wiley ad Sos I. 6. Shore, H. (). Geeral Cotrol Charts for Attributes, IIE Trasatios,, 49-6. 7. Sudamai Ramaswamy, A.R. (). Statistial Aalysis of M / M / Model, Researh Highlights,, 7-.