Skip to content

4.1.4 矩阵的算术运算

1. 矩阵的相等

两个矩阵 A=(aμν)B=(bμν) 相等,如果它们有相同的大小,并且对应元素相等:

(4.20)A=B,当 aμν=bμν(μ=1,,m;ν=1,,n).

2. 加法和减法

两个矩阵仅当大小相同时才可能相加减. 两个矩阵的和 (差) 由它们的对应元素相加 (相减) 得到:

(4.21a)A±B=(aμν)±(bμν)=(aμν±bμν).

(137214)+(350214)=(427408) .

对于矩阵的加法, 交换律和结合律成立:

a) 交换律 A+B=B+A ;(4.21b)

b) 结合律 (A+B)+C=A+(B+C) .(4.21c)

3. 数乘矩阵

实数或复数 α 与大小为(m, n)的矩阵 A 相乘,就是用 αA 的每个元素:

(4.22a)αA=α(aμν)=(αaμν).

3(137014)=(39210312).

由 (4.22a) 可知, 我们显然可以提出矩阵每个元素都含有的常数因子. 对于矩阵与标量的乘法, 乘法交换律、结合律及分配律都成立:、及

a) 交换律 αA=Aα ;(4.22b)

b) 结合律 α(βA)=(αβ)A ;(4.22c)

c) 分配律 (α±β)A=αA±βA;α(A±B)=αA±αB .(4.22d)

4. 数除以矩阵

用标量 γ0 除矩阵即用 α=1/γ 乘矩阵.

5. 两矩阵相乘

(1)矩阵乘积 两个矩阵 AB 的乘积 AB 仅当左边因子 A 的行数等于右边因子 B 的列数时才可计算. 如果 A 是大小为(m, n)的矩阵,那么矩阵 B 的大小必须是(n, p),并且乘积 AB 是大小为(m, p)的矩阵 C=(cμν) . 元素 cμν 等于左边因子 A 的第 μ 行与右边因子 B 的第 μ 列的标量积:

(4.23)AB=(ν=1naμνbνλ)=(cμλ)=C(μ=1,2,,m;λ=1,2,,p).

A=(137214101),B=(325103) ,依照公式 (4.23),积矩阵 C 的元素 c22=2211+43=15 .

(2) 矩阵乘积的不相等 即使乘积 ABBA 都存在,通常 ABBA ,也就是说,乘法交换律一般不成立. 如果等式 AB=BA 成立,那么我们称矩阵 AB 是可换的或互相交换的.

(3) 法尔克 (Falk) 格式 矩阵乘法 AB=C 可以借助法尔克格式实现 (图 4.1). 积矩阵 C 的元素 cμν 恰好出现在 A 的第 μ 行与 B 的第 λ 列的交点.

应用法尔克格式,矩阵 A(3,3)B(3,2) 相乘如图 4.2 所示.

019363af-d8ae-7006-ac42-15a9aafbc2ce_5_696_902_249_180_0.jpg

019363af-d8ae-7006-ac42-15a9aafbc2ce_5_681_1134_280_173_0.jpg

(4) 复元素矩阵 K1K2 的乘法 对于两个复元素矩阵的乘法,可以应用依据公式 (4.2b) 得到的实部和虚部分解式: K1=A1+iB1,K2=A2+iB2 ,这里 A1,B1,A2,B2 是实矩阵. 作此分解后,相乘得到一些矩阵之和,其中各加项是实矩阵之积.

I(A+iB)(AiB)=A2+B2+i(BAAB) (关于矩阵的幂,参见第 370 页 4.1.5,8.). 当然,在做这些矩阵的乘法时,应当考虑到交换律对于乘法一般不成立, 即矩阵 AB 通常是不可互相交换的.

6. 两个向量的标量积和并积

如果将向量 ab 分别看作 1 行矩阵和 1 列矩阵,那么依据矩阵乘法法则,存在两种可能将它们相乘.

如果 a 的大小是 (1,n),b 的大小是(n,1),那么它们的乘积的大小是(1,1),即是一个数. 它称作两个向量的标量积. 反之,如果 a 有大小 (n,1),b 有大小(1, m), 那么它们的乘积的大小是(n, m),即是一个矩阵. 这个矩阵称为这两个向量的并积.

(1) 两个向量的标量积 行向量 aT=(a1,a2,,an) 与列向量 b=(b1 , b2,,bn)T (两者都有 n 个元素) 的标量积定义为数,

(4.24)aTb=bTa=a1b1+a2b2++anbn=μ=1naμbμ.

