Ergodic approximations via matrix regularization approach