Polynomial filtering for linear discrete time non-Gaussian systems