Gibbs sampler for computing and propagating large covariance matrices   总被引：1，自引：1，他引：0  

Gundlich  B.  Koch  K.-R.  Kusche  J. 《Journal of Geodesy》2003,77(9):514-528

The use of sampling-based Monte Carlo methods for the computation and propagation of large covariance matrices in geodetic applications is investigated. In particular, the so-called Gibbs sampler, and its use in deriving covariance matrices by Monte Carlo integration, and in linear and nonlinear error propagation studies, is discussed. Modifications of this technique are given which improve in efficiency in situations where estimated parameters are highly correlated and normal matrices appear as ill-conditioned. This is a situation frequently encountered in satellite gravity field modelling. A synthetic experiment, where covariance matrices for spherical harmonic coefficients are estimated and propagated to geoid height covariance matrices, is described. In this case, the generated samples correspond to random realizations of errors of a gravity field model. AcknowledgementsThe authors are indebted to Pieter Visser and Pavel Ditmar for providing simulation output that was used in the GOCE error generation experiments. Furthermore, the NASA/NIMA/OSU team is acknowledged for providing public ftp access to the EGM96 error covariance matrix. The two anonymous reviewers are thanked for their valuable comments.  相似文献   

Confidence regions for GPS baselines by Bayesian statistics     

B. Gundlich  K.-R. Koch 《Journal of Geodesy》2002,76(1):55-62

 The global positioning system (GPS) model is distinctive in the way that the unknown parameters are not only real-valued, the baseline coordinates, but also integers, the phase ambiguities. The GPS model therefore leads to a mixed integer–real-valued estimation problem. Common solutions are the float solution, which ignores the ambiguities being integers, or the fixed solution, where the ambiguities are estimated as integers and then are fixed. Confidence regions, so-called HPD (highest posterior density) regions, for the GPS baselines are derived by Bayesian statistics. They take care of the integer character of the phase ambiguities but still consider them as unknown parameters. Estimating these confidence regions leads to a numerical integration problem which is solved by Monte Carlo methods. This is computationally expensive so that approximations of the confidence regions are also developed. In an example it is shown that for a high confidence level the confidence region consists of more than one region. Received: 1 February 2001 / Accepted: 18 July 2001  相似文献