警告:矩阵是奇异的工作精度:究竟什么生成这个警告?

6视图(30天)
我有一个与订购1458 x1458矩阵。这个结果测试模拟使用全套的麦克斯韦方程和混合有限元素。矩阵构造在FEniCS 2016.2和出口到matlab进行分析。这个矩阵构造使用自然边界条件在此设置(这个问题是齐次纽曼条件)并生成一个(理论上)奇异矩阵二维零空间(理论上证实)。当我分析这个矩阵,我发现:
(1)二维零空间如预期的那样,(运行计算和检查出最小奇异值我发现他们两个报告为零)。
(2)然而,矩阵是坏脾气的严重缺乏,rcond O (1 e-22)报告;的绝对值eig函数的输出,报道O (1 e-17)的最小值,即:,绝对最小的特征值。
一个发票在这个矩阵,这个矩阵或一个反斜杠解决不产生“奇异工作精度”警告,我想知道为什么不。即。,什么特征/属性矩阵的生成警告使用发票或反斜杠。

答案(1)

克里斯汀Tobler
克里斯汀Tobler 2021年3月24日
生成的警告是基于估计的条件(具体来说,1-norm互惠条件,如 //www.tatmou.com/help/matlab/ref/rcond.html)。
对密度矩阵,这使用LAPACK的条件估计量是相当准确的。
为稀疏矩阵,更启发式的选择是由估计这种情况数量,由于性能的原因。所以很可能不会注意到当一个矩阵是严重的。
另一个选择是使用conde函数来估计一个稀疏矩阵的条件数。
在实践中,它可能是更好的应用反斜杠,计算残余,然后决定在此基础上,如果矩阵是奇异的导致输出的准确性的问题。

类别

找到更多的在调整和重塑矩阵帮助中心文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!