This function computes the coefficient matrix Ck, which is a matrix of constants that allows us to obtain E[pλ(W)Wr], where r+∣λ∣=k and W∼Wmβ(n,Σ).
qk_coeff(k, alpha =2)
Arguments
k: The order of the Ck matrix
alpha: The type of Wishart distribution (α=2/β):
1/2: Quaternion Wishart
1: Complex Wishart
2: Real Wishart (default)
Returns
Ck, a matrix that allows us to obtain E[pλ(W)Wr], where r+∣λ∣=k and W∼Wmβ(n,Σ). The matrix is represented as a 3-dimensional array where each slice along the third dimension represents a coefficient matrix of the polynomial in descending powers of n.
Examples
# Example 1:qk_coeff(2)# For real Wishart distribution with k = 2# Example 2:qk_coeff(3,1)# For complex Wishart distribution with k = 3# Example 3:qk_coeff(2,1/2)# For quaternion Wishart distribution with k = 2