1

我有一个大小为 1*n ( X=cell(1, n)) 的单元格数组。此单元格数组（X）中的每个单元格都有不同的大小。所以，为了清楚起见，这里有一个例子：

X = {X1, X2, X3} (n=3)
X1 = [1 3 5 9]
X2 = [1 2 3 6 7]
X3 = [1 7 8 9]

我有一个长度为 m 的向量。例如：

Y = [0 3 6 9] (m=4)

我想以这种方式根据此向量划分单元格数组中的单元格：

X1-1=[1] (because 1>=Y(1) & 1<Y(2)), X1-2=[3 5] (because 3 and 5>=Y(2) & 1<Y(3)), X1-3=[9] (because 9>=Y(3) & 9<Y(4))

等等 X2 和 X3 ......最后我有：

Z = {X1-1 X2-1 X3-1
     X1-2 X2-2 X3-2
     X1-3 X2-3 X3-3}

这是我的for循环代码：

Z = cell(size(X,2), size(Y,2));
for i = 1:size(X,2)
    T = cell2mat(X(i));
    for j = 1:size(Y,2)-1
        idx = (T > Y(j) & T < Y(j+1));
        Z{i,j} = {T(idx)'};
    end
end
Z = cellfun(@(z) cell2mat(z), Z, 'UniformOutput', false);

我想在不使用 for 循环的情况下获得相同的结果。任何矢量化想法？

4

1 回答 1

0

您可以将两个for循环替换为 onecellfun和 one arrayfun。但不一定更快：

Z = cellfun(@(Xi) arrayfun(@(k) Xi(Y(k)<=Xi&Xi<Y(k+1)), 1:numel(Y)-1, 'uni', 0).', X, 'uni', 0);

想要的结果存储在Z{1}{1}（对应你的X1-1）、Z{1}{2}（对应X1-2）等。

示例：我正在使用您的数据，但更改Y为[0 3 6 9.1]（请注意，您Y = [0 3 6 9]的X1-3示例中的 不正确）：

>> Z{1}{1}

ans =

     1

>> Z{1}{2}

ans =

     3     5

>> Z{1}{3}

ans =

     9

>> Z{2}{1}

ans =

     1     2
于 2013-11-30T00:40:31.863 回答