matlab isosurface函数绘制隐函数曲面

转自http://blog.sina.com.cn/s/blog_6163bdeb0100pi9b.html

最近使用matlab，看到一些小技巧，总结至此

1.如xlabel('\fontname{times New Toman}\fontsize{40}你好')，可以更快捷、更方便的设置显示内容的属性，当然也可以直接在figure中设置

2.isosurface函数可以画隐函数的等值曲线或曲面

如画笛卡尔心形函数，其方程为（x^2 + (9/4)y^2 + z^2 - 1)^3 - x^2z^3 - (9/80)y^2z^3 = 0，可参考网址http://luo.bo/1924/

由于是隐函数，不好用一般的方法画出其图形来，这就可以用isosurface来实现，代码如下

figure
[x,y,z]=meshgrid(linspace(-3,3));                            %做出网格meshgrid                                    
p=(x.^2+(9/4)*y.^2+z.^2-1).^3-x.^2.*z.^3-(9/80)*y.^2.*z.^3;  %实现结果的表达
isosurface(x,y,z,p,0,x);
axis equal;
axis off;
view(0,45);          %视角的控制
colormap([1 0 0]);   %绘图颜色红色
brighten(0.5);       %增亮
camlight right;      %光源位置
lighting phong;      %光照模式

代码参考http://blog.sina.com.cn/s/blog_4b013fb10100llxh.html

 

另外这个函数还有一个好处是，可以得到等值面的顶点和面，然后可以直接调用patch画出来

fv = isosurface(x,y,z,p,0);
patch(fv)
这一点很好，可以将得到的顶点和面保存到文件，用opengl画出来，可以更好的渲染图形，这点已经实现。

 

下面简单翻译一下help里的解释吧：

isosurface函数用于由体积数据中提取等值数据。
使用方法如下

fv = isosurface(X,Y,Z,V,isovalue)
fv = isosurface(V,isovalue)
fvc = isosurface(...,colors)

[f,v] = isosurface(...)
[f,v,c] = isosurface(...)
isosurface(...)

参数意义为，从由X,Y,Z构造的体积V数据中提取由isovalue指定的等值数据，返回结果fv为一个结构体，包含了等值面的顶点和面（顶点的次序），这些参数可以直接传给patch命令画出图形来。

另外，可以在最后加上colors，返回值中也会得到每个面的颜色插值，这样能够使你用不同于计算等值面的数据控制等值面的颜色映射。

如果不设置返回值，就会自动在当前坐标系里用计算得到的顶点和面创建三维patch对象。