I am trying to find properties or constraints on a $(p \times n) $matrix $U$ such that upon left multiplying a real symmetric square positive semi-definite matrix $V$ with $U$ the resulting matrix $W$ of dimensions $(p \times n) $ is still real positive semi-definite matrix.
In other words i am trying to establish that there is $U$ such that
$$W = U V$$
where $W$ and $V$ are real positive semi-definite matrices.
If $W=V$ then obviously $U=I$ fits the equation.
I have a suspicion that this should be possible as by definition of real symmetric positive semi-definite matrix, $V$ should have non-negative eigenvalues, which implies that it should have a right inverse which I suspect would also be positive semi-definite matrix, so if $W$ (some positive semi-definite matrix) is multiplied with the right inverse we should be able to compute $U$.
Any help would be much appreciated.
edit
I suspect if my reasoning about computing $U$ above is correct then it should also be positive semi-definite matrix. So perhaps a way to construct $U$ may be by applying Gram-Schmidt process and finding new eigenvectors for W.
Best Answer
I am unsure how we should properly resolve this question.
As we have discussed at length in the comments, only square matrices can be positive semidefinite. Therefore, if $p\neq n$, there is no way that the product matrix $W=UV$ can be positive semidefinite, because it will also be non-square.
For grins, let's assume that $p=n$. Under what conditions is $W=UV$ positive semidefinite? This requires that $x^HUVx\geq 0$ for all complex vectors $x$. Alternatively, this is true if and only if $UV+VU^H=Q$ where $Q=W+W^H$ has nonnegative eigenvalues (and is therefore PSD itself). If $U$, $V$ are real, then you can relax the Hermitian transposes to real tranposes, and consider only reall vectors $x$.
Now for one special case, I know the answer. If $V$ is positive definite---i.e., not just PSD but nonsingular---and we require $W$ to be positive definite as well, then the answer follows from Lyapunov's theorem applied to linear systems:
What about the more relaxed cases? That is, what if $V$ is only positive semidefinite? What if $W$ is only required to be positive semidefinite? I am afraid I do not know. I'm sure people who study Lyapunov's theorem for linear systems in some depth know...