The simple answer is that the solution of Schrödinger's equation (in abstract notation):
[tex]i \frac{d}{dt} |\Psi> = H|\Psi>[/tex]
can be written as:
[tex]|\Psi(t)> = e^{-iHt} |\Psi(0)>[/tex]
The operator [itex]e^{-iHt}[/itex] is a unitary operator, and is called the time evolution operator, since it takes a state at time t' to time t+t'. Differentiating the last equation (treating H as an ordinary number) shows it is a solution.
More generally, any symmetry of the system is represented by a unitary operator on its Hilbert space (here the symmetry is time-translation invariance, ie, the outcome of an experiment is independent of when it is run). This is because we expect a symmetry to have no effect on the transition probabilities between various states, which means it should preserve the inner product on the Hilbert space, and this is precisely what a unitary operator does (you can think of unitary operators as the generalization to Hilbert spaces of orthogonal transformations (ie, rotations)). Then it is just a mathematical fact that unitary operators are of the form [itex]e^{iH}[/itex], where H is a Hermitian operator. In one (complex) dimension, this is just the statement that complex numbers of unit magnitude (1D unitary operators) are all of the form [itex]e^{it}[/itex], where t is a real number (1D Hermitian operator).