Conditional cumsum - how to create?

5 次查看(过去 30 天)
This is probably easy, but my brain isn't working today...
How can you do the following operation in a vectorized way? I'd think it should be possible with some combination of cumsum, diff & logical indexing:
input = rand(10,1);
output = zeros(size(input);
output(1) = input(1);
for ind = 2:numel(input)
dif = input(ind) - input(ind-1);
if dif < 0
output(ind) = output(ind-1) + dif;
else
output(ind) = output(ind-1);
end
end
  2 个评论
the cyclist
the cyclist 2013-4-2
It would be useful if you also described conceptually what you are trying to do.
Eric Sampson
Eric Sampson 2013-4-2
Basically it's a copy of the input, but anytime that the original increases from one val to the next, the output should be hold constant. Sort of like a copy that can only go down :)

请先登录,再进行评论。

采纳的回答

Roger Stafford
Roger Stafford 2013-4-2
Try this.
outp = cumsum([inp(1);min(diff(inp),0)]);

更多回答(1 个)

Matt Tearle
Matt Tearle 2013-4-2
There may be better ways, but this works:
d = [true;diff(input)<0];
idx = find(d);
output = input(idx(cumsum(d)));
When the array is large enough, there's a pretty decent speedup (~50x)
  5 个评论
Sean de Wolski
Sean de Wolski 2013-4-2
As you know, the DWIM Toolbox still hasn't been released to the public.
Eric Sampson
Eric Sampson 2013-4-2
Sean, wasn't Loren or Steve supposed to be on that? Slackers.

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 MATLAB 的更多信息

产品

Community Treasure Hunt

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

Start Hunting!

Translated by