学习视频地址
文章目录
- 1.01方程组的几何解释
- 2.02矩阵消元
- 3.03乘法和逆矩阵
- 乘法
- 逆
- 4.04矩阵A的LU分解
- 5.05转置,置换,向量空间
- 置换
- 转置
- 向量空间
- 6.06列空间和零空间
- 7.07求解Ax=0:主变量,特解
1.01方程组的几何解释
对于二元方程组:row picture
方程转化为矩阵 * 向量 = 向量的形式,并且二维三维可以用图像几何表示出来(解析几何)
用column picture进行表示(线性组合):
x [ 2 − 1 ] + y [ − 1 2 ] = [ 0 3 ] x\begin{bmatrix}2 \\ -1 \\ \end{bmatrix}+y\begin{bmatrix}-1 \\ 2 \\ \end{bmatrix}= \begin{bmatrix}0 \\ 3 \\ \end{bmatrix} \\ x[2−1]+y[−12]=[03]
x=1,y=2,上述是用向量的角度进行求解。
任意两个不共线的非零向量能够表示平面中的任意向量
对于三元方程组:row picture
column picture(向量表示):
x=0,y=0,z=1
A x = b Ax=b Ax=b
是否对于任意b都有上式可解 = 列的线性组合是否能覆盖整个三维(n维)空间
非奇异矩阵/可逆矩阵可达到上面的要求(可逆矩阵被称为非奇异矩阵,因为可逆矩阵的行列式不为0,在数学中如果一个矩阵的行列式为0,则称这个矩阵为奇异矩阵。)
是3个向量都处于同一平面,即维度减一,所以涉及第三维度的向量b无法构造。
对于n个向量也是如果最后消解只有n-1维那也无解(一种平面只能覆盖n维空间的某n-1维平面,最后的求解也只能在这n-1维平面展开)。
方程组的矩阵形式:
A x = b ( A x 看做 A 各列的线性组合 ) Ax=b(Ax看做A各列的线性组合) Ax=b(Ax看做A各列的线性组合)
2.02矩阵消元
x + 2 y + z = 2 3 x + 8 y + z = 12 4 y + z = 2 A = [ 1 2 1 3 8 1 0 4 1 ] x+2y+z=2\\ 3x+8y+z=12\\ 4y+z=2\\ A=\begin{bmatrix}1 & 2 & 1 \\3 & 8 & 1 \\0 & 4 &1 \end{bmatrix} x+2y+z=23x+8y+z=124y+z=2A= 130284111
消元
增广矩阵:Ab
回代(也就是加入b进行求解)
方程从上往下消元,从下往上回代求解
消元矩阵:
由上式可知,我们把左边的矩阵(初等矩阵)当成3个1*3的行向量,右边的结果的第一行=左边第一行的第一个数 * 中间第一行 + 左边第一行的第二个数 * 中间第二行 + 左边第一行的第三个数 *中间第三行,所以结果为:1,0,0,以此类推第二行为-3,1,0(-3 * 中间第一行+1 * 中间第二行),第三行为0,0,1
E 21 表示初等矩阵第二行第一个位置发生了变换:变为 0 E_{21}表示初等矩阵第二行第一个位置发生了变换:变为0 E21表示初等矩阵第二行第一个位置发生了变换:变为0
E 32 ( E 21 A ) = U 相当于 ( E 32 E 21 ) A = U 结合律 E_{32}(E_{21}A)=U \\ 相当于(E_{32}E_{21})A=U\\结合律 E32(E21A)=U相当于(E32E21)A=U结合律
置换矩阵P
行变换左乘,列变换右乘
矩阵中基本没有交换律
逆矩阵:从U变回A
我们要取消这次消元,我们就希望E’EA=A,所以让E’E等于单位矩阵即可,E’是逆矩阵
变换后的矩阵EA,如何变回A?以行变换为例(左侧乘),就是找一个矩阵X使得XEA为A,此时若XE为单位阵,则XEA为A
E − 1 E = I E^{-1}E=I E−1E=I
3.03乘法和逆矩阵
乘法
方法一:行与列之间的点积
C 34 = ( r o w 3 o f A ) ⋅ ( c o l u m n 4 o f B ) C_{34}=(row3ofA) ·(column4ofB) C34=(row3ofA)⋅(column4ofB)
A是m * n,B是n * P,C=AB大小为m * P
方法二:把B看做是P个单独的列向量,像上一节课一样对每个列向量单独乘积最后再合并:
方法三:把A看做m个单独的行向量,向上节课一样对每个行向量单独乘积最后再合并:
方法四:A中的列 * B中的行
矩阵分块:把每个块都看成单独的元素运用方法一就得到下图
逆
逆矩阵(非奇异矩阵) * 原矩阵 = I(单位矩阵)
奇异矩阵(非逆)可以找到一个向量x 使得Ax=0,也就是A中的向量组共面
Ax=0 有非零解,说明A中有的列没有贡献,所以不可逆
反证法:逆乘原得到单位矩阵,单位矩阵乘任一矩阵得原矩阵,所以得到X=0,但是X≠0,所以不成立
如何求逆矩阵:对于下列矩阵本质上就求两个二元一次方程组
同时求解两个方程:增广矩阵形式
高斯-若尔当消元法:
原理可理解为,假定有B矩阵,BA=I,那么右乘A逆得BI=A逆,如此可以得到逆矩阵
相当于对AI同时乘以A逆,得到I和A逆
A I = = > A ∗ A ′ I A ′ = = > I A ′ AI==>A*A^{\prime}I{\,}A^{\prime}==>IA^{\prime} AI==>A∗A′IA′==>IA′
4.04矩阵A的LU分解
两矩阵相乘,且它们的逆均已知,那么AB的逆是什么 B − 1 A − 1 B^{-1} A^{-1} B−1A−1
A转置乘以A逆的转置等于单位E,定理A转置乘以A转置的逆等于单位E,所以A逆转置等于A转置的逆
如果还是不明白,因为我们已经知道A * A逆=A逆 * A=E,所以下面那个式子还可以写成,A转置的逆*A的转置=E,逆矩阵又是唯一的,所以A逆的转置=A转置的逆
这个由“行”组成的系数矩阵的每一行应该和之前A-1代表的系数矩阵的“列”相对应,因此就是(A-1)T*AT
将原本的A*A-1=I中左侧的A中的所有列看作用于线性组合的列向量,A-1中的每一列代表每次线性组合的系数组成
A A − 1 = I ( A − 1 ) T A T = I ( A − 1 ) T = ( A T ) − 1 AA^{-1}=I\\ (A^{-1})^TA^T=I\\ (A^{-1})^T=(A^{T})^{-1} AA−1=I(A−1)TAT=I(A−1)T=(AT)−1
把-4倍的第一行和1倍的第二行相加,放到第二行
E 21 = [ 1 0 − 4 1 ] E_{21}=\begin{bmatrix}1 & 0 \\-4 & 1 \end{bmatrix} E21=[1−401]
L = E 21 − 1 L=E_{21}^{-1} L=E21−1
E 32 E 31 E 21 A = U ( n o r o w e x c h a n g e ) A = E 21 − 1 E 31 − 1 E 32 − 1 U E_{32}E_{31}E_{21}A=U(no{\,\,} row{\,}exchange)\\ A=E_{21}^{-1}E_{31}^{-1}E_{32}^{-1}U E32E31E21A=U(norowexchange)A=E21−1E31−1E32−1U
[ 1 0 0 2 1 0 0 5 1 ] = L ( l e f t o f U ) E A = U A = L U \begin{bmatrix}1 &0 & 0\\2 & 1&0\\0&5&1 \end{bmatrix}=L(left{\,}of{\,}U)\\ EA=U A=LU 120015001 =L(leftofU)EA=UA=LU
原因:然而L是恢复为original矩阵的逆过程,也就是从下向上处理的,同样使用较上的行来处理较下的行,然而自下而上不会更改较上的行。
原因:E是从上向下处理的,使用较上的行来处理较下的行,然而较上的行不一定是original的,有可能是经过处理修改的。
-
L中的矩阵相乘的顺序很好,从左往右从上往下
-
初等矩阵都是单位阵来的,对角线都是1
类似于时间复杂度如果n=100,对于n * n的矩阵A需要进行多少次操作(每个元素一次乘法+一次减法进行小消元)变成下三角矩阵?
首先把矩阵留下第一行第一列的主元,第n行第一列为0,需要大约100^2次运算
同理,把矩阵留下第二行第二列的主元,第n行第二列为0,需要大约99^2次运算…
总操作数是:大约是 1 / 3 ∗ n 3 1/3*n^3 1/3∗n3的总操作次数
a l l o p r a t i o n s = n 2 + ( n − 1 ) 2 + . . . 1 2 = n ( n + 1 )( 2 n + 1 ) / 6 all{\,}oprations=n^2+(n-1)^2+...1^2=n(n+1)(2n+1)/6 alloprations=n2+(n−1)2+...12=n(n+1)(2n+1)/6
对于增广矩阵右侧的b,需要n^2次操作次数
n * n的矩阵有 n ! n! n!种置换矩阵
P − 1 = P T P^{-1}=P^{T} P−1=PT
5.05转置,置换,向量空间
置换
permutations P:execute row exchanges
A=LU的分解其实就是L是A行变换的过程U是行变换的结果
如果要进行 行变换A=LU 就要变成PA=LU
置换矩阵是行重新排列了的单位矩阵
n * n的矩阵有n!种置换矩阵,就是各行重新排列后所有可能的数目
P都是可逆的,并且满足以下式子
P − 1 = P T P T P = I P^{-1}=P^{T}\\ P^{T}P=I P−1=PTPTP=I
转置
R^TR 得出来的矩阵都是对称的
( R T R ) T = R T R T T = R T R (R^TR)^T=R^{T}R^{TT}=R^TR (RTR)T=RTRTT=RTR
向量空间
向量空间对数乘和加法两种运算(线性组合)是封闭的
子空间内必须包含0向量
二维空间的子空间是过原点的一条直线
这一讲的核心思想是,通过某些向量构成一个向量组成的空间如果这些空间属于R3,那么它们构成的空间也在R3内,关键是,对其进行线性组合后仍然在子空间内
子空间:过原点、点(1,2,4)、点(3,3,1)的平面
如何构造矩阵列空间:取出每一列,线性组合,所有的线性组合就能构成列空间
6.06列空间和零空间
子空间S和T的交集是一个子空间
交集必含有零向量,线与线的交集是原点,线与面的交集要么是原点要么是线本身,面与面的交集要么是一条过原点的线要么是面本身
S和T的交集取向量v和w,则v,w既属于S也属于T,而S,T中的v,w是满足“加和数乘”规则的
三个列向量的线性组合无法充满整个四维空间
那么理解:四维空间里的三个向量最多只能组成四维空间里的三维空间,而要让B有解,就必须充满四维空间才行,因为B是四维空间里任意一个向量
问题:什么b让方程组有解:Ax=b
只有b是各列的线性组合时,Ax=b才有解
对于A,去掉某一列得到同样的列空间:去掉列三,因为列一和列二的线性组合能得到列三,也就是说列三没有贡献,列一列二为主列。(同理也可以去掉列一)
A的零空间是指是的Ax=0成立的向量空间
C ( A ) i n R 4 , x = [ x 1 x 2 x 3 ] i n R 3 C(A) in R^4, x=\begin{bmatrix}x_1 \\ x_2\\x_3 \end{bmatrix}in R^3 C(A)inR4,x= x1x2x3 inR3
Ax=0包括 c [ 1 1 − 1 ] c\begin{bmatrix}1 \\ 1 \\-1 \end{bmatrix} c 11−1 的一条直线
如果b不为零,Ax=b中x的解不构成子空间(0不是解),它是不过原点的一条直线。这条直线是 [c+1, c, -c],过(1,0,0)