本篇文章介绍四个应用:
A
的逆,求解
1、计算
A−1
。
2×2
矩阵展示了伴随矩阵如何表示
A−1
:
伴随矩阵还需乘以行列式,当
detA
非零时,
A
是可逆的。
a,b
分别和
C11,C12
相乘得到
ad−bc
,这是
detA
的代数余子式展开式,也是我们需要的线索:
A−1
就是用伴随矩阵除以
detA
我们的目标是验证这个公式,另外由此我们还看出为何
ACT=(detA)I
:
第一列
C11,…,C1n
乘以
a11,…,a1n
得到对角元素
detA
,
A
的每一行乘以它的代数余子式都得到
那么关键问题是:为什么对角线外的元素都是零呢?如果我们将第一行的元素
a1j
和第二行的代数余子式
C2j
进行组合,为什么结果是零呢?
答案就是:我们在计算一个新矩阵
B
的行列式。
乘法
ACT=(detA)I
直接得出
A−1
,记住去掉矩阵
A
第
例1:和矩阵的逆是微分矩阵:
因为代数余子式包含 (−1)i+j ,所以会有符号。
2、求解 Ax=b 。 x=A−1b 仅仅就是 CTb 除以 detA ,这里介绍一个非常出名的求 (x1,…,xn) 的方法:
克拉姆法则:
x=A−1b
的第
j
个元素是比值
证明:对
detB
按第
j
列展开成代数余子式的形式,因为代数余子式不考虑该列,所以
然后用它除以 detA 就得到 xj , x 的每个元素都是两个行列式的比值。
例2:下面方程组
的解
x1
需要将第一列换成0,6,
x2
需要将第二列换成0,6:
分母总是 detA ,如果有1000个方程,那么根据克莱姆法则,需要1001个行列式。
3、盒子的体积。当所有角都是直角——此时边互相垂直,行列式和体积的联系就非常紧密,此时的盒子的形状是矩形的,从而体积就是边长的乘积: volumn=ℓ1ℓ2⋯ℓn 。
当盒子的边长是
A
的行时,我们想从
ℓi
是行(边)的长度,因为行是正交的,所以对角线外的元素都是零。利用乘法和转置法则:
这个方程的平方根说明了行列式等于体积。而 detA 的符号暗示了是使用右手坐标(也就是通常的 x−y−z ),还是左手坐标(像 y−x−z )。
如果角度不是
90∘
,那么体积就不是长度的乘积。在平面上(如图1),平行四边形的“体积”等于底边
ℓ
乘以高
h
,向量
图1
对于 n 维的情况,将盒子变成矩形需要很长时间,但是思路就是如此,如果我们从每行减去它在前面行所生成空间的投影,体积和行列式都不变,通过格拉姆-施密特过程产生正交的行,这时体积=行列式,所以对于开始的行这个等式同样成立。
现在我们介绍完了体积和行列式之间的关系,但是很有必要回过头来看看最简单的情况,我们知道
这些行列式给出了体积——或者面积(因为这是在二维空间),如图2所示,平行四边形底和高都是1,所以面积也是1。
图2
4、主元公式。最后我们讨论一下消元,此时考虑不需要行交换的情况。观察的要点是第
k
个主元完全由
第一个主元只取决于第一行和第一列,第二个主元
(ad−bc)/a
值取决于
2×2
的子矩阵
A2
,
A
的其余部分不会影响到它。实际上,不仅仅是主元,整个左上角矩阵的
从上例我们看出,前两行和两列就是子矩阵 A2 的分解,如果没有行交换的前提下,这条规则恒成立。
4、如果
A
分解为
证明的关键在于进行其他行的消元之前,先将左上角的确定下来,或者使用块乘法规则:
将最后的矩阵和
A
进行比较,我们发现
前
k
个主元的乘积是
对于上面的例子,第二个主元就是比值
(ad−bc)/a
,也就是
A2
的行列式除以
A1
的行列式(约定
A0=1
,这样的话第一个主元就是
a/1=a
)。将所有主元相乘就回到我们的行列式:
根据方程(5)我们可以得出一个结论:当 detAk 都不为零时,主元也都不为零。
5、当且仅当子矩阵 A1,A2,…,An 都是非奇异时,消元法不需要行交换(所以 P=I,A=LU )。
现在考虑一个问题,行交换的次数会一直是奇数或偶数吗?如果是,那么根据2行列式要么为+1,要么为-1。
考虑 (3,2,1) ,我们可以只交换3,1就得到正序列 (1,2,3) ,交换3与2,然后3与1,最后2与1同样得到正序列,对于这两种方式,交换的次数都是奇数,这里得出的论断就是:偶数次的交换不可能从 (3,2,1) 得出正序列。
这里给出它的一个证明。观察置换中的每对数,我们用
N
表示大数在前的对数。这样的话
当相邻发生交换时,
N
要么+1,要么-1。所有交换都可以用奇数次的邻域交换实现,考虑一个例子,我们要交换第一和第四个元素,可以通过五次交换完成:
我们需要
ℓ−k
次邻域交换将
k
位置上的元素放到