A model for simulating the transient electronic absorption spectra of donor acceptor dyads undergoing ultrafast intramolecular charge transfer in solution has been developed. It is based on the stochastic multichannel transition approach and includes the reorganization of high frequency intramolecular modes (treated quantum ically) and of low frequency intramolecular and solvent modes (described classically). The relaxation of the slow modes is assumed to be exponential with time constants taken from experiments. The excited-state dynamics is obtained by simulating the population distribution of each quantum state after optical excitation and upon electronic and vibrational transitions. This model was used to simulate the transient electronic absorption spectra measured previously with a pyrylium phenolate in acetonitrile. A very good agreement between the simulated and measured spectra was obtained assuming a three-level model including the ground state, the optically excited state, and a dark state with large charge-transfer character and a substantially different geometry relative to that of the optically excited state. The merit of this approach to disentangle the contributions of both population changes and relaxation processes to the ultrafast spectral dynamics will be discussed.