Modeling multivariate ultra-high-frequency financial data by Monte Carlo simulation methods