数值分析课件ch_第1页
数值分析课件ch_第2页
数值分析课件ch_第3页
数值分析课件ch_第4页
免费预览已结束，剩余7页可下载查看

付费下载

下载本文档

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

文档简介

1、PAGE PAGE 11第二章 解线性方程组的直接法解线性方程组或写成矩阵式 其中 Gauss消去法（矩阵行变换法）第k次消元公式计算中，中间结果不必保留，进行一次变换后原来存放的单元存放，的单元存放。因此，我们得到Gauss消去法的算法：循环：n-1何时可行？即第k步 Gauss消去法可实行，易见充要条件是 若A的各阶顺序主子式 ，则有： 消元过程可进行到 。因此，可以用Gauss消去法解线性方程组的充要条件是系数矩阵的各阶顺序主子式不为0。最后得到是上三角阵与同解 解只需递推（回代过程）计算量 第k步消元计算用(n-k)次除法，算诸用乘法和次加减法， 对相加，可得消元过程共需 右端 矩阵的

2、三角分解（用矩阵乘法分解的观点看Gauss消去法）对A作行变换相当于左乘初等矩阵，例如其中类似的讨论易知：令 单位下三角阵定理：，则A可表示为A=LU L：单位下三角阵，U上三角阵，且分解唯一。证明：存在性上已证，下证唯一性，首先设A非奇异，若A有两种分解法：，则上式左边为单位下三角阵，右边为上三角阵，从而只能是单位阵，可得：.分解唯一。若A奇异，A的秩为n-1，将A写成分块矩阵如果A=LU分解好了，则解Ax=b相当于解Ly=b，递推即可（回代过程）L、U放在一个矩阵单元中，单元放Gauss主元素消去法Gauss消去法的缺点 遇到无法进行，很小，可能溢出，误差大解决办法： 进行行列交换，避免，

3、也不使其太小完全主元消去法思想：在做第k步消元时，交换行列，取中绝对值最大者为优点：计算精度高，稳定性好缺点：比较工作量大，条件转多，通用性小列主元消去法在做第k步消元时，交换行，取中，绝对值最大值为循环：k=1,2,n-11)按k列选主元，如果转3）2)主元行行与k行交换， j=k,n 3)消元过程 对固定的k同Gauss消去法回代（同Gauss消去法）列主元是最常用的方法Gauss-Jordan消去法（无回代过程的Gauss消去法）思想：第k次消元将第k列元素除主对角线上的其他元（对角线上、下）均化成0算法：1.循环 k= 1,n2. 3. 输出 无需回代，但消元过程计算量增加总的计算量（

4、乘除法）达，仅用于求逆矩阵线性代数中的行初等变换法求逆矩阵，A为n阶可逆距阵可选到列主元，也可不选（顺序主子式），但不能用全面主元下面用矩阵运算来描述列主元消去法用表示单位阵交换i行j行所得到显然所以以n=4为例（Gauss列主元消去法）由于例如相当于把中与单位阵相同的两列交换，再交换两行，仍为一单位下三角阵仅第二列与单位阵不同现在令，与L均为单位下三角阵置换方阵，I通过一系列交换行得到，每列恰有一个1.（n-1）个零定理：若存在置换阵P，使P不唯一，分解也不唯一（可能选，部分选，不选主元P=I）Gauss主元消去法， 平方根法有限元法解结构问题时，常要解Ax=b，其中A对称正定（常见的情况） 由唯一性 下三角矩阵定理（Cholesky分解）设A为对称正定，则存在下三角阵L使，若限定L的对角元0,则L是唯一的。给定

人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

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

最新文档

评论

0/150

提交评论