实习线性方程组AXB的数值解法

《实习线性方程组AXB的数值解法》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第3章线性方程组AX=B的数值解法7/18/2021华南师范大学数学科学学院谢骊玲引言在自然科学和工程技术中很多问题的解决常常归结为解线性代数方程组。例如电学中的网络问题，船体数学放样中建立三次样条函数问题，用最小二乘法求实验数据的曲线拟合问题，解非线性方程组问题，用差分法或者有限元法解常微分方程，偏微分方程边值问题等都导致求解线性方程组，而且后面几种情况常常归结为求解大型线性方程组。线性代数方面的计算方法就是研究求解线性方程组的一些数值解法与研究计算矩阵的特征值及特征向量的数值方法。7/18/2021华南师范大学数学科学学院谢骊玲线性方程组求解问

2、题考虑线性方程组Ax=b其中A是一个(n×n)的非奇异矩阵,x是要求解的n维未知向量,b是n维常向量7/18/2021华南师范大学数学科学学院谢骊玲线性方程组的解的存在性和唯一性定理3.4设A是N×N方阵，下列命题等价：给定任意N×1矩阵B，线性方程组AX=B有唯一解矩阵A是非奇异的（即A-1存在）方程组AX=0有唯一解X=0det(A)≠07/18/2021华南师范大学数学科学学院谢骊玲线性方程组的解最常见的求线性方程组Ax=b的解的方法是在方程组两侧同乘以矩阵A的逆Gram法则：Ax=b7/18/2021华南师范大学数学科学学院谢骊玲线性方程组

3、的解（续1）求逆运算和行列式计算由于运算量大，实际求解过程中基本不使用，仅作为理论上的定性讨论克莱姆法则在理论上有着重大意义，但在实际应用中存在很大的困难，在线性代数中，为解决这一困难给出了高斯消元法还有三角分解法和迭代求解法7/18/2021华南师范大学数学科学学院谢骊玲解法分类关于线性方程组的数值解法一般有两类直接法：若在计算过程中没有舍入误差，经过有限步算术运算，可求得方程组的精确解的方法迭代法：用某种极限过程去逐步逼近线性方程组精确解的方法迭代法具有占存储单元少，程序设计简单，原始系数矩阵在迭代过程中不变等优点，但存在收敛性及收敛速度等问题

4、7/18/2021华南师范大学数学科学学院谢骊玲3.3上三角线性方程组定义3.2N×N矩阵A=[aij]中的元素满足对所有i>j，有aij=0，则称矩阵A为上三角矩阵；如果A中的元素满足对所有i

5、件不满足，则可能无解或有无穷解联系定理3.4，可知要条件akk≠0成立才能保证方程组存在唯一解7/18/2021华南师范大学数学科学学院谢骊玲3.3上三角线性方程组（续2）求解上三角线性方程组的回代算法最后7/18/2021华南师范大学数学科学学院谢骊玲上三角线性方程组的求解基本算法：7/18/2021华南师范大学数学科学学院谢骊玲上三角线性方程组的求解（续1）7/18/2021华南师范大学数学科学学院谢骊玲3.4高斯消去法和选主元求解有N个方程和N个未知数的一般方程组AX=B的一般做法：构造一个等价的上三角方程组UX=Y，并利用回代法求解如果两个

6、N×N线性方程组的解相同，则称二者等价对一个给定方程组进行初等变换，不会改变它的解7/18/2021华南师范大学数学科学学院谢骊玲3.4高斯消去法和选主元（续1）考虑一个简单的例子：求解第二个方程，得第二个方程减去第一个方程除以3再乘以4得到的新方程，得到新的方程组：回代到第一个方程，得7/18/2021华南师范大学数学科学学院谢骊玲3.4高斯消去法和选主元（续2）考虑包含n个未知数的方程组or作如下行变换之后方程组的解向量x不变对调方程组的两行用非零常数乘以方程组的某一行将方程组的某一行乘以一个非零常数，再加到另一行上通过对增广矩阵[A

7、B]进行

8、如上的行变换求解7/18/2021华南师范大学数学科学学院谢骊玲3.4高斯消去法和选主元（续3）7/18/2021华南师范大学数学科学学院谢骊玲3.4高斯消去法和选主元（续4）7/18/2021华南师范大学数学科学学院谢骊玲3.4高斯消去法和选主元（续5）7/18/2021华南师范大学数学科学学院谢骊玲3.4高斯消去法和选主元（续6）利用3.3节的回代法求解上述上三角方程组7/18/2021华南师范大学数学科学学院谢骊玲3.4高斯消去法和选主元（续7）消去过程7/18/2021华南师范大学数学科学学院谢骊玲3.4高斯消去法和选主元（续8）回代过程7

9、/18/2021华南师范大学数学科学学院谢骊玲3.4高斯消去法和选主元（续9）上述消去过程中，如果akk=0，则不能使用第