We consider $n$ discrete random variables $Y_1, \ldots, Y_n$ with random vector $Y$ and covariance matrix $\Sigma_Y$. We denote by $P$ the change of basis matrix from the canonical basis of $\mathcal{M}_{n,1}(\mathbb{R})$ to an orthonormal basis formed by eigenvectors of $\Sigma_Y$. We define the discrete random variable $X = P^\top Y = \left(\begin{array}{c} X_1 \\ \vdots \\ X_n \end{array}\right)$.
Prove that the total variance of $X$ is equal to that of $Y$.