閾値までベクトル値を合計してから、もう一度やり直してください

Question

ベクトルがa = [1 3 4 2 1 5 6 3 2]です。今度は、cumsumがaの新しいベクトル 'b'を作成したいが、閾値に達した後、cumsumをリセットして再びしきい値に達するまで再開する必要があるので、新しいベクトルは次のようになります。

b = [1 4 4 2 3 5 6 3 5] 

Answer 1

1つの方法はループを使用することです。最初の累積合計csを作成し、csの要素がしきい値thより大きい場合は、残りの要素の累積合計の要素をaに置き換えます。

aの一部の要素がthより大きい可能性があるため、これらの要素も削除しない限り、このループは無限になります。

a = [1 3 4 2 1 5 6 3 2]; 
th = 5; 
cs = cumsum(a); 
while any(cs>th & cs~=a) % if 'cs' has values larger that 'th', 
         % and there are any values smaller than th left in 'a' 
    % sum all the values in 'a' that are after 'cs' reached 'th', 
    % excluding values that are larger then 'th' 
    cs(cs>th & cs~=a) = cumsum(a(cs>th & cs~=a)); 
end 

Answer 2

累積合計を計算し、条件に従ったインデックス値を置き換えます。

a = [1 3 4 2 1 5 6 3 2] ; 

b = [1 4 4 2 3 5 6 3 5] ; 

iwant = a ; 
a_sum = cumsum(a) ; 
iwant(a_sum<5) = a_sum(a_sum<5) ; 

Answer 3

閾値とcumsumの引数を正規化し、その後、cumsum明のGroupWiseを行うことができますaccumarray、のグループ化indizesを得ることができますフローリングで：：ここで

は whileループを備えたシンプルなソリューションです

t = 5; 
a = [1 3 4 2 1 5 6 3 2]; 

%// cumulative sum of normalized vector a 
n = cumsum(a/t); 
%// subs for accumarray 
subs = floor(n) + 1; 
%// cumsum of every group 
aout = accumarray(subs(:), (1:numel(subs)).', [], @(x) {cumsum(a(x))}); 
%// gather results; 
b = [aout{:}] 

Answer 4

元のベクトルを乗算すると累積合計を返す疎行列を作成できます。私はこのソリューションと他のソリューションとのタイミングを合わせることはできませんでしたが、私はこれが大規模な配列の中で最も速くなると強く思っています。

% Original data 
a = [1 3 4 2 1 5 6 3 2]; 
% Threshold 
th = 5; 
% Cumulative sum corrected by threshold 
b = cumsum(a)/th; 
% Group indices to be summed by checking for equality, 
% rounded down, between each cumsum value and its next value. We add one to 
% prevent NaNs from occuring in the next step. 
c = cumsum(floor(b) ~= floor([0,b(1:end-1)]))+1; 
% Build the sparse matrix, remove all values that are in the upper 
% triangle. 
S = tril(sparse(c.'./c == 1)); 
% In case you use matlab 2016a or older: 
% S = tril(sparse(bsxfun(@rdivide,c.',c) == 1)); 
% Matrix multiplication to create o. 
o = S*a.';