matlab计算

大海绵1990

大海绵1990 发表于 2014-7-5 11:11
应该是每五天累计的，还有，是三个金币不是5个金币，我刚来这里不久，没多少金币，不好意思啊

就是元胞数组里面是17个站点的数据，每个站点的数据是这样的，第一列是年，第二列是月，第三列是日，第四列是降水量，例如
1960   1     1      0
1960   1     2      3
。。。。
2011   12   31     0
我想知道怎么一次性把元胞数组里面的所有站点又有年份的数据都算出来

二爷名声在外

大海绵1990 发表于 2014-7-5 11:11
应该是每五天累计的，还有，是三个金币不是5个金币，我刚来这里不久，没多少金币，不好意思啊

要是五天累计那做一个循环不就出来了
for i=1:length-4
sum_5day=sum(i:1+4,4);
end

大海绵1990

沙颖凯 发表于 2014-7-5 10:30
http://bbs.06climate.com/forum.php?mod=viewthread&tid=24560

原来帮一个人写过类似的，是统计台站每 ...

for n = 1:17
    data = prec_all{n,1};
data(data<10)=0;
data(data>=10)=1;
% assume that your data contains 50 years and started from 1961
TDim=52; startYear=1960;
Result=zeros([TDim,1]);
% count rainfall days by a for loop
for i=1:TDim
    initial=datenum(startYear+i-1, 1, 1)-datenum(startYear, 1, 1)+1;
    final=datenum(startYear+i-1, 12, 31)-datenum(startYear, 1, 1)+1;
    Result(i)=sum(data(initial:final, 4));
end
end
这个用来算大于10mm的日降水日数挺好的，被我改了一下，可以出来最后一个站点的数据，怎么改可以出来我的17个站点的数据，求赐教啊！

rceclx

大海绵1990 发表于 2014-7-5 11:30
for n = 1:17
    data = prec_all{n,1};
data(data=10)=1;

你在第二个循环外，增加1个存放站点结果的向量即可。

大海绵1990

rceclx 发表于 2014-7-5 11:35
你在第二个循环外，增加1个存放站点结果的向量即可。

恩，没错的，谢谢！可是出来的结果里面有NAN怎么办

沙颖凯

大海绵1990 发表于 2014-7-5 14:02
恩，没错的，谢谢！可是出来的结果里面有NAN怎么办

有NaN应该是你原来的数据里有NaN吧，把sum(data(initial:final, 4));改成nansum(data(initial:final, 4));这样可以忽略原始数据中NaN的影响，或者你用个方法把原始数据里的缺测订正了也行。

大海绵1990

沙颖凯 发表于 2014-7-5 22:08
有NaN应该是你原来的数据里有NaN吧，把sum(data(initial:final, 4));改成nansum(data(initial:final, 4)) ...

恩，是这样的，已经解决了，谢谢你，那个95百分位我不知道该怎么计算，是一个极端降水指标，R95T，原定义是这样的，极端降水事件( 日降水量 1961- 1990年期间95%分位点) 的降水量之和占该年总降水量的百分率

皮特

http://bbs.06climate.com/forum.php?mod=viewthread&tid=24560

温文尔雅

有无大神回复程序呢？

jiefeng16

{:biggrin:}{:biggrin:}