7

我正在使用来自http://www.cis.upenn.edu/~jshi/software/Ncut_9.zip的标准化剪切包 (在 Windows 7 上)

这曾经与 Matlab2010a 一起工作得很好。但是我已经升级到 Matlab2013a(32 位学生版),现在我收到以下错误:

Error using arpackc
Expect 2 output arguments

Error in eigs_new (line 240)
        arpackc( aupdfun, ido, ...

Error in ncut (line 83)
[vbar,s,convergence] =
eigs_new(@mex_w_times_x_symmetric,size(P,1),nbEigenValues,'LA',options,tril(P));

Error in ncutW (line 9)
[NcutEigenvectors,NcutEigenvalues] =
ncut(W,nbcluster);

Error in NcutImage (line 18)
[NcutDiscrete,NcutEigenvectors,NcutEigenvalues]
= ncutW(W,nbSegments);

Error in demoNcutImage (line 25)
[SegLabel,NcutDiscrete,NcutEigenvectors,NcutEigenvalues,W,imageEdges]=
NcutImage(I,nbSegments);

很明显in中的new_eigs()函数ncuts与最新Matlab中的ARPACK版本不兼容。

  • 有人知道解决方法吗?
  • Normalized 使用 Matlabeigs()函数的修改版本。为什么我不能只使用 Matlab 的内置eigs()而不是eigs_new()(它似乎确实有效)。
4

3 回答 3

7

有两种解决方案:

  1. 使用 Matlab 的eigs()函数而不是eigs_new()规范化切割包中提供的函数。我想eigs_new()是为了解决与以前版本的 Matlab 的一些兼容性问题,现在它本身就引起了问题。

  2. 修改eigs_new(). 在有调用arpackc()的地方eigs_new(),像在工具箱版本中那样调用eigs()。所以,例如:

    arpackc( aupdfun, ido, ...
     bmat, intconvert(n), whch, nev, tol, resid, ncv, ...
        v, ldv, iparam, ipntr, workd, workl, lworkl, info );
    

变成

    [ido, info] = arpackc( aupdfun, ido, bmat, intconvert(n), whch, ...
        nev, tol, resid, ncv, v, ldv, iparam, ipntr, workd, workl, ...
        lworkl, info );

只有两个地方必须更改它才能使其工作(尽管还有其他调用arpackc()实际上并没有被执行)。

我仍然不确定为什么eigs_new()存在。

于 2013-05-08T15:10:38.860 回答
5

在 matlab2014a 64bit windows 中,你只需改变

[vbar,s,convergence] = eigs_new(@mex_w_times_x_symmetric,size(P,1),nbEigenValues,'LA',options,tril(P));

在 ncut.m 的第 81 行中

[vbar,s,convergence] = eigs(@mex_w_times_x_symmetric,size(P,1),nbEigenValues,'LA',options,tril(P));

然后 eigs_new.m 没用,不在乎。因为在 2014 年(或超过 2009 年) arpackc 已被 eigs 取代。您可以阅读 eigs 的帮助以获得更多帮助。

于 2014-12-22T03:59:24.767 回答
0

我曾经使用过这段代码。对于 32-64 位问题,我做了一些更改。我也把这eigs_new条线改成

[vbar,s,convergence] = eigs2(@mex_w_times_x_symmetric,size(P,1),nbEigenValues,'LA',options,tril(P)); 

应该管用...

于 2013-05-08T14:33:32.547 回答