The numerical technique based on two-dimensional block pulse functions (2D-BPFs) is proposed for solving the time fractional convection diffusion equations with variable coefficients(FCDEs).We introduce the block pulse operational matrices of the fractional order differentiation.Furthermore,we translate the original equation into a Sylvester equation by the proposed method.Finally,some numerical examples are given and numerical results are shown to demonstrate the accuracy and reliability of the above-mentioned algorithm.