Matlab 从开始/结束索引列表创建矢量化数组
我有一个两列矩阵Matlab 从开始/结束索引列表创建矢量化数组,matlab,vectorization,Matlab,Vectorization,我有一个两列矩阵M,其中包含一组区间的开始/结束索引: startInd EndInd 1 3 6 10 12 12 15 16 如何生成所有区间指数的向量: v = [1 2 3 6 7 8 9 10 12 15 16]; 我正在使用循环执行上述操作,但我想知道是否有更优雅的矢量化解决方案 v = []; for i=1:size(M,1) v = [v M(i,1):M(i,2)]; end 也许有一个更好
M
,其中包含一组区间的开始/结束索引:
startInd EndInd
1 3
6 10
12 12
15 16
如何生成所有区间指数的向量:
v = [1 2 3 6 7 8 9 10 12 15 16];
我正在使用循环执行上述操作,但我想知道是否有更优雅的矢量化解决方案
v = [];
for i=1:size(M,1)
v = [v M(i,1):M(i,2)];
end
也许有一个更好的解决方案我没有看到,但这里有一个版本使用
非常奇怪的解决方案是创建临时字符串并使用EVAL。也可以是一个班轮
tmp = cellstr(strcat(num2str(M(:,1)),{':'},num2str(M(:,2)),{' '}));
v = eval(['[' cell2mat(tmp') ']']);
我知道它在大型矩阵上可能不起作用。只是为了好玩。这里有一个向量化的解决方案,我想用它来解决这个特殊的问题,使用函数: 我没有IMFILL,但在我的机器上,这种方法比其他建议更快,我认为由于使用了find,它将优于IMFILL方法
如果将M设置为转置(并且我们调整arrayfun的第三个和第四个参数),它可以变得更快。是的,我知道有一个简单的解决方案。非常优雅+1正是我想要的。。感谢您,如果IDX可以很大,那么OPs解决方案实际上会更快,例如我的机器上的idxs=[1 3;6 11;12 12;25000 26000]。对于idxs=[1 3;6 11;12 12 12;25000 26000]等输入,此解决方案比gnovice的更快,但OP的简单For循环速度更快。对于OP的示例输入,这比gnovice的慢,也比OP的For循环慢。
tmp = cellstr(strcat(num2str(M(:,1)),{':'},num2str(M(:,2)),{' '}));
v = eval(['[' cell2mat(tmp') ']']);
v = zeros(1, max(endInd)+1); % An array of zeroes
v(startInd) = 1; % Place 1 at the starts of the intervals
v(endInd+1) = v(endInd+1)-1; % Add -1 one index after the ends of the intervals
v = find(cumsum(v)); % Perform a cumulative sum and find the nonzero entries
cell2mat(arrayfun(@colon,M(:,1)',M(:,2)','UniformOutput',false))