MATLAB双Y轴时间序列图：解决plotyy与datetick日期显示难题

1. 当时间轴遇上双Y轴：一个MATLAB绘图中的经典难题

如果你用MATLAB做过数据分析，尤其是处理那些带有时间序列、并且需要同时展示两种不同量纲数据（比如股价和成交量、温度和湿度）的图表，那么plotyy这个函数你大概率用过或者至少听说过。它曾经是MATLAB里绘制双Y轴图的“标准答案”，简单几行代码就能让两组数据共享同一个X轴，却拥有各自独立的Y轴刻度，非常直观。然而，当你兴冲冲地把时间数据作为X轴，准备用datetick函数把那些枯燥的数字转换成“2023-01-01”这样友好的日期格式时，问题就来了：你会发现datetick要么只对一个坐标轴生效，另一个Y轴对应的X轴刻度纹丝不动；要么就是转换后刻度标签错位、重叠，甚至直接报错。这感觉就像你组装好了一个精密的仪器，却发现两个表盘的指针无法同步——图表的核心功能还在，但可读性大打折扣，在报告或论文里根本拿不出手。

这个“plotyy+datetick”的组合难题，困扰过无数从学生到工程师的MATLAB用户。它本质上不是bug，而是MATLAB图形系统底层架构演进过程中留下的一个“历史兼容性”问题。plotyy函数在创建图形时，实际上在底层生成了两个重叠的坐标轴对象，一个用于左Y轴的数据，另一个用于右Y轴的数据。而datetick函数在设计上，一次只能作用于一个指定的坐标轴对象。更棘手的是，从MATLAB R2014b版本开始，图形系统转向了基于面向对象的全新架构，虽然带来了更强大的功能，但也让一些基于旧版句柄图形（Handle Graphics）的“技巧”变得不再稳定。网络上流传的解决方案五花八门，有的能用，有的已经失效，还有的会引发新的问题。本文将彻底拆解这个问题的根源，并提供一套经过实测、兼容新旧版本MATLAB的可靠解决方案，让你能轻松绘制出刻度清晰、日期格式正确的双Y轴时间序列图。

2. 理解问题的核心：plotyy的底层结构与datetick的工作机制

要解决问题，必须先理解问题是如何产生的。我们不能停留在“这个函数不好用”的层面，而要深入到MATLAB图形对象的层面去看。

2.1plotyy函数到底创建了什么？

当你调用[ax, h1, h2] = plotyy(x1, y1, x2, y2)时，MATLAB在幕后做了以下几件事：

1. 创建主坐标轴：首先，它创建一个我们可视为主画布的坐标轴对象。这个坐标轴的XAxisLocation设置为'bottom'，YAxisLocation设置为'left'。第一条数据(x1, y1)就被绘制在这个坐标轴上。ax(1)返回的就是这个左Y轴坐标轴的句柄。
2. 创建从属坐标轴：接着，它创建一个新的、与主坐标轴位置和大小完全重叠的坐标轴对象。但这个新坐标轴有一个关键设置：YAxisLocation被设置为'right'，并且它的Color属性被设为'none'（透明），这样我们才能透过它看到背后的主坐标轴。第二条数据(x2, y2)被绘制在这个透明的右Y轴坐标轴上。ax(2)返回的就是这个右Y轴坐标轴的句柄。
3. 链接X轴范围：为了保证两个数据序列共享相同的X轴视野，plotyy会将这两个坐标轴的XLim属性（X轴显示范围）强制同步。无论你缩放或平移哪一个，另一个都会跟着变化。

这里有一个至关重要的细节：虽然两个坐标轴共享“视觉上”的一个X轴，但在MATLAB的图形对象树中，这是两个独立的坐标轴对象，各自拥有自己的一套属性，包括XAxis（X轴对象）、XTick（刻度位置）、XTickLabel（刻度标签）等。这就是问题的根源所在。

2.2datetick函数如何工作？

