Customer-Oriented Finite Perturbation Analysis for Queueing Networks

We consider queueing networks for which the performance measureJ ( θ ) depends on a parameter θ, which can be a service time parameter or a buffer size, and we are interested in sensitivity analysis of J ( θ ) with respect to θ. We introduce a new method, called customer-oriented finite perturbation analysis (CFPA), which predicts J ( θ + Δ ) for an arbitrary, finite perturbation Δ from a simulation experiment at θ. CFPA can estimate the entire performance function (by using a finite number of chosen points and fitting a least-squares approximating polynomial to the observation) within one simulation experiment. We obtain CFPA by reformulating finite perturbation analysis (FPA) for customers. The main difference between FPA and CFPA is that the former calculates the sensitivities of timing epochs of events, such as external arrivals or service time completions, while the latter yields sensitivities of departure epochs of customers. We give sufficient conditions for unbiasedness of CFPA. Numerical examples show the efficiency of the method. In particular, we address sensitivity analysis with respect to buffer sizes and thereby give a solution to the problem for which perturbation analysis was originally built.

[1]  Onno Boxma,et al.  Sojourn times in queueing networks , 1989 .

[2]  Yu-Chi Ho,et al.  On derivative estimation of single-server queues via structural infinitesimal perturbation analysis , 1995, Discret. Event Dyn. Syst..

[3]  Jianqiang Hu Convexity of sample path performance and strong consistency of infinitesimal perturbation analysis estimates , 1992 .

[4]  Leyuan Shi Variance property of discontinuous perturbation analysis , 1996, Winter Simulation Conference.

[5]  Laurent Massoulié,et al.  Maximal coupling Rare Perturbation Analysis with a random horizon , 1995, Discret. Event Dyn. Syst..

[6]  Pierre Brémaud,et al.  On the pathwise computation of derivatives with respect to the rate of a point process: The phantom RPA method , 1992, Queueing Syst. Theory Appl..

[7]  Christos G. Cassandras,et al.  A new approach to the analysis of discrete event dynamic systems , 1983, Autom..

[8]  Paul Bratley,et al.  A guide to simulation , 1983 .

[9]  Hugo Lucca,et al.  A New Approach to Compute Performance Sensitivities of Stochastic Discrete Event Dynamic Systems (DEDS) , 1993 .

[10]  L. Dai,et al.  Optimizing discrete event dynamic systems via the gradient surface method , 1991, [1991] Proceedings of the 30th IEEE Conference on Decision and Control.

[11]  Gang Bao,et al.  First and Second Derivative Estimators for Closed Jackson-Like Queueing Networks Using Perturbation Analysis Techniques , 1997, Discret. Event Dyn. Syst..

[12]  Bernd Heidergott,et al.  Sensitivity analysis of a manufacturing workstation using perturbation analysis techniques , 1995 .

[13]  Christos G. Cassandras,et al.  Perturbation analytic methodologies for design and optimization of communication networks , 1988, IEEE J. Sel. Areas Commun..

[14]  Nikolai Krivulin Unbiased estimates for gradients of stochastic network performance measures , 1993 .

[15]  Paul Glasserman,et al.  Gradient Estimation Via Perturbation Analysis , 1990 .

[16]  Dinah W. Cheng,et al.  Tandem queues with general blocking: A unified model and comparison results , 1993, Discret. Event Dyn. Syst..

[17]  Pirooz Vakili,et al.  Uniformization based sensitivity estimation for a class of discrete-event systems , 1994, Discret. Event Dyn. Syst..

[18]  Y. Ho,et al.  Smoothed (conditional) perturbation analysis of discrete event dynamical systems , 1987 .

[19]  R. H. Myers,et al.  Response Surface Methodology: Process and Product Optimization Using Designed Experiments , 1995 .

[20]  Y. Ho,et al.  Extensions of infinitesimal perturbation analysis , 1988 .

[21]  Y. Ho,et al.  Structural infinitesimal perturbation analysis (SIPA) for derivative estimation of discrete-event dynamic systems , 1995, IEEE Trans. Autom. Control..

[22]  Y. C. Ho,et al.  A New Approach to Determine Parameter Sensitivities of Transfer Lines , 1983 .

[23]  V. Nollau Kushner, H. J./Clark, D. S., Stochastic Approximation Methods for Constrained and Unconstrained Systems. (Applied Mathematical Sciences 26). Berlin‐Heidelberg‐New York, Springer‐Verlag 1978. X, 261 S., 4 Abb., DM 26,40. US $ 13.20 , 1980 .

[24]  Linus Schrage,et al.  A guide to simulation , 1983 .

[25]  Yu-Chi Ho,et al.  Functional Estimation with Respect to a Threshold Parametervia Dynamic Split-and-Merge , 1997 .

[26]  R. Suri,et al.  Perturbation analysis: the state of the art and research issues explained via the GI/G/1 queue , 1989, Proc. IEEE.

[27]  Michael C. Fu,et al.  Optimization via simulation: A review , 1994, Ann. Oper. Res..

[28]  M. Johnson,et al.  Sensitivity analysis of serial transfer lines using finite perturbation analysis , 1989 .

[29]  George E. P. Box,et al.  Empirical Model‐Building and Response Surfaces , 1988 .

[30]  Xi-Ren Cao,et al.  The phantom customer and marked customer methods for optimization of closed queueing networks with blocking and general service times , 1983, SIGMETRICS '83.

[31]  Hans Daduna,et al.  Delay time distributions and adjusted transfer rates for Jackson networks , 1993 .

[32]  Bernd Heidergott Infinitesimal perturbation analysis for queueing networks with general service time distributions , 1999, Queueing Syst. Theory Appl..

[33]  Harold J. Kushner,et al.  wchastic. approximation methods for constrained and unconstrained systems , 1978 .

[34]  P. Glynn A GSMP formalism for discrete event systems , 1989, Proc. IEEE.

[35]  Leyuan Shi Discontinuous perturbation analysis of discrete-event dynamic systems , 1996, IEEE Trans. Autom. Control..

[36]  Xi-Ren Cao,et al.  First-Order Perturbation Analysis of a Simple Multi-Class Finite Source Queue , 1987, Perform. Evaluation.

[37]  Jean-Marc Lasgouttes,et al.  Stationary IPA estimates for nonsmooth G/G/1/∞ functionals via palm inversion and level-crossing analysis , 1993, Discret. Event Dyn. Syst..

[38]  P. Brémaud,et al.  Virtual customers in sensitivity and light traffic analysis via Campbell's formula for point processes , 1993, Advances in Applied Probability.

[39]  Jean-Marc Lasgouttes,et al.  Stationary IPA estimates for non-smooth functions of the GI/G/1/infini workload , 1992 .

[40]  P. Meyer,et al.  Probabilities and potential C , 1978 .

[41]  Paul Bratley,et al.  A guide to simulation (2nd ed.) , 1986 .

[42]  Pierre L'Ecuyer,et al.  Functional Estimation with Respect to a Threshold Parameter via Dynamic Split-and-Merge , 1997, Discret. Event Dyn. Syst..

[43]  Yu-Chi Ho,et al.  A gradient technique for general buffer storage design in a production line , 1979 .

[44]  George Ch. Pflug,et al.  Optimization of Stochastic Models , 1996 .

[45]  Christos G. Cassandras,et al.  Infinitesimal and finite perturbation analysis for queueing networks , 1982, 1982 21st IEEE Conference on Decision and Control.

[46]  Martin I. Reiman,et al.  Simterpolations: estimating an entire queueing function from a single sample path , 1987, WSC '87.