我有一个二维单元格,其中每个元素都是a)空或b)一个长度不一的向量,其值范围从0到2.我想得到某个值出现的单元格元素的索引或甚至更好,每个出现一定值的“完整”指数。
我目前正在开发一种基于药剂的疾病传播模型,这样做是为了找到感染药物的位置。
提前致谢。
答案 0 :(得分:5)
我将如何做到这一点:
% some example data
A = { [], [], [3 4 5]
[4 8 ], [], [0 2 3 0 1] };
p = 4; % value of interest
% Finding the indices:
% -------------------------
% use cellfun to find indices
I = cellfun(@(x) find(x==p), A, 'UniformOutput', false);
% check again for empties
% (just for consistency; you may skip this step)
I(cellfun('isempty', I)) = {[]};
调用此方法1。
循环也是可能的:
I = cell(size(A));
for ii = 1:numel(I)
I{ii} = find(A{ii} == p);
end
I(cellfun('isempty',I)) = {[]};
调用此方法2。
比较两种速度方法如下:
tic; for ii = 1:1e3, [method1], end; toc
tic; for ii = 1:1e3, [method2], end; toc
给出
Elapsed time is 0.483969 seconds. % method1
Elapsed time is 0.047126 seconds. % method2
在Matlab R2010b / 32bit w / Intel Core i3-2310M@2.10GHz w / Ubuntu 11.10 / 2.6.38-13。这主要是由于JIT on循环(以及如何实现cellfun
和匿名函数的实现,mumblemumble ..)
无论如何,简而言之,使用循环:它更易读,并且比矢量化解决方案快一个数量级。