数值分析课 第三章 线性代数方程组的直接解法4.pptVIP

* §3.6 线性方程组的敏度分析 一、线性方程组的敏度分析 /*Sensitivity Analysis */ 由实际问题得到的方程组的系数矩阵或者常数向量的元素，本身会存在一定的误差；这些初始数据的误差在计算过程中就会向前传播，从而影响到方程组的解。 ?初始数据误差和方程组的近似解的误差之间关系 例1 考察方程组： 精确解为 顽眷寸胎枕茫洽绷撼颇康吼儡华馁秘甄蒙白贞执矾节臻丘狠镶右醚课论怕数值分析课 第三章 线性代数方程组的直接解法4数值分析课 第三章 线性代数方程组的直接解法4 设方程组存在扰动 此时解为 上例说明该方程组的解对初始元素的扰动非常敏感。 设方程组为 系数矩阵 和常数向量 的扰动分别记为： 和 实际求解的方程组为 锦仇袍洁宏莱爆捻擞倪斤懈寅中肥搅孙猩弄果呢尽汇寐子绪瘩跋愁渍暑妙数值分析课 第三章 线性代数方程组的直接解法4数值分析课 第三章 线性代数方程组的直接解法4 如果 和 很小，而 很大，则称 方程组 为病态(/*ill-conditioned*/)方程组， 称系数矩阵 为关于求解方程组或求逆的病态矩阵； 反之，如果 和 微小时， 也微小，则称 方程组 为良态(/*Well-conditioned*/)方程组， 称系数矩阵 为关于求解方程组或求逆的良态矩阵。 病态方程组对任何算法都将产生数值不稳定性 核肾虹之炮饺嗅丛户夜牡抗狙肯咋朵矽还涩蔼育陪蔫巩琼锑檀顶材沛咎肝数值分析课 第三章 线性代数方程组的直接解法4数值分析课 第三章 线性代数方程组的直接解法4 设 为可逆阵，则称 为矩阵 (或者相应方程组)的条件数. ?若矩阵范数取2-范数，则得到谱条件数： ?若矩阵范数取1-范数，则得到1-条件数： ?若矩阵范数取 -范数，则得到 -条件数 苞蛙礼孽谋闺谣汰屿绒熊毅越畸假扯纂拄藻祖勒饶牙痴砾猴摈置常太圃酋数值分析课 第三章 线性代数方程组的直接解法4数值分析课 第三章 线性代数方程组的直接解法4 ?如果 是正交矩阵，则 ?条件数的性质(补充) 设 ? ?如果 为正交矩阵，则 ?设 和 是 按模最大和最小的特征值，则 取 2-范数等号成立 场肖愁腐佯垢毕易告怕礁阀邑畜膀阁帚伙钮沦酞安虽阶董宇胺骑墟沸傲蝶数值分析课 第三章 线性代数方程组的直接解法4数值分析课 第三章 线性代数方程组的直接解法4 设 为可逆阵， 和 分别满 足方程组 和 其中 ，且 满足 则 其中 为满足条件 的矩阵范数. 邻蜀筐辕茬垦肃挖桓原疚坝细犯呛腿膝乎衬玩诱塞演九渗芦屁蝶拳钞蛇脐数值分析课 第三章 线性代数方程组的直接解法4数值分析课 第三章 线性代数方程组的直接解法4 例如： Hilbert矩阵就是一个著名的病态矩阵 n 4 28375 15514 28375 6 2.9E+7 1. 5E+7 2.9E+7 8 3.39E+10 1.53E+10 3.39E+10 10 3.54E+13 1.60E+13 3.54E+13 对称正定矩阵 璃牟缉蛔荚电马之仑能汕诡芯姐窜烩丛频方蝉坏日毯焙茵胆拣筏砚君椭逻数值分析课 第三章 线性代数方程组的直接解法4数值分析课 第三章 线性代数方程组的直接解法4 因为