Python数据分析从小白到专家
上QQ阅读APP看本书,新人免费读10天
设备和账号都新为新人

4.1.3 矩阵及矩阵的运算

在讨论“矩阵”的概念之前,读者需要先了解“线性方程组”的概念。线性方程组分为“n元非齐次线性方程组”和“n元齐次线性方程组”两种。n元非齐次线性方程组是指由n个未知数、m个方程组成的方程组,其中,aij是第i个方程的第j个未知数的系数,bi是第i个方程的常数项,i=1,2,…,mj=1,2,…,n。当常数项b1b2,…,bm不全为0时,该方程组被称为“n元非齐次线性方程组”,示例如下:

反之,当b1b2,…,bm全为0时,该方程组被称为“n元齐次线性方程组”,示例如下:

n元线性方程组”往往会被简称为“线性方程组”或“方程组”。解的概念也在4.1.1节中提到过。既然知道了矩阵可以用来求解方程组,那么接下来就来介绍矩阵的定义。

m×n个数aiji=1,2,…,mj=1,2,…,n)排成的mn列的数表加上中括号后得到如下式子:

上式被称为mn列矩阵,简称为“m×n矩阵”。

若组成矩阵的元素全是实数,则称为“实矩阵”;若全为复数,则称为“复矩阵”。行和列都为n的矩阵称为“n阶方阵”。只有一行的矩阵称为“行矩阵”或“行向量”,示例如下:

A=[a1 a2an]

只有一列的矩阵被称为“列矩阵”或“列向量”,示例如下:

下面介绍几种不同类型的矩阵的概念,分别为“系数矩阵”、“未知数矩阵”、“常数项矩阵”和“增广矩阵”,示例如下:

式中,A是“系数矩阵”,x是“未知数矩阵”,b是“常数项矩阵”,B是“增广矩阵”。

紧接着介绍对角阵,构成矩阵Λ的主对角线的数字全不为0,而其他的元素全为0,则称为“对角阵”,示例如下:

而主对角线的数字全为1的对角阵称为“单位阵”,记为E,示例如下:

不仅是线性方程的求解需要借助矩阵运算,众所周知,图像在计算机上显示为二维的像素点阵,想要将图像旋转,也需借助矩阵运算,实际上是先将二维的像素点阵转化为2阶矩阵,再通过矩阵乘法计算。在此之前,读者还需要了解矩阵的相关运算法则。

先来看矩阵的加法运算,设有两个m×n的矩阵A=(aij)和B=(bij),那么它们的和记为

只有当矩阵A、矩阵B和矩阵C都是mn列时,AB才可以相加,且需要符合以下两条规则。

(1)A+B=B+A

(2)(A+B)+C=A+(B+C)。

同时规定矩阵的减法运算为A-B=A+(-B),即将矩阵B中的值全部加负号和矩阵A相加求和。下面来介绍数与矩阵的乘法以及矩阵与矩阵的乘法。

λ与矩阵A的乘积记作λAAλ,式子如下:

数与矩阵的乘法满足下面的规则(假设ABm×n的矩阵,λμ为数)。

(1)(μλ)A=μ(λA)。

(2)(μ+λ)A=μA+λA

(3)μ(A+B)=μA+μB

矩阵相加和数乘矩阵统称为矩阵的线性运算。

这里以一个2×3矩阵与一个3×2矩阵相乘为例,介绍矩阵与矩阵相乘的用法,示例如下:

由此可以得出结论:设A=(aij)是一个m×s矩阵,B=(bij)是一个s×n矩阵,那么矩阵A与矩阵B相乘的结果必定是一个m×n矩阵,记为C=(cij),并且满足

需要注意的是,矩阵乘法不满足交换律,即ABBA,当然这只是在通常情况下。若满足AB=BA这一条件,就称“矩阵A与矩阵B是可交换的”,而这只是一个巧合。矩阵乘法虽然不满足交换律,但是满足结合律和分配律,即

(1)(AB)C=A(BC);

(2)λ(AB)=(λA)B;

(3)A(B+C)=AB+AC;

(4)EmAm×n=Am×n=Am×nEm(可简记为EA=A=AEE是单位矩阵)。

还有一个概念是矩阵的转置,矩阵D的转置矩阵记为DT。示例如下:

它的转置矩阵为

矩阵的转置作为一种单独的运算,遵循如下的运算法则。

(1)(AT)T=A

(2)(A+B)T=AT+BT

(3)(λA)T=λAT

(4)(AB)T=BTAT

BT=B,则称矩阵B是“对称矩阵”,简称“对称阵”,它的特点是元素以对角线为对称轴对应相等。

下面介绍方阵的行列式:由n阶方阵A的元素所构成的行列式(各元素的位置不变),称为方阵A的行列式,记为detA或|A|。

需要注意的是,方阵与行列式是两个不同的概念,n阶方阵是n2个数按一定方式排列成的数表,而n阶行列式则是这些数构成的数表A按一定运算法则得到的一个数。方阵的行列式的运算法则如下。

(1)|AT|=|A|。

(2)|λA|=λn|A|。

(3)|AB|=|A||B|。

下面介绍“伴随矩阵”。行列式|A|的各个元素的代数余子式Aij所构成的如下矩阵:

称为矩阵A的“伴随矩阵”,简称“伴随阵”,记为A*

同样地,根据伴随矩阵的性质,其满足如下公式:

A*A=|A|E

“逆矩阵”又被称为“逆阵”。假设矩阵A可逆,那么矩阵A只有一个逆矩阵,且满足如下式子:

AB=BA=E

其中,E是单位阵,A是示例矩阵,B是矩阵A的逆矩阵,也可记为A-1,即B=A-1。逆矩阵具有如下性质:

其中,A*是矩阵A的伴随矩阵;A-1是矩阵A的逆矩阵;|A|是矩阵A的行列式。这里有一个关于逆矩阵的重要推论:若AB=EBA=E,则B=A-1,这也印证了在开始讲解逆矩阵时的说法。

克拉默法则是用来求nn元线性方程组的。这里使用本节开头所示的线性方程组,具体如下:

如果上式的系数矩阵行列式不等于0,那么齐次方程组有唯一解,如下:

其中,Ajj=1,2,…,n)是把系数矩阵A中第j列的元素用方程组右端的常数项代替后得到的n阶矩阵,如下: