MATLAB三维绘图拓展

（一）MATLAB三维绘图功能的叙述

• MATLAB 的三维绘图功能主要用于可视化复杂的三维数据，帮助用户直观地理解数据的分布、趋势和关系。

• 通过不同的绘图函数（如 plot3、mesh、surf 等），用户可以绘制各种类型的三维图形，包括曲线、曲面、散点图和柱状图。

• 通过调整视角、坐标轴范围、着色方式等，用户可以优化图形的显示效果，使其更符合分析需求。

• 三维绘图在科学研究、工程设计、数据分析等领域有广泛应用，例如流体动力学模拟、地形可视化、数据拟合等。

（二）三维绘图拓展

1.三维观察视角的设置

三维图的空间观察点的视角方位关系如图所示。

三个视图的代码指令：（1）正视图：view(0,0),即沿y轴观察；（2）侧视图：view(90,0),即沿x轴观察；（3）俯视图：view(0,90),即沿z轴观察。用户可以使用view(a,b)进行设定任意的观察角。

示例：绘制参数方程：的三维曲线及其不同视角图。

clc
clear
close all
 
t = 0:0.05:6*pi;  %设置t的采样点向量
%T = 0:0.1:4*pi;  %设置函数变量，采样数>40帧x2
x = exp(-0.2*t).*cos(t);
y = exp(-0.2*t).*sin(t);
z = t;
subplot(2,2,1);
plot3(x,y,z);
grid on;
title('三维视图');
xlabel('X axis');
ylabel('Y axis');
zlabel('Z axis');
subplot(2,2,2);
plot3(x,y,z);
grid on;
view(0,0);%正视图
title('三维正视图');
xlabel('X axis');
ylabel('Y axis');
zlabel('Z axis');
subplot(2,2,3);
plot3(x,y,z);
grid on;
view(90,0);%侧视图
title('三维侧视图');
xlabel('X axis');
ylabel('Y axis');
zlabel('Z axis');
subplot(2,2,4);
plot3(x,y,z);
grid on;
view(0,90);%俯视图
title('三维俯视图');
xlabel('X axis');
ylabel('Y axis');
zlabel('Z axis');
 

结果：

2.三维图的表面修饰

shading flat :对每一个网格块用相同的颜色着色且不显示网格

shading interp :采用插值法的光滑表面图

shading faceted :绘制表面图

曲面着色处理：colormap 色图名 或colormap (色图名），常见的色图名如下图：

示例：

clc
clear
close all
 
v = -2:0.1:2;
w = -2:0.1:2;%设定x-y平面上的取值区间
[x,y] = meshgrid(v,w);%生成矩阵数据[x,y]
z = (x.^2-2*x).*exp(-x.^2-x.*y-y.^2);%计算函数对应的曲面矩阵z
figure(1);
mesh(x,y,z);%绘制表面网格图
colormap summer;%给网格着夏色
figure(2);
mesh(x,y,z);%绘制表面网格图
colormap hot;%给网格着火色
figure(3);
surf(x,y,z);%绘制表面图
colormap copper;%给表面着古铜色
figure(4);
surf(x,y,z);%绘制表面图
colormap prism;%给表面着棱镜色

结果：

3.三维动画图的绘制

与二维动画图的绘制类似。

moviein(n):创建具有列的矩阵，可以存放n帧画面，每一列对应一帧图像。

getframe:将plot()等图形命令函数绘制的第i幅图作为第i帧画面以列的形式保存在矩阵中。

movie(m,k,fps):影片播放函数，以每秒fps帧放映矩阵m中的画面k遍。

示例：绘制曲面函数：绕z轴旋转的模型，分别绘制以曲面光照方位为变化量的三维动画图和以曲面函数的幅度为变化量的三维动画图。

clc
clear
close all

v = -8:0.25:8;
w = v;%二维平面绘图网格取样点
[x,y] = meshgrid(v,w);%生成网格
R = sqrt(x.^2+y.^2)+eps;%计算曲面函数的分子项，其中用eps定义最小误差避免R为0
Z = sin(R)./R;%计算曲面函数
M = moviein(16);%欲产生一个存储16帧图的矩阵M
for j = 1:16
    surfl(Z,[15*j,8]);%产生不同光照方位角的帧图，光源方位15*j随帧图序号变化
    axis tight;%以有数据区域紧凑坐标
    colormap copper;%曲面颜色设置为古铜色
    shading flat;%曲面同一色修饰
    M(j) = getframe;%保存帧图
end
movie(M,5,10);%以每秒10帧播放5次

clc
clear
close all

v = -8:0.25:8;
w = v;%二维平面绘图网格取样点
[x,y] = meshgrid(v,w);%生成网格
R = sqrt(x.^2+y.^2)+eps;%计算曲面函数的分子项，其中用eps定义最小误差避免R为0
Z = sin(R)./R;%计算曲面函数
M = moviein(16);%欲产生一个存储16帧图的矩阵M
figure(1);
for j = 1:16
    Z = sin(0.15*j*R)./R;
    surf(v,w,Z);%产生三维表面图
    colormap copper;%曲面颜色设置为古铜色
    shading flat;%曲面同一色修饰
    M(j) = getframe;%保存帧图
end
movie(M,5,9);%以每秒9帧播放5次
figure(2);%只有第一个窗口的三维动画播放完毕，才能播放第二个三维动画
for j = 1:16
    Z = sin(0.15*j*R)./R;
    surf(v,w,Z);%产生三维表面图
    colormap hot;%曲面颜色设置为火色
    shading flat;%曲面同一色修饰
    M(j) = getframe;%保存帧图
end
movie(M,5,9);%以每秒9帧播放5次

结果：