工程矩阵理论(5章-范数及矩阵函数).pptVIP

请双面打印/复印(节约纸张) 272365083@ 主讲: 张小向 第五章 范数及矩阵函数 第五章 范数及矩阵函数 §5.4 矩阵函数 例1 已知A = , B = , 求eA, eB, eA+B. 1 0 0 0 0 1 0 0 解: 注意到A2 = A, 因而对于任意的正整数k, 有Ak = A. 于是 eA = I + A + … + + … A m! = + ? 1 0 0 1 ? m=1 1/m! 0 0 0 = . e 0 0 1 第五章 范数及矩阵函数 §5.4 矩阵函数 注意对于任意的正整数k 1, 有Bk = O. 于是 eB = I + B = + 1 0 0 1 0 1 0 0 = . 1 1 0 1 最后, 对于任意的正整数k, 有 故 (A+B)k = A+B = , 1 1 0 0 eA+B = I + (A+B) … + + … A+B m! = + ? 1 0 0 1 ? m=1 1/m! 1/m! 0 0 = . e e?1 0 1 第五章 范数及矩阵函数 §5.4 矩阵函数 注: eA = , e 0 0 1 在上面的例子中, eB = , 1 1 0 1 eA+B = . e e?1 0 1 eAeB = e 0 0 1 1 1 0 1 = . e e 0 1 eBeA = 1 1 0 1 e 0 0 1 = . e 1 0 1 可见, 对于矩阵指数, eA+B, eAeB, eBeA三者未必相等. 第五章 范数及矩阵函数 §5.4 矩阵函数 二. 利用Jordan标准形计算矩阵函数 定理5.4.1 若A ? n?n, 且?(A) R, 则 f(A) = Pdiag[ f(J1) , …, f(Js)]P?1, 设 f(x)在{x ? | |x| R}内解析. P?1AP = J = diag(J1, …, Js), 其中 ri?ri Ji = , ?i 1 0 … 0 0 ?i 1 … 0 0 0 0 … 1 0 0 0 … ?i … … … … 第五章 范数及矩阵函数 §5.4 矩阵函数 其中f(Ji) = ri?ri f(?i) … … … … f(?i) f(?i) f ?(?i) f ?(?i) f ( )(?i) (ri?1)! ri?1 f ( )(?i) (ri?2)! ri?2 … 0 0 0 0 0 … … , i = 1, 2, …, s, 第五章 范数及矩阵函数 §5.4 矩阵函数 证明: 设 f(x) = ? amxm, |x| R. ? m=0 记 ? amxm = Sk(x), m=0 k = Pdiag[Sk(J1), …, Sk(Js)]P?1. ? amAm = Sk(A) = PSk(J)P?1 m=0 k 则Sk(x) ? f(x) (k?+?), 而且 第五章 范数及矩阵函数 §5.3 两个收敛定理 因为Sk(Ji) = ri?ri Sk(?i) … … … … Sk(?i) Sk(?i) Sk?(?i) Sk?(?i) Sk( )(?i) (ri?1)! ri?1 Sk( )(?i) (ri?2)! ri?2 … 0 0 0 0 0 … … , 其中i = 1, 2, …, s, 而且幂级数在收敛域内可以逐项求导, 所以当?(A) R时, {Sk(Ji)}? f(Ji) (k?+?). 对Sk(A) = Pdiag[Sk(J1), …, Sk(Js)]P?1两边取极限… 第五章 范数及