The problem of compressive detection of random subspace signals is studied. We consider signals modeled as $\mathbf{s} = \mathbf{H} \mathbf{x}$ where $\mathbf{H}$ is an $N \times K$ matrix with $K \le N$ and $\mathbf{x} \sim \mathcal{N}(\mathbf{0}_{K,1},σ_x^2 \mathbf{I}_K)$. We say that signal $\mathbf{s}$ lies in or leans toward a subspace if the largest eigenvalue of $\mathbf{H} \mathbf{H}^T$ is strictly greater than its smallest eigenvalue. We first design a measurement matrix $\mathbfΦ=[\mathbfΦ_s^T,\mathbfΦ_o^T]^T$ comprising of two sub-matrices $\mathbfΦ_s$ and $\mathbfΦ_o$ where $\mathbfΦ_s$ projects the signals to the strongest left-singular vectors, i.e., the left-singular vectors corresponding to the largest singular values, of subspace matrix $\mathbf{H}$ and $\mathbfΦ_o$ projects it to the weakest left-singular vectors. We then propose two detectors which work based on the difference in energies of the samples measured by two sub-matrices $\mathbfΦ_s$ and $\mathbfΦ_o$ and prove their optimality. Simplified versions of the proposed detectors for the case when the variance of noise is known are also provided. Furthermore, we study the performance of the detector when measurements are imprecise and show how imprecision can be compensated by employing more measurement devices. The problem is then re-formulated for the case when the signal lies in the union of a finite number of linear subspaces instead of a single linear subspace. Finally, we study the performance of the proposed methods by simulation examples.