博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
线性代数(MIT)
阅读量:3953 次
发布时间:2019-05-24

本文共 6476 字,大约阅读时间需要 21 分钟。

线性代数

矩阵消元

矩阵行变化, 相当于左乘矩阵.

如下矩阵乘法为对第一行进行行变化 r o w 1 : = r o w 1 + 2 ∗ r o w 2 − r o w 3 row1 := row1 + 2*row2 - row3 row1:=row1+2row2row3.

[ 1 2 − 1 0 1 0 0 0 1 ] [ − . r o w 1 − − . r o w 2 − − . r o w 3 − ] \begin{bmatrix} 1 & 2 & -1\\ 0 & 1 & 0\\ 0 & 0 & 1\\ \end{bmatrix} \begin{bmatrix} -. & row1 & - \\ -. & row2 & - \\ -. & row3 & - \\ \end{bmatrix} 100210101...row1row2row3

矩阵乘法的几种理解

A n × m ⋅ B m × p = C n × p A_{n \times m}\cdot B_{m \times p} = C_{n \times p} An×mBm×p=Cn×p

从单个元素的角度

C i , j = ∑ k = 1 m A i , k ⋅ B k , j C_{i,j} = \sum_{k=1}^m{A_{i,k} \cdot B_{k,j}} Ci,j=k=1mAi,kBk,j

列角度

[ ∣ ∣ ⋯ c o l 1 c o l 2 ⋯ ∣ ∣ ⋯ ] ⋅ [ ∣ ⋯ c o l 1 ⋯ ∣ ⋯ ] = [ ∣ ⋯ c o l 1 ⋯ ∣ ⋯ ] \begin{bmatrix} | & | & \cdots \\ col1 & col2 & \cdots \\ | & | & \cdots \\ \end{bmatrix} \cdot \begin{bmatrix} | & \cdots \\ col1 & \cdots \\ | & \cdots \\ \end{bmatrix}= \begin{bmatrix} | & \cdots \\ col1 & \cdots \\ | & \cdots \\ \end{bmatrix} col1col2col1=col1

c o l 1   o f   C col1\ of\ C col1 of C A A A的每一列的线性组合, 每一列的权重为 c o l 1   o f   B col1\ of\ B col1 of B.

行角度

[ − . r o w 1 − ⋯ ⋯ ⋯ ⋮ ⋮ ⋮ ] ⋅ [ − . r o w 1 − − . r o w 2 − ⋮ ⋮ ⋮ ] = [ − . r o w 1 − ⋯ ⋯ ⋯ ⋮ ⋮ ⋮ ] \begin{bmatrix} -. & row1 & - \\ \cdots & \cdots & \cdots \\ \vdots & \vdots & \vdots \\ \end{bmatrix} \cdot \begin{bmatrix} -. & row1 & - \\ -. & row2 & - \\ \vdots & \vdots & \vdots \\ \end{bmatrix}= \begin{bmatrix} -. & row1 & - \\ \cdots & \cdots & \cdots \\ \vdots & \vdots & \vdots \\ \end{bmatrix} .row1..row1row2=.row1

r o w 1   o f   C row1\ of\ C row1 of C B B B的每一行的线性组合, 每一行的权重为 r o w 1   o f   A row1\ of\ A row1 of A.

列乘行的角度

A = [ c o l 1 c o l 2 ⋯ ] , B = [ r o w 1 r o w 2 ⋮ ] A ⋅ B = ∑ i = 1 m c o l i ⋅ r o w i A = \begin{bmatrix} col1 & col2 & \cdots \end{bmatrix}, B = \begin{bmatrix} row1 \\ row2 \\ \vdots \end{bmatrix}\\ A \cdot B = \sum_{i=1}^m{col_i \cdot row_i} A=[col1col2],B=row1row2AB=i=1mcolirowi

例子

c o l 1 = [ 2 3 4 ] , r o w 1 = [ 1 6 ] c o l 1 ⋅ r o w 1 = [ 2 3 4 ] ⋅ [ 1 6 ] = [ 2 12 3 18 4 24 ] col1 = \begin{bmatrix}2 \\ 3 \\ 4 \end{bmatrix},row1=\begin{bmatrix}1 & 6\end{bmatrix}\\ col1 \cdot row1 = \begin{bmatrix}2 \\ 3 \\ 4 \end{bmatrix} \cdot \begin{bmatrix}1 & 6\end{bmatrix}= \begin{bmatrix} 2 & 12\\ 3 & 18\\ 4 & 24\\ \end{bmatrix} col1=234,row1=[16]col1row1=234[16]=234121824

分块角度

以上三种角度都可以用分块角度很直观地看出.

矩阵的逆

可逆

若存在 x ≠ 0 ,   A ⋅ x = 0 x \neq 0,\ A \cdot x = 0 x=0, Ax=0, 则 A A A是奇异矩阵, 不可逆.

证明:

