本实验取材于中南大学《MATLAB与科学计算》
今天学习了三维曲面绘制,三维曲面绘制需要用相应的函数是,但是步骤总体包括:
- 平面网格数据的生成
- 绘制三维曲面的mesh函数和surf函数
- fmesh函数和fsurf函数
一、平面网络数据的生成
(1)利用矩阵运算生成
>> x=2:6;
>> y=(3:8)';
>> X=ones(size(y))*x;
>> Y=y*ones(size(x));
>>
(2)利用meshgrid函数生成
[X,Y]=meshgrid(x,y);
其中,参数x,y为向量,存储网格点坐标的X,Y为矩阵
>> x=2:6;
>> y=(3:8)';
>> [X,Y]=meshgrid(x,y);
例子1:绘制空间曲线。
>> x=2:6;
>> y=(3:8)';
>> [X,Y]=meshgrid(x,y);
>> Z=randn(size(X));
>> plot3(X,Y,Z);
>> grid on;
>>
图片
二、绘制三维曲面的函数
mesh(x,y,z,c)
surf(x,y,z,c)
其中,x、y是网格坐标矩阵,z是网格点上的高度矩阵,c用于指定在不同高度下的曲面颜色,
>> t=-2:0.2:2;
>> [X,Y]=meshgrid(t);
>> Z=X.*exp(-X.^2-Y.^2);
>> subplot(1,3,1);
>> mesh(X,Y,Z);
>> subplot(1,3,2);
>> surf(X,Y,Z);
>> subplot(1,3,3);
>> plot3(X,Y,Z);
>> grid on
>>
mesh函数和surf函数的其他调用格式:
mesh(z,c);
surf(z,c);
当x,y省略时,z矩阵的第2维下表当作x轴坐标,z矩阵的第1维下标当作y轴坐标,
>> t=1:5;
>> z=[0.5*t;2*t;3*t];
>> mesh(z);
其他曲面函数公式:
- 带等高线的三维网格曲面函数meshc
- 带底座的三维网格曲面函数meshz
- 具有等高线的曲面函数surfc
- 具有光照效果的曲面函数surfl
>> [x,y]=meshgrid(0:0.1:2,1:0.1:3);
>> z=(x-1).^2+(y-1).^2-1;
>> subplot(2,2,1);
>> meshc(x,y,z);title('meshc(x,y,z)')
>> subplot(2,2,2);
>> meshz(x,y,z);title('meshz(x,y,z)')
>> subplot(2,2,3);
>> surfc(x,y,z);title('surfc(x,y,z)')
>> subplot(2,2,4);
>> surfl(x,y,z);title('surfl(x,y,z)')
>>
三、标准三维曲面
(1)sphere函数
[x,y,z]=sphere(n)
(2)cylinder函数
[x,y,z]=cylinder(R,n)
例子:用cylinder函数分别绘制柱面、花瓶和圆锥面
>> subplot(1,3,1);
>> [x,y,z]=cylinder;
>> surf(x,y,z);
>> subplot(1,3,2);
>> t=linspace(0,2*pi,40);
>> [x,y,z]=cylinder(2+cos(t),30);
>> surf(x,y,z);
>> subplot(1,3,3);
>> [x,y,z]=cylinder(0:0.2:2,30);
>> surf(x,y,z);
>>
例子:用cylinder函数绘制两个相互垂直且直径相等的圆柱面的相交图形
>> [x,y,z]=cylinder(1,60);
>> z=[-1*z(2,:);z(2,:)];
>> surf(x,y,z);
>> hold on
>> surf(y,z,x);
>> axis equal
(3)peaks函数
多峰函数:
peaks函数的调用格式:
peaks(n)
peaks(V)
peaks(x,y)
peaks
fsurf函数和fmesh函数
fsurf(funx,funy,funxz,uvlims)
fmesh(funx,funy,funz,uvlims)
其中,funx,funy、funz代表定义曲面x、y、z坐标的函数,通常采用函数句柄的形式。uvlims为funx,funy和funz的自变量的取值范围,用4元向量[umin,umax,vmin,vmax]描述默认为[-5,5,-5,5].