The fractal dimension is an important parameter of detecting and characterizing chaos produced from a dynamical system. According to the relationship between time scale and multidimensional super-body’s volume, a method to compute the fractal dimension of multivariate time series is proposed. Numerical simulations show that the length of the time series is shorter, the method is more efficient and doesn’t need reconstructing phase space so as to avoid the infection from parameters determination such as delay time and embedding dimension. Furthermore, it can resist noise to a certain degree and can be used as an effective method to calculate the fractal dimension of multivariate time series in many areas.