[工学]研究生矩阵理论课后答案第6-7章.pptVIP

第六章 矩阵函数 ?矩阵函数一般定义:矩阵函数是从Cm?n到Cu?v的一个对应规则f:Cm?n?Cu?v,使对每个x?Cm?n,都对应于唯一的f(x)?Cu?v. 例如:det:Cn?n?C1?1,?A?Cn?n,det(A)?C1?1; f:Cn?n?Cn?n,?A?Cn?n,f(A)=2A-E?Cn?n. ?矩阵函数的概念十分广泛,其应用也相应地十分广泛. ?我们仅限于讨论从Cn?n到自身的函数 f:Cn?n?Cn?n. 特别更限于最简单的矩阵多项式函数和由矩阵幂级数定义的矩阵函数. 矩阵函数举例 例1:我们知道幂级数ex=?k=0?(1/k!)xk的收敛半径R=?,故可以定义Cn?n到自身的矩阵函数: ?A?Cn?n, eA = ?k=0?(1/k!)Ak (绝对收敛) 例2:ln(1+x)=?k=1?((-1)k/k)xk,收敛半径R=1,可以定义Cn?n的单位球‖A‖1内的一个矩阵函数: ?A?Cn?n,ln(E+A)=?k=1?((-1)k/k)Ak (绝对收敛) 例3:对任意多项式p(x)=?k=0makxk,可以定义Cn?n的一个矩阵函数: p:Cn?n?Cn?n,?A?Cn?n,p(A)=?k=0makAk.(A0=E) 例如,p(x)=x2-2,p(A)=A2-2E. 当A= 时,有 p(A)= 矩阵多项式 定义:对任意复系数多项式p(x)=?k=0makxk,可以定义矩阵函数: p:Cn?n?Cn?n,?A?Cn?n,p(A)=?k=0makAk. 并称它为A的一个矩阵多项式. 由定义不难直接推出矩阵多项式的下列性质: 矩阵多项式的性质 ①:与A(同阶)可交换性 p(A)A=Ap(A). 证:不难看出上式两边都等于?k=0makAk+1. ②:经相似变换 B=PAP-1 ? p(B)=Pp(A)P-1 证:p(B)=?k=0mak(PAP-1)k =?k=0makPAkP-1=P(?k=0makAk)P-1=Pp(A)P-1 ③:直和 若A=diag(J1,…Jr), 每个Ji为方阵,则 p(A)=?k=0makAk=diag(?k=0makJ1k,…,?k=0makJrk) =diag(p(J1),…,p(Jr)). ④:p(A)的Jordan表示：A=Pdiag(J1,…,Jr)P-1 ? p(A)=Pdiag(p(J1),…,p(Jr))P-1 diag(J1,…,Jr)为A的Jordan标准形. 矩阵多项式的Jordan表示应用举例 ?定理6.1.1:对任意方阵A和多项式p(x), σ(A)={λ1,…,λn} ? σ(p(A))={p(λ1),…,p(λn)} A?=?i?, ??0 ? p(A)?=p(λi)? 证: 由p(A)的Jordan表示直接推出,因p(A)的特 征值集等于σ(p(J1))∪…∪σ(p(Jr)), 并且p(Ji)是对角元为p(?i)上三角矩阵. p(A)?=?k=0makAk? =?k=0mak?ik? =p(?i)? 对任意方阵A和多项式p(x), σ(A)={λ1,…,λn} ? σ(p(A))={p(λ1),…,p(λn)} 例6.1.1:A=Pdiag(J1,J2)P-1, J1= , J2=(2), A= p(J1)= ,p(J2)=(p(2)). p(A)=Pdiag(p(J1), p(J2))P-1 (p(x)可以是任意多 = P P-1 项式,如?4?2?3+??1) 不变因子 ?定理2.1.4(p.76):任何满秩的n阶?-方阵A(λ)都等价于一个对角矩阵:(P,Q为满秩?-方阵) PAQ=diag(d1(λ),…,dn(λ)), d1(λ),…,dn(λ)为?的首一多项式,满足 di(λ)|di+1(λ), i=1,…,n 称为 A(λ)的不变因子. 特别,A(λ)=?E-A的不变因子称为A的不变因子. 行列式因子 ?定义2.1.7(p.82):满秩的n阶?-方阵A(λ)的全部k阶子式的首一的最大公因式称为A(λ)的k阶行列式因子,记为:Dk(λ). ?定理2.1.6(p.83):A(λ)的不变因子被其行列式因子唯一决定: d1(λ)=D1(λ), d2(λ)=D2(λ)/D1(λ), …, dn(λ)=Dn(λ)/Dn-1(λ). 不变因子决定初等因子 ? A(?)=?E-A的所有不变因子分解为一次因式