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

133 阅读 0 评论 88 点赞

我是靠谱客的博主快乐曲奇，最近开发中收集的这篇文章主要介绍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)
%显示小刻度的值