]> Generate Jacobian for unknown coupling coefficients

Generate Jacobian for Unknown Coupling Coefficients

genSDJacobian() generates the Jacobian for the (presumably) unknown source-detector coupling coefficients. See the code for fitSD() and newFitBackground() for examples of its use. Most users will want to call these functions directly, however, instead of calculating the Jacobians directly.

Function Summary

Syntax: J = genSDJacobian(SD, Medium, MeasList, data, Method);
Inputs: SD SD structure
Medium Medium structure
MeasList Measurement List that corresponds to experimental data
data Experimental data to fit against
Method How the Jacobian should be calculated, "Born" or "Rytov"
Outputs: J The Jacobian

Detailed Descriptions

The Jacobian (and the system of linear equations that it represents) is naturally rank deficient. One column of the matrix should be set to 0 before inverting.

For more information see:

David A. Boas, Thomas Gaudette, and Simon R. Arridge, "Simultaneous imaging and optode calibration with diffuse optical tomography," Optics Express, 8(5), pg. 263, 2001.