您所在位置：主题：解大型线性方程组遇到的问题
loveyjc
[专家分：0] 发布于 2011-09-01 22:45:00
用CVF编写的小型有限元程序，遇到解线性方程组，如AX=B，A为刚度矩阵，B为荷载向量，X为需要求解的节点位移。
进行计算，用高斯消去法，解出X没问题。然后根据X的值来，调整刚度矩阵A。由于A中的元素很大，对A中元素调整的数差了1E-5个量级，这时再用调整后的A矩阵，求X，所求结果与上次基本没变化。
比如 10000 10001
A 20000 25000 调整后的 A= 20000 25001
15000 25000 35000 15000 25000 35000
其中A为对称矩阵，类似这样的问题叫什么问题，可以得到理想的精度解吗？
肯定是高斯消去法不能处理的，程序中所有的数据均为双精度。
进行计算，用高斯消去法，解出X没问题。然后根据X的值来，调整刚度矩阵A。由于A中的元素很大，对A中元素调整的数差了1E-5个量级，这时再用调整后的A矩阵，求X，所求结果与上次基本没变化。
比如 10000 10001
A 20000 25000 调整后的 A= 20000 25001
15000 25000 35000 15000 25000 35000
其中A为对称矩阵，类似这样的问题叫什么问题，可以得到理想的精度解吗？
肯定是高斯消去法不能处理的，程序中所有的数据均为双精度。