Modelling Volatility Surface Dynamics With Matrix Autoregressive Models