... This is done by incorporation of the inverse weighting matrices and into the quadratic error term . The Toeplitz matrix determines a linear shift invariant regularization function: . ... Application of matrix differentiation rules and some matrix algebra show that the diagonal elements of the optimal diagonal matrices must be equal to the diagonal elements of the previously determined optimal general weighting matrices so that . ... is estimated from the measured data. ... Therefore we estimate: . ...