用很小的输入求解矩阵方程

Question

我在Matlab中编程，在我的程序中，我需要求解一个系统Ax=b，其中A是m乘m的平方矩阵，具有非常小的条目。如果m增加，则A的条目变小。

A是一个稀疏矩阵，所以我用Matlab的函数sparse将这个矩阵重命名为B=sparse(A)。

我可以使用一种特殊的方法来求解Bx=b而不出现下溢或数值错误吗？

Answer 1

如@Parag S.Chandakkar所述，您可以求助于vpasolve。首先，让我们构建一个条件较差的矩阵系统

n = 10;
A = hilb(n)
b = randn(n, 1);

解决该系统的常用方法是调用

x = Ab;

而CCD_ 12的是

xs = sym('x', [n 1]);
sol = vpasolve(A*xs==b);