Note: This content is accessible to all versions of every browser. However, this browser does not seem to support current Web standards, preventing the display of our site's design details.


Simulation of stochastic systems via polynomial chaos expansions and convex optimization


L. Fagiano, M. Khammash

Physical Review E., no. 036702

Polynomial Chaos Expansions represent a powerful tool to simulate stochastic models of dynamical systems. Yet, deriving the expansionís coefficients for complex systems might require a significant and non-trivial manipulation of the model, or the computation of large numbers of simulation runs, rendering the approach too time consuming and impracticable for applications with more than a handful of random variables. We introduce a novel computationally tractable technique for computing the coefficients of polynomial chaos expansions. The approach exploits a regularization technique with a particular choice of weighting matrices, which allow to take into account the specific features of Polynomial Chaos expansions. The method, completely based on convex optimization, can be applied to problems with a large number of random variables and uses a modest number of Monte Carlo simulations, while avoiding model manipulations. Additional information on the stochastic process, when available, can be also incorporated in the approach by means of convex constraints. We show the effectiveness of the proposed technique in three applications in diverse fields, including the analysis of a nonlinear electric circuit, a chaotic model of organizational behavior, finally a chemical oscillator.


Type of Publication:


No Files for download available.
% Autogenerated BibTeX entry
@Article { FagKha:2012:IFA_4445,
    author={L. Fagiano and M. Khammash},
    title={{Simulation of stochastic systems via polynomial chaos
	  expansions and convex optimization}},
    journal={Physical Review E.},
Permanent link