Objective: Due to the non-linearity of numerous biomedical signals, non-linear analysis of multi-channel time series, notably multivariate multiscale entropy (mvMSE), has been extensively used in biomedical signal processing. However, mvMSE has three drawbacks: 1) mvMSE values are either undefined or unreliable for short signals; 2) mvMSE is not fast enough for real-time applications; and 3) the computation of mvMSE for signals with a large number of channels requires the storage of a huge number of elements. Methods: To deal with these problems and improve the stability of mvMSE, we introduce multivariate multiscale dispersion entropy (MDE -mvMDE), as an extension of our recently developed MDE, to quantify the complexity of multivariate time series. Results: We assess mvMDE, in comparison with mvMSE and multivariate multiscale fuzzy entropy (mvMFE), on correlated and uncorrelated multichannel noise signals, bivariate autoregressive processes, and three biomedical datasets. The results show that mvMDE takes into account dependencies in patterns across both the time and spatial domains. The mvMDE, mvMSE, and mvMFE methods are consistent in that they lead to similar conclusions about the underlying physiological conditions. However, the proposed mvMDE discriminates various physiological states of the biomedical recordings better than mvMSE and mvMFE. In addition, for both the short and long time series, the mvMDE-based results are noticeably more stable than the mvMSE-and mvMFE-based ones. Conclusion: For short multivariate time series, mvMDE, unlike mvMSE, does not result in undefined values. Furthermore, mvMDE is noticeably faster than mvMFE and mvMSE and also needs to store a considerably smaller number of elements. Significance: Due to its ability to detect different kinds of dynamics of multivariate signals, mvMDE has great potential to analyse various physiological signals.