三维分组针状图可以看作是三维分组散点图的升级，能够直观地展示各组分、各元素的位置、对比情况。
由于Matlab中未收录三维分组针状图的绘制函数，因此需要大家自行设法解决。
本文使用自制的groupedstem3小工具进行三维分组针状图的绘制，先来看一下成品效果：
特别提示：本期内容『数据+代码』已上传资源群中，加群的朋友请自行下载。有需要的朋友可以关注同名公号【阿昆的科研日常】，后台回复关键词【全家桶】查看加入方式。
1. 数据准备
此部分主要是读取原始数据，并初始化绘图参数。
% 读取数据load data.mat% 初始化绘图参数X = x;Y = y;Z = z;Clusters = clusters; % 分组信息
2. 颜色定义
作图不配色就好比做菜不放盐，总让人感觉少些味道。
但颜色搭配比较考验个人审美，需要多加尝试。
这里直接使用TheColor配色工具中的SCI权威配色库：
%% 颜色定义map = TheColor('sci',635);C1 = map([1,9,5],1:3);C2 = map([2,10,6],1:3);
3. 三维分组针状图绘制
调用‘groupedstem3’函数，绘制初始三维分组针状图。
st = groupedstem3(X,Y,Z,Clusters,C1,C2);view(-43,74)hTitle = title('Grouped Scatter chart');hXLabel = xlabel('Xaxis');hYLabel = ylabel('Yaxis');
4. 细节优化
为了插图的美观，对坐标轴细节等进行美化：
% 坐标区调整set(gca, 'Box', 'on', ... % 边框'LineWidth', 1, 'GridLineStyle', '-',... % 坐标轴线宽'XGrid', 'on', 'YGrid', 'on','ZGrid', 'on', ... % 网格'TickDir', 'out', 'TickLength', [.015 .015])% LegendhLegend = legend([st(1),st(2),st(3)], ...'Cluster1','Cluster2', 'Cluster3',...'Location', 'northwest');% 字体和字号set(gca, 'FontName', 'Arial', 'FontSize', 11)set([hLegend,hXLabel,hYLabel], 'FontSize', 11, 'FontName', 'Arial')set(hTitle, 'FontSize', 12, 'FontWeight' , 'bold')% 背景颜色set(gcf,'Color',[1 1 1])
设置完毕后，以期刊所需分辨率、格式输出图片。
%% 图片输出figW = figureWidth;figH = figureHeight;set(figureHandle,'PaperUnits',figureUnits);set(figureHandle,'PaperPosition',[0 0 figW figH]);fileout = 'test';print(figureHandle,[fileout,'.png'],'-r300','-dpng');
以上。
