如何告诉倍频程索引是一个矩阵，而不是一个向量

Question

当索引具有可变尺寸矩阵的向量时，如何指示八度索引是矩阵而不是向量？

例如，在下面的代码中，在第四次迭代中，Octave将索引视为行向量，并且将索引操作的返回值转置为列向量（以匹配被索引的向量，而不是所使用的向量）。用于索引）。

当我跑步时：

v = rand(16,1);

t = magic(4);

f = @(m)(sum(m, 1));

for i = 4:-1:1
    s = t(1:i,:);
    f(v(s))
endfor

我得到：

ans =

   1.47780   2.28879   1.29786   2.98981

ans =

   1.24705   1.31940   0.87484   2.18276

ans =

   0.89387   0.55288   0.50312   1.61950

ans =  1.9294

看到在前三个迭代中，答案是行向量，但是在最后一个迭代中，答案是单例值。

我该如何告诉octave变量s应该是矩阵，而不是向量，并且索引操作的结果应始终与s形状相同？

Answer 1

解决此“功能”的一种方法是这样的：

v = rand(16,1);
t = magic(4);
f = @(m)(sum(m, 1));

for i = 4:-1:1
    w(1:i,:) = v(t(1:i,:));
    f(w(1:i,:))
end

这迫使v的输出以正确的方向存储在w中。 请注意，我通过w(1:i,:)入函数f而不仅仅是w因为你的反向索引for循环。

这是一种令人惊讶的常见模式，并且是重复设置矩阵的列等于行向量或矩阵的行等于列向量的有用方法，而无需执行转置甚至不知道向量的方向（谁知道Matlab / Octave JIT在幕后做了什么？

Answer 2

如何告诉八度，变量s应该是矩阵，而不是向量

在Matlab / Octave中，矩阵和矢量是同一回事。 向量只是一个1xm或mx1矩阵。

总和函数维选项仅在两个维都大于1时才有效。

只是使vector成为特例。

for i = 4:-1:1
    s = t(1:i,:);
    if i~=1
        f(v(s))
    else
        v(s)'
    end
end