datetick(axis_handle, 'x', date_format)函数的核心任务，是将指定坐标轴axis_handle的X轴上那些以“序列日期数”（如738522）或“日期时间向量”表示的刻度位置（XTick），转换成人类可读的日期字符串（如'01-Jan-2023'），并用这些字符串替换原有的刻度标签（XTickLabel）。

它的工作流程可以简化为：

1. 获取目标坐标轴当前的XTick值（一组数字）。
2. 将这组数字通过datestr或datetime函数（取决于MATLAB版本和输入）转换为日期字符串。
3. 将坐标轴的XTickLabel属性设置为这组字符串。
4. 一个重要的副作用：为了防止图形被后续的plot、xlim等操作自动重置刻度，datetick通常会将坐标轴的XTickMode和XTickLabelMode属性设置为'manual'（手动模式）。

2.3 冲突的诞生

现在，我们把两者结合起来看：

1. 你用plotyy画图，得到了两个句柄ax(1)和ax(2)。
2. 你调用datetick(ax(1), 'x', 'yyyy-mm-dd')。这个命令成功了，它修改了左Y轴坐标轴（ax(1)）的XTickLabel。
3. 但是，右Y轴坐标轴（ax(2)）的XTickLabel并没有被改变！因为它是一个独立的对象。所以，从右Y轴视角看到的X轴刻度，仍然是原始的序列日期数。
4. 更糟糕的是，由于两个坐标轴的XLim是链接的，当你用鼠标拖动或缩放图表时，MATLAB可能会尝试自动调整刻度（XTickMode为'auto'时），这可能会覆盖掉datetick在ax(1)上设置的XTickLabel，导致日期标签消失，变回数字。

所以，用户看到的现象就是：只有一个X轴显示了日期，或者日期标签时有时无，图形表现不稳定、不专业。

3. 解决方案一：手动同步坐标轴属性（最可靠的基础方法）

这是最根本、兼容性最好的方法。其核心思想是：既然plotyy给了我们两个坐标轴对象的句柄，那我们就对这两个句柄都执行一遍datetick操作，并手动确保它们的所有相关属性保持一致。

3.1 基础同步步骤

假设我们有以下数据：

% 生成示例数据：2023年1月1日到1月10日 x = datetime(2023, 1, 1):days(1):datetime(2023, 1, 10); y1 = rand(1, 10)*10 + 20; % 温度数据 y2 = rand(1, 10)*100 + 500; % 湿度数据 % 使用plotyy绘图 [ax, h1, h2] = plotyy(x, y1, x, y2); ylabel(ax(1), '温度 (°C)'); ylabel(ax(2), '湿度'); title('温度与湿度时间序列');

此时，两个坐标轴的X轴显示的都是日期时间对象，但刻度标签可能不是我们想要的格式。我们开始手动同步：

% 步骤1：对两个坐标轴分别应用datetick，使用相同的格式 datetick(ax(1), 'x', 'mm/dd', 'keepticks'); % ‘keepticks’ 保持当前刻度位置 datetick(ax(2), 'x', 'mm/dd', 'keepticks'); % 步骤2：确保两个坐标轴的X轴范围、刻度位置、刻度标签模式完全一致 % 获取ax(1)设置好的属性 xLimits = get(ax(1), 'XLim'); xTicks = get(ax(1), 'XTick'); xTickLabels = get(ax(1), 'XTickLabel'); % 将这些属性同步到ax(2) set(ax(2), ... 'XLim', xLimits, ... 'XTick', xTicks, ... 'XTickLabel', xTickLabels, ... 'XTickMode', 'manual', ... % 设置为手动模式，防止自动调整 'XTickLabelMode', 'manual'); % 步骤3：同样，将ax(2)的模式也同步给ax(1)以确保万无一失（虽然ax(1)已被datetick设为manual） set(ax(1), 'XTickMode', 'manual', 'XTickLabelMode', 'manual');

关键点解释：

