MATLAB分时段分类汇总代码

Posted 2021-11-10 Wild Iris

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

• 待处理数据

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
为了方便，转置为列向量