八度音阶 nmf_bpas 错误：垂直尺寸不匹配（8x1 与 1x400）

Question

我在八度音阶中遇到了非负矩阵分解的问题。我试图从 Emg 数据中估计协同效应，但 octave 只能让我为两个或更多协同效应做这个，但不能为一个。我能够使用以下代码重现该问题。 nmf_bpas 来自 octave-forge 的线性代数 pkg。

V=rand(4, 20);
k=1; k2=2;
[W, H, Iter, HIS]=nmf_bpas(V,k);

使用此输入，我收到以下错误：

error: vertical dimensions mismatch (4x1 vs 1x20)
error: called from
    nmf_bpas>getStopCriterion at line 373 column 19
    nmf_bpas at line 266 column 26

当我如下定义 k>1 时，它会起作用

[W2, H2, Iter2, HIS2]=nmf_bpas(V,k2);

有了这个输入，它工作得很好，输出矩阵是 W2 (4x2) 和 H2 (2x20)。

同样的问题是，当我尝试使用其他 nnmf 代码并为矩阵 (4x1) 和 (1x20) 指定输入数据 Winit 和 Hinit 时（例如：[W, H] = nmf_pg (V, Winit, Hinit, tol, timelimit, maxiter))

在 matlab 中它与 nnmf 函数一起使用。

我会很高兴得到帮助

Answer 1

这似乎是nmf_bpas 中的一个错误。

据我所知，错误在第 373 行。更改

pGrad = [gradW(gradW<0|W>0); gradH(gradH<0|H>0)];

到

pGrad = [gradW(gradW<0|W>0); gradH(gradH<0|H>0)(:)];

它会起作用的。