计算方法第5章1、解线性方程组的直接方法(5.15.4)_第1页
计算方法第5章1、解线性方程组的直接方法(5.15.4)_第2页
计算方法第5章1、解线性方程组的直接方法(5.15.4)_第3页
计算方法第5章1、解线性方程组的直接方法(5.15.4)_第4页
计算方法第5章1、解线性方程组的直接方法(5.15.4)_第5页
已阅读5页,还剩35页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

计算方法解线性方程组的直接方法1/4/20241第五章

解线性方程组的直接方法

§5.1引言解线性方程组的两类方法:直接法:经过有限次运算后可求得方程组精确解的方法(不计舍入误差)迭代法:从解的某个近似值出发,通过构造一个无穷序列去逼近精确解的方法。(一般有限步内得不到精确解)1/4/20242n阶线性方程组1/4/202431/4/20244§5.2高斯消去法转化为同解的方程组1/4/20245§5.2.1高斯消去法计算过程1/4/202461/4/202471/4/202481/4/202491/4/202410系数矩阵与常数项:1/4/202411回代过程:1/4/2024121/4/202413消去第一列的n-1个系数要计算n*(n-1)

个乘法。§5.2.2高斯消去法计算量1/4/202414每一步消去过程相当于左乘初等变换矩阵Lk§5.2.3矩阵的三角分解1/4/2024151/4/202416i+1行

i+1行依次递推1/4/202417定理7(矩阵的LU分解)设A为n阶矩阵,如果A的顺序主子式

Di≠0(i=1,2,…,n-1),则A可分解为一个单位下三角矩阵L和一个上三角矩阵U的乘积,且这种分解是

唯一的。1/4/2024181/4/202419§5.3高斯主元素消去法为避免此种情况的发生,可通过交换方程的次序,选取绝对值大的元素作主元。§5.3.1

列主元素消去法1/4/202420选取或称此方法为全主元素高斯消去法称此方法为列主元素高斯消去法1/4/20242131/4/2024221/4/202423定理8(列主元素的三角分解定理)如果A为非奇异矩阵,则存在排列矩阵P使

PA=LU

其中L为单位下三角阵,U为上三角阵。1/4/202424§5.3.2高斯—若当消去法1/4/2024251/4/202426本章作业

P17671/4/202427§5.4矩阵三角分解法

§5.4.1直接三角分解法将高斯消去法改写为紧凑形式,可以直接从矩阵A的元素得到计算L,U元素的递推公式,而不需要任何中间步骤,这就是直接三角分解法。由于A=LU,求解Ax=b的问题就等价于求解两个三角形方程组①Ly=b,求y;②

Ux=y,求x.1/4/2024281、不选主元的三角分解法A=LU其中L为单位下三角阵,U为上三角阵(4.1)1/4/202429一、直接计算A的LU分解(例)1/4/2024301/4/202431二、一般计算公式1/4/202432三、LU分解求解线性方程组1/4/202433矩阵A的直接分解法称为杜利特尔(Doolittle)分解1/4/202434例1:将方程组的系数矩阵A作LU分解,并求方程组的解1/4/202435解LU分解的紧凑格式为1/4/202436推出:1/4/202437由Ux=y,即用回代法解得即为线性方程组的解1/4/2024382、选主元的三角分解法采用与列主

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论