我一直被告知,在MATLAB中几乎所有的for循环都可以省略,而且它们通常会减慢进程有办法吗?:
我有一个单元格数组(tsCelltsCell存储长度可变的时间数组我想为所有时间数组(InterSection)找到一个相交时间数组:

InterSection = tsCell{1}.time
for i = 2:length{tsCell};
    InterSection = intersect(InterSection,tsCell{i}.time);
end

最佳答案

这里有一个使用uniqueaccumarray的矢量化方法,假设输入单元格数组的每个单元格中没有重复项-

[~,~,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

10-06 06:25