We investigate the charged particle spectra produced in the heavy-ion collisions at nine centralities from different systems, i.e., Pb+Pb at $sqrt{s_{NN}}=2.76$ TeV and 5.02 TeV as well as Xe+Xe at $sqrt{s_{NN}}=5.44$ TeV, at Large Hadron Collider (LHC) using one empirical formula inspired by the solution of the Fokker-Planck equation, dubbed as the generalized Fokker-Planck solution (GFPS). Our results show that the GFPS can reproduce the experimental particle spectrum up to transverse momentum $p_T$ about 45 GeV/c with the maximum discrepancy 30% covering 10 orders of magnitude. The discrepancy between the data and the results from the GFPS decreases to 15% when the maximum of the charged particle transverse momentum is cut to 20 GeV/c. We confirmed that the Tsallis distribution derived from the non-extensive statistics, which can reproduce the particle spectra produced in small collision systems, such as p+p, up to few hundreds GeV/c, can only apply to systematically study the particle spectra up to 8 GeV/c in A+A collisions at LHC, as pointed out in the study of identified particle spectra in Pb+Pb collisions at $sqrt{s_{NN}}=2.76$ TeV. The possible explanation why GFPS functions well is also discussed.