• 'keepticks'参数：这个参数非常重要。它告诉datetick函数：“使用坐标轴上当前的刻度位置（XTick），只帮我转换标签，不要尝试计算新的、‘更美观’的刻度位置。” 这能避免两个坐标轴因datetick内部算法差异而产生不同的刻度位置。
• 手动模式（'manual'）：将XTickMode和XTickLabelMode设置为'manual'，是锁定当前刻度设置的关键。这能防止MATLAB在图形重置、调整大小或与其他绘图函数交互时，自动将刻度改回去。
• 属性同步：我们以ax(1)为“主”，将其设置好的XLim，XTick，XTickLabel直接复制给ax(2)。这样就保证了两个坐标轴在视觉上完全一致。

3.2 处理日期时间（datetime）对象与序列日期数

在上面的例子中，我们的x是datetime数组。从MATLAB R2014b开始，对datetime和duration对象的绘图支持越来越好。当你直接用datetime数组绘图时，plotyy可能会自动处理X轴标签，但格式可能不理想。此时，datetick可能不是最佳选择，因为datetick主要针对序列日期数。

更好的做法是直接设置坐标轴的XTickLabelFormat属性：

% 如果x是datetime数组，可以不用datetick [ax, h1, h2] = plotyy(x, y1, x, y2); % 直接设置日期显示格式 ax(1).XTickLabelFormat = 'MMM dd'; ax(2).XTickLabelFormat = 'MMM dd'; % 同样需要同步其他属性以确保一致 ax(2).XLim = ax(1).XLim; ax(2).XTick = ax(1).XTick;

如果你的原始数据是序列日期数（如从Excel读取的日期数字），那么datetick就是必需品。此时，上述3.1节的同步方法完全适用。

注意：在较新的MATLAB版本中，推荐使用点表示法（ax(1).XTickLabelFormat）来设置对象属性，这比旧的set/get函数更直观、高效。但为了兼容旧版本，本文也保留了set/get的写法。

4. 解决方案二：使用yyaxis函数（现代MATLAB的推荐方案）

从MATLAB R2016a开始，官方引入了yyaxis函数来替代plotyy。yyaxis采用了更现代的图形对象模型，它只创建一个坐标轴对象，但为其添加了两个独立的Y轴。这从根本上解决了双坐标轴属性不同步的问题。

4.1yyaxis的基本用法

% 准备数据 x = datetime(2023, 1, 1):days(1):datetime(2023, 1, 15); y_left = sin(0:14) * 5 + 10; y_right = cos(0:14) * 100 + 200; % 创建图形和坐标轴 figure; % 激活左侧Y轴并绘图 yyaxis left; plot(x, y_left, 'b-o', 'LineWidth', 1.5); ylabel('左侧指标 (单位A)'); ylim([0 20]); % 激活右侧Y轴并绘图 yyaxis right; plot(x, y_right, 'r--s', 'LineWidth', 1.5); ylabel('右侧指标 (单位B)'); ylim([100 300]); % 设置标题和X轴标签 title('使用yyaxis绘制双Y轴时间序列图'); xlabel('日期'); % 现在，轻松设置X轴日期格式！因为只有一个坐标轴对象。 % 方法A：使用datetick (如果x是序列日期数) % datetick('x', 'yyyy-mm-dd', 'keeplimits'); % 方法B：直接设置XTickLabelFormat (如果x是datetime，推荐) ax = gca; % 获取当前坐标轴句柄，现在只有一个！ ax.XTickLabelFormat = 'MM/dd'; % 例如设置为“月/日”格式 % 可以进一步美化 grid on; legend('左侧数据', '右侧数据', 'Location', 'best');

4.2yyaxis的优势与注意事项

核心优势：

1. 单一坐标轴对象：所有X轴属性（XTick，XTickLabel，XLim，XScale）都是唯一的。你只需要设置一次，左右两侧的数据都共享它。彻底解决了plotyy的同步难题。
2. 更简洁的API：通过yyaxis left和yyaxis right切换，逻辑清晰。
3. 更好的兼容性：与MATLAB新的图形系统结合更紧密，对datetime、categorical等新型数据类型的支持更好。

