第3章线性代数方程组直接解法_第1页
第3章线性代数方程组直接解法_第2页
第3章线性代数方程组直接解法_第3页
第3章线性代数方程组直接解法_第4页
第3章线性代数方程组直接解法_第5页
已阅读5页，还剩13页未读， 继续免费阅读

下载本文档

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

文档简介

§5直接法的误差分析/*ErrorAnalysisofDirectMethod*/一、扰动方程组的误差界/*ErrorBoundofPerturbedSystems*/

由实际问题得到的方程组的系数矩阵或者常数向量的元素，本身会存在一定的误差；这些初始数据的误差在计算过程中就会向前传播，从而影响到方程组的解。初始数据误差和方程组的近似解的误差之间关系。例12考察方程组：精确解为设方程组存在扰动精确解为上例说明该方程组的解对初始元素的扰动非常敏感。设方程组为系数矩阵和常数向量的扰动分别记为：和实际求解的方程组为

如果和很小，而很大，则称方程组为病态(/*ill-conditioned*/)方程组，称系数矩阵为关于求解方程组或求逆的病态矩阵；反之，如果和微小时，也微小，则称方程组为良态(/*Well-conditioned*/)方程组，称系数矩阵为关于求解方程组或求逆的良态矩阵。病态方程组对任何算法都将产生数值不稳定性

设为可逆阵，为一种从属矩阵范数，则称为矩阵的条件数。

若矩阵范数取2-范数，则得到谱条件数：

若矩阵范数取1-范数，则得到条件数：

若矩阵范数取

-范数，则得到条件数：例如：Hilbert矩阵就是一个著名的病态矩阵n46810

283752.9E+73.39E+103.54E+13155141.5E+71.53E+101.60E+13283752.9E+73.39E+103.54E+13对称正定矩阵

如果是正交矩阵，则

条件数的性质（）设

正交矩阵：

如果为正交矩阵，则Lemma设，在正交变换下，谱范数和F-范数保持不变，即设则有：证明：（1）同理可证（2）对设

设和是按模最大和最小的特征值，则取2-范数等号成立

设为可逆阵，和分别满足方程组和其中，且适当小，使得则其中所用的是任意向量范数和从属于它的矩阵范数证明：设，如果则可逆，且且由条件所以可逆，且（摄动定理）方程组有唯一解：Corollary(推论)

如果，则

如果，则在上述定理的条件下，仅由或者引起的解的相对误差界可大致认为是原始扰动的相对误差的倍数，倍数近似等于条件数。

如果，则证明：结论得证实际相对误差界与条件数的关系设是方程组的精确解，是近似解，则有，称之为剩余（或残余）向量证明：定理给出了相对误差的估计范围（与条件数有关）结论得证二、病态方程组的解法

常用的几种判定方程组为病态的经验方法

当相对来说很小（大），或者矩阵的某些行（列）近似线性相关时，可能为病态；

如果采用选主元消去法求解方程组时，在消元过程中出现很小的主元，可能为病态；

解方程组得到了一个很大的解，或者特征值相差大数量级，可能为病态；

当系数矩阵的元素间数量级相差很大，且无一定规则时，可能为病态。n较大时，条件数的计算困难。

求解病态方程

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

温馨提示

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

最新文档

评论

0/150

提交评论