假设存在 B ⋅ A = I B \cdot A = I BA=I, 则有

A ⋅ x = 0 B ⋅ A x = 0 I ⋅ x = 0 x = 0 A\cdot x = 0\\ B \cdot Ax = 0\\ I \cdot x = 0\\ x = 0 Ax=0BAx=0Ix=0x=0
矛盾. 故 B 不 存 在 B不存在 B.

逆的计算

Gauss-Jordan消元法对 [ A I ] \begin{bmatrix}A & I\end{bmatrix} [AI]进行消元, 最后100左边消为 I I I, 右边即为 A − 1 A^{-1} A1. 由于矩阵行变化相当于左乘矩阵, 故

∵ E ⋅ A = I ∴ E ⋅ I = A − 1 E ⋅ [ A I ] = [ I A − 1 ] \because E \cdot A = I\\ \therefore E \cdot I = A^{-1}\\ E \cdot \begin{bmatrix}A & I\end{bmatrix} = \begin{bmatrix}I & A^{-1}\end{bmatrix} EA=IEI=A1E[AI]=[IA1]

LU分解

把A进行行变化, 得到上三角矩阵U.(不使用行交换).

A 3 × 3 A_{3\times 3} A3×3例子

E 2 , 1 ⋅ E 3 , 1 ⋅ E 3 , 2 ⋅ A = U A = E 3 , 2 − 1 ⋅ E 3 , 1 − 1 ⋅ E 2 , 1 − 1 ⋅ U = L U E_{2,1}\cdot E_{3,1}\cdot E_{3,2}\cdot A = U\\ A = E_{3,2}^{-1}\cdot E_{3,1}^{-1}\cdot E_{2,1}^{-1}\cdot U = LU E2,1E3,1E3,2A=UA=E3,21E3,11E2,11U=LU

转置、置换、向量空间

置换矩阵

置换矩阵(permutations)是指可以让矩阵进行 行交换(row exchanges)的矩阵。 n × n n\times n n×n的置换矩阵就是通过单位矩阵行交换的得到,个数为 n ! n! n!。他们的逆矩阵都在这个群里。进一步来说,逆矩阵相当于把做的交换逆回去,故有

P − 1 = P T P^{-1} = P^{T} P1=PT

所以,当需要行交换的时候, LU分解应该写成

P A = L U PA=LU PA=LU

转置

对于任意矩阵 R R R R T R R^{T}R RTR是对称矩阵恒成立。即 ( R T R ) T = R T R (R^TR)^T=R^TR (RTR)T=RTR.

向量空间

满足数乘和加法封闭性。最小的向量空间只有零向量。

列空间和零空间

两个子空间 S S S T T T的交集也是一个子空间。因为 v , w ∈ S ∩ T v,w\in S\cap T v,wST,所以 v + w ∈ S , v + w ∈ T v+w\in S, v+w\in T v+wS,v+wT,即 v + w ∈ S ∩ T v+w\in S\cap T v+wST

列向量空间

A = [ 1 1 2 2 1 3 3 1 4 4 1 5 ] A=\begin{bmatrix}1 & 1 & 2 \\ 2 & 1 & 3 \\ 3 & 1 & 4 \\ 4 & 1 & 5\end{bmatrix} A=123411112345

A的列向量子空间是4维空间中的子空间,不能覆盖整个4维空间。用方程组来说,就是 A x = b Ax=b Ax=b并不是对于所有的 b b b都有解, x x x是各个列向量线性组合的权重。有解的 b b b是各列的线性组合,即 b b b 属于 A A A的列空间 C ( A ) C(A) C(A)

A ⋅ x = [ 1 1 2 2 1 3 3 1 4 4 1 5 ] ⋅ [ x 1 x 2 x 3 ] = [ b 1 b 2 b 3 b 4 ] A \cdot x = \begin{bmatrix}1 & 1 & 2 \\ 2 & 1 & 3 \\ 3 & 1 & 4 \\ 4 & 1 & 5\end{bmatrix}\cdot \begin{bmatrix}x_1 \\ x_2 \\ x_3 \end{bmatrix}= \begin{bmatrix} b_1\\ b_2\\ b_3\\ b_4 \end{bmatrix} Ax=123411112345x1x2x3=b1b2b3b4

零空间

A x = 0 Ax=0 Ax=0的所有解 x ∈ R 4 x\in R^4 xR4组成零空间 N ( A ) N(A) N(A)。至少包含零向量。

用上面的例子,零空间的向量形如以下形式,是一条直线。

[ c c − c ] \begin{bmatrix} c\\ c\\ -c\\ \end{bmatrix} ccc
证明为什么 A x = 0 Ax=0 Ax=0得到的所有向量一定构成一个向量空间

x , x ∗ 为 两 个 解 x,x^*为两个解 x,x,则有

