fem_master_full.m

用来实现三维阻抗及光学断层成像重建的matlab程序

function [E,pp] = fem_master_full(vtx,simp,mat,gnd_ind,elec,zc,sym);
%function [E,pp] = fem_master_full(vtx,simp,mat,gnd_ind,elec,zc,sym);
%
%Builds up the system matrix based on the complete electrode model. E is not 
%yet permuted. To permute E -> E(pp,pp) as in fwd_solver.
%
%
%
%E       = The full rank system matrix based on the 3D complete electrode model
%vtx     = The vertices matrix
%simp    = The simplices matrix
%mat     = The conductivity vector
%gnd_ind = The index of the ground node
%elec    = The bounary electrodes matrix
%zc      = The contact impedance vector, satisfying size(elec,1) = length(zc)
%sym     = Column permutation of E, either '{y}' to opt or '{n}' to avoid.       

%if isreal(mat)==0   
%  mat = conj(mat);
%end
     
   [Ef] = bld_master_full(vtx,simp,mat,elec,zc); 
   
   [E] = ref_master(Ef,gnd_ind);  
   
 if sym == '{n}'
    pp = 1:size(E,1);
end
 if sym == '{y}'
    pp = symmmd(E);
end
   

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This is part of the EIDORS suite.
% Copyright (c) N. Polydorides 2001
% Copying permitted under terms of GNU GPL
% See enclosed file gpl.html for details.
% EIDORS 3D version 1.0
% MATLAB version 5.3 R11
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%