occam {PEIP}

Occam inversion

Description

Occam's inversion

Usage

occam(afun, ajac, L, d, m0, delta)

Arguments

afun

character, function handle that computes the forward problem

ajac

character, function handle that computes the Jacobian of the forward problem

L

regularization matrix

d

data that should be fit

m0

guess at the model

delta

cutoff to use for the discrepancy principle portion

Value

vector, model found

Note

This is a simple brute force way to do the line search. Much more sophisticated methods are available. Note: we've restricted the line search to the range from 1.0e-20 to 1. This seems to work well in practice, but might need to be adjusted for a particular problem.

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

bayes