文摘
A Crank-Nicolson-type difference scheme is presented for the spatial variable coefficient subdiffusion equation with Riemann-Liouville fractional derivative. The truncation errors in temporal and spatial directions are analyzed rigorously. At each time level, it results in a linear system in which the coefficient matrix is tridiagonal and strictly diagonally dominant, so it can be solved by the Thomas algorithm. The unconditional stability and convergence of the scheme are proved in the discrete \(L_{2}\) norm by the energy method. The convergence order is \(\min \{2-\frac{\alpha}{2}, 1+\alpha \}\) in the temporal direction and two in the spatial one. Finally, numerical examples are presented to verify the efficiency of our method.