INV_KALMAN_DIFF inverse solver for difference EIT img= inv_kalman_diff( inv_model, data1, data2) img => output image inv_model => inverse model struct data1 => differential data at earlier time data2 => differential data at later time if inv_model.fwd_model.stimulation(:).delta_time exists and is non_zero, then the kalman filter will be applied to each data measurement separately Note that the classic Kalman filter assumes that the time step between each measurement is constant (ie as part of the state update eqn). inv_kalman_diff cannot work with non-constant time steps if inv_model.inv_kalman_diff.keep_K_k1 = 1 then img outputs img.inv_kalman_diff.K_k1 = K_k1 this can be used to estimate noise properties
0001 function img= inv_kalman_diff( inv_model, varargin) 0002 % INV_KALMAN_DIFF inverse solver for difference EIT 0003 % img= inv_kalman_diff( inv_model, data1, data2) 0004 % 0005 % img => output image 0006 % inv_model => inverse model struct 0007 % data1 => differential data at earlier time 0008 % data2 => differential data at later time 0009 % 0010 % if inv_model.fwd_model.stimulation(:).delta_time 0011 % exists and is non_zero, then the kalman filter will 0012 % be applied to each data measurement separately 0013 % 0014 % Note that the classic Kalman filter assumes that the 0015 % time step between each measurement is constant 0016 % (ie as part of the state update eqn). inv_kalman_diff 0017 % cannot work with non-constant time steps 0018 % 0019 % if inv_model.inv_kalman_diff.keep_K_k1 = 1 0020 % then img outputs img.inv_kalman_diff.K_k1 = K_k1 0021 % this can be used to estimate noise properties 0022 0023 % (C) 2005 Andy Adler. License: GPL version 2 or version 3 0024 % $Id: inv_kalman_diff.m 5701 2018-03-19 15:39:36Z aadler $ 0025 0026 warning('EIDORS:deprecated','INV_KALMAN_DIFF is deprecated as of 08-Jun-2012. Use INV_SOLVE_DIFF_KALMAN instead.'); 0027 if isfield(inv_model,'inv_kalman_diff'); 0028 inv_model.inv_solve_diff_kalman = inv_model.inv_kalman_diff; 0029 end 0030 img = inv_solve_diff_kalman(inv_model, varargin{:});