背景介绍
Matlab提供了一系列绘图函数,常见的包括绘制2D曲线的plot函数、绘制2D隐函数曲线的ezplot函数、绘制3D曲面的mesh和surf函数、绘制3D显函数曲面的ezmesh和ezsurf函数。值得注意的是,ez系列的绘图函数里只有ezplot是绘制隐函数曲线的,ezmesh和ezsurf都是画显函数曲面的(不要被ez的名字误解了)。遗憾的是,matlab里并没有提供直接绘制3D隐函数曲面的函数。本帖的目的就是归纳总结几种方便易用的绘制隐函数曲面的办法。
问题描述
如何绘制 3 元方程 f(x, y,z) = 0 确立的隐函数曲面 z = g(x,y) ?其中,方程 f(x, y,z) = 0 无法求解 z 关于 x 、 y 的表达式,即 g(x, y) 的显式表达式无法获取。
准备工作——基础函数介绍
为了解决上述问题,我们需要先 对几个重要的图形函数 isosurface 、 patch 、 isonormals 取得初步的了解,如果您已经对这三个函数很熟悉,可以直接跳过这一步。
l. isosurface 等值面函数
调用格式: fv = isosurface(X,Y,Z,V,isovalue)
作用:返回某个等值面(由 isovalue 指定)的表面( faces )和顶点( vertices )数据,存放在结构体 fv 中( fv 由 vertices 、 faces 两个域构成)。如果是画隐函数 v = f(x,y,z) = 0 的三维图形,那么等值面的数值为 isovalue = 0 。
2. patch函数
调用格式: patch(X,Y,C) 以平面坐标 (X, Y) 为顶点,构造平面多边形,