Monday, May 14, 2012

4:15 pm

## Applied Mathematics Colloquium

Implicit Sampling, with Applications to Data Assimilation

Alexandre Chorin, Professor, Mathematics, UC Berkeley

It is often necessary to sample a given probability density (pdf), i.e., use a computer to construct
a sequence of independent random vectors x_i, i=1,2,..., whose histogram converges to the given pdf.
This can be difficult because the sample space can be huge, and more important, because
the portion of the space where the density is significant can be very small, so that
one may miss it by an ill-designed sampling scheme. Indeed, Markov-chain Monte Carlo (MCMC), the most
widely used sampling scheme, can be thought of as a search algorithm, where one starts at an
arbitrary point and one advances step-by-step towards the high probability region of the space.
This can be expensive, in particular because one is typically interested in independent samples,
while the chain has a memory.

I will present an alternative, in which samples are found by solving an algebraic equation with a random input rather than by following a chain; each sample is independent of the previous samples. I will explain the construction in the context of numerical integration, and then apply it to filtering and data assimilation.

I will present an alternative, in which samples are found by solving an algebraic equation with a random input rather than by following a chain; each sample is independent of the previous samples. I will explain the construction in the context of numerical integration, and then apply it to filtering and data assimilation.

Contact Sydney Garstang sydney@caltech.edu at x4555

For more information see http://www.acm.caltech.edu

Event Sponsors: