Abstract To jointly track frequency drift (or fine frequency offset) and channel state variation for orthogonal frequency division multiplexing (OFDM) communications over mobile wireless channels, the maximum-likelihood (ML) estimation techniques are commonly adopted. A major difficulty arises from the highly nonlinear nature of the log-likelihood function which renders local extrema or multiple solutions. The use of mathematical approximation coupled with an adaptive iteration algorithm is a viable approach to ease problem solving. The approximation methods used in existing works are usually of the first-order level. In this work, we devise a high-order approximation algorithm to improve the tracking performance. As shall be seen, the problem amounts to finding the roots of an approximate high-degree polynomial equation derived from the log-likelihood function. To this end, we triangularize the companion matrix constructed from the high-degree polynomial using iterative Gram-Schmidt QR transformations, an eigenvalue problem in matrix computation, coupled with another iterative correction process to produce a frequency offset estimate to good accuracy. It is found that a high-order approximation algorithm formulated can achieve much better tracking performance in terms of estimation accuracy, tracking range, and error rate results than various other tracking algorithms appearing in the literature.