MATLAB 时间变量 和 常用时间处理函数 datetime datestr datevec datenum

% 月份缩写
mon = {‘Jan’;‘Feb’;‘Mar’;‘Apr’;‘May’;‘Jun’;‘Jul’;‘Aug’;‘Sep’;‘Oct’;‘Nov’;‘Dec’};

几个年份：
日历年：当年1-12月
水文年：10月1日-次年9月30日 USGS station
气候年：4月1日-次年3月31日

按照季度划分时间段

% 在气候研究中

春季spring 3-5月
夏季summer 6-8月
秋季autumn 9-11月
冬季winter 12-次年2月

春季天数：31+30+31=92
夏季天数：30+31+31=92
秋季天数：30+31+30=91
冬季天数：31+31+28/29=90/91
每个季节的天数略有差异，因此采取补足天数的方法使分析时间段长度一致为92天

春季：[3,1] -[5,31]
夏季：[6,1] -[8,31]
秋季：[9,1]-[12,1] 【11月30日之后补一天】
冬季：[11,30]-[3,1]【平年补两天，前补11月30日，后补3月1日】或
[11,30]-[2,29]【闰年补一天，前补11月30日】

MATLAB的小时数举例：
t_spr = datetime([2005,3,1,0,0,0]):hours(1):datetime([2005,5,31,23,0,0]);  
t_sum = datetime([2005,6,1,0,0,0]):hours(1):datetime([2005,8,31,23,0,0]);  
t_fal = datetime([2005,9,1,0,0,0]):hours(1):datetime([2005,12,1,23,0,0]);  
t_win = datetime([2005,11,30,0,0,0]):hours(1):datetime([2006,3,1,23,0,0]);    % 平年
t_win = datetime([2005,11,30,0,0,0]):hours(1):datetime([2006,2,29,23,0,0]);   % 闰年

MATLAB的小时数举例：
t_spr = datetime([YR,3,1,0,0,0]):hours(1):datetime([YR,5,31,23,0,0]);  
t_sum = datetime([YR,6,1,0,0,0]):hours(1):datetime([YR,8,31,23,0,0]);  
t_fal = datetime([YR,9,1,0,0,0]):hours(1):datetime([YR,12,1,23,0,0]);  
t_win = datetime([YR,11,30,0,0,0]):hours(1):datetime([YR+1,3,1,23,0,0]);     % 平年
t_win = datetime([YR,11,30,0,0,0]):hours(1):datetime([YR+1,2,29,23,0,0]);   % 闰年

% 上述变量规律在于 前面的时间后三个为0 0 0，后面的时间后三个为23 0 0 即一天的开始（0点）到当天的23点（共24小时，24时就转到下一天的0时）。
% 这些数据的年份都是一样的 除了冬季 末尾年份需要加一，且需判断加一的时间年份是否为闰年。把2005替换为Yr即可便利的生成不同年份的时间段

按照月份划分时间段

%以上是按照季度划分的时间段，也可以按照月份划分时间段。由于月份天数的差异 28 29 30 31 因此分析时间长度也需要补足。
%1 3 5 7 8 10 12 这七个月为31天 不需要改动
% 4 6 9 11 这四个月为30天 往前 或往后各补一天 都可以 ，甚至4月前补1天，6月后补1天都可以，若考虑月份所属季节，则补同季节的月份 4月前补后补都可以 6月后补 9月后补 11月前补
% 2月 平年 28天 前补2天 后补1天
% 2月 闰年 29天 前补1天 后补1天
% 2月 结束天数均为 3月1日23时

t_JAN = datetime([YR,1,1,0,0,0]):hours(1):datetime([YR,1,31,23,0,0]); % 31 DAY

t_FEB= datetime([YR,1,30,0,0,0]):hours(1):datetime([YR,3,1,23,0,0]);  % 28 DAY          **2月 平年 28天  前补2天 后补1天**
or
t_FEB= datetime([YR,1,31,0,0,0]):hours(1):datetime([YR,3,1,23,0,0]);  % 29 DAY          **2月 闰年 29天  前补1天 后补1天**

t_MAR = datetime([YR,3,1,0,0,0]):hours(1):datetime([YR,3,31,23,0,0]);   % 31 DAY

