Numerical evaluation of the Lambert W function and application to generation of generalized Gaussian noise with exponent 1/2

  • Authors:
  • F. Chapeau-Blondeau;A. Monir

  • Affiliations:
  • Lab. d'Ingenierie des Systemes Automatises, Univ. d'Angers, Angers;-

  • Venue:
  • IEEE Transactions on Signal Processing
  • Year:
  • 2002

Quantified Score

Hi-index 35.68

Visualization

Abstract

We address the problem of synthesizing a generalized Gaussian noise with exponent 1/2 by means of a nonlinear memoryless transformation applied to a uniform noise. We show that this transformation is expressable in terms of a special function known under the name of the Lambert W function. We review the main methods for numerical evaluation of the relevant branch of the (multivalued) Lambert W function with controlled accuracy and complement them with an original rational function approximation. Based on these methods, synthesis of the generalized Gaussian noise can be performed with arbitrary accuracy. We construct a simple and fast evaluation algorithm with prescribed accuracy, which is especially suited for Monte Carlo simulation requiring large numbers of realizations of the generalized Gaussian noise.