如何在Matlab中计算多年内同一天数据的3D矩阵平均值?

首先,我们去掉2月29日的所有数据,因为这些天都在数据中间,不会出现在每年的数据中,这会影响平均值:

Feb29=60+365*[1:4:32];
mean_Feb29=mean(GPH(:,:,Feb29),3);    % A matrix of 95x38 with the mean of all February 29th

GPH(:,:,Feb29)=[];           % omit Feb 29th from the data
Last_Jan_1=GPH(:,:,end);     % for Jan 1st you have additional data set, of year 2011
GPH(:,:,end)=[];            % omit Jan 1st 2011

re_GPH=reshape(GPH,95,38,365,[]);
av_GPH=mean(re_GPH,4);    

re_GPH 是一个95x38x365的矩阵,其中第三维中的每个切片是一年中一天的平均值,从1月1日开始,以此类推。

如果希望包含最后一个1月1日(2011年1月1日),请在上一个代码后运行此行:

av_GPH(:,:,1)=mean(cat(3,av_GPH(:,:,1),Last_Jan_1),3);

为了便于了解每个日期对应的切片数字,您可以制作一个包含一年中所有日期的数组:

t1 = datetime(2011,1,1,'format','MMMMd')
t2 = datetime(2011,12,31,'format','MMMMd')
t3=t1:t2;

现在,例如:

t3(156)=

  datetime

   June5

所以 av_GPH(:,:,156) 是6月5日的平均值。

对于您的评论,如果您想从平均值中减去每天:

sub_GPH=GPH-repmat(av_GPH,1,1,32);

对于2月29日,在从数据中删除它们之前,您需要这样做(上面的第3行):

sub_GPH_Feb_29=GPH(:,:,Feb29)-repmat(mean_Feb29,1,1,8);