Numerical study of excitonic orders in strongly correlated systems