Sampling Uniformly From the Set of Positive Definite Matrices With Trace Constraint

View Researcher's Other Codes

Disclaimer: The provided code links for this paper are external links. Science Nest has no responsibility for the accuracy, legality or content of these links. Also, by downloading this code(s), you agree to comply with the terms of use as set out by the author(s) of the code(s).

Please contact us in case of a broken link from here

Authors Martin Mittelbach, B. Matthiesen, E. Jorswieck
Journal/Conference Name IEEE Transactions on Signal Processing
Paper Category
Paper Abstract We derive a parameterization of positive definite matrices using the Cholesky decomposition in combination with hyperspherical coordinates. Based on the parameterization we develop a simple and efficient method to randomly generate positive definite matrices with constant or bounded trace according to a uniform distribution. Further, we present an efficient implementation using the inversion method and either rejection sampling or transforming a beta distribution. The matrix parameterization might be of independent interest, whereas the random sampling algorithm finds applications in Monte Carlo simulations, testing of algorithms, and performance studies. With the help of an abstract example we describe how the sampling method can be used to approximate the optimum in a difficult, e.g., nonconvex, optimization problem for which no solution or efficient global optimization algorithm is known. In this paper we consider real as well as complex matrices.
Date of publication 2012
Code Programming Language Matlab

Copyright Researcher 2022