ch03-线性方程组的数值解法-迭代法省公开课一等奖全国示范课微课金奖课件.pptx

第三章线性方程组数值解——迭代法第1页

3.4向量与矩阵范数求解线性方程组误差是不可防止，因：Ax=b中，A、b往往是前面计算结果误差分析中：e=|x-x*|ε来度量真值与近似值之间误差？？怎样预计一个向量（方程组解向量）、矩阵误差？？？怎样预计向量、矩阵大小、“距离”？？第2页

向量范数定义：Rn空间向量x范数||·||是一个实数，且满足非负性：齐次性：三角不等式：则称||·||是向量x一个范数第3页

性质：‖x‖≠0时，‖-x‖=‖x‖｜‖x‖-‖y‖｜≤‖x-y‖向量范数‖x‖是Rn上向量x连续函数.第4页

常见向量范数：1-范数：2-范数：∞-范数：p-范数：示例：求x=(1，2，3)T三种范数第5页

定理：向量x，有：范数等价：设‖·‖a和‖·‖b是Rn上任意两种范数，若存在常数C1、C20，使得C1‖x‖a≤‖x‖b｜≤C2‖x‖a，则称‖·‖a和‖·‖b等价。||x||∞≤||x||2≤||x||1≤n1/2||x||2≤n||x||∞推论：Rn上一切范数都等价第6页

向量间距离：||x-y||为向量x、y间距离，故对线性方程组近似解x*与准确解x间误差可用二者间差距（距离）（即范数||x-x*||）进行刻画向量收敛性：Rn中向量序列{x(k)}，即：x(0)、x(1)、…x(k),……，其中x(k)=（x1(k),x2(k)x3(k)….xn(k))T，对全部分量xi(k)，都有：则第7页

定理：对任意一个向量范数‖·‖而言，向量序列｛x(k)｝收敛于向量x*充要条件是因为向量各种范数是等价，故只要向量在一个范数下收敛，则另一个范数依然收敛。针对不一样问题可选不一样范数进行讨论第8页

矩阵范数定义：对任意A∈Rn?n，||A||表示按照一定法则确定一个实数，且满足：非负性：齐次性：三角不等式：则称||A||是矩阵A范数如：定义||A||=，满足上述三个条件，故也可作为矩阵范数第9页

矩阵算子范数：A∈Rn?n,x∈Rn其中，||A||是全部非0向量x中，使得比值取得最大或最小上界。则||A||为向量范数导出矩阵算子范数算子范数满足相容性条件：‖Ax‖≤‖A‖‖x‖等价定义：全部非0且其范数为1向量x集合中，||Ax||最大值即为矩阵算子范数第10页

由向量范数||·||p，诱导出关于矩阵A?Rn?n范数称为矩阵p-范数惯用矩阵算子范数（行和范数）（列和范数）（谱范数）表示最大特征值第11页

例：求矩阵3种范数第12页

矩阵谱半径：对于Rn×n上矩阵A，设其特征值为λ1,λ2,…λn，称即：特征值模最大值为矩阵谱半径性质：Rn×n上矩阵A，则：ρ(λ)≤||A||（||A||为任一个范数)若矩阵A对某个算子范数满足||A||1，则必有:I±A可逆第13页

1)证实：证实：设λ是矩阵A任一特征值，其对应特征向量为x,则有Ax=λx，故||λx||=|λ|||x||=||Ax||≤||A||||x||所以任意|λ|≤||A||，故ρ(λ)≤||A||2）证I±A可逆，用反证法若I±A不可逆，则|I±A|=0，故方程组（I±A)x=0有非零解向量x，即存在x≠0(向量），使得：Ax=±x成立所以：||x||=||Ax||≤||A||||x||,又因为x非0向量，故||x||≠0所以||A||≥1（矛盾）第14页

3)证实：因：(I±A)(I±A)-1=(I±A)-1±A(I±A)-1=I所以：(I±A)-1=I±A(I±A)-1故：||(I±A)-1||=||I±A(I±A)-1||≤||I||+||A(I±A)-1||≤1+||A||||(I±A)-1||故：第15页

3.6线性方程组数值解误差分析求解Ax=b时,A和b误差对解x有何影响？常向量b有δb扰动，b?b+δb?方程近似解为x+δx即：即：近似解相对误差是常向量相对误差||A||||A-1||倍条件数第16页

常向量不变，而系数矩阵A有δA扰动时:条件数第17页

结论：当||A||||A-1||较小时，方程组常向量、系数矩阵扰动，不会引发解向量巨大改变，此时称该方程在是良态反之，当||A||||A-1||较大时，b、A微小扰动，会引发解向量巨大改变（扩大||A||||A-1||倍），则该方程组是病态，对应该系数矩阵也是病态条件数