This study concentrates on time fractional convection–diffusion equations (TFCDEs) with variable coefficients and their numerical solutions. Caputo derivative is used to calculate the time fractional order derivatives. In order to give an approximate solution to the TFCDE, an effective approach is proposed utilizing Fibonacci wavelet and block pulse functions. The Fibonacci wavelets operational matrices of fractional order integration are constructed. By combining the collocation technique, they are used to simplify the fractional model to a collection of algebraic equations. The suggested approach is quite practical for resolving issues of this nature. The comparison and analysis with other approaches demonstrate the effectiveness and precision of the suggested approach.