【高级绘图】画一个SCI论文里的柱状图(文末送书)

硬件电子工程师

点击上方蓝字和“好玩的matlab”一起快乐的玩耍吧
# K# Z: t1 M1 H3 q8 @6 R9 ^

dzbh43o50sz6408890423.jpg

好玩的matlab
! ?) w# Q4 B; \; k" c$ W, W带你解锁不一样的matlab新玩法
. T" ~/ h2 N7 i8 V* B
! O5 {# G3 F7 M' N喜欢此推文的小伙伴们记得点赞+关注+分享！
相信大家都绘制过柱状图，可是大家绘制的柱状图真的好看嘛~现在来看看SCI论文里面的柱状图是什么样的。以下是来自science2020[1]论文里的配图，接下来手把手带大家用MATLAB来复现一下这两个柱状图。

7 g5 ]9 B9 E$ \; R

lozlvuflv0s6408890523.jpg

# \$ H! e% f5 V



9 E: N: ?7 P# X( Y: M4 a6 w; E

pis3pzas24t6408890623.gif

效果

左边是复现的图，右边是原图。虽然不是一模一样，但大部分元素还是复现出来了。


) j' ~% e6 [$ G2 W# g

& C9 n1 |% H  _ wfevc1focxw6408890723.png	& T; j4 u6 J! o- v' m- m tjbsdb2x0vg6408890823.png
hzu1qungcmo6408890923.png	8 s4 U8 [4 X3 @0 D9 ~% E 4nsca450ygt6408891023.png ) C3 B' _5 A& l! D' q

. V' N& Z. _: M; n3 j
; s% w/ ~9 d3 h

- {  y) z5 ?" H. e0 `" m' M
绘制要素

, x- A' S& H% |" G# u

第一个绘图要素：颜色
0 X3 C) n& K& v0 J3 G" I4 L看到这张图，首先映入眼帘的是啥呢？当然是颜色了！颜色配的好，作图就成功了一大半。首先要获得图片里的颜色~

第二个绘图要素：数据
3 \# q, C) ?7 x' ]. L作图数据是必不可少的，但你要画图的画，你肯定已经有数据了。当然小编没有这篇SCI论文里的数据，所以就仿造了一点数据。

