The model is derived in a natural way from the theory of coupled modes. It is based on an expansion of the Hamiltonian in terms of elementary transition operators, including direct rotation-vibration coupling with phonons. The treatment is limited to three types of phonons: ( I = K = 0), S (I = K = 1) and (I = K = 2). The basis of the operators, acting on the ground state is truncated by an inclusion of a reasonable number of phonon states. In the framework of this approximation one may evaluate the matrix elements of the model Hamiltonian and diagonalize it by standard numerical methods to fit the experimental spectrum. The well known picture of band hybridization is obtained as a special case of the model under consideration.