Stochastic Optimization by Simulation: Some Experiments with a Simple Steady-State Queue

Abstract : New approaches like perturbation analysis and the likelihood ratio method have been proposed recently to estimate the gradient of a performance measure with respect to some continuous parameters in a dynamic stochastic system. In this paper, we experiment the use of these estimators in stochastic approximation algorithms, to perform so-called single-run optimizations. We also compare them to finite difference estimators, with and without common random numbers. The experiments are done on a simple M/M/1 queue. The performance measure involves the average system time per customer, and the optimal solution is easy to compute analytically, which facilitates the evaluation of the algorithms. We also demonstrate some properties of the algorithms. In particular, we show that using perturbation analysis, the single-run optimization converges to the optimum even with a fixed (and small) number of ends of service per iteration, while under the same conditions, the algorithm that uses the finite difference estimators converges to the wrong answer.

[1]  H. Kesten Accelerated Stochastic Approximation , 1958 .

[2]  J. Talavage,et al.  Optimization of stochastic simulation models , 1980 .

[3]  H. Kushner,et al.  An Invariant Measure Approach to the Convergence of Stochastic Approximations with State Dependent Noise. , 1984 .

[4]  M. Metivier,et al.  Applications of a Kushner and Clark lemma to general classes of stochastic algorithms , 1984, IEEE Trans. Inf. Theory.

[5]  Peter W. Glynn,et al.  Proceedings of Ihe 1986 Winter Simulation , 2022 .

[6]  Reuven Y. Rubinstein,et al.  Monte Carlo Optimization, Simulation and Sensitivity of Queueing Networks , 1986 .

[7]  R Y Rubinstein,et al.  The score function approach for sensitivity analysis of computer simulation models , 1986 .

[8]  Peter W. Glynn,et al.  Optimization of stochastic systems , 1986, WSC '86.

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

[10]  Rajan Suri,et al.  Infinitesimal perturbation analysis for general discrete event systems , 1987, JACM.

[11]  Peter W. Glynn,et al.  Likelilood ratio gradient estimation: an overview , 1987, WSC '87.

[12]  Marc S. Meketon,et al.  Optimization in simulation: a survey of recent results , 1987, WSC '87.

[13]  M. Métivier,et al.  Théorèmes de convergence presque sure pour une classe d'algorithmes stochastiques à pas décroissant , 1987 .

[14]  Xi-Ren Cao,et al.  Convergence properties of infinitesimal perturbation analysis , 1988 .

[15]  R. Suri,et al.  Perturbation analysis gives strongly consistent sensitivity estimates for the M/G/ 1 queue , 1988 .

[16]  Rajan Suri,et al.  Single Run Optimization of Discrete Event Simulations—An Empirical Study Using the M/M/l Queue , 1989 .