makenode2nd.m

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

function [Node]=MakeNode2nd(Element,Nodelist,g);

%MakeNode2nd Computes the Node data structure for a 2D mesh having quadratic triangles
% Function [Node]=MakeNode2nd(Element,g);
% computes the Node data for MeshData.
% Node is a structure including all the nodal coordinates and
% for each node there is information to which nodes (NodeConnection) 
% and elements (ElementConnection) the node is
% connected.  
%
% INPUT
%
% Element = Element structure, see MakeElement.m
% Nodelist = list of all the faces, see MakeElement.m
% g = coordinates of the nodes
%
% OUTPUT
%
% Node = node structure


[rg,cg]=size(g);
msE=max(size(Element));
for ii=1:rg
 ElementConnection=[];
 Node(ii).Coordinate=[g(ii,:)];
  for jj=1:msE
   if find(Element(jj).Topology==ii)
    ElementConnection=[ElementConnection,jj];
   end
  end
  Node(ii).ElementConnection=ElementConnection;
 Nc=[];
 [I,J]=find(Nodelist==ii);
 apu1=Nodelist(I(find(J==1)),:); %Corner node
 apu2=Nodelist(I(find(J==2)),:); %Middle node
 apu3=Nodelist(I(find(J==3)),:); %Corner Node
 apu=[apu1(:,2);apu2(:,1);apu2(:,3);apu3(:,2)];
 Nc=unique(apu);
  Node(ii).NodeConnection=Nc;
end