perturbation based method for approximating impuse response and ergodic distribution moments