matlab中添加网格命令,Matlab网格划分

之前转载了一篇博客http://blog.sina.com.cn/s/blog_6163bdeb0102dvay.html，讲Matlab网格划分程序Distmesh，看了看程序，感觉程序写得有很多值得学的地方，所以又自己重新看了一看，加了一些注释，最后再总结一下学到的东西吧。

源代码的地址已经改变，如下http://people.sc.fsu.edu/~jburkardt/m_src/distmesh/distmesh.html。程序给出了20个划分的例子，文件名为p1_demo.m~p20_demo.m，直接运行程序可以看到划分的动画效果，每个例子基本都是设置一些参数，然后调用distmesh_2d函数进行网格划分，最后得到划分的节点p和各三角形t，最后将划分的三角形绘制出来。划分结果如下

p1_demo p14_demo

p5_demo p18_demo

进一步加注释的程序(需要学习的地方用颜色标记)：

function [ p, t ] = distmesh_2d ( fd, fh, h0, box,

iteration_max, pfix, varargin )

%% DISTMESH_2D is a 2D mesh generator using distance

functions.

% Example:

% Uniform

Mesh on Unit Circle:

% fd = inline('sqrt(sum(p.^2,2))-1','p');

% [p,t] = distmesh_2d(fd,@huniform,0.2,[-1,-1;1,1],100,[]);

% Rectangle

with circular hole, refined at circle boundary:

% fd =

inline('ddiff(drectangle(p,-1,1,-1,1),dcircle(p,0,0,0.5))','p');