注意事项与技巧：

• 版本要求：确保你的MATLAB版本在R2016a或以上。
• 获取坐标轴句柄：在yyaxis绘图后，使用gca或ax = gca;获取的句柄就是那个唯一的坐标轴对象。你可以通过ax.YAxis(1)和ax.YAxis(2)来分别访问左、右Y轴对象，进行更精细的控制（如修改颜色）。
• 属性继承：在yyaxis left模式下设置的某些线条属性（如ColorOrder），可能会影响后续在yyaxis right模式下绘图的颜色。建议在每次plot时显式指定颜色。
• hold on的使用：在yyaxis模式下使用hold on需要小心。最好在第一次plot前就声明hold on，以确保所有图形元素被添加到正确的坐标轴上下文中。

强烈建议：如果你的项目不受限于旧版MATLAB，应优先使用yyaxis替代plotyy。它不仅解决了日期显示问题，也是MATLAB官方推动的现代化绘图方式，在未来会获得更好的支持和维护。

5. 解决方案三：封装实用函数与高级美化技巧

即使你决定继续使用plotyy（例如为了兼容旧代码），我们也可以将同步逻辑封装成一个函数，方便复用。同时，无论使用哪种方法，图表的美化都至关重要。

5.1 创建plotyy_with_datetick工具函数

你可以将以下代码保存为plotyy_with_datetick.m文件：

function [ax, h1, h2] = plotyy_with_datetick(x1, y1, x2, y2, date_format, varargin) % PLOTYY_WITH_DATETICK 增强版plotyy，自动同步双Y轴图的日期刻度。 % [AX, H1, H2] = PLOTYY_WITH_DATETICK(X1, Y1, X2, Y2, DATE_FORMAT) % 参数： % X1, Y1 - 左侧Y轴的数据。 % X2, Y2 - 右侧Y轴的数据。 % DATE_FORMAT - 传递给datetick的日期格式字符串，如'yyyy-mm-dd'。 % 如果为''或[]，则跳过datetick处理。 % VARARGIN - 可选的，传递给原始plotyy函数的参数对。 % % 返回句柄与标准plotyy相同。 % % 示例： % x = today-30:today; % [ax, h1, h2] = plotyy_with_datetick(x, rand(1,31), x, rand(1,31)*100, 'mm/dd'); % ylabel(ax(1), '数据A'); % ylabel(ax(2), '数据B'); % 调用原始plotyy if nargin > 5 [ax, h1, h2] = plotyy(x1, y1, x2, y2, varargin{:}); else [ax, h1, h2] = plotyy(x1, y1, x2, y2); end % 如果提供了日期格式，则处理日期刻度 if nargin >= 5 && ~isempty(date_format) % 确保输入是数值（序列日期数）以供datetick使用 % 如果输入是datetime，转换为序列日期数 if isdatetime(x1) x1_numeric = datenum(x1); else x1_numeric = x1; end % 设置X轴范围为数据范围，避免空白 x_combined = [x1_numeric(:); (isdatetime(x2) ? datenum(x2(:)) : x2(:))]; x_min = min(x_combined); x_max = max(x_combined); buffer = (x_max - x_min) * 0.02; % 增加2%的边距 set(ax, 'XLim', [x_min-buffer, x_max+buffer]); % 对两个坐标轴应用datetick，使用'keepticks'保持刻度一致 datetick(ax(1), 'x', date_format, 'keepticks', 'keeplimits'); datetick(ax(2), 'x', date_format, 'keepticks', 'keeplimits'); % 强制同步所有X轴相关属性，并设置为手动模式 sync_xaxis_properties(ax); end end function sync_xaxis_properties(ax) % 同步两个坐标轴的X轴属性 xLim = get(ax(1), 'XLim'); xTick = get(ax(1), 'XTick'); xTickLabel = get(ax(1), 'XTickLabel'); set(ax(2), ... 'XLim', xLim, ... 'XTick', xTick, ... 'XTickLabel', xTickLabel, ... 'XTickMode', 'manual', ... 'XTickLabelMode', 'manual'); set(ax(1), 'XTickMode', 'manual', 'XTickLabelMode', 'manual'); end

