矩阵的基本运算.docx

矩阵的基本运算 （摘自：华东师范大学数学系； / ） § 3.1加和减 § 3.2矩阵乘法 § 3.2.1矩阵的普通乘法 § 3.2.2 矩阵的Kronecker乘法 § 3.3矩阵除法 § 3.4矩阵乘方 § 3.5矩阵的超越函数 § 3.6数组运算 § 3.6.1数组的加和减 § 3.6.2数组的乘和除 § 3.6.3数组乘方 § 3.7矩阵函数 § 3.7.1三角分解 § 3.7.2正交变换 § 3.7.3奇异值分解 § 3.7.4特征值分解 § 3.7.5 秩 § 3.1加和减 如矩阵A和B的维数相同，贝U A+B与A-B表示矩阵A与B的和与差.如果矩 阵A和B的维数不匹配,Matlab 阵A和B的维数不匹配, Matlab会给出相应的错误提示信息.如: A= B= 1 2 3 1 4 7 4 5 6 2 5 8 7 8 0 3 6 0 C =A+B返回： C = 2 6 10 6 10 14 10 14 0 如果运算对象是个标量（即 1 x 1 矩阵） x= -1 y=x-1= -2 0 -1 2 1 ,可和其它矩阵进行加减运算.例如: § 3.2矩阵乘法 Matlab中的矩阵乘法有通常意义上的矩阵乘法， 也有Kronecker乘法，以下 分别介绍. § 3.2.1矩阵的普通乘法 矩阵乘法用“* ”符号表示，当A矩阵列数与B矩阵的行数相等时，二者可 以进行乘法运算，否则是错误的.计算方法和线性代数中所介绍的完全相同. 如：A=[1 2 ; 3 4]; B=[5 6 ; 7 8]; C=A*B 结果为 Ix5+2x78J = Ix5+2x7 8J =^3x5+4x7 Ix6+2x8? 3x6+4x8j r19羽 ⑴50, 即Matlab返回: C = 19 22 43 50 如果A或B是标量，则A*B返回标量A （或B）乘上矩阵B （或A）的每一个 元素所得的矩阵. § 322 矩阵的Kronecker乘法 对nx m阶矩阵A和px q阶矩阵B, A和B的Kronecher乘法运算可定义 为： 疋anB …令沪 TOC \o 1-5 \h \z C^A?B=弓 3? . * 9 V 9 9 ■ ■ * ? 由上面的式子可以看出，Kronecker乘积A_ B表示矩阵A的所有元素与 B之间的乘积组合而成的较大的矩阵， B- A则完全类似.A】B和B- A均为np x mq矩阵，但一般情况下 A〔B^B〔A.和普通矩阵的乘法不同，Kronecker乘 法并不要求两个被乘矩阵满足任何维数匹配方面的要求. Kronecker乘法的 Matlab命令为C=kron（A,B），例如给定两个矩阵 A和B： 『】2、 3空 A= 3 4 B= 2 4 ? 则由以下命令可以求出 A和B的Kronecker乘积C: A=[1 2; 3 4]; B=[1 3 2; 2 4 6]; C=kro n(A,B) C = TOC \o 1-5 \h \z 3 2 2 6 4 4 6 4 8 12 9 6 4 12 8 6 12 18 8 16 24 作为比较，可以计算B和A的Kronecker乘积D,可以看出C D是不同的: A=[1 2; 3 4]; B=[1 3 2; 2 4 6]; D=kron(B,A) D = 1 2 3 6 2 4 3 4 9 12 6 8 2 4 4 8 6 12 6 8 12 16 18 24 § 3.3 矩阵除法 在Matlab中有两种矩阵除法符号：即左除和“/”即右除?如果 A 矩阵是非奇异方阵，则 A\B是A的逆矩阵乘B,即inv(A)*B ;而B/A是B乘A 的逆矩阵，即 B*inv(A) ．具体计算时可不用逆矩阵而直接计算． 通常： x=A\B就是A*x=B的解； x=B/A就是x*A=B的解. 当B与A矩阵行数相等可进行左除.如果A是方阵，用高斯消元法分解因数.解 方程：A*x(:, j)=B(:, j) ，式中的(:,j) 表示B矩阵的第j列，返回的结果x 具有与B矩阵相同的阶数，如果A是奇异矩阵将给出警告信息. 如果A矩阵不是方阵，可由以列为基准的 Householder正交分解法分解，这种分解法可以解决在最小二乘法中的欠定方程或超定方程，结果是 mK n的x矩 阵.m是A矩阵的列数，n是B矩阵的列数.每个矩阵的列向量最多有 k个非零 元素，k是A的有效秩. 右除 B/A 可由 B/A=(A\B) 左除来实现. §3.4 矩阵乘方 AAP意思是A的P次方.如果A是一个方阵，P是一个大于1的整数，贝U AAP 表示A的P次幕，即A自乘P次.如果P不是整数，计算涉及到特征值和特征向 量的问题，如已经求得： [V,D]=eig(A) ，则： AAp=v*D.Ap/v (注：这里的.A表示数组乘方，