我使用有限元法(FEM)来解决一些问题。刚度矩阵K是一个非常大的稀疏矩阵,我想解决这个等式:
U(freedofs,:) = K(freedofs,freedofs) \ F(freedofs,:)
,
freedofs
意味着自由节点。
有时我会在MATLAB中收到警告
MATLAB:singularMatrix
或MATLAB:nearlySingularMatrix
,
我计算的是condest(K(freedofs,freedofs))
值为Inf,但我仍然得到U的答案,这是非常合理的,我不知道是否有任何错误? U可靠的答案