病态矩阵求解 matlab,病态矩阵的例子.ppt

病态矩阵的例子

从一个数值结果谈起 从这里至少可以看出三点，一理论上的满秩不能保证实际计算的满秩，有时甚至相差很大；二数学软件只能解决常态的问题，难以解决大量存在的有特殊性态的问题；三要研究非方阵、非满秩的问题。 * 以为基底，构造不超过n 次最佳平方逼近多项式时，其法方程的系数矩阵会出现所谓Hilbert矩阵。

n=5

>> A=hilb(5)

A =

1.0000 0.5000 0.3333 0.2500 0.2000

0.5000 0.3333 0.2500 0.2000 0.1667

0.3333 0.2500 0.2000 0.1667 0.1429

0.2500 0.2000 0.1667 0.1429 0.1250

0.2000 0.1667 0.1429 0.1250 0.1111

>> Eigenvalue=eig(A)

Eigenvalue =

0.0000

0.0003

0.0114

0.2085

1.5671

>> cond(A)

ans =

4.7661e+005

>> b=A*ones(5,1);

>> x=inv(A)*b;

>> x

solution =x

1.0000

1.0000

1.0000

1.0000

1.0000

>> norm(ones(5,1)-x)

ans &#