Die Matrix \( A \) ist hermitesch und somit auch normal. Der Spektralsatz sichert dir die Existenz einer unitären Matrix \( U \) mit
$$ A = U^*DU $$
Wobei \( D \) eine Diagonalmatrix ist, die die Eigenwerte der Matrix \( A \) auf der Diagonalen besitzt.
Überlege dir mal warum \( A \) nur reelle Eigenwerte hat. Wenn du das überprüft hast, kannst du \( D \) in der Form
$$ D = \begin{pmatrix} \lambda_1 & & & \\ & \lambda_2 & &\\ &&&\ddots&\\ &&&&\lambda_n \end{pmatrix}$$
mit \( \lambda_1, ..., \lambda_n \in \mathbb{R} \) darstellen. Sind die Eigenwerte alle nicht-negativ? Wenn ja könnten wir \( D \) einfach zerlegen:
$$ D = \begin{pmatrix} \lambda_1 & & & \\ & \lambda_2 & &\\ &&&\ddots&\\ &&&&\lambda_n \end{pmatrix}= \begin{pmatrix} \sqrt{\lambda_1} & & & \\ & \sqrt{\lambda_2} & &\\ &&&\ddots&\\ &&&&\sqrt{\lambda_n} \end{pmatrix} ^ 2 =: \left(\sqrt{D}\right)^2 $$
Also
$$ A = U^* \left(\sqrt{D}\right)^2 U = (U^* \sqrt{D})(\sqrt{D}U) $$
Wenn du jetzt noch zeigst, dass
$$\sqrt{D} = \sqrt{D}^* $$
folgt doch direkt:
$$ A = (U^* \sqrt{D})(\sqrt{D}U) = (U^* \sqrt{D}^*)(\sqrt{D}U) = (\sqrt{D}U)^*(\sqrt{D}U) $$
Setze jetzt \( M:= (\sqrt{D}U)^* \), dann stehts schon da.