Statistical Inference for High-Dimensional Matrix-Variate Factor Models