a(k ik
1)
) /a a(k)
ij
(k) kk
(i lik
k a(k)
kj
(i
1,...,n) k 1,...,n;
j
k
1,...,n, n
1)
该Gauss消去法为顺序高斯消去法
第7页/共87页
Gauss
for k 1, 2, , n 1
for i k 1, k 2, , n
Cramer法则:
xi
Di D
i 1, 2,
,n
所需乘除法的运算量大约为(n+1)!+n
n=20时,每秒1亿次运算速度的计算机要算30多万年!
直接法
在没有舍入误差的情况下,经过有限次 运算可以得到方程组的精确解的方法。
第2页/共87页
§3.1 Gauss消去与矩阵LU分解
属于解方程的直接法
一 Gauss消去 1 直接法的关键思想
ln,k
1
第26页/共87页
A L1L2 Ln2 Ln1U LU L为单位下三角
1
l21 1 l31 l32 1
L l41 l42 l43 1
u11 u12 ... u1n
U
u22 ... u2n ...
1 ln1 ln2 ln3 lnk lnn1
unn
A LU 矩阵分解为单位下三角 和上三角矩阵的乘积
aii
第13页/共87页
例:在8位制计算机上解方程组
109
x1
x2
1
x1 x2 2
要求用Gauss消去法计算。
解:l21 a21 / a11 109 8个
x1 x2 1
a22 1 l21 1 0.0 ...01109 109 109