Orthogonal, Mean-Invariant Matrix Generation Between Python and Matlab

Question

I needed a program that produces a random orthogonal matrix, which also preserves the mean and covariance of the columns of a matrix. In other words, the matrix should have