Use repelem, which is the built-in efficient solution for repeating elements:
x = [linspace(1,4,4)' linspace(5,20,4)']
I = [2 3 2 1];
y = repelem(x, I, 1)
A speed comparison:
x = rand(1e4, 2);
I = repmat(1:10, 1, 1000);
tic
for k = 1:100
C = cell(length(I),1) ;
for i = 1:length(C)
C{i} = repmat(x(i,:),I(i),1) ;
end
C = cell2mat(C);
end
toc
tic
for k = 1:100 % Slightly modified for I is a row vector:
idx = arrayfun(@(x,y)x*ones(1,y),1:numel(I),I,'UniformOutput',false);
C = x([idx{:}],:);
end
toc
tic
for k = 1:100
C = repelem(x, I, 1);
end
toc
% R2018b:
% Elapsed time is 3.676185 seconds. loop
% Elapsed time is 6.848004 seconds. arrayfun
% Elapsed time is 0.038391 seconds. repelem
