PERTURB_JACOBIAN: J= perturb_jacobian( fwd_model, img) Calculate Jacobian Matrix, based on small perturbations in the forward model. This will tend to be slow, but should be best used to 'sanity check' other code J = Jacobian matrix fwd_model = forward model fwd_model.perturb_jacobian.delta - delta perturbation to use img = image background for jacobian calc
0001 function J= perturb_jacobian( varargin ) 0002 % PERTURB_JACOBIAN: J= perturb_jacobian( fwd_model, img) 0003 % Calculate Jacobian Matrix, based on small perturbations 0004 % in the forward model. This will tend to be slow, but 0005 % should be best used to 'sanity check' other code 0006 % 0007 % J = Jacobian matrix 0008 % fwd_model = forward model 0009 % fwd_model.perturb_jacobian.delta - delta perturbation to use 0010 % img = image background for jacobian calc 0011 0012 % (C) 2006 Andy Adler. License: GPL version 2 or version 3 0013 % $Id: perturb_jacobian.m 3289 2012-07-01 10:40:31Z aadler $ 0014 0015 warning('EIDORS:deprecated','PERTURB_JACOBIAN is deprecated as of 08-Jun-2012. Use JACOBIAN_PERTURB instead.'); 0016 0017 if isfield(inv_model,'perturb_jacobian'); 0018 inv_model.jacobian_perturb = inv_model.perturb_jacobian; 0019 end 0020 0021 J = jacobian_perturb(varargin{:});