1 function A = random_positive_definite_matrix(integerN)
3 % Generate a random, symmetric positive-definite (SPD) matrix.
5 % Since all (SPD) matrices are diagonalizable and have positive
6 % eigenvalues, it seems likely that we can generate an SPD matrix by
7 % combining random orthogonal matrices with a diagonal matrix of
10 % I have no proof/evidence that this approach is sound.
14 % - ``integerN`` -- The dimension of the resulting matrix.
18 % - ``A`` -- A symmetric, positive definite matrix.
20 U = random_orthogonal_matrix(integerN);
21 d = unifrnd(eps, realmax, 1, integerN);