Based on the methods of complex function, conformal mapping, and multipolar coordinate system, dynamic response of an elliptical inclusion embedded in an anisotropic half space is investigated. In order to find the solution of SH waves, the governing equation is transferred into its normalized form. Then, the scattering wave induced by the inclusion and the standing wave in the inclusion is deduced. Different incident wave angles and the corresponding anisotropy of the half space are considered to obtain the reflected waves. The elliptical inclusion is transferred into a unit circle by conformal mapping method, and then the undetermined coefficients in scattering wave and standing wave are solved by using the continuous condition at the boundary of the inclusion. Subsequently, the dynamic stress concentration factor (DSCF) around the inclusion is calculated and analyzed. Numerical results demonstrate that the distribution of the DSCF is mainly influenced by the incident wave angle and the incident wave number. It is affected by anisotropic parameters as well.