Based on the Huygens–Fresnel integral and ABCD matrix, the propagation equation for the Lorentz–Gauss vortex beam (LGVB) in a gradient-index medium (GRIN) is rederived. The evolution of the intensity and phase distributions of an LGVB through a GRIN medium are numerically calculated as a function of the gradient-index parameter with changes in the incident beam parameters. The results showed that the propagation path and intensity distributions changed periodically with increasing propagation distance. In contrast, phase distributions change at multiples of π/β or 2π/β, depending on whether the M values are odd or even, respectively. At the same time, the parameters of the gradient index determine the periodic values of the Lorentz–Gauss vortex beams during propagation, and as β increased, the period of evolution decreased. The Lorentz–Gauss vortex beam propagating through the gradient index will develop from a square beam to a Gaussian vortex beam more quickly with an increase of w0x=w0y. In addition, the topological charge affects the size of the dark spot at the center of the beam and the size of the beam, causing the phase distributions to change periodically in the medium. This study is beneficial for laser optics and optical communications.