第2讲 线性方程组的数值方法.ppt

文档介绍：线性代数方程组的数值解法(1)
Gauss 消去法
(Demos in Matlab: airfoil in 2D)
线性代数方程组的数值解法
直接法:Gauss 消去法,SuperLU
迭代法:定常迭代(Jacobi, GS, SOR, SSOR)
Krylov 子空间方法(CG, MINRES , GMRES, QMR, BiCGStab)
The Landscape of Ax=b Solvers
Direct
A = LU
Iterative
y’= Ay
Non-
symmetric
Symmetric
positive
definite
More Robust
Less Storage
More Robust
More General
刘徽(约220-280)
Gaussian elimination, which first appeared in the text Nine Chapters on the Mathematical Art written in 200 BC, was used by Gauss in his work which studied the orbit of the asteroid Pallas. Using observations of Pallas taken between 1803 and 1809, Gauss obtained a system of six linear equations in six unknowns. Gauss gave a systematic method for solving such equations which is precisely Gaussian elimination on the coefficient matrix. (The MacTutor History of Mathematics, -)
Gauss(1777-1855)
今有上禾三秉,中禾二秉,下禾一秉,实三十九斗;上禾二秉,中禾三秉,下禾一秉,实三十四斗;上禾一秉,中禾二秉,下禾三秉,实二十六斗。问上、中、下禾实一秉各几何?答曰:上禾一秉九斗四分斗之一。中禾一秉四斗四分斗之一。下禾一秉二斗四分斗之三。-------《九章算术》
一个两千年前的例子
Basic idea: Add multiples of each row to later rows to make A upper triangular
一个两千年前的例子(2)
Solving linear equations is not trivial.
Forsythe (1952)