(1)采样间隔为1h，求天平均的时间序列

 time_vec=datevec(time);

 [s,ii,jj]=unique(time_vec(:,1:3),'rows');

 data_daily=accumarray(jj,(1:numel(jj))',[],@(x)nanmean(data_hourly(x,:),1));

(2)采样间隔为3min，求天平均的时间序列

  %先求小时平均

 time_vec=datevec(time);

 [s,ii,jj]=unique(time_vec(:,1:4),'rows');

 data_hourly=accumarray(jj,(1:numel(jj))',[],@(x)nanmean(data_Min(x,:),1));

 再用上述方法求天平均。

(3)采样间隔为1h，求2h平均的时间序列

 flag=0;

 for i = 1:2:length(time_hourly)

   flag = flag + 1;

   time_2hourly=time_hourly(i,1);

   data_2hourly(flag,1)=nanmean(data_hourly(i:i+1,1),1);

 end