倍频程 nmf_bpas 错误：垂直尺寸不匹配8x1 vs 1x400

问题描述

我在八度音阶中遇到非负矩阵分解问题。我尝试从 Emg-data 估计协同效应，但 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)(:)];

它会起作用。