做线性方程组计算的时候,想找一些病态矩阵来测试算法的优劣性。希尔伯特矩阵太难搞了,于是就想找一个没那么病态的,而且方便记忆和使用的。
这里我会分享一些碰到的病态矩阵的例子
A=(61861861⋱⋱⋱86186186)A=\begin{pmatrix} 6 & 1 & & &&\\ 8 & 6 & 1 & &&\\ & 8 & 6 & 1 & \\ & &\ddots &\ddots &\ddots\\ & & &8 &6 &1\\ & & & &8 &6 &1\\ & & & & &8 &6 \\ \end{pmatrix}A=⎝⎜⎜⎜⎜⎜⎜⎜⎜⎛6816816⋱1⋱8⋱6816816⎠⎟⎟⎟⎟⎟⎟⎟⎟⎞
供参考的线性方程组:
(61861861⋱⋱⋱86186186)(111⋮111)=(71515⋮151514)\begin{pmatrix}
6 & 1 & & &&\\
8 & 6 & 1 & &&\\
& 8 & 6 & 1 & \\
& &\ddots &\ddots &\ddots\\
& & &8 &6 &1\\
& & & &8 &6 &1\\
& & & & &8 &6 \\
\end{pmatrix}
\begin{pmatrix}
1\\
1\\
1\\
\vdots\\
1\\
1\\
1\\
\end{pmatrix}
=\begin{pmatrix}
7\\
15\\
15\\
\vdots\\
15\\
15\\
14\\
\end{pmatrix}
⎝⎜⎜⎜⎜⎜⎜⎜⎜⎛6816816⋱1⋱8⋱6816816⎠⎟⎟⎟⎟⎟⎟⎟⎟⎞⎝⎜⎜⎜⎜⎜⎜⎜⎜⎜⎛111⋮111⎠⎟⎟⎟⎟⎟⎟⎟⎟⎟⎞=⎝⎜⎜⎜⎜⎜⎜⎜⎜⎜⎛71515⋮151514⎠⎟⎟⎟⎟⎟⎟⎟⎟⎟⎞
2022年4月28日13:44:58
再分享一个不能用Gauss_Seidel,jacobi迭代法求解的例子
2022年6月13日18:07:27