因此在某些輸入數(shù)據(jù)上使用MATLAB的svds()函數(shù)：

[U, S, V, flag] = svds(data, nSVDs, 'L')

我注意到從運(yùn)行到使用相同的數(shù)據(jù)運(yùn)行，我會(huì)從運(yùn)行到運(yùn)行得到截然不同的輸出SVD大小。當(dāng)我檢查'flag'是否已設(shè)置時(shí)，我發(fā)現(xiàn)它是，表明SVD沒有收斂。我的常規(guī)系統(tǒng)是，如果真的需要收斂，我會(huì)做這樣的事情：

flag = 1

svdOpts = struct('tol', 1e-10, 'maxit', 600, 'disp', 0);

while flag:

if svdOpts.maxit > 1e6

error('There''s a real problem here.')

end

[U, S, V, flag] = svds(data, nSVDs, 'L', svdOpts)

svdOpts.maxit = svdOpts.maxit*2

end

但是從我所知道的，當(dāng)你使用'L'作為第三個(gè)參數(shù)時(shí)，第四個(gè)參數(shù)被忽略了，這意味著我只需要處理它沒有收斂的事實(shí)？我甚至不確定如何使用'sigma'參數(shù)代替'L'參數(shù)。我也試過減少計(jì)算的SVD數(shù)量無濟(jì)于事。對(duì)此事的任何幫助都將不勝感激。

修改

在跟進(jìn)下面的評(píng)論時(shí)，我發(fā)現(xiàn)問題與我構(gòu)建數(shù)據(jù)矩陣的方式有關(guān)。原來我不小心顛倒了一個(gè)矩陣，輸入的大小(4000x1)而不是(20x200)，這是拒絕收斂的。 這不是問題

第二次修改

如果有人仍然關(guān)注這個(gè)，我實(shí)際上是錯(cuò)的，我的數(shù)據(jù)只是收斂了，因?yàn)槲野演斎脲e(cuò)誤地縮放了。這是一個(gè)在我生成數(shù)據(jù)時(shí)生成數(shù)據(jù)的程序：

% Generate data for SVD failure to converge

% Kernel functions

data_fun1 = @(t, tau)(exp(-t*(1./tau)));

t = linspace(0, 10, 26)';

tau1 = logspace(-1, log10(5), 150);

k1 = data_fun1(t, tau1);

gamma = 4257;

n = 6;

tau = 0.075;

A = -(2/3)*(2*pi*gamma)^2*n*tau.^3;

data_fun2 = @(V, t)exp(A*t*(V.^2));

V = linspace(0, 0.4, 29);

tau3 = logspace(-1, log10(5), 150)';

tau3 = tau3*1e-5;

k2 = data_fun2(V, tau3)';

svdOpts = struct('tol', 1e-10, 'maxit', 1e3, 'disp', 0);

svdOpts2 = svdOpts;

flag = 1;

while flag

if svdOpts2.maxit > 1e8

break

end

[U1, S1, V1, flag] = svds(k1, length(t), 'L', svdOpts);

svdOpts2.maxit = svdOpts2.maxit * 2;

end

flag

% flag == 0

flag = 1;

while flag

if svdOpts2.maxit > 1e8

break

end

[U2, S2, V2, flag] = svds(k2, length(V), 'L', svdOpts);

svdOpts2.maxit = svdOpts2.maxit * 2;

end

flag

% flag == 1

我也試過讓它跑到svdOpts2.maxit＆gt; 1e9，但整個(gè)周末都跑了，從未超過4.096e8。任何建議將不勝感激。

總結(jié)

以上是生活随笔為你收集整理的matlab 收敛,MATLAB：svds（）没有收敛的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。