In this paper,we propose a method for theoretically simulating the response functions and then the 2D,third-order,non-linear,photon response spectra obtained from photon echo experiments of some complex molecular system.Here,we supposed that the exciton modes of the molecular system are coupled to their environments which are described with damped harmonic oscillators,and the Hamiltonian of the system-bath can be expressed into the normally ordered form.