针对你这种跨44年7月的逐6小时数据，用MATLAB的datetime类型来处理时间轴是最省心高效的，直接上实操步骤：

1. 数据读取与时间格式标准化 #

首先读取你的数据文件（假设是txt/csv格式），把时间列转换成MATLAB原生的datetime对象，这是后续自定义轴标签的核心基础：

% 读取数据，假设第一列是时间点（数值型小时数/字符串格式都可）
data = readmatrix('your_data_file.csv');
time_raw = data(:,1);
param_data = data(:,2);

% 分两种常见情况处理时间：
% 情况1：如果time_raw是从1958-07-01 00:00开始的累计小时数
time_vec = datetime(1958,7,1,0,0,0) + hours(time_raw);

% 情况2：如果time_raw是字符串格式（比如'1958-07-01 00:00'）
% time_vec = datetime(time_raw, 'InputFormat', 'yyyy-MM-dd HH:mm');

2. 绘制曲线并自定义X轴刻度 #

直接用plot绘制数据，然后手动锁定X轴刻度为每年7月的起始日期，这样X轴只会显示44个年份标签，完全符合你的需求：

figure('Position', [100 100 1200 600]); % 设宽画布避免标签拥挤
plot(time_vec, param_data, 'LineWidth', 1);

% 生成每年7月1日的datetime数组作为刻度锚点
year_list = 1958:2001;
tick_dates = datetime(year_list, 7, 1);

% 设置X轴刻度和对应标签
xticks(tick_dates);
xticklabels(year_list); % 直接显示年份，简洁清晰

% 旋转X轴标签避免重叠，调整字体提升可读性
xtickangle(45);
set(gca, 'FontName', 'Arial', 'FontSize', 10);

3. 可视化细节优化 #

最后加上必要标注和美化，让图表更专业：

xlabel('年份（7月）', 'FontSize', 12);
ylabel('参数值', 'FontSize', 12);
title('1958-2001年7月逐6小时参数时间序列', 'FontSize', 14, 'FontWeight', 'bold');
grid on;

方案优势 #

• 用datetime原生处理时间，避免手动计算时间点对应关系，减少出错概率
• 自定义刻度直接定位到每年7月，X轴显示干净，无需额外分组处理数据
• 代码逻辑清晰，后续调整刻度（比如显示7月中旬）或添加其他元素都很灵活

内容的提问来源于stack exchange，提问作者Chayan Roychoudhury