A x = 0 ,   A x ∗ = 0 Ax=0,\ Ax^*=0\\ Ax=0, Ax=0
A ( x + x ∗ ) = A x + A x ∗ = 0 + 0 = 0 A ( k x ) = k A x = 0 A(x + x^*) =Ax+Ax^*=0+0=0\\ A(kx)=kAx=0 A(x+x)=Ax+Ax=0+0=0A(kx)=kAx=0
故向量加法和数乘都在向量空间内,故符合向量空间的定义。

A x = 0 Ax=0 Ax=0的主元和特解

A x = 0 Ax=0 Ax=0 的解不会因对 A A A矩阵进行初等行变化而改变。可以从两个角度得到这个结论。

  1. 进行行变化时,只是对A矩阵每一列内部的行之间进行变化,每列的变化相同。列与列之间的关系不变。
  2. 进行初等行变化相当于左乘一个矩阵 E E E,而 E A x = 0 EAx=0 EAx=0仍然成立。

矩阵的秩就是矩阵进行行变化之后的主元的个数。

A = [ 1 2 2 2 2 4 6 8 3 6 8 10 ] = > U = [ 1 2 2 2 0 0 2 4 0 0 0 0 ] = > R = [ 1 2 0 − 2 0 0 1 2 0 0 0 0 ] A=\begin{bmatrix}1 & 2 & 2 & 2\\ 2 & 4 & 6 & 8 \\ 3 & 6 & 8 & 10\end{bmatrix}\\ =>\\ U=\begin{bmatrix}1 & 2 & 2 & 2\\ 0& 0& 2 & 4 \\ 0 & 0& 0 & 0\end{bmatrix}\\ =>\\ R=\begin{bmatrix}1 & 2 & 0 & -2\\ 0& 0& 1 & 2 \\ 0 & 0& 0 & 0\end{bmatrix}\\ A=1232462682810=>U=100200220240=>R=100200010220
A A A的秩 r = 2 r=2 r=2 U U U的第一列和第三列为主列(pivot columns),第二列和第四列为自由列(free columns)。 A x = 0 Ax=0 Ax=0的自由元的个数为 n − r = 4 − 2 = 2 n-r=4-2=2 nr=42=2。零空间就是让自由元逐一等于1得到的向量之间的线性组合。
x = c ⋅ [ − 2 1 0 0 ] + d ⋅ [ 2 0 − 2 1 ] x=c\cdot \begin{bmatrix} -2\\ 1\\ 0\\ 0 \end{bmatrix} +d \cdot \begin{bmatrix} 2\\ 0\\ -2\\ 1 \end{bmatrix} x=c2100+d2021
R R R矩阵中的主列凑在一起,把自由列组合在一起,就可以得到rref形式的 R R R矩阵。
I = [ 1 0 0 1 ] , F = [ 2 − 2 0 2 ] R = [ [ 1 0 0 1 ] [ 2 − 2 0 2 ] 0 0 0 0 ] I=\begin{bmatrix}1&0\\0&1\end{bmatrix},F=\begin{bmatrix}2&-2\\0&2\end{bmatrix} R=\begin{bmatrix}\begin{bmatrix}1&0\\0&1\end{bmatrix}&\begin{bmatrix}2&-2\\0&2\end{bmatrix}\\ \begin{matrix}0&0\end{matrix}&\begin{matrix}0&0\end{matrix}\end{bmatrix} I=[1001],F=[2022]R=[1001]00[2022]00

在这里插入图片描述

则可以通过解 R N = 0 RN=0 RN=0来得到列向量为 A x = 0 Ax=0 Ax=0 n − r n-r nr个特解,即得到了零基。

N = [ − F I ] N=\begin{bmatrix} -F\\ I\\ \end{bmatrix} N=[FI]

转载地址:http://fwuzi.baihongyu.com/

你可能感兴趣的文章
十张图看懂未来大数据世界
查看>>
“揭秘”大数据的10个神话!
查看>>
《中国区块链行业发展报告2018》全文发布!
查看>>
高盛发布区块链报告:从理论到实践(中文版)
查看>>
用Python从零开始创建区块链
查看>>
使用 Charles 抓取 app 数据包
查看>>
Linux time subsystem 详解(1) ----概述
查看>>
大牛很通俗地介绍《信号与系统》
查看>>
执行程序(例如UltraEdit)在WIN7下添加到右键菜单
查看>>
flash and root your Nexus10
查看>>
深入学习Make命令和Makefile(上)(2)
查看>>
深入学习Make命令和Makefile(下)(2)
查看>>
10大玄机为你揭开炒土豆丝爽脆的秘密——尖椒土豆丝
查看>>
grep与正则表达式
查看>>
git patch 使用
查看>>
如何进行Linux Kernel 开发
查看>>
技术人攻略访谈二十九:平行世界守护者
查看>>
制作initramfs/initrd镜像
查看>>
浅析busybox查找命令和调用相应命令函数的实现流程框架
查看>>
利用linux dd和tr命令生成特定的数据
查看>>