如何解决倍频程 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)(:)];
它会起作用。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。