这个函数自动化了属性同步的过程，你只需要关心数据和日期格式。

5.2 图表美化与实战技巧

一个专业的图表不仅功能正确，还要美观易读。以下是一些针对双Y轴时间序列图的美化技巧：

1. 区分左右数据：

   • 颜色：这是最直接的区分方式。通常，左Y轴数据和其坐标轴标签、刻度线使用一种颜色（如蓝色），右Y轴使用对比色（如红色）。

   % 使用plotyy时 [ax, h1, h2] = plotyy(x, y1, x, y2); set(h1, 'Color', 'b', 'LineWidth', 2); set(h2, 'Color', 'r', 'LineStyle', '--', 'LineWidth', 2); set(ax(1), 'YColor', 'b'); % 左Y轴颜色 set(ax(2), 'YColor', 'r'); % 右Y轴颜色 % 使用yyaxis时更简单 yyaxis left; plot(x, y1, 'b-', 'LineWidth', 2); yyaxis right; plot(x, y2, 'r--', 'LineWidth', 2); ax = gca; ax.YAxis(1).Color = 'b'; ax.YAxis(2).Color = 'r';

2. 优化日期刻度密度： 时间序列数据点过多时，自动生成的日期刻度可能会过于密集，导致标签重叠。需要手动调整。

   % 假设x是datetime数组，我们想每3天显示一个标签 ax = gca; % 或 ax(1) for plotyy % 设置刻度位置 ax.XTick = x(1:3:end); % 从第1个开始，每3个取一个 % 确保刻度标签格式 ax.XTickLabelFormat = 'MMM-dd'; % 如果需要旋转标签防止重叠 ax.XTickLabelRotation = 45;

3. 添加图例和网格：

   legend([h1, h2], {'温度 (°C)', '湿度'}, 'Location', 'northwest'); grid(ax(1), 'on'); % 对于plotyy，通常对主坐标轴加网格即可 % 或对于yyaxis/gca grid on;

4. 处理非均匀时间戳： 如果时间数据不是等间隔的（例如股票交易数据，没有周末），datetick或自动刻度可能会在无数据的位置也生成标签。此时，最好根据实际数据点的时间来手动指定XTick。

   % x_dates 是datetime或序列日期数 ax.XTick = x_dates; % 在每个数据点位置显示刻度 ax.XTickLabel = datestr(x_dates, 'mm/dd'); % 手动生成标签 ax.XTickLabelRotation = 90; % 垂直旋转以节省空间

6. 疑难排查与版本兼容性指南

即使按照上述方法操作，在实际使用中仍可能遇到一些“怪现象”。这里汇总一些常见问题及排查思路。

6.1 日期标签在交互后消失或还原为数字

现象：用datetick设置好日期标签后，一旦用鼠标缩放、平移图形，或者在该图形窗口内绘制新的图形，日期标签就变回了原始数字。

根因：坐标轴的XTickMode和XTickLabelMode属性被重置为'auto'。当图形交互或更新时，MATLAB的自动刻度计算机制被触发，覆盖了手动设置的标签。

解决方案：

• 彻底方案：如前所述，在应用datetick后，务必将这两个属性设置为'manual'。

  set(ax, 'XTickMode', 'manual', 'XTickLabelMode', 'manual');

• 检查时机：确保这段设置代码在所有可能改变坐标轴状态的操作（如plot，xlim，ylim，hold on/off）之后执行。

6.2 左右X轴刻度位置对不齐

