Abstract

We propose a new approach to weighting initial parameter misfits in a least squares optimization problem for linear parameter estimation. Parameter misfit weights are found by solving an optimization problem which ensures the penalty function has the properties of a χ ² random variable with n degrees of freedom, where n is the number of data. This approach differs from others in that weights found by the proposed algorithm vary along a diagonal matrix rather than remain constant. In addition, it is assumed that data and parameters are random, but not necessarily normally distributed.

The proposed algorithm successfully solved three benchmark problems, one with discontinuous solutions. Solutions from a more idealized discontinuous problem show that the algorithm can successfully weight initial parameter misfits even though the two-norm typically smoothes solutions. For all test problems sample solutions show that results from the proposed algorithm can be better than those found using the L-curve and generalized cross-validation. In the cases where the parameter estimates are not as accurate, their corresponding standard deviations or error bounds correctly identify their uncertainty.