This paper presents a mathematical model for photo-excited carrier decay in a semiconductor. Due to the carrier trapping states and recombination centers in the bandgap, the carrier decay process is defined by the system of nonlinear differential equations. The system of nonlinear ordinary differential equations is approximated by linearized backward Euler scheme. Some a priori estimates of the discrete solution are obtained and the convergence of the linearized backward Euler method is proved. The identifiability analysis of the parameters of deep centers is performed and the fitting of the model to experimental data is done by using the genetic optimization algorithm. Results of numerical experiments are presented.