A model based on two coupled generalized Langevin equations is proposed to investigate the trans-stilbene photoisomerization dynamics. In this model, a system which has two independent coordinates is considered and these two system coordinates are coupled to the same harmonic bath. The direct coupling between the system coordinates is assumed negligible and these two coordinates influence each other through the frictional coupling mediated by solvent molecules. From the Hamiltonian which is equivalent to the coupled generalized Langevin equations, we obtain the transition state theory rate constants of the stilbene photoisomerization. The rates obtained from this model are compared to experimental results in n-alkane solvents.