4-3矩阵函数.ppt

第三节 矩阵函数 矩阵函数的幂级数定义 说明：在矩阵分析中， 定理 如果AB=BA，则 矩阵函数的Jordan标准形表示—矩阵函数值的求法一 J= diag (J1, J2, …, Js )，且 A=TJT -1 其中 例1 求 例2 求矩阵函数eAx 故矩阵函数为 纯量函数在矩阵谱上的值 设A的最小多项式为 定理 设A是n阶矩阵， p1(?)、 p2(?)是两个多项式，则 p1(A)= p2(A）当且仅当p1(?)、 p2(?)在矩阵A谱上的值一致。即 矩阵函数的多项式定义 设A的最小多项式为 例3、Jordan 块矩阵的函数值 设J是n阶Jordan 块矩阵，主对角元为?0?， f (z)是在?0上有n -1阶导数的函数，求f (J) 由于 所以 矩阵函数的性质 (1) A f (A) = f (A)A (2) (f+g)(A) = f (A)+g(A) (3)若B=T-1 AT，则f (B)=T-1 f (A)T 矩阵函数值的求法二--待定系数法 由矩阵函数的定义，只要函数f (z) 在A的谱上有意义，就可以利用矩阵多项式 p(A)来表示矩阵函数f (A) 。 注 若矩阵A的最小多项式的次数为m，则待定系数多项式的次数最大为m-1，记为 例5 例6 解 A的最小多项式为(z-1)(z-2)2 求解线性方程组得 求 则可设 由 得方程组 解之得： 则 设A的待定系数多项式为 对于 由 * * 主要内容： 1·利用幂级数定义矩阵函数 2·利用多项式定义矩阵函数 2·矩阵函数值的计算方法-- Jordan 形方法和待定系数法 定义：设解析函数 可展开成幂级数，即 设矩阵A的谱半径为ρ(A),幂级数收敛半径为R；则当 ρ(A)R时，对应的矩阵幂级数是绝对收敛的，定义： 例1 下列函数展开成幂级数在整个复平面C上都是收敛的。则对于任何n阶矩阵A，定义矩阵函数为 矩阵指数函数 矩阵三角函数 而对于 其对应的矩阵函数为 其对应的矩阵函数为 矩阵函数恒等式 利用矩阵指数函数、三角函数的幂级数定义式可以证明. 矩阵三角函数恒等式 设AB=BA，则 e A+B= e A e B = e B e A 一般不再成立（P131） 反例 e A+B= e A e B = e B e A 推论： （3）设m为正整数，则 f (A)= Tdiag ( f (J1 ),f (J2 ), …,f (Js ))T -1 设A是n阶矩阵， f (z)是一个可导函数， 则 若A的Jordan 标准形为 其中 设 则 解 矩阵A的特征多项式为 对应于 的特征向量为 而对应于 的线性无关的特征向量只有一个为 则A的Jordan标准形为 从而，矩阵A的标准形及相似变换矩阵为 再利用 ，则对应于 的一个广义特征向 量X3满足： 求得 p(z)是复数域上在?i上有d i -1阶导数的函数， 定义函数p(z)在矩阵A的谱上的值为 说明： A的最小多项式在矩阵A的谱上的值都为0. f (z)是在?i上有d i -1阶导数的函数，p(z)是任意多项式，如果f (z)、 p(z)在矩阵A的谱上的值一致，即 则定义矩阵函数f(A) 为 解 J的最小多项式为?( ?) = (?- ?0?)n ，令 则f (z)、 p(z)在矩阵A谱上的值一致，故 (4)若A= diag (A1,A2, …,As)，则 f (A)=diag ( f (A1 ),f (A2 ), …,f (As )) 说明：利用性质（4）及Jordan 块矩阵的函数值，可求出Jordan 形矩阵的函数值 例4:设 求 解：J为Jordan 形矩阵，即 设A是n阶矩阵， A的特征多项式为 的互异零点为 相应的重数为 则有 设 其中 是次数低于n的 多项式，于是可由 确定出 ，再利用 可得： 因此，n阶矩阵A的函数f（A）最多可用A的m-1次多项式表示。应用 f (z)、 p(z)在矩阵A谱上的值一致，可以求出m个待定系数