AA_CALC_SYSTEM_MAT: SS= aa_calc_system_mat( fwd_model, img) Calc system matrix for Andy Adler's EIT code fwd_model = forward model img = image background for system matrix calc s_mat.E = CC' * SS * conductivites * CC; where: SS = Unconnected system Matrix CC = Connectivity Matrix
0001 function s_mat= aa_calc_system_mat( fwd_model, img) 0002 % AA_CALC_SYSTEM_MAT: SS= aa_calc_system_mat( fwd_model, img) 0003 % Calc system matrix for Andy Adler's EIT code 0004 % fwd_model = forward model 0005 % img = image background for system matrix calc 0006 % s_mat.E = CC' * SS * conductivites * CC; 0007 % where: 0008 % SS = Unconnected system Matrix 0009 % CC = Connectivity Matrix 0010 0011 % (C) 2005 Andy Adler. License: GPL version 2 or version 3 0012 % $Id: aa_calc_system_mat.m 3087 2012-06-07 18:11:42Z aadler $ 0013 0014 warning('EIDORS:deprecated','AA_CALC_SYSTEM_MAT is deprecated as of 07-Jun-2012. Use SYSTEM_MAT_1ST_ORDER instead.'); 0015 0016 FC= system_mat_fields( fwd_model); 0017 lFC= size(FC,1); 0018 0019 elem_sigma = kron( img.elem_data(:), ones(elem_dim(fwd_model),1) ); 0020 0021 ES= ones(lFC,1); 0022 ES(1:length(elem_sigma))= elem_sigma; 0023 ES= spdiags(ES,0,lFC,lFC); 0024 0025 s_mat.E= FC' * ES * FC; 0026