Matlab网格划分

最新推荐文章于 2024-08-09 09:33:25 发布

weixin_30783913

最新推荐文章于 2024-08-09 09:33:25 发布

阅读量5.4k

点赞数 1

文章标签： matlab c/c++

原文链接：http://www.cnblogs.com/haoyul/p/5909427.html

版权

本文介绍了Matlab中的网格划分程序Distmesh，详细解析了源代码，并通过示例展示了其在不同场景下的应用。程序通过设置参数，如距离函数、网格大小和最大迭代次数，进行二维网格划分。此外，文章还探讨了Matlab中用于三角形划分的delaunayn函数和图形输出trimesh函数的使用，以及如何利用数学函数的思想来提高编程效率。

摘要由CSDN通过智能技术生成

之前转载了一篇博客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:
%