【MATLAB】分时段分类汇总代码

在微气象观测中,经常会需要对短时间内相同时段的数据做分类汇总。但是由于Excel并不支持隔项分类汇总功能。笔者写了一段小代码来实现这个功能。

  • 待处理数据

【MATLAB】分时段分类汇总代码_第1张图片
5月到9月每天相同时段(0-23时)做汇总平均值

  • 代码
S=table2array(data);

导入名为data的数据,table转array是所有通过“导入数据”功能进行数据处理的必要步骤

for t=0:23
hour=S(:,3)';
O3=S(:,4)';

hour-数据中第三列为小时(分类项),并转置为行向量
O3-数据中第四列为臭氧浓度(汇总项),并转置为行向量

ind=find(hour==t);
O3hour=O3(ind);
O3hour_mean=mean(O3hour)

利用find函数,查找与t相同的hour所在索引,ind同样为数组
此处[参考文章],但事实上,find函数在参数==0时也可以查找0值(https://blog.csdn.net/Marvelous_Morty/article/details/96432373?utm_source=app&app_version=4.16.0&code=app_1562916241&uLinkId=usr1mkqgl919blen)
O3hour为根据索引ind数组查找O3中对应的值,并利用mean函数输出一个时段的平均值

result(t+1)=O3hour_mean;
end
result_=result'
xlswrite('station_ozone_hour_mean.xls',result_);

因为MATLAB数组是从1开始,而我们的t是0-23,所以此处需要(t+1)将各个时段的O3mean写入数组result
为了方便,转置为列向量

你可能感兴趣的:(MATLAB记录,matlab,分类)