两个单元阵列之间的叉积

Question

我有以下两个向量场：

>> orient

orient = 

    [1x3 double]    [1x3 double]    [1x3 double]
    [1x3 double]    [1x3 double]    [1x3 double]
    [1x3 double]    [1x3 double]    [1x3 double]

>> distance

distance = 

    [1x3 double]    [1x3 double]    [1x3 double]
    [1x3 double]    [1x3 double]    [1x3 double]
    [1x3 double]    [1x3 double]    [1x3 double]

我需要取成对元素的叉积，即

b = (cross(orient{1,1},distance{1,1}) + cross(orient{1,2},distance{1,2})..... and so on

然后重塑以匹配距离和方向的尺寸。

我可以在不使用 for 循环的情况下执行此操作吗？

如果我有

orient{1,1} = 

[1x3 double]

distance = 

        [1x3 double]    [1x3 double]    [1x3 double]
        [1x3 double]    [1x3 double]    [1x3 double]
        [1x3 double]    [1x3 double]    [1x3 double]

我该怎么做

sum1 = (cross(orient{1,1},distance{1,1}) + cross(orient{1,1},distance{1,2}) +...) sum2 = (cross(orient{1,2},distance{1,1}) + cross(orient{1,2},distance{1,2}) +...)

其中每个“和”只是单个 orient 元素的迭代，与所有距离元素相交，然后将这些叉积​​相加。然后我会：

mastersum = sum1 sum2 sum3
            sum4 sum5 sum6
            sum6 sum8 sum9

在哪里

sum1 = 

[1x3 double]

我只是把它放在一个令人困惑的地方吗？

Answer 1

您需要使用cellfun 来遍历没有for-循环的元胞数组。

对于两个向量场（两个元胞数组），你应该这样做：

crosses = cellfun(@(u, v)cross(u, v)', orient, distance, 'UniformOutput', 0);
b = sum(cell2mat({crosses{:}})', 1)  %# Summing all vectors in all cells

来自orient 的单个单元格的类似过程，比如orient{1, 2}，将是：

u = orient{1, 2};
crosses = cellfun(@(v)cross(u, v)', distance, 'UniformOutput', 0);
b = sum(cell2mat({crosses{:}})', 1)  %# This command remains the same

要从 orient 获取所有向量的结果而不使用 for 循环，请改为：

b_func = @(u)sum(cell2mat(cellfun(@(v)cross(u,v)', {distance{:}}, 'Un', 0))', 1);
U = cellfun(b_func, orient, 'UniformOutput', 0)

现在U 也是一个元胞数组（与orient 具有相同的维度）：U{1, 1} 具有orient{1, 1} 的交叉总和，U{1, 2} 具有orient{1, 2} 的交叉总和，等等...