在《Matlab論文插圖繪制模板第56期—曲面圖（Surf）》中，我分享過曲面圖的繪制模板。

然而，有的時(shí)候，需要在一張圖上繪制兩個(gè)及以上的曲面圖，且每個(gè)曲面圖使用不同的配色方案。

在Matlab中，一張圖上只支持一種colormap/colorbar，所以想要繪制兩個(gè)及以上的曲面圖，需要大家自行設(shè)法解決。

本文利用freezeColors工具（John Iversen, MathWork, 2023），以及我自己制作的colorar_k2工具，進(jìn)行雙曲面的繪制，先來看一下成品效果：

特別提示：本期內(nèi)容『數(shù)據(jù)+代碼』已上傳資源群中，加群的朋友請(qǐng)自行下載。有需要的朋友可以關(guān)注同名公號(hào)【阿昆的科研日?！?#xff0c;后臺(tái)回復(fù)關(guān)鍵詞【繪圖桶】查看加入方式。

---

1. 數(shù)據(jù)準(zhǔn)備

此部分主要是讀取原始數(shù)據(jù)并初始化繪圖參數(shù)。

% 讀取數(shù)據(jù)load data.mat% 初始化繪圖參數(shù)% 曲面1x1 = X;y1 = Y;z1 = Z1;% 曲面2x2 = X;y2 = Y;z2 = Z2;

2. 顏色定義

作圖不配色就好比做菜不放鹽，總讓人感覺少些味道。

但顏色搭配比較考驗(yàn)個(gè)人審美，需要多加嘗試。

這里直接使用TheColor配色工具中的SCI權(quán)威配色庫：

%% 顏色定義map1 = TheColor('sci',2068);% map1 = flipud(map1);map2 = TheColor('sci',2073);

3. 雙曲面圖繪制

調(diào)用‘surf’和‘freezColors’命令，繪制初始雙曲面圖。

ax = gca;% 曲面1繪制s = surf(X,Y,Z1,'EdgeColor','none');caxis([min(Z1(:)) max(Z1(:))]);colormap(map1)freezeColors; hold on% 曲面2繪制s2 = surf(X,Y,Z2,'EdgeColor','none');caxis([min(Z2(:)) max(Z2(:))]);colormap(map2)freezeColors; % 標(biāo)題、標(biāo)簽、視角hTitle = title('DoubleSurface Plot');hXLabel = xlabel('x');hYLabel = ylabel('y');hZLabel = zlabel('z');view(-35,30)

4. 細(xì)節(jié)優(yōu)化

為了插圖的美觀，利用colorbar_k2工具添加顏色條，并對(duì)圖形細(xì)節(jié)等進(jìn)行美化：

% 添加顏色條colorbar_k2('right',Z1,map1,Z2,map2)% 坐標(biāo)區(qū)調(diào)整axes(ax)axis tightset(gca, 'Box', 'off', ...                                                          % 邊框         'LineWidth', 1, 'GridLineStyle', '-',...                                   % 坐標(biāo)軸線寬         'XGrid', 'on', 'YGrid', 'on', 'ZGrid', 'on',...                            % 網(wǎng)格         'TickDir', 'out', 'TickLength', [.01 .01], ...                             % 刻度         'XColor', [.1 .1 .1],  'YColor', [.1 .1 .1],'ZColor', [.1 .1 .1],...       % 坐標(biāo)軸顏色         'zlim',[0 700])% 字體和字號(hào)set(gca, 'FontName', 'Arial', 'FontSize', 11)set([hXLabel,hYLabel,hZLabel], 'FontName',  'Arial', 'FontSize', 11)set(hTitle, 'FontSize', 12, 'FontWeight' , 'bold')% 背景顏色set(gcf,'Color',[1 1 1])

設(shè)置完畢后，以期刊所需分辨率、格式輸出圖片。

%% 圖片輸出figW = figureWidth;figH = figureHeight;set(figureHandle,'PaperUnits',figureUnits);set(figureHandle,'PaperPosition',[0?0?figW figH]);fileout = 'test';print(figureHandle,[fileout,'.png'],'-r300','-dpng');

以上。