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

在微气象观测中,经常会需要对短时间内相同时段的数据做分类汇总。但是由于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
为了方便,转置为列向量

  • 5
    点赞
  • 12
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值