removevolt.m

来自「实现对电磁层析模型的建立」· M 代码 · 共 38 行

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 


% M. Vauhkonen 6.9.1999, 
% University of Kuopio, Department of Applied Physics, PO Box 1627,
% FIN-70211 Kuopio, Finland, email: [e-mail]

  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;