fem_master_full.sci

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

function [E,pp]=fem_master_full(vtx,%simp,mat,gnd_ind,elec,zc,sym)
E=[];pp=[];
//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}' then
  pp = 1:size(E,1);
end
if sym=='{y}' then
   
  //!! Unknown function symmmd ,the original calling sequence is used
  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
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%