Optimal discrete-time Prony series fitting method for viscoelastic materials