Implicit Sampling, with Application to Data Assimilation |
| |
Authors: | Alexandre J. CHORIN Matthias MORZFELD Xuemin TU |
| |
Affiliation: | 1. Department of Mathematics, University of California, Berkeley, CA, 94720, USA; Lawrence Berkeley National Laboratory, CA, 94720, USA 2. Lawrence Berkeley National Laboratory, CA, 94720, USA 3. Department of Mathematics, University of Kansas, Lawrence, KS, 66045, USA |
| |
Abstract: | There are many computational tasks, in which it is necessary to sample a given probability density function (or pdf for short), i.e., to 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 importantly, 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, Markovchain Monte Carlo, 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. The authors present an alternative, in which samples are found by solving an algebraic equation with a random right-hand side rather than by following a chain; each sample is independent of the previous samples. The construction in the context of numerical integration is explained, and then it is applied to data assimilation. |
| |
Keywords: | Importance sampling Bayesian estimation Particle filter Implicit filter Data assimilation |
本文献已被 CNKI 维普 万方数据 SpringerLink 等数据库收录! |
| 点击此处可从《数学年刊B辑(英文版)》浏览原始摘要信息 |
|
点击此处可从《数学年刊B辑(英文版)》下载全文 |
|