<p>This paper presents a higher-order conservative finite difference scheme for solving the coupled nonlinear Schrödinger-Boussinesq equations. The scheme adopts a novel sixth-order spatial difference operator, along with the Crank-Nicolson method and Richardson extrapolation technique for temporal discretization, yielding sixth-order accuracy in space and fourth-order accuracy in time. Subsequently, the discrete mass conservation, energy conservation, and convergence properties of this scheme are analyzed one by one. Furthermore, a fast algorithm is designed to improve the computational efficiency of the proposed method. Numerical simulations are conducted to validate the theoretical findings, which demonstrate the scheme’s high accuracy and reliability.</p>