Question

所以我有这段代码（在MATLAB中）

% Define vectorfield
 g1=[5,0,0];
 g2=[0,3,0];
 g3=[0,0,4];

% Define on-off
 u=[0;1];

% Define set to make field symmetric
symm=[1;-1];
k=1;

%% Generate possible combinations of vector fields
for a=1:length(u)
 for b=1:length(symm)
   for c=1:length(u)
     for d=1:length(symm)
       for e=1:length(u)
         for f=1:length(symm)
           allvecfields(k,:)=u(a).*symm(b).*g1+u(c).*symm(d).*g2+u(e).*symm(f).*g3;
           k=k+1;
         end
       end
     end
   end
 end
end

realfields=transpose(unique(allvecfields,'rows'));

realfields的每一列都是g的唯一正，负或零组合。我需要一些帮助来概括这一点。也就是说，每个g的大小可以是n，g的数量可以是m。代码仍应返回g的所有唯一可能组合。我有一种感觉，必须使用递归，但到目前为止我的所有尝试都失败了。

同样allvecfields(k,:)仅表示第k行，所有列。即使你的答案包含C / C ++或Java代码（没有任何特殊功能），对我来说也没问题。我将把它翻译成MATLAB。

我查看了combvec和allcomb文件，但他们没有做我需要的事情。例如transpose(unique(combvec(g1,g2,g3,-g1,-g2,-g3)','rows'))返回6x63矩阵而不是我想要的3x27。做

vals=transpose(unique(combvec(g1,g2,g3)','rows'));
vals=[vals transpose(unique(combvec(-g1,g2,g3)','rows'))];
vals=[vals transpose(unique(combvec(g1,-g2,g3)','rows'))];
vals=[vals transpose(unique(combvec(g1,g2,-g3)','rows'))];
vals=[vals transpose(unique(combvec(g1,-g2,-g3)','rows'))];
vals=[vals transpose(unique(combvec(-g1,g2,-g3)','rows'))];
vals=[vals transpose(unique(combvec(-g1,-g2,g3)','rows'))];
vals=[vals transpose(unique(combvec(-g1,-g2,-g3)','rows'))];
vals=unique(vals','rows');

是否给出我想要的东西，但这无助于概括。

编辑：修复了第一个代码块最后一行的错误。这种情况下所需的输出非常大（27列）但如果我们只有g1和g2那么输出将是：

realfields =

-5    -5    -5     0     0     0     5     5     5
-3     0     3    -3     0     3    -3     0     3
 0     0     0     0     0     0     0     0     0

编辑：根据评论中的建议，我已经能够将上述代码重写为，

u=[-1,0,1];
k=1;
for a=1:length(u)
  for b=1:length(u)
    for c=1:length(u)
        uMat(k,:)=[u(a) u(b) u(c)];
        k=k+1;
    end
  end
end

g1=[5,0,0];
g2=[0,3,0];
g3=[0,0,4];
gMat=[g1' g2' g3'];

for a=1:size(uMat,1)
    allvecfields(k,:)=sum(bsxfun(@times,gMat,uMat(a,:)),2);
end

realfields=transpose(unique(allvecfields,'rows'))

我认为这稍微优雅一点，但考虑到uMat中的列数，我仍然坚持如何动态生成gMat。我无法相信这个功能并不存在。任何帮助，将不胜感激。

Answer 1

这是解决方案

g1=[5,0,0];
g2=[0,3,0];
g3=[0,0,4];

gMat=[g1' g2' g3'];

m=size(gMat',1);
allvecfields = zeros(3^m,m);

for k = 1:3^m
 allvecfields(k,:) = double(dec2base(k-1,3,m)-'1');
end

realfields = (allvecfields*gMat')'

致信：Roger Stafford这个巧妙的解决方案。

MATLAB：向量总和的所有组合

1 个答案: