Variational approximation for heteroscedastic linear models and matching pursuit algorithms