简介MATLAB用户如果需要快速复现论文级图表这份压缩包提供了一套覆盖面很广的绘图代码集。涵盖分组柱状图、热图差异气泡图、分层聚类分析图、组合泰勒图、半小提琴图、桑基图与弦图等18种复刻案例既有基础图表配色与布局调整也有树状图环形热图、截断误差棒柱状图、弦图末端弧形块单独上色等进阶技巧同时还包含灰色背景叠加、NaN图例地图、旋转相关系数热图等细节处理。压缩包大小约59MB各个案例相对独立便于按需调用可参考其中的坐标轴、误差棒、聚类分支与图例设计嵌入自己的数据中。已有179人学习适合希望在较短时间内掌握多种MATLAB绘图方法并迁移到科研写作中的中高阶用户。对于需要绘制出版级图表的科研人员而言这套合集能显著减少搜索与调参时间。1. 复刻了一批数据图形代码先弄清这批资源能解决什么问题实际处理论文配图或课题汇报图时最消耗时间的往往不是算法本身而是「同一组数据要换五种图形表达」这件事。收到的这批 MATLAB 绘图资源共包含 18 个复刻项目覆盖了数据分析中最常遭遇的图形场景分组柱状图、热图与差异气泡图、分层聚类分析图、组合泰勒图、折线图叠加误差棒与散点抖动、带树状图的环形热图、堆叠柱状图加哑铃图、桑基图、弦图、半小提琴图等。适合正在做科研绘图、实验数据可视化或者需要在毕业设计里一次性输出多种成熟图形的读者。这批代码的价值不在于某一张图画得多漂亮而在于它搭好了一套「拿到原始数据后按模板改路径就能出图」的框架。如果之前画图靠 Excel 或 Origin 手动摆数据这套资源能直接帮你把大批量数据的图形产出流程接进 MATLAB 脚本里。2. MATLAAB 复刻绘图的核心机制坐标区、对象句柄与批量导出要真正利用这批复刻项目需要先理解它们共同依赖的 MATLAB 绘图底座。每张复刻图本质上是几个基础对象按顺序叠加的结果figure管理画布axes决定坐标范围bar、heatmap、scatter等函数生成图层所有图层都以句柄handle的形式存在。所谓「复刻」在技术上就是拿到目标图片然后拆解它包含了哪几个图层再按顺序构造相同对象。% 建立一个基本画布和坐标区 fig figure(Color, [1 1 1], Units, centimeters, Position, [2 2 16 10]); ax axes(Parent, fig); hold(ax, on); % 后续图层都叠加在当前坐标区 grid(ax, on); % 生成示例数据两组处理 三个时间点 x 1:3; y1 [2.5 3.1 2.8; 2.9 3.4 3.0]; err1 [0.2 0.3 0.25; 0.25 0.2 0.3]; b bar(ax, x, y1, BarWidth, 0.7);这段代码是统计型复刻项目分组柱状图、误差棒、堆叠柱状图的公用骨架。bar函数在x1:3的三组位置上生成两条柱子y1的第 1 行对应第一组处理第 2 行对应第二组处理。hold(ax,on)让后续的errorbar、scatter都落在同一坐标区内而不是新开一个axes——这是复刻「组合图形」时最容易出问题的地方漏掉hold on会导致后面的图层直接覆盖掉之前的图。Position单位指定为centimeters方便在导出的 Word 或 PPT 里控制图片宽度一般两栏论文插图控制在 12 cm 到 16 cm 之间。2.1 读入原始数据时应保持的统一数据结构这批复刻项目里每张图的数据组织方式不完全相同但多数项目约定了一个隐式规则输入数据以矩阵或table形式存放行是样本/组别列是特征/时间点。比如分组柱状图中 y1 的行数等于柱子分组数heatmap 的行列需要配RowLabels和ColumnLabels。建议在读取任意项目前先把数据整理为「特征列 分组列」的宽表结构再按需 splitapply 分组。% 推荐的数据读取方式从 Excel 读取后转矩阵 dataTable readtable(实验数据.xlsx); matrixData dataTable{:, 2:end}; groupNames dataTable.Properties.VariableNames(2:end); sampleNames string(dataTable.组别);读入后先用size(matrixData)确认维度再检查是否存在NaN。在热图与聚类分析这类项目里NaN并不是被忽略而是会被当作缺失值参与计算直接影响距离矩阵的实际取值。至于从 Excel 读入的字符串变cell数组的常见问题后续章节会专门给出处理方案。2.2 分组柱状图复刻的细节从条形布局误差棒「复刻一分组柱状图」以及「复刻十一截断含误差棒分组柱状图」都涉及一个共同埋点误差棒的位置和柱子的中心位置必须严格一一对应否则误差棒会整体偏移一列。MATLAB 的bar返回的句柄b的XEndPoints属性直接给出每个条形的中心 x 坐标直接用这个坐标画误差棒就不会错。% 为第一组处理添加误差棒 errorbar(b(1).XEndPoints, y1(1,:), err1(1,:), k, ... LineStyle, none, LineWidth, 1.2, CapSize, 6);这里b(1).XEndPoints是长度为 3 的向量正好对应第一组柱子的中心坐标第三个参数给出误差的大小。CapSize, 6控制误差线上下横线宽度数值过小时横线短到看不见偏大则抢柱子视觉焦点。误差棒这种看起来不起眼的元素在截断图复刻十一里还要额外处理纵轴截断和截断标记推荐做法是把纵轴范围上限设为实际最大值的 1.3 倍再画一条双斜线标记截断点。3. 热图系列复刻的进阶逻辑从基础热图到嵌套聚类树热图类是这批资源里最多的部分也是物理含义最丰富的一类共涉及分组热图、热图差异气泡图、分层聚类分析图、带树状图的环形热图、旋转相关系数热图等至少 5 种变体。要高效复用这些代码必须弄清它们在heatmap和clustergram两个函数上的分叉点。3.1 差异气泡图与热图的叠加方式「复刻七热图差异气泡图」在基础热图上再叠加一层气泡来同时表征两个变量热图颜色表示表达量/数值大小气泡半径表示差异倍数或 p 值。实现思路不是把气泡画在热图上方而是使用散点层透明叠加的手法本质上是绘制完热图后再用scatter在每个单元格中心绘制透明圆点。% 假设 heat 为 6x8 矩阵, pval 为同尺寸显著性矩阵 h heatmap(heat); h.Colormap parula; h.CellLabelFormat %.1f; % 覆盖一层差异气泡 figure; imagesc(heat); colormap(parula); hold on; [X, Y] meshgrid(1:size(heat,2), 1:size(heat,1)); scatter(X(:), Y(:), pval(:)*120, [0 0 0], filled, MarkerFaceAlpha, 0.35);这里的scatter第四个参数控制点的大小pval(:)需要是经过 min-max 归一化后的 0~1 之间的数值数组乘以 120 映射到像素级。实际复刻时若气泡过密互相遮挡需要先对 p 值做负对数转换再归一化。这部分的最终效果对坐标系的整数对齐非常敏感务必确认meshgrid的行列顺序与imagesc的显示方向一致否则气泡会整体错位到相邻格子里。3.2 分层聚类分析图的簇间距与标尺问题「复刻三分层聚类分析图」与「复刻十五环形聚类树状图」的核心差异在于一个使用clustergram输出矩形树状图另一个改用dendrogram重排图形并绘制环形布局。基于clustergram的版本是最成熟、最适合论文的因为它能自动完成归一化、特征筛选以及两侧的联动排列。但要注意它的两个参数Standardize设为row表示按行 z-score 标准化这决定了距离计算的实际语义Linkage指定聚类算法完整选项包括ward、average、complete。默认ward方差最小化适合大多数基因表达数据和实验测量数据但如果数据噪声较大比如重复间差异明显改average往往更稳。% 最小可用的聚类热图 cg clustergram(matrixData, ... RowLabels, sampleNames, ... ColumnLabels, featureNames, ... Standardize, row, ... Linkage, average, ... Cluster, 3);Cluster, 3强制聚为 3 簇此时左侧树状图会用不同颜色区分样本簇。注意这个参数在版本较新的 MATLAB 里要求输入数值向量或all等字符串如果提示非法直接删掉它改为在图上点击树状图分支来手动调整。环形聚类树状图复刻十五的技术难度较大需要先拿到聚类顺序再计算每个叶子的角度用polaraxes重新排布。但它在可视化形式上更适合组别数量少、只需要展示层次关系的场景。3.3 环形热图能做什么分组环形热图复刻六是一种接收率很高的「展示型图形」每个样本组沿圆环分布热图铺在环形带内常用在物种丰度、基因模块或传感器组的空间分布类数据上。核心代码段是利用极坐标绘制矩形块% ringHeatmap: 按角度和半径绘制热力块 thetaStart linspace(0, 2*pi, nGroups1); rMin 3; rMax 4; for g 1:nGroups for f 1:nFeatures theta1 thetaStart(g1); theta0 thetaStart(g); % 用 patch 绘制一个环形热力块 patch([theta0 theta1 theta1 theta0], ... [rMin rMin rMax rMax], heatVal(g,f), ... EdgeColor, none); end endpatch同时接受极坐标值括号里的前两个数组决定顶角角度范围后两个数组决定径向范围每次循环画一个环带。这类图形应用场景有限不建议在没有原始数据支撑硬复刻的前提下去做否则会呈现很强的「为了画图而画图」的感觉。4. 组合类图形与特殊坐标系泰勒图、桑基图、弦图的实现要点这批资源里还包含几类专业细分领域才用得到的图泰勒图复刻九、桑基图复刻十二/十四、弦图复刻四/十六。大多数普通数据分析场景用不到但一旦需要常规绘图函数里没有现成方案复刻价值很高。它们的共同点是需要在polaraxes或自定义几何坐标上完成大量人工映射。4.1 组合泰勒图的指标映射泰勒图用于评价模型的综合表现笛卡尔坐标系下同时编码相关系数、标准差比和中心均方根误差三个统计量。MATLAB 没有内置泰勒图函数复刻九的思路是先建立极坐标系再把散点线性映射上去最后叠加参考等值线。% 极坐标模板泰勒图 pax polaraxes; pax.ThetaZeroLocation right; % 0 度在右对应相关系数 1 pax.ThetaDir clockwise; % 顺时针增大角度 rlim(pax, [0 2.5]); % 半径表示标准差比 % 极坐标转换公式相关系数 r 映射角度标准差的归一化比值映射半径 corrVal [0.96 0.84 0.71]; stdRatio [1.05 1.82 1.33]; thetaRad acos(corrVal); % 角度 arccos(r) polarscatter(pax, thetaRad, stdRatio, 60, filled);acos(corrVal)把相关系数换算成极坐标中的角度相关系数不能大于 1否则会报错或者出现 NaNstdRatio是模式标准差与观测标准差的比值它和相关系数分别控制点的两个坐标。复刻组合泰勒图时通常会叠加两种不同模型的精度点使用不同颜色区分。实际使用时还要处理一个隐藏细节rlim(pax, [0 2.5])的半径范围必须大于所有stdRatio的上限否则点会跑到坐标轴外。半径上限取所有模型比值最大值的 1.2 倍更合理如果观察或模型数据的标准差接近 0比如归一化后泰勒图会失真此时应退回用散点图展示原始变量。4.2 桑基图与弦图的共同复用基础桑基图复刻十二/十四表达流量在节点间的分配关系。MATLAB 没有官方桑基图组件复刻代码实现方式基本为用patch或fill按从上到下的顺序绘制多层色带节点之间用梯形色带连接。左侧对齐默认与右侧对齐复刻十四的区别仅仅在于流量块的起点 x 坐标是固定还是按累计值反向排布代码主体逻辑相同。% 构建一个最小桑基图层三条流从 A 到 B/C/D nodes {A, B, C, D}; flow [0.4 0.35 0.25]; xLeft 0; xRight 1; yBase 0: cumsum(flow); % 从上到下累积 for k 1:3 yBottom yBase(k); yTop yBase(k1); patch([xLeft xLeft xRight xRight], ... [yBottom yTop yTop yBottom], ... colorMap(k,:), FaceAlpha, 0.75, EdgeColor, none); endcumsum(flow)是核心它把各流量占比换算成「从上往下累积的边界」这样每条流的厚度占比跟它的流量成比例。想做出右侧对齐的桑基图复刻十四只需要把 x 坐标轴颠倒从右向左绘制或用flip旋转所有几何对象的坐标。桑基图调试时的多数异常现象来自flow总和不为 1如 0.40.350.251.0 正常但如果从 Excel 读入的是百分比还要除以 100。弦图复刻四/十六是复刻难度最高的图纸类型需要将各端点分布在圆周上再绘制贝塞尔曲线连接不同端点并在弦末端加强弧形色标。复刻十六涉及一个复杂技巧给每段弦的末端单独上色本质上是把贝塞尔曲线的终点改成不同颜色值而不是统一颜色。这个功能在patch一对多构建弧形物体时需要给「每段弧」单独设置CData否则 MATLAB 默认会对所有对象应用同一个颜色FaceColor。% 弦图末端的单独弧色利用 patch 的每对象属性 hold on; for ii 1:numChords p patch(xx{ii}, yy{ii}, EdgeColor, chordColors(ii,:), ... FaceColor, none, LineWidth, 1.5); end这里采用了循环而不是向量化一次性绘制因为不同弦的颜色不同循环赋值更直观可维护。如果使用plot3绘制贝塞尔曲线再叠加线段则需要把每一条线段的Color属性单独赋一个 RGB 行向量这也是复刻十六中「末端单独上色」的实现方式。弦图初学时本质上就是在做几何绘图与数据统计无关。核心工作量在于贝塞尔控制点计算和标签防重叠角度调整。5. MATLAB 复刻绘图时最容易被误解的几个角落位置标签与图像导出所有复刻代码调试到最后排查的其实不是图形构造而是坐标信息如何在各个对象间共享。这里汇总几个本人复刻这批图纸时踩过的坑算是长期做 MATLAB 图形处理沉淀下来的一手经验。5.1NaN不做清洗的后果「复刻十三含 NaN 图例地图绘制」专门演示了一个常见现象地图上某个测站的数据是缺失值复刻时如果不处理对应的图例项会变成空白块或者拉出一条不存在的数据线。地图绘制里处理NaN的规范做法是validIdx ~isnan(dataValue); geoscatter(lat(validIdx), lon(validIdx), dataValue(validIdx)*100, filled);用逻辑索引一次性过滤掉所有含NaN的行或坐标而不是在画完图之后手动删图例项。如果需要在图例中保留「缺失」这一类别可以单独构造一个虚拟数据点并手动添加一个图例条目名为无数据。5.2 中文标签与字体环境这批资源涉及「标签旋转」和各类图例的内容特别多。在导出 PDF 或 EMF 时中文标签乱码是一种极其普遍的「水土不服」现象。处理方法是首行统一指定字体set(groot, defaultAxesFontName, 宋体); set(groot, defaultTextFontName, 宋体); set(groot, defaultLegendFontName, 宋体);如果电脑装了 SimHei黑体或微软雅黑可以直接替换。字体名设置有误时 MATLAB 才会自动替换为默认字体多数情况下只需把中文系统的字体名填对即可避免后续问题。标签旋转复刻四/十在heatmap里用h.XDisplayLabels或h.CellLabelFormat控制手动旋转需要用rotate配合位置属性。在散点图或弦图中想要旋转标签更直接的办法是text函数配合Rotation参数角度为 0 到 90 整数。5.3 导出清晰度与尺寸控制论文投稿对图分辨率和尺寸有硬性要求exportgraphics是当前推荐的导出函数不再建议使用已废弃的print加-r300的组合。完整导出语句建议按以下风格组织exportgraphics(gcf, output.png, Resolution, 600, ... BackgroundColor, white, ContentType, auto);Resolution参数直接控制 DPI600 是绝大多数期刊的硬性标准如果写的是ContentType, vector导出 PDF/EPS 格式还能保留矢量信息文字和线条在放大后不模糊。要检查当前图是不是软件渲染而非硬件渲染可以避免黑白或透明背景带来的问题。另外很多图形叠加灰色背景的实现复刻二「灰色背景」不是指figure的Color设置为灰而是要设置坐标区内底色最常见实现是set(ax,Color,[0.92 0.92 0.92])然后再叠加网格。这样能让点层更显眼且不会让整页变灰。6. 用三个指标判断复刻成果尺寸、对比度与重叠度这批资源本身不带统一的校验脚本使用半个月后的实际心得是复刻是否算成功不应只看「看起来像不像」而应从三个可量化维度去验证。第一个指标是输出图形的物理尺寸Position的宽高比与预期版本相差超过 15% 时文字和线条会在缩放时出现明显变形。第二个指标是色彩对比度热图类复刻项目中默认 parula 色带在灰阶打印下效果不佳如果目标期刊是黑白印刷建议提前用colormap(gray)模拟。第三个指标是散点/气泡的重叠度散点抖动复刻二的实际作用在于让重叠点分离用jitter或scatter增加随机偏移后一定要检查等价于sum(duplicate)的重复点计数是否明显下降否则抖动失效。最后一个细节建议给每个复刻项目外加一个exportgraphics的自动化批处理脚本把全部图形一次导出并统一重命名这样整理成支持材料或会议 PPT 时能省去大量手动导出操作。本文还有配套的精品资源点击获取