Matlab网格划分.doc

《Matlab网格划分.doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、Matlab网格划分http://www-math.mit.edu/~persson/mesh/（已删除）Distmesh是一个matlab语言写的网格划分软件。源文件可以从上面的网址获取。这里按行讲解各个算例。p01_demo:概算例是一个单位圆（半径为1）的网格划分，划分后的网格为：以下逐行讲解该算例：functionp01_demo(iteration_max,h)%  Parameters:%  Input,integerITERATION_MAX,themaximumnumberofiterationsthatDISTMESH%    shouldtake.  (The

2、programmighttakefeweriterationsifitdetectsconvergence.)%  Input,realh,themeshspacingparameter.%这里有两个输入参数，一个是ITERATION_MAX,迭代的最大次数。%另一个是h,网格划分的大小。0

3、x^2+y^2)-1;fh,scalededgelengthfunctionh(x,y).也就是网格大小的函数。h0也就是h,网格的大小realBOX(2,2),theboundingbox[xmin,ymin;xmax,ymax].最大外围矩形范围。本例中为[0,0;1,1]ITERATION_MAX,themaximumnumberofiterations.realPFIX(NFIX,2),thefixednodepositions.网格中需要固定的点坐标，也就是一定需要出现在网格中的点。输出参数：realP(N,2),thenodepositions.  网格点的x,y坐标

4、integerT(NT,3),thetriangleindices.输出网格任一一个三角形的三个顶点。第一步：  [x,y]=meshgrid(box(1,1):h0           :box(2,1),...                        box(1,2):h0*sqrt(3)/2:box(2,2));根据h0，网格的大小，先把能涵盖欲划分区域的最大矩形划分为结构网格。然后把偶数行的点整体向右平移半格， x(2:2:end,:)=x(2:2:end,:)+h0/2;效果如下：第二步：根据fd的函数定义，移除边界外的点。 p=p(feval_r(fd,p,va

5、rargin{:})<=geps,:);varagin为fd，fh的附加参数，这里为空。geps=0.001*h0;也就是保留了到边界的距离以外0.001*h0以内的点。 根据网格密度函数fh，每个点上产生一个0-1随机数，判断是否小于r0/max(r0)大于的话，改点被删除。 p=[pfix;p(rand(size(p,1),1)

6、veafew%  duplicates.  Getridofanyyouspot.%  p=unique(p,'rows');  N=size(p,1);这个时候产生的网格如下：第三步：迭代  pold=inf;%第一次迭代前设置旧的点的坐标为无穷while(iteration

7、节点      t=delaunayn(p);  %利用delauny算法，生成三角形网格      triangulation_count=triangulation_count+1;      pmid=(p(t(:,1),:)+p(t(:,2),:)+p(t(:,3),:))/3;  %计算三角形的重心。      t=t(feval_r(fd,pmid,varargin{:})<=-geps,:);%移除重心在边界外部的三角形%    4.Describeeachbarb