octave nmf_bpas error: vertical dimensions mismatch (8x1 vs 1x400)

octave nmf_bpas error: vertical dimensions mismatch (8x1 vs 1x400)

我在八度音程中遇到非负矩阵分解问题。我尝试从 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 函数一起工作。

我很乐意提供帮助

这似乎是 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)(:)];

它会起作用。