Skip to content

第三讲:乘法和逆矩阵

上一讲大概介绍了矩阵乘法和逆矩阵,本讲就来做进一步说明。

矩阵乘法

  • 行列内积:有m×n矩阵An×p矩阵BA的总列数必须与B的总行数相等),两矩阵相乘有AB=CC是一个m×p矩阵,对于C矩阵中的第i行第j列元素cij,有:

    cij=rowicolumnj=k=inaikbkj

    其中aikA矩阵的第i行第k列元素,bkjB矩阵的第k行第j列元素。

    可以看出cij其实是A矩阵第i行点乘B矩阵第j[rowi][columnj]=[cij]

  • 整列相乘:上一讲我们知道了如何计算矩阵乘以向量,而整列相乘就是使用这种线性组合的思想:

    [Acol1Acol2Acoln][b1jb2jbnj]=[(b1jAcol1+b2jAcol2++bnjAcoln)]

    上面的运算为B的第j个列向量右乘矩阵A,求得的结果就是C矩阵的第j列,即C的第j列是A的列向量以B的第j列作为系数所求得的线性组合,Cj=b1jAcol1+b2jAcol2++bnjAcoln

  • 整行相乘:同样的,也是利用行向量线性组合的思想:

    [ai1ai2ain][Brow1Brow2Brown]=[(ai1Brow1+ai2Brow2++ainBrown)]

    上面的运算为A的第i个行向量左乘矩阵B,求得的结果就是C矩阵的第i行,即C的第i行是B的行向量以A的第i行作为系数所求的的线性组合,Ci=ai1Brow1+ai2Brow2++ainBrown

  • 列乘以行:用A矩阵的列乘以B矩阵的行,得到的矩阵相加即可:

    [Acol1Acol2Acoln][Brow1Brow2Brown]=Acol1Brow1+Acol2Brow2++AcolnBrown

    注意,AcoliBrowi是一个m×1向量乘以一个1×p向量,其结果是一个m×p矩阵,而所有的m×p矩阵之和就是计算结果。

  • 分块乘法:[A1A2A3A4][B1B2B3B4]=[A1B1+A2B3A1B2+A2B4A3B1+A4B3A3B2+A4B4]

    在分块合适的情况下,可以简化运算。

逆(方阵)

首先,并不是所有的方阵都有逆;而如果逆存在,则有A1A=I=AA1。教授这里提前剧透,对于方阵,左逆和右逆是相等的,但是对于非方阵(长方形矩阵),其左逆不等于右逆。

对于这些有逆的矩阵,我们称其为可逆的或非奇异的。我们先来看看奇异矩阵(不可逆的):A=[1236],在后面将要学习的行列式中,会发现这个矩阵的行列式为0

观察这个方阵,我们如果用另一个矩阵乘A,则得到的结果矩阵中的每一列应该都是[12]的倍数,所以我们不可能从AB的乘积中得到单位矩阵I

另一种判定方法,如果存在非零向量x,使得Ax=0,则矩阵A不可逆。我们来用上面的矩阵为例:[1236][31]=[00]

证明:如果对于非零的x仍有Ax=0,而A有逆A1,则A1Ax=0,即x=0,与题设矛盾,得证。

现在来看看什么矩阵有逆,设A=[1327],我们来求A1[1327][abcd]=[1001],使用列向量线性组合的思想,我们可以说A乘以A1的第j列,能够得到I的第j列,这时我会得到一个关于列的方程组。

接下来介绍高斯-若尔当(Gauss-Jordan)方法,该方法可以一次处理所有的方程:

  • 这个方程组为{[1327][ab]=[10][1327][cd]=[01],我们想要同时解这两个方程;

  • 构造这样一个矩阵[13102701],接下来用消元法将左侧变为单位矩阵;

  • [13102701]row22row1[13100121]row13row2[10730121]

  • 于是,我们就将矩阵从[AI]变为[IA1]

而高斯-若尔当法的本质是使用消元矩阵E,对A进行操作,E[AI],利用一步步消元有EA=I,进而得到[IE],其实这个消元矩阵E就是A1,而高斯-若尔当法中的I只是负责记录消元的每一步操作,待消元完成,逆矩阵就自然出现了。

本站没有备案,因为不需要备案