I read this problem from Kostrikin's algebra problem book. Here is my solution:
By complex spectral theorem, we know the normal operator $A$ can be diagonalised with respect to an orthonormal basis consisting of $A$'s eigenvectors, therefore we have \begin{equation*} \mathcal{M}(A)= \begin{bmatrix} \lambda_{1} && 0\\ & \ddots &\\ 0 && \lambda_{n} \end{bmatrix} \end{equation*} Given $AB=BA$, we have \begin{equation*} AB_{ij}=\sum_{k=1}^{n}a_{ik}b_{kj}=\lambda_{i}b_{ij} \end{equation*} \begin{equation*} BA_{ij}=\sum_{k=1}^{n}b_{ik}a_{kj}=\lambda_{j}b_{ij} \end{equation*} Then we yield $\lambda_{i}b_{ij}=\lambda_{j}b_{ij}$, which just implies $\lambda_{i}=\lambda_{j}$. We conclude that $A$ is just a scalar multiple of the identity matrix. Since identity matrix commutes with any matrix, we conclude that $AB^{*}=B^{*}A$. Q.E.D
Is my idea correct?