removejacob.m

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

function Jh=RemoveJacob(J,L);

%RemoveJacob Removes the derivatives corresponding to the current  carrying electrodes
% Function Jh=RemoveJacob(J,L);
% removes all the derivatives that correspond to the
% current carrying electrodes.
% Change the voltages accordingly, see RemoveVolt.m.
%
% INPUT
% J = Jacobian
% L = number of electrodes
% OUTPUT
% Jh = new Jacobian

[rJ,cJ]=size(J);
J=reshape(J,L,rJ/L*cJ);
kk=1;p=1;
for jj=1:cJ
 for ii=1:L
  J(:,p)=[J(kk:L,p);J(1:kk-1,p)];
 kk=kk+1;p=p+1;
end
kk=1;
end
J([1 2 L],:)=[];
Jh=reshape(J,(L-3)*L,cJ);