线性方程组直接解法
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
3-8
W
Y
Gauss消元法的基本步骤2（4阶）
可以检查，分别以li1乘第一个方程加到第i个方程上 可以完成第一次消元，得同解方程组：
a1(11) x1 a1(12) x2 a1(13) x3 a1(14) x4 b1(1)
a2(22) x2 a2(23) x3 a2(24) x4 b2(2) (3- 4)(a) a3(22) x2 a3(23) x3 a3(24) x4 b3(2)
3-9
W
Y
Gauss消元法的基本步骤3（4阶）
第二步: 消x2 ,首先找到乘数
li 2
a(2) i2
a(2) 22
,i
3,4
程。然后按方程相反顺序求解上三角形方程组，得到原
方程组的解，此过程称为回代过程。
我们的目的，是要总结归纳出一般情况下的n阶线性方程 组的消元公式和回代求解公式，从而得到求解n阶线性方
程组的能顺利在计算机上实现的行之有效的算法。
3-7
W
a1(11) x1 a1G(12)axu2ss消a1元(13) 法x3 的 基a1(14本) x步4 骤b11(（1)4阶）
个方程中,并以
a3(11) a1(11)
,
a4(11) a1(11)
分别
乘第一个方程加到第三 ,第四个方程上消 x1,这些乘数实际上可记为
l21
a2(11) a1(11)
, l31
a3(11) a1(11)
, l41
a4(11) a1(11)
或记为li1
ai(11) a1(11)
百度文库
（i 2,3,4）
a
(1) 43
x3
a4(14) x4
b4(1)
统一加上标 ,并简记为 A(1) x b(1) , A(1) A,b(1) b,首先消元 :
第一步
: 找乘数, 假定a1(11)
0, 要消第二个方程中
x1 , 可以
a2(11) a1(11)
为乘数,
而以
a2(11) a1(11)
乘第一个方程加到第二
a4(22) x2 a4(23) x3 a4(24) x4 b4(2)
变化以后的方 程组系数及右 边的常数项可 总结出如下的 计算公式：
完成第一次消元之后的 方程组记为：
A(2) x = b (2)
ai(j2)
a (1) ij
li1a1(1j)
bi(2) bi(1) li1b1(1)
i, j 2,3,4
12x1 3x2 3x3 15 18x1 3x2 x3 15
(3- 3)
解：消去x1，进行第一次消元：首先找乘数，以 -12乘第一个方程加到第二个方程，以18乘第一个 方程加到第三个方程上可得同解方程组：
x115xx22
x3 6 9x3 57
(3- 3)(a)
21x2 17x3 93
其矩阵形式为： an1x1 an2 x2 ann xn bn
Ax=b
(3-2) 其中： a11 a12 a1n
x1
b1
A
a21 a22 an1an2
a2n ann
，x
x2 ：
xn
，b
b2 ∶
3-3 bn
W
Y
线性方程组的概念（续）
如果线性方程组Ax = b的系数行列式不为零， 即det(A) 0，则该方程组有唯一解。
第三章
线性方程组
直接解法
3-1
W
Y
第三章目录
§1. Gauus 消元法 §2. 主元素法
2.1 引入主元素法的必要性 2.2 列主元素法 2.3 全主元素法 2.4 解三对角方程组的追赶法 §3. 矩阵分解法 3.1 Gauss消去法的矩阵形式 3.2 矩阵的三角分解 3.3 直接三角分解法 §4. 平方根法与改进的平方根法 §5. 矩阵求逆
数，即不可避免地存在着舍入误差的影响，因 而即使是确解法，也只能求到近似解。 直接法在求解中小型线性方程组（≤100个），特别是 系数矩阵为稠密型时，是常用的、非常好的方法。
3-5
W
Y
§1 Gauss消元法
Gauss消元法是最基本的一种方法，下例说明其基本思想:
例1 解线性方程组： x1 x2 x3 6
求解Ax = b，曾经学过高斯（Gauss）消元法， 克莱姆（Cramer）法则，矩阵变换法等，但已远 远满足不了实际运算的需要，主要体现两个方面： 一是运算的快速和准确，其次是方程组的个数增 大时的计算问题。如何建立能在计算机上可以实 现的有效而实用的解法，具有极其重要的意义， 我们也曾指出过，Cramer法则在理论上是绝对正 确的，但当n较大时，在实际计算中却不能用。
§6.方程组的性态和条件数 3-2
W
Y
方程在组科的学求研解究问和题工是程基技线本术性的中，方所常提程见出组的的的，计很概算多问念问题题中如，插线值性函
数，最小二乘数据拟合，构造求解微分方程的差分格式等， 都包含了解线性方程组问题，因此，线性方程组的解法在 数值计算中占有较重要的地位。
设n阶线性方程组：a11x1 a12 x2 a1n xn b1 a2 1x1a22 x2 a2n xn b2 (3-1)
Y
求为aa解32((能1111步))更xx骤11清，楚并aa地32((且1122得))很xx到22容算易法aa地，32((11可33))下xx推33面广以至aa4一阶32((1144))般线xx的性44 n方阶程bb线32((组11性))为方(例3程总-组4结。)
a
(1) 41
x1
a4(12) x2
3-6
W
Y
例1（续）
再消一次元得：
二次消元后将方程化为
倒三角形式，然后进行
x115x2x29 xx33
6 57
22 5
x3
66 5
(3- 3)(b)
回代容易解出： x3 = 3, x2 = 2, x1 = 1。
上述Gauss消元法的基本思想是：先逐次消去变量，将
方程组化成同解的上三角形方程组，此过程称为消元过
3-4
W
Y
线性方程组的数值解法
解线性方程组的数值方法大致分为两类： 1. 直接法：指假设计算过程中不产生含入误差，经过有
限步四则运算可求得方程组准确解的方法。
2. 迭代法：从给定的方程组的一个近似值出发，构造某 种算法逐步将其准确化，一般不能在有限步内得到准确解。
这一章介绍计算机上常用的直接法，它们都是以Gauss 消元法为基本方法，即先将线性方程组化为等价的三角形 方程组，然后求解。 请注意：由于在计算中某些数据实际上只能用有限位小
W
Y
Gauss消元法的基本步骤2（4阶）
可以检查，分别以li1乘第一个方程加到第i个方程上 可以完成第一次消元，得同解方程组：
a1(11) x1 a1(12) x2 a1(13) x3 a1(14) x4 b1(1)
a2(22) x2 a2(23) x3 a2(24) x4 b2(2) (3- 4)(a) a3(22) x2 a3(23) x3 a3(24) x4 b3(2)
3-9
W
Y
Gauss消元法的基本步骤3（4阶）
第二步: 消x2 ,首先找到乘数
li 2
a(2) i2
a(2) 22
,i
3,4
程。然后按方程相反顺序求解上三角形方程组，得到原
方程组的解，此过程称为回代过程。
我们的目的，是要总结归纳出一般情况下的n阶线性方程 组的消元公式和回代求解公式，从而得到求解n阶线性方
程组的能顺利在计算机上实现的行之有效的算法。
3-7
W
a1(11) x1 a1G(12)axu2ss消a1元(13) 法x3 的 基a1(14本) x步4 骤b11(（1)4阶）
个方程中,并以
a3(11) a1(11)
,
a4(11) a1(11)
分别
乘第一个方程加到第三 ,第四个方程上消 x1,这些乘数实际上可记为
l21
a2(11) a1(11)
, l31
a3(11) a1(11)
, l41
a4(11) a1(11)
或记为li1
ai(11) a1(11)
百度文库
（i 2,3,4）
a
(1) 43
x3
a4(14) x4
b4(1)
统一加上标 ,并简记为 A(1) x b(1) , A(1) A,b(1) b,首先消元 :
第一步
: 找乘数, 假定a1(11)
0, 要消第二个方程中
x1 , 可以
a2(11) a1(11)
为乘数,
而以
a2(11) a1(11)
乘第一个方程加到第二
a4(22) x2 a4(23) x3 a4(24) x4 b4(2)
变化以后的方 程组系数及右 边的常数项可 总结出如下的 计算公式：
完成第一次消元之后的 方程组记为：
A(2) x = b (2)
ai(j2)
a (1) ij
li1a1(1j)
bi(2) bi(1) li1b1(1)
i, j 2,3,4
12x1 3x2 3x3 15 18x1 3x2 x3 15
(3- 3)
解：消去x1，进行第一次消元：首先找乘数，以 -12乘第一个方程加到第二个方程，以18乘第一个 方程加到第三个方程上可得同解方程组：
x115xx22
x3 6 9x3 57
(3- 3)(a)
21x2 17x3 93
其矩阵形式为： an1x1 an2 x2 ann xn bn
Ax=b
(3-2) 其中： a11 a12 a1n
x1
b1
A
a21 a22 an1an2
a2n ann
，x
x2 ：
xn
，b
b2 ∶
3-3 bn
W
Y
线性方程组的概念（续）
如果线性方程组Ax = b的系数行列式不为零， 即det(A) 0，则该方程组有唯一解。
第三章
线性方程组
直接解法
3-1
W
Y
第三章目录
§1. Gauus 消元法 §2. 主元素法
2.1 引入主元素法的必要性 2.2 列主元素法 2.3 全主元素法 2.4 解三对角方程组的追赶法 §3. 矩阵分解法 3.1 Gauss消去法的矩阵形式 3.2 矩阵的三角分解 3.3 直接三角分解法 §4. 平方根法与改进的平方根法 §5. 矩阵求逆
数，即不可避免地存在着舍入误差的影响，因 而即使是确解法，也只能求到近似解。 直接法在求解中小型线性方程组（≤100个），特别是 系数矩阵为稠密型时，是常用的、非常好的方法。
3-5
W
Y
§1 Gauss消元法
Gauss消元法是最基本的一种方法，下例说明其基本思想:
例1 解线性方程组： x1 x2 x3 6
求解Ax = b，曾经学过高斯（Gauss）消元法， 克莱姆（Cramer）法则，矩阵变换法等，但已远 远满足不了实际运算的需要，主要体现两个方面： 一是运算的快速和准确，其次是方程组的个数增 大时的计算问题。如何建立能在计算机上可以实 现的有效而实用的解法，具有极其重要的意义， 我们也曾指出过，Cramer法则在理论上是绝对正 确的，但当n较大时，在实际计算中却不能用。
§6.方程组的性态和条件数 3-2
W
Y
方程在组科的学求研解究问和题工是程基技线本术性的中，方所常提程见出组的的的，计很概算多问念问题题中如，插线值性函
数，最小二乘数据拟合，构造求解微分方程的差分格式等， 都包含了解线性方程组问题，因此，线性方程组的解法在 数值计算中占有较重要的地位。
设n阶线性方程组：a11x1 a12 x2 a1n xn b1 a2 1x1a22 x2 a2n xn b2 (3-1)
Y
求为aa解32((能1111步))更xx骤11清，楚并aa地32((且1122得))很xx到22容算易法aa地，32((11可33))下xx推33面广以至aa4一阶32((1144))般线xx的性44 n方阶程bb线32((组11性))为方(例3程总-组4结。)
a
(1) 41
x1
a4(12) x2
3-6
W
Y
例1（续）
再消一次元得：
二次消元后将方程化为
倒三角形式，然后进行
x115x2x29 xx33
6 57
22 5
x3
66 5
(3- 3)(b)
回代容易解出： x3 = 3, x2 = 2, x1 = 1。
上述Gauss消元法的基本思想是：先逐次消去变量，将
方程组化成同解的上三角形方程组，此过程称为消元过
3-4
W
Y
线性方程组的数值解法
解线性方程组的数值方法大致分为两类： 1. 直接法：指假设计算过程中不产生含入误差，经过有
限步四则运算可求得方程组准确解的方法。
2. 迭代法：从给定的方程组的一个近似值出发，构造某 种算法逐步将其准确化，一般不能在有限步内得到准确解。
这一章介绍计算机上常用的直接法，它们都是以Gauss 消元法为基本方法，即先将线性方程组化为等价的三角形 方程组，然后求解。 请注意：由于在计算中某些数据实际上只能用有限位小