Modelling matrix time series via a tensor CP-decomposition