f_hsplcigrid.m

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

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

% plots a given 2D FEM mesh.
%
% INPUT
% g = node coordinate matrix
% H = topology 
% E = elements under the electrodes (can also be empty)
% No = labe the No.of Element and node if No=='Y'
%By Huangsong CQU at 2003.2.16

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),'r');
    [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),'r');
  end
  if No=='Y'
      text(Ho(1),Ho(2),num2str(ii),'FontSize',6.5,'Color','black');  
      text(Hii(1,1),Hii(1,2),num2str(H(ii,1)),'FontSize',7,'Color','b');
      text(Hii(2,1),Hii(2,2),num2str(H(ii,2)),'FontSize',7,'Color','b');
      text(Hii(3,1),Hii(3,2),num2str(H(ii,3)),'FontSize',7,'Color','b');
  else
  %    disp('Not Labed')
  end
end