get_multi_meas.sci

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

function [voltage,ind,df]=get_multi_meas(protocol,elec,V,I,vtx,no_pl)
voltage=[];ind=[];df=[];
//function [voltage,ind,df] = get_multi_meas(protocol,elec,V,I,vtx,no_pl);
// 
//The function can be used in the occasions where plane current patterns 
//(adjacent or polar) are adopted for systems with more planes, i.e. set by
//set_multi_currents function. Only non-current carrying electrodes are 
//involved in the measurements.
// 
// 
// 
//protocol= The selected protocol '{op}' or '{ad}'
//elec    = The electrodes
//no_pl   = The number of planes
//I       = The RHS vectors, i.e., the current patterns padded with zeroes 
//          for the forward calculations
//V       = The forward solution
//voltage = The array of measureemnts according to the selected protocol
//ind     = The index of the measurements
 
 
if size(V,2)~=size(I,2) then
  error('Unmatched pattens');
end
 
 
no_el = size(elec,1);
el_pp = no_el/no_pl;
 
 
//! mtlb(end) can be replaced by end() or end whether end is an m-file or not
I = I(size(vtx,1)+1:mtlb(end),:);
//Lower chunk
 
//! mtlb(end) can be replaced by end() or end whether end is an m-file or not
Vm = V(size(vtx,1)+1:mtlb(end),:);
//Lower chunk 
 
 
voltage = [];
ind = [];
df = [];
 
d = size(I,2);
//Injections
 
for i = 1:d
  //for each injection
   
  cn = 0;
   
  L = [];
  fst = 0;
  kk = 0;
   
  for ej = 1:el_pp:no_el
     
    kk = ej-1;
     
    for ei = 1:el_pp-1
      if ei==1 then
        fst = ej;
      end
      L = [L;[kk+ei,kk+ei+1]];
    end
    L = [L;[fst,fst+16-1]];
  end
   
  LMP = [];
  //Exclude measurements engaging electrodes from different electrodes
   
  for j = 1:size(L,1)
     
    if ceil(L(j,1)/el_pp)==ceil(L(j,2)/el_pp) then
       
      LMP = [LMP;L(j,:)];
    end
  end
   
  el_in = matrix(find(I(:,i)~=0),1,-1);
  // The current carring electrodes
   
  LM = [];
   
  for k = 1:size(LMP,1)
     
    //!! Unknown function ismember ,the original calling sequence is used
    if ismember(LMP(k,:),el_in)==[0,0] then
      LM = [LM;LMP(k,:)];
      //Modified list current carring electrodes removed. 
    end
  end
   
  lf = size(LM,1);
  dd = lf/no_pl;
   
  LN = [];
   
  k = 0;
   
  for in = 1:dd:lf
    //   1    13    25    37
     
    k = k+1;
     
    th_LM = LM(in:in+dd-1,:);
    //size(th_LM) = (12 2)
     
    if protocol=='{op}' then
      for j = 1:size(th_LM,1)
         
        //!! Unknown function intersect ,the original calling sequence is used
        if th_LM(j,1)==(intersect(matrix(find(I(:,i)>0),1,-1),in:k*el_pp)+1) then
           
          //! mtlb(end) can be replaced by end() or end whether end is an m-file or not
          LN = [th_LM(j:mtlb(end),:);th_LM(1:j-1,:)];
          break
           ;
           
        end
      end
    end
     
    if protocol=='{ad}' then
      for j = 1:size(th_LM,1)
         
        //!! Unknown function intersect ,the original calling sequence is used
        if th_LM(j,1)==(intersect(matrix(find(I(:,i)<0),1,-1),in:k*el_pp)+1) then
           
          //! mtlb(end) can be replaced by end() or end whether end is an m-file or not
          LN = [th_LM(j:mtlb(end),:);th_LM(1:j-1,:)];
          break
           ;
           
        end
      end
    end
     
    for u = 1:size(LN,1)
      voltage = [voltage;Vm(LN(u,1),i)-Vm(LN(u,2),i)];
    end
     
     
    ind = [ind;LN];
    cn = cn+1;
     
  end
  //for in
   
  df = [df;cn];
end
//for injections
 
 
 
 
 
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
// 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
//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%