plcigrid.m

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

function []=plcigrid(g,H,E,co)

% plots a given 2D FEM mesh.
%
% INPUT
% g = node coordinate matrix
% H = topology 
% E = elements under the electrodes (can also be empty)
% co = colour of the mesh

nH=max(size(H));

axis('xy'),axis('square')
hold on

for ii=1:nH                        %对每个单元H
  Hii=g(H(ii,:),:);                %单元3个节点坐标
    Ho=sum(Hii)/3;                   %单元号标注位置  
  Hii=[Hii;Hii(1,:)];              %单元3个节点+第1节点坐标，4点3边
  if ~isempty(E)
    hHii=plot(Hii(:,1),Hii(:,2),co);
    [pp,cc]=find(ii==E);
    if ~isempty(find(ii==E)),
     hHii=fill(Hii(:,1),Hii(:,2),pp);
    end
  else
    hHii=plot(Hii(:,1),Hii(:,2),co);
  end
end