华中科技大学研究生矩阵论Matrix.ppt
第 6 章 矩阵的Kronecker积和Hadamard积,The Kronecker Product and Hadamard Product,概述：,主要内容：介绍Kronecker积和Hadamard积讨论 K-积，H-积的运算性质、之间的关系 K-积与矩阵乘积的关系 K-积，H-积的矩阵性质 K-积的矩阵等价与相似关系应用：求解矩阵方程 向量化算子重点：K-积及其应用,6.1 Kronecker积和Hadamard积的定义,定义6.1（P.136）设矩阵 A=aijmn和 B=bijst，则A和B的 Kronecker被定义为 AB：AB=aijBmsnt 设A=aijmn和 B=bijmn为同阶矩阵，则A和B的Hadamard被定义为 AB：AB=aijbijm n,6.1 K-积和H-积的定义,例题1 设，计算 AB，BA，I2B，AB，I2A,例题1 设，计算 AB，BA，I2B，AB，I2A,分块对角矩阵,对角矩阵,6.1 K-积和H-积的定义,例题2 设分块矩阵A=(Ast)，则 AB=(Ast B)特别地，若A=(A1,A2,An)，则 AB=(A1B,A2B,AnB),例题3 快速Walsh(Hadamard)变换 yN=HNxN，其中于是有,6.1 K-积和H-积的定义,例题2 设分块矩阵A=(Ast)，则 AB=(Ast B)特别地，若A=(A1,A2,An)，则 AB=(A1B,A2B,AnB),例题3 快速Walsh(Hadamard)变换 yN=HNxN，其中于是有,6.1 K-积和H-积的定义,K-积，H-积的基本结果：A和B中有一个为零矩阵，则 AB=0，AB=0 II=I，II=I 若A为对角矩阵，则AB为分块对角矩阵，AB为对角矩阵。K-积的基本性质 定理6.1（P.138）设以下矩阵使计算有意义，则(kA)B=A(kB)A(B+C)=AB+AC(AB)C=A(BC)(AB)H=AH BHAB BA,H-积的基本性质：设A，B为同阶矩阵，则 AB=BA(kA)B=A(kB)A(B+C)=AB+AC(AB)C=A(BC)(AB)H=AH BHKronecker和Hadamard的关系：定理6.3（P.139）AB 可由AB的元素构成。,K-积与矩阵乘法 定理6.2（P.138）设矩阵A，B，C，D使得下列运算有意义，则有(AB)(CD)=(AC)(BD)意义：建立Kronecker积和矩阵乘法的相互转换。特别情形：设 AFmm，B Fnn，则 AB=(ImA)(BIn)=(AIm)(InB)=(ImB)(AIn)=(AIn)(ImB)(AB)k=Ak Bk,(A1B1C1)(A2B2C2)=(A1A2)(B1B2)(C1C2),(A1B1)(A2B2)(A3B3)=(A1A2A3)(B1B2B3),6.2 Kronecker积和Hadamard积的性质,Kronecker积的矩阵性质定理6.4（P.140）设矩阵使下列运算有意义，则当A，B分别为可逆矩阵时，AB和BA均为可逆矩阵，而且有(AB)1=A1 B1当方阵AFmm，BFnn时，方阵ABFmnmn的行列式为|AB|=|BA|=|A|n|B|m若A，B是Hermite矩阵，则AB 和BA均是Hermite矩阵若A，B是酉矩阵，则AB和BA均是酉矩阵。,Kronecker与矩阵等价、相似关系定理6.5（P.141）设矩阵A，B，为等价矩阵，则(AI)等价于(BI)设方阵A相似与JA，方阵B相似于JB，则(AB)相似于(JAJB)K-积特征值和特征向量定理6.6（P.142）设AFmm 的特征值、特征向量分别是i，xi，B Fnn的特征值、特征向量分别是 j，yj，则(AB)的特征值是ij。特征向量是(xiyj)。(AIn)+(ImB)的特征值是i+j，特征向量是(xiyj),Kronecker和，记为AB,Kronecker与矩阵等价、相似关系推论 若A，B正定（半正定），则AB和AB均正定（半正定）；若A相似于JA，B相似于JB，则 AB 相似于 JAJB，AB 相似于 JAJB。更一般的结果：定理6.7（P.142）的特征值为,Kronecker积的矩阵函数性质定理6.8（P.143）设是f(z)解析函数，f(A)有意义，则 f(IA)=If(A)f(AI)=f(A)I特例：,定理的证明思路：利用定理5.12，矩阵函数可由多项式表示。也可以直接用极限性质证明。,SN(IA)=ISN(A),SN(AI)=SN(A)I,例题1 设 AFmn，BFst，证明 rank(A B)=rank(A)rank(B),例题2（P.144），设，求(AB)的特征值和特征向量 求(AI)+(IB)的特征值和特征向量,例题3：证明对任何方阵A,B,有,Hadamard积的性质定理6.9（Schur积定理）设A、B为同阶方阵。若A和B半正定（正定），则AB亦半正定（正定）。,证明思路：利用定理3.6，有,推出 AB可表示为,6.3 矩阵的向量化算子和K-积,向量化算子Vec:Fmn Fmn定义（P.143）设 A=aijmn,则 Vec(A)=(a11 a21 am1；a12 a22 am2；a1n a2n amn)T 性质：（P.146）Vec是线性算子，并保持线性关系不变：Vec(k1A+k2B)=k1Vec(A)+k2Vec(B)2.定理6.10（P.146）Vec(ABC)=(CT A)VecB 3.Vec(AX)=(I A)VecX4.Vec(XC)=(CTI)VecX,令 B=X,C=I,令 B=X,A=I,用向量化算子求解矩阵方程,思想：用Vec算子，结合Kronecker积将矩阵方程化为线性方程组求解。1、AFmm，BFnn，DFmn，AX+XB=D分析：AX+XB=D(IA+BTI)VecX=VecDG=(IA+BTI)，方程有惟一解的充要条件是G为可逆矩阵，即A和-B没有共同的特征值。例题1（P.147）,用向量化算子求解矩阵方程,2、A，XFnn，AX XA=kX分析：AX XA=kX(I AATI)VecX=kVecX H=(IA ATI)，方程(kI H)y=0 有非零解的充要条件是k为H的特征值，k=i j。例题2 求解矩阵方程 AX XA=2X,用向量化算子求解矩阵方程,3 A，B，D，XFnn，AXB=D分析：AXB=D(BTA)VecX=VecD L=BT A，方程有惟一解的充要条件是L为可逆矩阵.例题3 求解方程 A1XB1+A2XB2=D,例题4 设A Cmm，B Cnn，D Fmn，证明谱半径(A)(B)1 时方程：X=AXB+D 的解为,证,用向量化算子求解矩阵微分方程,A，B，XFnn，X(t)=AX(t)+X(t)B,X(t0)=C VecX(t)=(I A+BTI)VecX(t),VecX(t0)=VecC。,交换矩阵Kmn及其性质,定理6.11(1)(2)(3),定理6.12 设 则,定理6.13 设 则,复习选讲：,线性空间的表示线性变换与变换矩阵线性变换的确定方法相应变换矩阵的求法矩阵分解与空间分解准对角矩阵分解与不变子空间的分解可对角化矩阵的分解与特征子空间的分解幂等矩阵的空间分解JA，mA()，f()=|I-A|之间的关系A与f(A)在Jordan标准形上的关系f(A)的矩阵性质,复习选讲,正规矩阵的性质与应用向量范数与矩阵范数 向量的p范数 矩阵的F范数和p范数矩阵幂级数和矩阵函数 矩阵幂级数的收敛与矩阵函数的意义 矩阵幂级数的求和与矩阵函数的计算 矩阵函数与矩阵多项式,习题选讲,P150：9P31:1(3),17,P58：6,11，20P92：11，12，15，,
