 Solution

Let $A{v}_j=a_j{v}_j$. Commutation gives $A(B{v}_j)=B(A{v}_j)=a_{j}B{v}_j$, so $B$ preserves each eigenspace of $A$. Since the eigenvalues of $A$ are nondegenerate, each eigenspace is one-dimensional and $B{v}_j=b_j{v}_j$. The spectral theorem for a Hermitian matrix lets us choose these vectors as an orthonormal eigenbasis. With $U=(v_1,\dots ,v_n)$,This is simple spectrum gives simultaneous diagonalization. The common eigenvectors can be chosen up to their phases and ordering. In fact degeneracy does not prevent a common eigenbasis for commuting Hermitian matrices, but the one-dimensional argument already suffices here.

The convergent matrix exponential series respects similarity, so $e^A=U{e}^{{\Lambda }_A}{U}^{†}$ and likewise for $B$. Scalar exponentials multiply on each diagonal entry, givingThe essential hypothesis is commutation; Hermiticity supplies this convenient diagonalization proof.