makeregmatrix.m

MATLAB二维电阻抗断层成像算法!用于医学成像,里面包括有限元剖分正问题,及反问题的算法.并且附有网络剖分数据表!

function R=MakeRegmatrix(Element);

%MakeRegmatrix Computes a regularisation matrix which includes smoothness assumptions
% Function R=MakeRegmatrix(Element);
% computes a regularization matrix R which is a version of 2D difference. 
%
% INPUT
% Element = element structure
%
% OUTPUT
% R = regularisation matrix

msE=max(size(Element)); %剖分1的单元数
R=sparse(msE,msE);      %重构稀疏矩阵R
for ii=1:msE
 Inds=[Element(ii).Face{:,2}];
 R(ii,Inds(Inds~=0))=-1;
 R(ii,ii)=abs(sum(R(ii,:)));
end