matlab - 非常に小さな入力で行列方程式を解く

Question

私はMatlabでプログラミングしており、私のプログラムではシステムを解く必要があります。Ax=bここで、Aは非常に小さなエントリmをm持つ正方行列です。増加するとm、 のエントリはA小さくなります。

Aは疎行列なので、この行列の名前sparseを Matlab: の関数で変更しますB=sparse(A)。

特別な方法を使用して、Bx=bアンダーフローや数値エラーなしで解決できますか?

Answer 1

@Parag S. Chandakkar が指摘したように、に頼ることができますvpasolve。まず、条件の悪い行列でシステムを構築しましょう

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

このシステムを解決する通常の方法は、

x = A\b;

一方vpasolve、

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