现象：两个坐标轴的日期标签虽然都是日期格式，但刻度线（tick mark）的位置没有严格对齐，导致视觉上错位。

根因：plotyy创建的两个坐标轴初始的XTick可能因为数据范围或自动布局算法而有细微差异。datetick的'keepticks'参数是基于各自当前的XTick进行转换的。

解决方案：

1. 强制统一XTick：在调用datetick之前，先计算一个合理的刻度位置，然后同时设置给两个坐标轴。

   % 计算一个合适的刻度向量，例如在数据范围内等间距取5个点 x_numeric = datenum(x); % 如果x是datetime tick_locations = linspace(min(x_numeric), max(x_numeric), 5); set(ax(1), 'XTick', tick_locations); set(ax(2), 'XTick', tick_locations); % 然后再应用datetick datetick(ax(1), 'x', 'yyyy-mm-dd', 'keepticks'); datetick(ax(2), 'x', 'yyyy-mm-dd', 'keepticks');

2. 使用linkaxes函数（谨慎）：linkaxes([ax(1), ax(2)], 'x')可以链接两个坐标轴的X轴范围，但通常plotyy已经做了这件事。链接主要有助于交互时保持同步，但不能保证初始XTick一致。

6.3 不同MATLAB版本的差异处理

• R2014b之前（旧版句柄图形）：plotyy是主流选择。本文的“手动同步属性”方案（第3节）在此版本上完全有效。注意使用set和get函数。
• R2014b ~ R2015b（过渡期）：MATLAB引入了新的图形系统，但plotyy仍可用。可能会出现一些属性设置行为的变化。如果遇到问题，尝试在设置属性后加上drawnow命令强制刷新图形。
• R2016a及以后（现代图形系统）：强烈推荐使用yyaxis。如果必须使用plotyy，其行为可能与旧版略有不同，但上述同步方法通常仍有效。你可以通过graphicsversion函数检查图形系统版本。
• 检查函数是否存在：在编写兼容性代码时，可以先判断yyaxis是否存在。

  if exist('yyaxis', 'builtin') || exist('yyaxis', 'file') % 使用yyaxis figure; yyaxis left; ... else % 使用plotyy及兼容方案 [ax, h1, h2] = plotyy(...); % ... 应用同步代码 end

6.4 从plotyy迁移到yyaxis的代码改写示例

假设有一段旧的plotyy代码：

% 旧代码 (使用plotyy) figure; [ax, h1, h2] = plotyy(time_serial, data1, time_serial, data2, 'plot', 'stem'); set(get(ax(1),'Ylabel'),'String','数据1'); set(get(ax(2),'Ylabel'),'String','数据2'); title('旧式双轴图'); datetick(ax(1), 'x', 'HH:MM'); % ... 需要手动同步ax(2)的代码

可以改写为：

% 新代码 (使用yyaxis， 假设MATLAB >= R2016a) figure; % 左侧轴 yyaxis left; if strcmpi(plot_style2, 'stem') % 处理原来的'stem'参数 stem(time_serial, data1, 'LineWidth', 1.5); else plot(time_serial, data1, 'LineWidth', 1.5); end ylabel('数据1'); % 右侧轴 yyaxis right; if strcmpi(plot_style2, 'stem') % 假设第二个图也是stem stem(time_serial, data2, 'LineWidth', 1.5); else plot(time_serial, data2, 'LineWidth', 1.5); end ylabel('数据2'); title('新式双轴图 (使用yyaxis)'); % 统一设置X轴日期格式 ax = gca; if isnumeric(time_serial) % 如果是序列日期数 datetick('x', 'HH:MM', 'keeplimits'); elseif isdatetime(time_serial) % 如果是datetime ax.XTickLabelFormat = 'HH:mm'; end grid on;

改写后，代码更简洁，且彻底避免了日期刻度不同步的问题。关键在于理解yyaxis left/right的切换逻辑，以及所有X轴属性现在由唯一的坐标轴对象gca控制。