乘法交换律对于标量积一般不成立,因此必须准确保持 aTb 的顺序. 如果交换乘法顺序,那么乘积 baT 将是并积.

(2) 两个向量的并积或张量积 n 维行向量 a=(a1,a2,,an)Tm 维列向量 bT=(b1,b2,,bm) 的并积定义为下列大小为(n, m)的矩阵:

(4.25)abT=(a1b1a1b2a1bma2b1a2b2a2bmanb1anb2anbm).

在此乘法交换律一般也不成立.

(3) 关于两个向量的向量积概念的提示 在多向量或交错张量的范围中有所谓外积, 其三维形式就是熟知的向量积或叉积 (参见第 247 页 3.5.1.5, 2. 及其后). 本书不讨论高秩多向量的外积.

7. 矩阵的秩

(1)定义 矩阵 A 中线性无关的列向量的最大个数 r 等于线性无关的行向量的最大个数. 这个数 r 称为矩阵的秩,并且记作 rank(A)=r .

(2)关于矩阵的秩的一些说明

a) 因为在 m 维向量空间中不存在多于 m 个线性无关的行向量或列向量 (参见第 490 页 5.3.8.2),所以大小为(m, n)的矩阵 A 的秩不可能大于 mn 中的较小者:

(4.26a)rank(A(m,n))=rmin{m,n}.

b) 方阵 A(n,n) 称为正则矩阵,如果

(4.26b)rank(A(n,n))=r=n.

大小为(n, n)的方阵是正则的,当且仅当它的行列式不为零,即 detA0 (参见第 373 页 4.2.2,3.) . 不然它是奇异方阵.

c) 因此,对于奇异方阵 A(n,n) ,即 detA=0 ,有

(4.26c)rank(A(n,n))=r<n.

d) 零矩阵0的秩等于零:

(4.26d)rank(0)=0.

e) 矩阵的和及积的秩满足关系式

(4.26e)|rank(A)rank(B)|rank(A+B)rank(A)+rank(B),(4.26f)rank(AB)min{rank(A),rank(B)}.

(3) 确定秩的法则 初等变换不改变矩阵的秩. 其中初等变换是

a) 交换两列或两行.

b) 用数乘一行或一列.

c) 将一行加到另一行或将一列加到另一列.

为了确定它们的秩, 可以通过适当的行的线性组合将每个矩阵变换为这种形式: 在第 μ 行中至少最初 μ1 个元素等于零 (μ=2,3,,m) (高斯算法原理,参见第 417 页 4.5.2.4). 在变换得到的矩阵中非零向量的行向量的个数等于矩阵的秩 r .

8. 逆矩阵

对于正则矩阵 A=(aμν) 总存在 (对于乘法的) 逆矩阵 A1 ,即

(4.27a)AA1=A1A=I.

A1=(βμν) 的元素是

(4.27b)βμν=AνμdetA

其中 Aνμ 是属于矩阵 A 的元素 aνμ 的余因子 (参见第 372 页 4.2.1,1.). 为实际计算 A1 ,可应用第 373 页 4.2.2,2. 所给出的方法. 在大小为(2,2)的矩阵的情形中, 有

(4.28)A=(abcd),A1=1adbc(dbca).

注 为什么不定义矩阵的除法, 而是应用逆矩阵进行计算? 这与除法不能唯一地定义的事实相关. 方程

BX1=A 与 X2B=A (其中 B 正则) 

的解

(4.29)X1=B1A 与 X2=AB1

一般是不同的.

9. 正交矩阵

如果对于方阵 A 关系式

(4.30)AT=A1 或 AAT=ATA=I

成立, 则称它是一个正交矩阵, 这就是说, 它的任一行与另一行的转置的标量积, 或任一列的转置与另一列的标量积都为零, 同时任一行与它自己的转置或任一列的转置与它自己的标量积都等于 1 .

正交矩阵有下列性质:

a) 正交矩阵的转置和逆也是正交矩阵; 并且其行列式

(4.31)detA=±1.

b) 正交矩阵的积也是正交矩阵.

旋转矩阵 D 也是正交矩阵,它被用来刻画坐标系的旋转,其元素是新轴向的方向余弦 (参见第 284 页 3.5.3.3, 2.).

10. 酉矩阵

如果对于复元素矩阵 A 关系式

(4.32)(A)T=A1 或 A(A)T=(A)TA=I

成立, 则称它是一个酉矩阵. 在实数情形, 酉矩阵与正交矩阵相同.

version 1.24.0