UpdateFemMatrix Assembles the system matrix for EIT Function [A]=UpdateFemMatrix(Agrad,Kb,M,S,sigma); updates the system matrix A given a new conductivity vector sigma. INPUT Agrad = the gradient part of the system matrix (see FemMatrix.m) Kb,M and S = other blocks of the system matrix (see FemMatrix.m) sigma = conductivity (or admittivity) vector OUTPUT A = the updated system matrix
0001 function [A]=UpdateFemMatrix(Agrad,Kb,M,S,sigma); 0002 0003 %UpdateFemMatrix Assembles the system matrix for EIT 0004 % Function [A]=UpdateFemMatrix(Agrad,Kb,M,S,sigma); 0005 % updates the system matrix A given a new conductivity 0006 % vector sigma. 0007 % 0008 % INPUT 0009 % 0010 % Agrad = the gradient part of the system matrix (see FemMatrix.m) 0011 % Kb,M and S = other blocks of the system matrix (see FemMatrix.m) 0012 % sigma = conductivity (or admittivity) vector 0013 % 0014 % OUTPUT 0015 % 0016 % A = the updated system matrix 0017 0018 % M. Vauhkonen 23.3.2000, University of Kuopio, Department of Applied Physics, 0019 % PO Box 1627, FIN-70211, Kuopio, Finland, email:Marko.Vauhkonen@uku.fi 0020 0021 rAgrad=sqrt(size(Agrad,1)); 0022 A=[reshape(Agrad*sigma,rAgrad,rAgrad)+Kb,M;M.',S]; 0023 0024 0025 0026 0027