目录
一、功能说明
二、plot3函数说明
三、代码实现
四、代码优化
一、功能说明
近期在实验中采集了多组二维样本数据,每组数据都是随时间变化的曲线(如下图1.1),现在需要将所有样本都放在1张三维图中显示出来,从而直观的对比各组数据的差异。最终使用MATLAB的plot3函数实现了此功能,现将此方法记录下来。
图1.1 单组二维样本曲线 图1.2 最终效果:多组二维曲线堆叠成三维图
二、plot3函数说明
MATLAB官方使用说明:
plot3(X,Y,Z)
绘制三维空间中的坐标。
要绘制由线段连接的一组坐标,请将 X
、Y
、Z
指定为相同长度的向量。
要在同一组坐标轴上绘制多组坐标,请将 X
、Y
或 Z
中的至少一个指定为矩阵,其他指定为向量。
在绘制图象时我是这样理解:(x,y,z)对应着空间中的每个点,在绘图时,先将Y轴全部视为某一常数值,然后再将X轴视为自变量,Z轴视为因变量,这样的话,plot3(x,常数值,z)实际上与plot(x,y)函数的使用方式是一样的了。下面这张图是将Y向量全部赋值常数3(即第3组),X向量取值1:100,Z向量取X对应的样本值。
图2 三维图像(以取第3组为例)Q:为什么Y轴的值是个向量而非一个常数值呢?
这是因为第3组数据平面上的点不止一个,确切的说是由length(x)这么多个,也就是有length(x)个这样的点(x,常数值,z),所以在上面的情况中,Y轴的值是个长度为length(x),数值全为3的向量。
图3 标注出Y轴坐标值:全为3如果实在不能理解上述说明,可以先尝试将下面的代码运行起来,一边调整参数,一边理解使用方法。
三、代码实现
有了前面的基础,现在可以这样编写代码:
首先,生成一个样本集Data用于后续绘图展示(移植代码时,要根据自身数据大小来定义Data尺寸,并为其赋值),这个样本集共有10组样本,每组样本有100个数据,分别存储在Data的每一行中,所以Data初始值暂时定义为zeros(20,100)。[!行与列不要弄混!]
clc;clf;close;clear;%% 产生数据源 Data% 待显示数据的大小:20组数据样本,每组样本100个数据Data = zeros(20,100); % 为待显示数据(Data)赋值,自定义曲线t = 1:100;for i = 1 :10 Noise = randn(1,length(t))*0.1; %生成随机噪声 Data(i,:) = [ones(1,20)*2 ones(1,60)*4 ones(1,20)*2] + Noise; end
然后对Data中的数据进行绘图(核心代码):
x为1:100的向量;y为常数值向量Group_num*ones(size(x)),Group_num的值对应组别的序号;z是x对应的因变量,每次取Data矩阵的第Group_num行即可。
%% 绘制三维图像figure(1);for Group_Num = 1:10 x = 1:100; % 数据长度为100,将全部显示出来 y = Group_Num * ones(size(x)); % 共10组数据(组别号:1:10) z = Data(Group_Num, 1:100); % 每组数据的曲线值 plot3(x,y,z); % 绘图 hold on; % 保持endgrid on; title(['多组样本值变化曲线']);xlabel("时间");ylabel("组别");zlabel("曲线值");
运行结果如下所示:
图4 三维图像(多组二维图像叠加)四、代码优化
实验中有时只需要观察部分数据,比如只看第4组到第8组的数据,或是只看曲线中间部分的数据,所以根据需求可以优化一下代码。
限制X轴范围(设置起止数据点)限制Y轴范围(选取样本集组别)限制Z轴范围(设置数据值范围)优化后的代码全部附在下面(可以直接编译运行):
clc;clf;close;clear;%% 产生数据源 Data% 待显示数据的大小:20组数据样本,每组样本100个数据Data = zeros(20,100); % 为待显示数据(Data)赋值,自定义曲线t = 1:100;for i = 1 :10 Noise = randn(1,length(t))*0.1; %生成随机噪声 Data(i,:) = [ones(1,20)*2 ones(1,60)*4 ones(1,20)*2] + Noise; end%% 绘图代码优化%---3D图中显示数据的起始点---Range = [20,90]; %---导入数据量的组别范围---Number = [3,9] ;%---三维样本值变化曲线---%figure(2);for Group_Num = Number(1) : Number(2) x = Range(1) : Range(2); y = Group_Num*ones(size(x)); z = Data(Group_Num, Range(1):Range(2) ); plot3(x,y,z); hold on;endgrid on; %---3D图中数据值显示范围---zlim([1,5]); title(['多组样本值变化曲线']);xlabel("时间");ylabel("组别");zlabel("样本值");
下图为选取第3组至第9组样本,每组只显示第30到90数据点的结果。
图5 X,Y,Z限制后图像
如果想固定观察角度,可以在代码结尾处添加view(a,b)函数进行调整:
view(az,el)
为当前坐标区设置照相机视线的方位角和仰角。
将两个参数都设置为0,则是相当于从正前方观察,运行结果如下图所示。
图6 view(0,0)视角观察~------------~
~---Over---~
~------------~