bound2.m

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

function int=bound2(g);

%bound2 Computes the boundary integral of the product of two linear basis function in 2D FEM
% Function int=bound2(g) calculates the boundary integral
% of the product of two basis function from g(1,:) to g(2,:).
%
% INPUT
%
% g = integration endpoints
%
% OUTPUT
%
% int = value of the integral



w=[1/2,1/2];
ip=[1/2-1/6*sqrt(3),1/2+1/6*sqrt(3)];
dJt=sqrt((g(2,1)-g(1,1))^2+(g(2,2)-g(1,2))^2); 
int=0;
 for ii=1:2
  S=[1-ip(ii);ip(ii)];
  int=int+w(ii)*S*S';
 end
int=int*dJt;