removevolt.m

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

function Vh=RemoveVolt(Uel,L);

%RemoveVolt Removes all the voltages measured on the current carrying electrodes
% Function Vh=RemoveVolt(Uel,L);
% removes all the voltages measured on the current carrying electrodes.
% Change the Jacobian accordingly, see RemoveJacob.m.
%
% INPUT
% Uel = voltages on the electrodes
% L = number of electrodes
% 
% OUTPUT
% Vh = new voltages 

  Uel=Uel(:);
  [rUel,cUel]=size(Uel);
  Uh=reshape(Uel,L,rUel/L);
  Vadj=Uh;
  kk=1;
   for ii=1:L
    Vadj(:,ii)=[Vadj(kk:L,ii);Vadj(1:kk-1,ii)];
    kk=kk+1;
   end
  Vadj([1 2 L],:)=[];
  Vh=Vadj;