t_APR= datetime([YR,4,1,0,0,0]):hours(1):datetime([YR,5,1,23,0,0]);  % 30 DAY  4月后补一天

t_MAY = datetime([YR,5,1,0,0,0]):hours(1):datetime([YR,5,31,23,0,0]);   % 31 DAY

t_JUN = datetime([YR,6,1,0,0,0]):hours(1):datetime([YR,7,1,23,0,0]);  % 30 DAY  6月后补一天

t_JUL = datetime([YR,7,1,0,0,0]):hours(1):datetime([YR,7,31,23,0,0]);   % 31 DAY

t_AUG = datetime([YR,8,1,0,0,0]):hours(1):datetime([YR,8,31,23,0,0]);   % 31 DAY

t_SEP = datetime([YR,9,1,0,0,0]):hours(1):datetime([YR,10,1,23,0,0]);  % 30 DAY  9月后补一天

t_OCT = datetime([YR,10,1,0,0,0]):hours(1):datetime([YR,10,31,23,0,0]);   % 31 DAY

t_NOV = datetime([YR,10,31,0,0,0]):hours(1):datetime([YR,11,30,23,0,0]);   % 30 DAY  11月前补一天

t_DEC= datetime([YR,12,1,0,0,0]):hours(1):datetime([YR,12,31,23,0,0]);   % 31 DAY

t_JAN = datetime([YR,1,1,0,0,0]):hours(1):datetime([YR,1,31,23,0,0])； % 31 DAY

t_FEB= datetime([YR,1,30,0,0,0]):hours(1):datetime([YR,3,1,23,0,0])；% 28 DAY 2月 平年 28天 前补2天 后补1天
or
t_FEB= datetime([YR,1,31,0,0,0]):hours(1):datetime([YR,3,1,23,0,0])；% 29 DAY 2月 闰年 29天 前补1天 后补1天

t_MAR = datetime([YR,3,1,0,0,0]):hours(1):datetime([YR,3,31,23,0,0])； % 31 DAY

t_ APR= datetime([YR,4,1,0,0,0]):hours(1):datetime([YR,5,1,23,0,0])；% 30 DAY 4月后补一天

t_MAY = datetime([YR,5,1,0,0,0]):hours(1):datetime([YR,5,31,23,0,0])； % 31 DAY

t_JUN = datetime([YR,6,1,0,0,0]):hours(1):datetime([YR,7,1,23,0,0])；% 30 DAY 6月后补一天

t_JUL = datetime([YR,7,1,0,0,0]):hours(1):datetime([YR,7,31,23,0,0])； % 31 DAY

t_AUG = datetime([YR,8,1,0,0,0]):hours(1):datetime([YR,8,31,23,0,0])； % 31 DAY

t_SEP = datetime([YR,9,1,0,0,0]):hours(1):datetime([YR,10,1,23,0,0])；% 30 DAY 9月后补一天

t_OCT = datetime([YR,10,1,0,0,0]):hours(1):datetime([YR,10,31,23,0,0])； % 31 DAY

t_NOV = datetime([YR,10,31,0,0,0]):hours(1):datetime([YR,11,30,23,0,0])； % 30 DAY 11月前补一天

t_DEC= datetime([YR,12,1,0,0,0]):hours(1):datetime([YR,12,31,23,0,0])； % 31 DAY

%% MATLAB 绘制 时间变量的例子和代码
http://blog.sina.com.cn/s/blog_6c640c790102xbhl.html

x = linspace(datenum('2013/1/1'),datenum('2014/1/1'),360);   %生成360个时间
y = sin(2*x+1);    % 生成时间序列y
plot(x,y);          % 画图
% 设计X轴坐标
N = 10;   %坐标轴上显示N个刻度
% x轴显示刻度的时间区间，以及区间数N
date_point = linspace(datenum('2013/1/1'),datenum('2014/1/1'),N);
set(gca,'xtick',date_point );  % x轴分m个小刻度
date_point_str = datestr(x,'yyyy-mm-dd HH:MM:SS');  %X轴刻度上对应的字符
set(gca,'xticklabel',date_point_str)  %显示小刻度的值