我有一个$D$路Tensor$H$，维度为$I \times I \times \dots \times I$（$D$ times），表示多项式的系数。为了更好地理解，请参阅这篇文章：https://math.stackexchange.com/questions/4718760/higher-order-tensor-computations
我试图将H中某个位置的元素和这个位置索引的所有可能的唯一排列求和，并将它们重新排列在向量$F$中。向量$F$首先包含2阶单项式的系数，然后是3阶等等。我怎样才能将这段代码扩展到任何高阶多项式，而不必编写许多嵌套的for循环？
我正在寻找可能的方法，使下面的代码工作的任何给定的$Order$。目前，代码适用于$Order = 5$。对于更高的顺序，我必须添加更多的嵌套for循环。但这不是最佳解决方案。

pos = 1; start = 2; Memory = 3; Order = 3; I = Memory+2; H = randn(I*ones(1,Order));
F = zeros(1,1000);

if Memory > 0
    if Order > 1
        for i = start:I
            for j = i:I
                indices = num2cell(unique(perms([j,i,ones(1,Order-2)]),'rows'));
                [F,pos] = calcF(indices,F,pos,H);
            end
        end
    end
    if Order > 2
        for k = start:I
            for i = k:I
                for j = i:I
                    indices = num2cell(unique(perms([j,i,k,ones(1,Order-3)]),'rows'));
                    [F,pos] = calcF(indices,F,pos,H);
                end
            end
        end
    end
    if Order > 3
        for m = start:I
            for k = m:I
                for i = k:I
                    for j = i:I
                        indices = num2cell(unique(perms([j,i,k,m,ones(1,Order-4)]),'rows'));
                        [F,pos] = calcF(indices,F,pos,H);
                    end
                end
            end
        end
    end
    if Order > 4
        for n = start:I
            for m = n:I
                for k = m:I
                    for i = k:I
                        for j = i:I
                            indices = num2cell(unique(perms([j,i,k,m,n,ones(1,Order-5)]),'rows'));
                            [F,pos] = calcF(indices,F,pos,H);
                        end
                    end
                end
            end
        end
    end
end

function [F,pos] = calcF(indices,F,pos,H)
    for l = 1:size(indices,1)
        F(pos) = F(pos) + H(indices{l,:});
    end
    pos = pos + 1;
end

字符串