updatefemmatrix.m

来自「MATLAB二维电阻抗断层成像算法!用于医学成像,里面包括有限元剖分正问题,及反」· M 代码 · 共 22 行

function [A]=UpdateFemMatrix(Agrad,Kb,M,S,sigma);

%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 

rAgrad=sqrt(size(Agrad,1));
A=[reshape(Agrad*sigma,rAgrad,rAgrad)+Kb,M;M.',S];