We consider the problem of estimating signal priors, which are used in
the -regularization and the resulting inverse problem has a
stabilized solution. From a Bayesian viewpoint, we first define a
multivariate
-Laplace density function
and then solve a maximum likelihood problem with an added
-norm penalty term. The problem as formulated is convex but
the memory requirements and the nonlinear non-smooth sub-gradient
equations are prohibitive for large-scale problems. We develop an
iterative algorithm to efficiently solve such large problems and
demonstrate the selected priors generally behave better than those
commonly used ones in the signal processing.