Abstract: The use of the finite-difference (FD) optimized coefficients algorithm to suppress numerical dispersion is a straightforward and efficient strategy in the field of FD forward modeling.For the current optimal algorithms, Remez exchange algorithm (REA) has a broader effective bandwidth than other methods, but suffers stronger numerical dispersion in low-wavenumber regions.To reduce the numerical dispersion in low-wavenumber region and maintain the characteristic of broad effective bandwidth, an optimized FD coefficients method coupled with REA and Lagrange algorithm (LA) was proposed.First, the FD coefficients are obtained using the REA, and the zero points of the curve are then derived using the dispersion relation curve of the FD coefficients.Subsequently, the zero points and dispersion relation curve are used to construct the restricted condition.Finally, the optimized FD coefficients are obtained using the LA to solve the L₂-norm objective function and the restricted condition.Theoretical analysis and numerical experiment reveal that the proposed method has a low dispersion error in low-wavenumber regions and wide effective bandwidth.Under a low space step, the numerical experiment of modified Marmousi model demonstrated that the proposed method has a lower numerical dispersion compared with that of the REA and least square methods.