A Higher Order Unscented Transform

06/24/2020
by   Deanna Easley, et al.
0

We develop a new approach for estimating the expected values of nonlinear functions applied to multivariate random variables with arbitrary distributions. Rather than assuming a particular distribution, we assume that we are only given the first four moments of the distribution. The goal is to summarize the distribution using a small number of quadrature nodes which are called σ-points. We achieve this by choosing nodes and weights in order to match the specified moments of the distribution. The classical scaled unscented transform (SUT) matches the mean and covariance of a distribution. In this paper, introduce the higher order unscented transform (HOUT) which also matches any given skewness and kurtosis tensors. It turns out that the key to matching the higher moments is the rank-1 tensor decomposition. While the minimal rank-1 decomposition is NP-complete, we present a practical algorithm for computing a non-minimal rank-1 decomposition and prove convergence in linear time. We then show how to combine the rank-1 decompositions of the moments in order to form the σ-points and weights of the HOUT. By passing the σ-points through a nonlinear function and applying our quadrature rule we can estimate the moments of the output distribution. We prove that the HOUT is exact on arbitrary polynomials up to fourth order. Finally, we numerically compare the HOUT to the SUT on nonlinear functions applied to non-Gaussian random variables including an application to forecasting and uncertainty quantification for chaotic dynamics.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
05/13/2022

Distribution Estimation for Probabilistic Loops

We present an algorithmic approach to estimate the value distributions o...
research
11/10/2019

Estimating Higher-Order Moments Using Symmetric Tensor Decomposition

We consider the problem of decomposing higher-order moment tensors, i.e....
research
09/02/2020

An Approximation Scheme for Multivariate Information based on Partial Information Decomposition

We consider an approximation scheme for multivariate information assumin...
research
03/23/2020

Moment State Dynamical Systems for Nonlinear Chance-Constrained Motion Planning

Chance-constrained motion planning requires uncertainty in dynamics to b...
research
08/16/2008

Higher Order Moments Generation by Mellin Transform for Compound Models of Clutter

The compound models of clutter statistics are found suitable to describe...
research
11/15/2018

A Schur transform for spatial stochastic processes

The variance, higher order moments, covariance, and joint moments or cum...
research
06/29/2021

Tensor decomposition of higher-order correlations by nonlinear Hebbian plasticity

Biological synaptic plasticity exhibits nonlinearities that are not acco...

Please sign up or login with your details

Forgot password? Click here to reset