我一直被告知,在MATLAB中几乎所有的for循环都可以省略,而且它们通常会减慢进程有办法吗?:
我有一个单元格数组(tsCell
)tsCell
存储长度可变的时间数组我想为所有时间数组(InterSection
)找到一个相交时间数组:
InterSection = tsCell{1}.time
for i = 2:length{tsCell};
InterSection = intersect(InterSection,tsCell{i}.time);
end
最佳答案
这里有一个使用unique
和accumarray
的矢量化方法,假设输入单元格数组的每个单元格中没有重复项-
[~,~,idx] = unique([tsCell_time{:}],'stable')
out = tsCell_time{1}(accumarray(idx,1) == length(tsCell_time))
样本运行-
>> tsCell_time = {[1 6 4 5],[4 7 1],[1 4 3],[4 3 1 7]};
>> InterSection = tsCell_time{1};
for i = 2:length(tsCell_time)
InterSection = intersect(InterSection,tsCell_time{i});
end
>> InterSection
InterSection =
1 4
>> [~,~,idx] = unique([tsCell_time{:}],'stable');
out = tsCell_time{1}(accumarray(idx,1) == length(tsCell_time));
>> out
out =
1 4