calculate the accumulation of a variable over time within a loop?

7 次查看(过去 30 天)
Hello,
I have evaporation as daily values (m). I have calculated the weighted fieldmean, which gives me a vector with length = no. of days. However, I would like to display these values as monthly sums.
How can I do this in a loop?
For example:
newsum(1) = nansum(data(1:30)) and so forth...
  2 个评论
dpb
dpb 2014-9-17
Do you have a corresponding date/day of month vector or does the series simply begin at a known date? And, of course, year to know if leap year or not...
Michael
Michael 2014-9-17
it starts in 1993 1st January and ends 2010 31st December, i'm not sure how to do it other than every 30 days..

请先登录,再进行评论。

采纳的回答

Sean de Wolski
Sean de Wolski 2014-9-17
Once you get the month data, the code will be almost identical to this answer:
Instead of @std, you'll want to pass in @sum
  6 个评论
Michael
Michael 2014-9-17
I have read the documents but I can't seem to find a way of increasing the interval to 3 or 6 months, is there a way to do this?
Michael
Michael 2014-9-17
I have realised it makes more sense to calculate 3-month sums..I know how to do this the long way but wondered if I can modify this function

请先登录,再进行评论。

更多回答(1 个)

dpb
dpb 2014-9-17
编辑:dpb 2014-9-17
... to calculate 3-month sums...
There it's simpler to use another Matlab "trick" of recognizing and utilizing internal storage order and Matlab's handling of arrays. Generate the vector of monthly sums as above then to group by threes reshape and sum over the reshaped array and recast back...
mnthsum=accumarray(... % the monthly sums as above
threemnth=sum(reshape(mnthsum,3,[])).'; % sum over columns of 3 row/each
The above relies on the linear storage order and that Matlab orders arrays/matrices in column-major order and the default action of functions such as sum is to operate over columns. This "trick" is useful for any such operation over a given evenly-divisible subset of any array.
ADDENDUM
If it's decided that the monthly sums aren't needed any longer, the above can be done in place, of course...
mnthsum=sum(reshape(mnthsum,3,[])).';

类别

Help CenterFile Exchange 中查找有关 Loops and Conditional Statements 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by