From the definition of eigenvectors of a (square) matrix $A$, I can obtain the result: $$ A=P\Omega P^{-1} $$ where $P$ is unitary matrix (columns are eigenvectors), and $\Omega $ is a diagonal matrix of eigenvalues.
But I am wondering a very different situation: is it possible to decompose a given matrix $A$ as: $$ A=M \Omega M^{-1} $$ where $M$ is NOT unitary (the columns are not orthogonal). In other words, how to find a non-orthogonal matrix $M$ making $A=M \Omega M^{-1}$ holds ($\Omega$ is diagonal), given an arbitrary square matrix $A$?
That $P$ is unitary is a relativly rare case! Always happens, when $A$ is "normal", which means $AA^H = A^H A$. It can occure more often, but I#m not sure about that.
One example to make a difference is by using $$ A = \begin{pmatrix}1 &2 \\ 0 & 1\end{pmatrix} \begin{pmatrix} 2 & 0 \\ 0 &5\end{pmatrix}\begin{pmatrix}1 &2 \\ 0 & 1\end{pmatrix}^{-1} = \begin{pmatrix}2&6\\0&5 \end{pmatrix} $$ The collumns of $P$ are not orthogonal.
Keep in mind, that it always depends on $A$ and it's eigenvalues and eigenvectors. If the vectors are orthogonal, you can find an orthonormal basis. If not, you can not. However you can just judge this by calculationg the eigenvalues.