我想以编程方式找出 A\b 何时失败(对于稀疏 A),以便我可以运行一些特定于问题的逻辑。使用反斜杠运算符
A\b
我将警告打印到控制台,但我想以编程方式了解这些条件(单数或几乎单数),以便我可以做一些特定于问题的事情。
对于密集系统,我可以
[soln, cond_recip] = linsolve(A,b);
if cond_recip < 1e-15, ..., end
但是 linsolve 不适用于稀疏矩阵,我不想使矩阵变密。
我想以编程方式找出 A\b 何时失败(对于稀疏 A),以便我可以运行一些特定于问题的逻辑。使用反斜杠运算符
A\b
我将警告打印到控制台,但我想以编程方式了解这些条件(单数或几乎单数),以便我可以做一些特定于问题的事情。
对于密集系统,我可以
[soln, cond_recip] = linsolve(A,b);
if cond_recip < 1e-15, ..., end
但是 linsolve 不适用于稀疏矩阵,我不想使矩阵变密。
尝试以下操作:
%# temporarily set warning to issue errors (maybe there are others?)
s = warning('error', 'MATLAB:nearlySingularMatrix'); %#ok<CTPCT>
try
x = magic(4)\[34; 34; 34; 34];
catch ME
disp(ME.message)
%#.. problem specific stuff..
end
%# restore warning state
warning(s);