Da \( A \) und B kommutieren und normal sind, sind sie simultan diagonalisierbar. Es gibt also eine unitäre Matrix \( S \) und Diagonalmatrizen \(D_1, D_2\) mit
\( A=S D_{1} S^{*}, \quad B=S D_{2} S^{*} \)
Damit ergibt sich
\( \begin{aligned} A B^{*} &=\mathbf{S D}_{1} \mathbf{S}^{*}\left(\mathbf{S D}_{2} \mathbf{S}^{*}\right)^{*} \\ &=\mathbf{S D}_{1} \mathbf{S}^{*} \mathbf{S D}_{2}^{*} \mathbf{S}^{*} \\ &=\mathbf{S D}_{1} \mathbf{D}_{2}^{*} \mathbf{S}^{*} \\ &=\mathbf{S D}_{2}^{*} \mathrm{D}_{1} \mathbf{S}^{*} \\ &=\mathbf{S D}_{2}^{*} \mathbf{S}^{*} \mathbf{S D}_{1} \mathbf{S}^{*} \\ &=\left(\mathbf{S D}_{2} \mathbf{S}^{*}\right)^{*} \mathbf{S D}_{1} \mathbf{S}^{*}=\mathbf{B}^{*} \mathbf{A} \end{aligned} \)