第三个绘图要素：框架
9 I: J* e) V( y2 ^4 S* X! \9 D3 P         所谓画图的框架，也可以称为画图的函数吧，明显可以看出这不是          一张简单的柱状图，它还叠加了误差图，散点图，需要用到三个                     MATLAB里面的绘图函数，bar,errorbar,scatter。知道了这三个要素咱            就可以开始画了。
% E  c  c7 M0 m! d7 a$ o
# i4 e5 i& a/ B- H; a

vyg4mxsvjdg6408891123.gif

11lw5ycyeay6408891223.png

1、取色
% f( I3 p8 U) [

jk31xezf3yf6408891323.png

还是使用和【图片配色】那篇推送一样的方法，还没有看过的可以去看一眼，要想复现出一样的图，肯定先得拿出这张图的颜色。

n02mv0l14pg6408891423.png

: |* L7 q5 E+ _* Y! ], S1 l那么问题来了，这张图有几种颜色呢？
4 X( O0 [$ l8 o- d4 P当然是7种颜色了，注意散点图和柱状图边缘颜色是不一样的奥。用以下代码取色图片就是上面的截图，运行之后在图片不同颜色点处打点，完了之后按enter,颜色就自动取出来了。

% 公众号：好玩的MATLABhua=imread('1.png');figure(1)imshow(hua);[c1,c2]=getpts; %选择颜色点c1=ceil(c1); %将选点的位置化为整数c2=ceil(c2);for i=1:length(c1)    color(i,:)=hua(c2(i),c1(i),:);endcolor=double(color);color=color/255;
以下就是取出的七种颜色了，建议把它记录下来，省得下一次麻烦。记录下来的颜色变量如下

color=[0.717647058823529,0.682352941176471,0.741176470588235;...    0.807843137254902,0.792156862745098,0.890196078431373;0.580392156862745,...    0.525490196078431,0.729411764705882;0.392156862745098,0.345098039215686,...    0.470588235294118;0.560784313725490,0.533333333333333,0.741176470588235;...    0.372549019607843,0.282352941176471,0.600000000000000;0.0156862745098039,...    0.0196078431372549,0.0156862745098039];


4 H0 O4 T  U7 E! h8 M* A! r
# @0 H, Q# J  H' x# U2、数据
, X, }0 z+ V  @" ]6 e8 f- [9 K
1 e$ O3 m) {, E+ o5 }" s1 s1 Q
这张看起来小小的柱状图，包含的数据信息还真不少呢。以下是根据原图自拟的数据，凑合着用。

% 公众号：好玩的MATLABy=[28,34,18];%柱状图纵坐标x=[1,2,3];     %柱状图横坐标%散点图的横纵坐标s1{1,1}=[0.9,0.95,0.95,1,1,1.05,1.1,1.1];s1{1,2}=[1.85,1.9,1.95,2,2,2,2,2.2];s1{1,3}=[2.7,2.7,2.95,3,3,3,3,3,3.2,3.2,3.2];s1{2,1}=[18,24,36,30,42,20,29,38];s1{2,2}=[26,22,24,30,38,42,48,50];s1{2,3}=[9,19,8,12,18,24,30,40,10,17,19];low=[4,3,3]; %误差图的上限high=[3,2,1]; %误差图的下限
9 [/ D( S( [; Z+ R

0 y1 I, Z' B6 N; J3 v/ H
3、画图
2 f( s" r9 m! d: l0 ]
- k& _, |. Z7 d( k: p. |
  _# h( Z' T2 g7 ^" ]1 初始化显示位置
0 g! k% `0 ~) V
看起来不起眼的一点实则非常重要，可以根据论文的大小来设置图片的大小，很大程度上可避免后期拖拉图片造成的字体变形，不清晰等问题。
% S$ [/ W/ E: d2 `4 @" V

%%初始化显示位置 [x,y,dx,dy] 从 x,y坐标开始，dx,dy为沿升的长度figure('Units', 'pixels', ...    'Position', [100 100 360 375]);
2 循环画图和配色

( e! Q4 `' \# V. w' q   设置柱状图的宽度让柱状图间隔大一点，注意柱状图的填充颜色和边缘颜色是需要区分开来的。其他改变画图的标记都在下面了~。

width=0.6; %柱状图宽度for  i=1:length(y)   %bar(i,y(i),width) width 柱状图宽占比  ，FaceColor 柱状图填充颜色 ，EdgeColor 柱状图边缘颜色   % LineWidth  边缘线宽    set(bar(i,y(i),width),'FaceColor',color(i,:),'EdgeColor',color(i+3,:),'LineWidth',2)   hold on      %散点图 大小，填充颜色   scatter(s1{1,i},s1{2,i},60,'filled','MarkerFaceColor',color(i+3,:))      %在柱状图 x,y 基础上 绘制误差 ,low为下误差，high为上误差，LineStyle 误差图样式，'Color' 误差图颜色     % 'LineWidth', 线宽,'CapSize',误差标注大小   errorbar(i, y(i), low(i), high(i), 'LineStyle', 'none', 'Color', color(i+3,:), 'LineWidth', 1.5,'CapSize',18);end
3 坐标和字体的设置
, [# ]" U  i" Z! ^4 ?
' `+ R/ R* J* ?6 V, h1 t& H" h坐标和字体的适当配置可以大大增加作图的颜值

%绘制虚线plot([0,3.5],[50,50],'--','Color',color(7,:),'LineWidth',2)ylabel('Time on Warm Floor(%)')ax = gca;%x轴标签ax.XTick = [1 2 3];set(gca,'XTickLabel',{['\color[rgb]{',num2str(color(4,:)) ,'}','Pre-Injection'],['\color[rgb]{',num2str(color(5,:)) ,'}','1 hr CFA'] , ['\color[rgb]{',num2str(color(6,:)) ,'}','1wk CFA']},"LineWidth",2);%y轴标签ax.YTick= 0:20:80;%y轴范围ax.YLim=[0,80];%使得只显示x,y轴两根线box off%坐标轴字体和大小set(gca,"FontName","Times New Roman","FontSize",12,"LineWidth",2)title("TPT Timecourse","FontName","Times New Roman","FontSize",12,"LineWidth",2);
运行代码，就可以得到这样一张图，还是还原了大部分原图的。

ygyhyfqzdzo6408891523.png

! T. w1 {2 ~1 f

$ K8 T6 L: |4 d! M: A
- s8 G  B1 A+ J8 @) J( z3 E
两幅图完整源码
; |- S+ _$ S/ C( V9 j5 \, ^
) J5 Z, b. n/ [# z  G0 R$ O* X7 f
8 ?: s. T2 i; ?1 m

clc;color=[0.717647058823529,0.682352941176471,0.741176470588235;...    0.807843137254902,0.792156862745098,0.890196078431373;0.580392156862745,...    0.525490196078431,0.729411764705882;0.392156862745098,0.345098039215686,...    0.470588235294118;0.560784313725490,0.533333333333333,0.741176470588235;...    0.372549019607843,0.282352941176471,0.600000000000000;0.0156862745098039,...    0.0196078431372549,0.0156862745098039];% 颜色1
& ^& e! \0 x6 h5 B& Y7 ~" p9 A6 Sy=[28,34,18];%柱状图纵坐标x=[1,2,3];     %柱状图横坐标%散点图的横纵坐标s1{1,1}=[0.9,0.95,0.95,1,1,1.05,1.1,1.1];s1{1,2}=[1.85,1.9,1.95,2,2,2,2,2.2];s1{1,3}=[2.7,2.7,2.95,3,3,3,3,3,3.2,3.2,3.2];s1{2,1}=[18,24,36,30,42,20,29,38];s1{2,2}=[26,22,24,30,38,42,48,50];s1{2,3}=[9,19,8,12,18,24,30,40,10,17,19];low=[4,3,3]; %误差图的上限high=[3,2,1]; %误差图的下限
%%初始化显示位置 [x,y,dx,dy] 从 x,y坐标开始，dx,dy为沿升的长度figure('Units', 'pixels', ...    'Position', [100 100 360 375]);
width=0.6; %柱状图宽度for  i=1:length(y)   %bar(i,y(i),width) width 柱状图宽占比  ，FaceColor 柱状图填充颜色 ，EdgeColor 柱状图边缘颜色   % LineWidth  边缘线宽    set(bar(i,y(i),width),'FaceColor',color(i,:),'EdgeColor',color(i+3,:),'LineWidth',2)   hold on      %散点图 大小，填充颜色   scatter(s1{1,i},s1{2,i},60,'filled','MarkerFaceColor',color(i+3,:))      %在柱状图 x,y 基础上 绘制误差 ,low为下误差，high为上误差，LineStyle 误差图样式，'Color' 误差图颜色     % 'LineWidth', 线宽,'CapSize',误差标注大小   errorbar(i, y(i), low(i), high(i), 'LineStyle', 'none', 'Color', color(i+3,:), 'LineWidth', 1.5,'CapSize',18);end
%绘制虚线plot([0,3.5],[50,50],'--','Color',color(7,:),'LineWidth',2)ylabel('Time on Warm Floor(%)')ax = gca;%x轴标签ax.XTick = [1 2 3];set(gca,'XTickLabel',{['\color[rgb]{',num2str(color(4,:)) ,'}','Pre-Injection'],['\color[rgb]{',num2str(color(5,:)) ,'}','1 hr CFA'] , ['\color[rgb]{',num2str(color(6,:)) ,'}','1wk CFA']},"LineWidth",2);%y轴标签ax.YTick= 0:20:80;%y轴范围ax.YLim=[0,80];%使得只显示x,y轴两根线box off%坐标轴字体和大小set(gca,"FontName","Times New Roman","FontSize",12,"LineWidth",2)title("TPT Timecourse","FontName","Times New Roman","FontSize",12,"LineWidth",2);
" i# S: A# z. }

uppzfaf5mtp6408891623.png

& f6 o, q2 z- K/ p3 E" ~4 y: i

clc;clear;close allcolor=[0.741176470588235,0.729411764705882,0.725490196078431;0.525490196078431,...    0.623529411764706,0.752941176470588;0.631372549019608,0.803921568627451,...    0.835294117647059;0.588235294117647,0.576470588235294,0.576470588235294;...    0.0745098039215686,0.407843137254902,0.607843137254902;0.454901960784314,...    0.737254901960784,0.776470588235294;0.0156862745098039,0.0196078431372549,0.0156862745098039];% 颜色1
%第二个图的y坐标y1=[18,48,38];x=[1,2,3];%第一图的散点图位置{1，}为横坐标，{2，}为纵坐标s2{1,1}=[0.7,0.75,1,1,1,1.05,1.25,1.25];s2{1,2}=[1.85,1.9,1.95,2,2,2.05,2.1,2.2];s2{1,3}=[2.7,2.75,2.95,3,3,3,3,3,3.2,3.2,3.2];s2{2,1}=[8,18,36,18,24,30,17,18];s2{2,2}=[49,51,50,68,52,51,52,50];s2{2,3}=[30,50,20,24,38,40,50,54,58,48,50];low=[4,3,3];high=[3,2,1];figure('Units', 'pixels', ...    'Position', [100 100 360 375]);%初始化显示位置marker={'^','o','s'};width=0.6; %柱状图宽度for  i=1:length(y1)   set(bar(i,y1(i),width),'FaceColor',color(i,:),'EdgeColor',color(i+3,:),'LineWidth',2)   hold on   scatter(s2{1,i},s2{2,i},70,'filled',marker{1,i},'MarkerFaceColor',color(i+3,:))   errorbar(i, y1(i), low(i), high(i), 'LineStyle', 'none', 'Color', color(i+3,:), 'LineWidth', 1.5,'CapSize',18);end%绘制虚线plot([0,3.5],[50,50],'--','Color',color(7,:),'LineWidth',2)ylabel('Time on Warm Floor(%)')ax = gca;ax.XTick = [1 2 3];set(gca,'XTickLabel',{['\color[rgb]{',num2str(color(4,:)) ,'}','CFA-con'],['\color[rgb]{',num2str(color(5,:)) ,'}', 'CFA-mon'] , ['\color[rgb]{',num2str(color(6,:)) ,'}',' CFA-Angle']},"LineWidth",2);% ax.XTickLabels ={'CFA-con', 'CFA-mon', ' CFA-Angle'};ax.YTick= 0:20:80;ax.YLim=[0,80];box off%坐标轴字体和大小set(gca,"FontName","Times New Roman","FontSize",12,"LineWidth",2)title({'TPT';'opQ' },"FontName","Times New Roman","FontSize",12,"LineWidth",2);

upkblbadsen6408891723.png

" L6 u4 m4 X5 ^# z+ {1 v& C
END


参考资料：
【1】图片来源（sciencemag.org）【10.1126/science.abe3040】

【送书活动】
包邮赠送2021年新书《MATLAB智能优化算法：从写代码到算法思想》共5本，这是一本MATLAB语言的新书籍，刚刚上架！由「 北京大学出版社」赞助提供 ，感兴趣的朋友推荐入手一本。
2 y6 L0 e' K2 o  O9 g

lkeuzrwadui6408891823.png

7 W. u0 f, ^- _! _& C
7 E: g1 P/ M, f% L8 ]# R( k【抽奖方式及满足条件】：
1 k- Z' G3 K7 P1.关注「好玩的MATLAB 」公众号
/ ^7 c9 m; L+ M" W. s2.给本文点【赞】+【在看】；
3.留言区评论点赞最多的前5名。
; |5 [+ D2 v- g8 F' P同时满足上述3个条件的读者朋友，包邮赠送一本：《MATLAB智能优化算法：从写代码到算法思想》
【开奖时间】 ：2022年3月3日夜晚8:00 (周四)
, Z& ^1 U- \0 F3 }【领奖方式】：在开奖时，后台将会回复点赞最多的前5名小伙伴的留言。后台回复的信息截图+姓名+地址+联系方式发送小编私人微信：idmatlab
- P1 _. {) s9 j( c9 b5 ^

当然，没有中奖但有需要的小伙伴可以通过下方链接进行购买：【本书亮点】
: f) u. h* U  F& m6 ?) U本书便以实战和快速教会智能优化算法初学者为主旨，通过9个常见的组合优化问题、5个经典的智能优化算法及4个新颖的智能优化算法，让读者全面、深入、透彻地理解智能优化算法求解问题时的算法设计思路及代码编写思路，进一步提高使用智能优化算法求解实际问题的实战能力。【内容简介】
2 a0 v* t8 K4 L' l; s! l本书以简单的组合优化问题作为MATLAB智能优化算法实战应用的切入点，逐步深入使用MATLAB编写更复杂的智能优化算法和求解更复杂的组合优化问题，让读者逐渐理解智能优化算法的实际求解过程。【目录一览表】
往下滑动查看目录详情（1）遗传算法求解0-1背包问题；
（2）变邻域搜索算法求解旅行商问题
（3）大规模邻域搜索算法求解旅行商问题
6 X2 }. ~7 H- T- x- Z3 N1 K1 ~（4）灰狼优化算法求解多旅行商问题
$ x! Z1 a0 ^5 N" J& @4 }1 ]（5）蚁群算法求解容量受限的车辆路径问题
（6）模拟退火算法求解同时取送货的车辆路径问题
（7）遗传算法求解带时间窗的车辆路径问题
（8）萤火虫算法求解订单分批问题
（9）头脑风暴优化算法求解带时间窗和同时取送货的车辆路径问题
2 {" l7 E- Z" Q, G（10）鲸鱼优化算法求解开放式车辆路径问题。
往期精彩回顾

twpsia5tttq6408891923.png

7 l; i! b. M4 E  x0 M推荐 | 【问题】分享一些小伙伴问的问题 推荐 | 【高级绘图】Matlab绘制阴影误差图 推荐 | 【高级绘图】MATLAB绘制3D雷达图推荐 | 【高级绘图】MATLAB怎么将图形局部放大
$ F0 t( F* u: F) P/ E; S5 h

luvp1c5c1oz6408892023.png

1 L; u) t! L; G% o3 a$ ~3 I

ag5jy13z0yz6408892123.jpg

$ j" }+ l. X5 F) x" n+ C
↓↓↓ 点击"阅读原文" 【加入QQ群】