第六篇 章 矩阵分析及其应用 矩阵理论课件.ppt

第六章 矩阵分析及其应用;虽然在微积分开端时期贝克莱将无穷小称为“上帝的幽灵”，进而导致“第二次数学危机”，直到柯西的“极限论”和戴德金等的“实数理论”的出现危机才算彻底解决。但微积分在近代社会的巨大作用我们早已深有体会，将微积分中的极限、导数、积分、级数等分析思想和方法应用于矩阵的研究，自然就在情理之中。;§1、矩阵序列与矩阵级数;一、矩阵序列的收敛性;根据矩阵序列收敛性的定义，可证明下列性质。;用矩阵的范数理论来研究矩阵序列的收敛性是最常用、最简洁的方法。;证明：;所以;由于向量是特殊的矩阵，因此我们有;联想到等比数列 收敛当且仅当 ，类似地，我们有;证明：;这里规定 时，;由于谱半径不易计算，联系到谱半径不超过任何一种矩阵范数，实际常用范数来判断矩阵是否是收敛矩阵。只有很难找到这样的范数，才计算出矩阵的所有特征值，进而得到谱半径。;二、矩阵级数;这是因为;定义9 中的矩阵级数 称为绝对收敛的，如果数项级数;定理10 中的矩阵级数 绝对收敛的充要条件是正项级数 收敛，这里矩阵范数是任意的。;证明：必要性。;所以正项级数 收敛。;证明：充分性。;定义11 中的矩阵级数;定理12 设幂级数的收敛半径为 ，则;证明：;其中;绝对收敛，故矩阵幂级数 绝对收敛。;最后讨论最特殊的诺伊曼(Neumann)级数,即;定理13 上的诺伊曼(Neumann)级数收敛的充要条件是 。并且诺伊曼(Neumann)级数收敛于;定理14 对 上满足 的相容矩阵范数 。如果 ，则有误差估计式;定理14的证明需要用到上一章的引理6，即：;证明：;两边取范数，并利用引理6，得;例 15 判断方阵幂级数;解：方阵 的谱半径满足;§2、矩阵函数及其计算;一、矩阵函数的定义及性质;在高等数学和复变函数中，有幂级数展开式：;相应地，我们有矩阵函数：;以及含参数的矩阵函数：;根据欧拉公式 ，可以推出：;遗憾的是，指数运算规则一般不成立：;可以验证 确实两两不等。;那么什么条件下指数运算规则成立呢？;推论 设 ，则;二、矩阵函数的计算;Jordan分解法;其中;特征值分解法;例 3 求矩阵函数 、 和 ，其中;解： 求得 的Jordan分解为;当 时 ，则;当 时;%ex601.m A=[-1 1 0 ;-4 3 0; 1 0 2]; expm(A) %调用expm函数 %expm uses the Padé approximation with scaling % and squaring.;%ex601.m（续） A=[-1 1 0 ;-4 3 0; 1 0 2]; syms t % 声明符号变量t S=expm(A*t)； S=simple(S) % 简化矩阵函数的结果;当 时;%ex601.m（续） A=[-1 1 0 ;-4 3 0; 1 0 2]; syms t % 声明符号变量t A1=sin(A*t) % 内置函数sin（A）给出错误结果;%ex601.m（续） A=[-1 1 0 ;-4 3 0; 1 0 2]; j=sqrt(-1); %虚数单位i A2=(expm(j*A*t) - expm((-1)*j*A*t))/(2*j) %利用Euler公式，调用函数expm A2=simple(A2) % 简化矩阵函数的结果 ;例 4 求矩阵函数 和 ，其中;解： 矩阵 的特征值为; 从而;%ex602.m A=[4 6 0;-3 -5 0;-3 -6 1]; syms t % 声明符号变量t S=expm(A*t)； S=simple(S) % 简化矩阵函数的结果;%ex602.m （续） A=[4 6 0;-3 -5 0;-3 -6 1]; j=