冲击图是一种特殊的堆叠柱状图。

与堆叠柱状图相比，冲击图添加了相邻柱子中相同组分之间的连线，可以更加清晰地表达各组分占比情况。

由于Matlab中未收录冲击图的绘制函数，因此需要大家自行解决。

本文使用自制的Fbarstacked小工具进行冲击图的绘制，先来看一下成品效果：

1. 数据准备

此部分主要是读取原始数据并初始化绘图参数。

% 读取数据

load data.mat

% 初始化绘图参数

y = A(:,1:6);

2. 颜色定义

作图不配色就好比做菜不放盐，总让人感觉少些味道。

但颜色搭配比较考验个人审美，需要多加尝试。

这里直接使用TheColor配色工具中的SCI权威配色库：

%% 颜色定义

map = TheColor('sci',2064,'map',10);

map = flipud(map);

C = map([1 2 3 6],1:3);

3. 冲击图绘制

调用‘Fbarstacked’命令，绘制初始冲击图。

GO = Fbarstacked(y,0.4,C,0.7);

hTitle = title('Filled stacked bar chart');

hXLabel = xlabel('Samples');

hYLabel = ylabel('RMSE (m)');

4. 细节优化

为了插图的美观，将初始冲击图坐标区基本参数进行调整：

% 坐标区基本属性调整

set(gca, 'Box', 'off', ... % 边框

'LineWidth', 1, ... % 线宽

'XGrid', 'off', 'YGrid', 'on', ... % 网格

'TickDir', 'out', 'TickLength', [.01 .01], ... % 刻度

'XMinorTick', 'off', 'YMinorTick', 'off', ... % 小刻度

'XColor', [.1 .1 .1], 'YColor', [.1 .1 .1]) % 坐标轴颜色

% 坐标轴刻度调整

set(gca, 'YTick', 0:0.3:1.5,...

'Ylim' , [0 1.5], ...

'Xlim' , [0.5 6.5], ...

'XTick', 1:7,...

'Xticklabel',{1:7},...

'Yticklabel',{0:0.3:1.5})

% legend

hLegend = legend([GO(1),GO(2),GO(3),GO(4)],...

{'A', 'B', 'C', 'D'});

% 字体和字号

set(gca, 'FontName', 'Arial', 'FontSize', 10)

set([hLegend, hXLabel, hYLabel], 'FontName', 'Arial', 'FontSize', 11)

set(hTitle, 'FontSize', 12, 'FontWeight' , 'bold')

% 背景颜色

set(gcf,'Color',[1 1 1])

然后，以期刊所需分辨率、格式输出图片。

%% 图片输出

exportgraphics(gca,'test.png','Resolution',300)

以上。

转自：“阿昆的科研日常”微信公众号

如有侵权，请联系本站删除！