# MatrixPropertyDistribution

MatrixPropertyDistribution[expr,xmdist]

represents the distribution of the matrix property expr where the matrix-valued random variable x follows the matrix distribution mdist.

MatrixPropertyDistribution[expr,{x1mdist1,x2mdist2,}]

represents the distribution where x1, x2, are independent and follow the matrix distributions mdist1, mdist2, .

# Details • MatrixPropertyDistribution is a transformation from the space of matrices to some property that is often of much lower dimension.
• MatrixPropertyDistribution is typically used to study properties of a distribution of matrices such as eigenvalues, singular values, determinant, norm or any property that can be computed.
• xdist can be entered as x dist dist or x[Distributed]dist.
• MatrixPropertyDistribution can be used with such functions as NProbability, NExpectation, and RandomVariate.

# Examples

open allclose all

## Basic Examples(3)

Approximate the mean of for Gaussian orthogonal matrix :

Draw a sample solution of a random linear system:

Estimate distribution of Log10 of condition number of a random matrix:

## Scope(3)

Define distribution of a scalar-valued function of matrix argument:

Approximate the mean of the function:

Define distribution of a vector-valued function of matrix argument:

Sample from the distribution:

Define distribution from a random matrix and a random vector:

Approximate quartiles of the distribution:

## Applications(4)

Sample determinant of matrix from GaussianOrthogonalMatrixDistribution:

Compare the histogram to the known PDF:

Estimate the spectral density of matrix from GaussianUnitaryMatrixDistribution:

The closed form is known to be the following:

For smaller matrices, there is a characteristic oscillatory pattern:

Compare the histogram of the sample to the known PDF:

The number of density maxima is equal to the matrix size:

In the limit of large matrices, the density converges to WignerSemicircleDistribution:

The zeros of the Riemann zeta function have been conjectured to be related to the eigenvalues of Hermitian operators and matrices. Compare the normalized spacing of the zeros to the normalized spacing of the bulk eigenvalues of samples from GaussianUnitaryMatrixDistribution:

Compare the histogram of normalized spacings to the known PDF:

Compare this PDF to the normalized spacing for the zeros of the zeta function in the critical line for a series of the zeros starting at the  zero (Odzlyko):

Check that they indeed are zeros:

Compare the histogram of normalized spacings to the known PDF for the random matrices:

Define distribution for scaled condition number of a WishartMatrixDistribution:

Sample the scaled condition number of a large matrix and check that it agrees with asymptotic closed-form distribution:

The asymptotic scaled condition number distribution has infinite mean:

Simulate whether LinearSolve determines the random matrix to be illconditioned:

Infer the probability that the random Wishart matrix is badly conditioned:

Use asymptotic distribution to infer the critical ratio of the largest and the smallest eigenvalues: