线性代数笔记

版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/mobius_strip/article/details/79313105

1 矩阵

1.1 矩阵的定义及运算

1.1.1 矩阵的定义

矩阵是一个数表,由 m x n 个数排成的 m 行 n 列的数表

[ a 11 a 12 a 1 n a 21 a 22 a 2 n a m 1 a m 2 a m n ]

称为一个m行n列的矩阵, a i j 为矩阵第i行j列的元素。
m=n时矩阵为方阵;只有 a i i 不全为0的方阵叫对角阵; a i i =1的对角阵为单位矩阵,记做I;元素都是0的矩阵为0矩阵。

1.1.2 线性方程组

{ 2 x 1 + x 2 x 3 = 1 x 1 + x 3 = 2

对应的系数矩阵为:
A = [ 2 1 1 1 0 1 ]

对应的曾广矩阵为:
B = [ 2 1 1 1 1 0 1 2 ]

1.1.3 矩阵的运算

运算 说明
同型矩阵 行列对应相等
矩阵相等 对应元素相同的同型矩阵
矩阵加法 同行矩阵才有定义,对应不元素相加
矩阵数乘 所有元素乘以对应的实数:kA = (k a i j
矩阵乘法 A m × t B t × n = C m × n = ( c i j ) m × n ,其中 c i j = k = 1 t a i k b k j
矩阵等价 A通过初等变换得到B
矩阵相似 A = M 1 B M

1.2 矩阵的两个图

2x2矩阵

{ 2 x y = 0 x + 2 y = 3 [ 2 1 1 2 ] [ x y ] = [ 0 3 ] x [ 2 1 ] + y [ 1 2 ] = [ 0 3 ] { x = 1 y = 2

1.2.1 行图(row picture)

通过交点可以判断方程是否有解
这里写图片描述

1.2.2 列图(column picture)

列向量的线性组合
这里写图片描述

1.2.3 行列运算

[ c o l 1 c o l 2 c o l 3 ] [ a b c ] = a c o l 1 + b c o l 2 + c c o l 3 [ a b c ] [ r o w 1 r o w 2 r o w 3 ] = a r o w 21 + b r o w 2 + c r o w 3

1.3 消元法与LU分解

1.3.1 矩阵的初等变换

ID 初等行(列)变换
1 交换两行(列)的位置
2 用一非零数程昱某一行的所有元素
3 把矩阵的某一行(列)的适当倍数加到另一行(列)上去

经过初等变换后,将矩阵转化成阶梯形式,如果曾广矩阵有非零主元则方程无解;如果阶梯矩阵对应的系数部分主对角含零,则有无穷解;

初等矩阵:对单位矩阵做一次初等变换得到的矩阵;行变换等于左乘初等矩阵,列变换等价于右乘初等矩阵

1.3.2 LU分解

高斯消元法:对曾广矩阵实施行初等变换化为行(简化)阶梯形

通过消元法将矩阵A转换成矩阵B(阶梯,上三角),对应行变换左乘的初等矩阵记为: E 1 , . . . , E k 则有:
E k . . E 1 A = B A = E 1 1 . . E k 1 B L = E 1 1 . . . E k 1 A = L U

补:matlab求矩阵A的LU分解:rref(A)

1.3.3 线性方程组的解

在行变换的基础上加上交换行,化简后的矩阵为:

[ I F 0 0 ]
其中I为单位矩阵,F是自由变元的矩阵
特解: A X b = b
通解: A ( X b + X n u l l ) = b X n u l l

1.4 Gram-Schmidt 正交化

q i T q j = { 0 i j 1 i = j

Q = [ q 1 , . . , q n ] , Q T Q = I

例: Q 1 = 1 2 [ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ] , Q 2 = [ c o s ( θ ) s i n ( θ ) a i n ( θ ) c o s ( θ ) ]

步骤 说明
单位化 q 1 = a 1 | | a 1 | | , q 2 = b 1 | | b 1 | | , . . .
正交化 a = q 1 , b = q 2 a T b a T a a , . . .
整理 Q = [ a b . . ] A = Q R R = Q T A
注1:Q可以不是方阵,如果是方阵则 Q 1 = Q T
注2:Fourier变换基于正交基 [ s i n ( t ) , c o s ( t ) , . . ] [ e t , e i t , . . ]

1.5 分块矩阵

矩阵可以分块计算

1.6 Markov矩阵

定义:所有元素都不小于0,每列和为1。
性质: λ 1 = 1 | λ i | < 1 ( i 1 )
用途:可以用于表示人口的迁移。


2 矩阵的逆

2.1 逆矩阵定义

A B = B A = I B A A 1

等价命题
A可逆
AX = 0只有零解
A与I等价
A可以表示成有限个初等矩阵的乘积

2.2 Gass-Jordan

solve more equations at once:
( A | I ) ( I | A 1 )

A A x = b x = b / A = A 1 b

2.3 伪逆(sudo inverse)

伪逆 条件1 条件2 形式 结论
left inverse r=n N ( A ) = 0 ( A T A ) 1 A T 0或1个解
right inverse r=m N ( A T ) = 0 A T ( A A T ) 1 多解,n-m个自由变量

A ( A T A ) 1 A T 是在列空间中的投影
A T ( A A T ) 1 A 是在行空间中的投影

2.4 奇异值分解(SVD)

任何矩阵A都可以分解成: A = U Σ V T 其中 Σ 为对角矩阵
A A T = U Σ Σ T U T = U D U T
A T A = V Σ T Σ V T = V W V T
D m × m = [ σ 1 2 0 0 . 0 . 0 . 0 0 σ k 2 . . . . 0 ] W n × n = [ σ 1 2 0 0 . 0 . 0 . 0 0 σ k 2 . . . . 0 ]
对应的逆为: A + = V Σ + U T Σ + = [ 1 σ 1 0 0 . 0 . 0 . 0 0 . 1 σ k . . . . 0 ]


3 线性空间

3.1 线性相关和张成

定义 说明 例子
向量空间 < V,+,*>是定义了加法和数乘的代数系统 过原点直线,过原点平面
线性空间 满足叠加性:加性+齐性 f ( a x + b y ) = a f ( x ) + b f ( y )
线性无关 不能互相表示的向量构成的集合 笛卡尔坐标轴
张成(span) 一组向量的全部线性组合的集合
基(base) 张成整个空间V的线性无关向量组称为V的基 (1,0,0,…),(0,1,0,…),..
维度(dimention) 基的长度(或基的个数)
注:0向量是所有向量空间的子集
注:向量空间的交是向量空间,并不一定是

3.2 四个空间

名称 说明
列空间(column space) 矩阵列向量张成的空间
零空间(null space) A x = 0 的X的解构成的空间,记 N ( A )
行空间(转置的列空间) 矩阵行向量张成的空间
转置零空间(左零空间) A T y = y T A = 0 ,记 N ( A T )

这里写图片描述
如图row space 与 N(A) 正交补,column space 与 N ( A T ) 正交补,row space与column space 构成双射,其他部分都是零空间。

定理:dim(V)= dim(N)+dim(值域)

3.3 线性映射

每个矩阵 A m × n 对应一个从 R m R n 的线性映射,我们一般想要对矩阵A进行化简(特征值构成的对角矩阵),得到不同的基下的表示(分解成独立的子空间),可以简化问题或者看到问题的本质。

注:线性映射到矩阵的映射也是线性映射

3.3.1 特征值和特征向量(解耦、对角化)

特征向量在矩阵A的变化下只会出现缩放,不会出现方向变化:
A x = λ x λ = 0 A x = 0 x N ( A )
求解: A x = λ x ( A λ I ) x = 0
x不为0时,有 | A λ I | = 0 x i A λ i I

注: { A x = α x B x = β x ( A + B ) x = ( α + β ) x

3.3.2 特殊的特征值

矩阵 特征值 例子
负对阵 λ [ 0 1 1 0 ] λ 1 = i λ 2 = i v 1 = [ i 1 ] v 2 = [ i 1 ]
正对称 λ
畸形 λ n [ 0 1 0 0 ] λ = 1 v = [ 1 0 ]

3.3.3 特征值性质

i = 1 n λ i = i = 1 n a i i
i = 1 n λ i = d e t A

3.3.4 QR法

利用QR分解迭代求解特征值:
A i = Q i R i
R i Q i = A i + 1
迭代数次得到对角线矩阵,其中值为特征值

3.3.5 矩阵的多项式

f ( A ) = a 1 A k + . . . + a 1 A + a 0 I
f ( λ ) = a 1 λ k + . . . + a 1 λ + a 0

性质:f(A)g(A) = g(A)f(A)

4 行列式

4.1 定义

d e t A = | a 11 a 12 . . . a 1 n a 21 a 22 . . . a 2 n . . . a n 1 a n 2 . . . a n n |

n = 1 d e t A = a 11
n 2 d e t A = a 11 A 11 + a 12 A 12 + . . . a 1 n A 1 n , A i j = ( 1 ) i + j M i j
A i j M i j

4.2 对比矩阵和行列式

行列式 矩阵
数表
D n A m × n
| | (),[]

4.3 行列式的性质

ID 性质 推论
1 d e t I = 1 ( ) d e t A = 0
2 行列式按任意一行展开,其值相等: d e t A = a i 1 A i 1 + . . . + a i n A i n d e t A 0 d e t A = 0
3 d e t A d e t A = 0
4
| a 11 . . a 1 n . . . b i 1 + c i 1 . . b i n + c i n . . . a n 1 . . a n n | = | a 11 . . a 1 n . . . b i 1 . . b i n . . . a n 1 . . a n n | + | a 11 . . a 1 n . . . c i 1 . . c i n . . . a n 1 . . a n n |
5 A k d e t A 1 = k d e t A k d e t A d e t A d e t A = 0
6 d e t ( A B ) = d e t A d e t B d e t ( A T ) = d e t A

4.4 伴随矩阵

A 1 = C T d e t A
C T 为伴随矩阵, c i j
验证:

[ a 11 a 12 . a 1 n a 21 a 22 . a 2 n . . . a n 1 a n 2 . a n n ] [ c 11 c 12 . c 1 n c 21 c 22 . c 2 n . . . c n 1 a n 2 . c n n ] = [ d e t A 0 . 0 0 d e t A . 0 . . . 0 0 . d e t A ] = d e t A I

4.5 集合意义

行列式的值为向量组维成的超体的体积,特征向量方向的缩放比为对应的特征值。


5 应用

5.1 电路

这里写图片描述
令矩阵A的行为边,列为节点,出度为-1,入度为1,x为顶点(电压),y为边(电流)。

5.1.1 电压

A x = [ 1 1 0 0 0 1 1 0 1 0 1 0 1 0 0 1 0 0 1 1 ] [ x 1 x 2 x 3 x 4 ] = C [ y 1 y 2 y 3 y 4 y 5 ]

左侧为压降矩阵,右侧为欧姆定理,C为电阻矩阵。

5.1.2 电流

A T y = [ 1 0 1 1 0 1 1 0 0 0 0 1 1 0 1 0 0 0 1 1 ] [ y 1 y 2 y 3 y 4 y 5 ] = [ 0 0 0 0 ]

基尔霍夫定律,节点电流和为0;加入电流源 f = A T y = A T C A X A T C A

5.2 投影(projection)

A x = b b A A x ^ = p p b a

5.2.1 一维投影

这里写图片描述
p = x a x x e p e
e p = a T ( b p ) = a T ( b x a ) = 0
x = a T b a T a p = a x = a a T b a T a = a a T a T a b = P b

5.2.2 高维投影

因为e与A正交,故e N ( A T )
定义: p = x 1 a 1 + . . + x k a k = A x
正交: a i T e i = a i T ( b i A x i ) = 0
p = A x = A ( A T A ) 1 A T b = P b

5.2.3 投影性质

投影矩阵(P)的性质
P T = P
P 2 = P
I = P + E , b p = e = E b

5.3 最小二乘

数据拟合,求解到已知点距离和最小的目标曲线。
这里写图片描述

5.3.1 利用矩阵求解

设直线为C+Dx=y,带入三点坐标:

[ 1 1 1 2 1 3 ] [ C D ] = [ 1 2 2 ]

A T A x ^ = A T b x ^
这里写图片描述

5.3.2 利用导数求解

m i n | | A x b | | 2 = | | e | | 2 = e 1 2 + e 2 2 + e 3 2
f ( C , D ) = ( C + D 1 ) 2 + ( C + 2 D 2 ) 2 + ( C + 3 D 2 ) 2
f C = 0 6 C + 4 D = 0
f D = 0 3 C + 6 D = 0

微分方程 5.4

A x = λ x A S = Λ S A = S 1 Λ S
S Λ

5.4.1 通解

方程 通解
U k + 1 = A U k U k = C 1 λ k x x + . .
d u d t = A u u ( t ) = C 1 e λ 1 t + . .

5.4.2 稳定性

线代控制理论用矩阵表示系统,如上面第二种微分方程,指数中的虚部代表震荡,实部影响稳定性。

λ 稳定性
所有 λ i 1 李雅普诺夫稳定(不发散)
所有 λ i = 1 渐进稳定(随着时间收敛到0)
存在 λ i > 1 不稳定

5.5 基变换和图像压缩

图片(512x512 pixel) J P E G 64block(8x8 pixel) 系数c c ^ 0

变换基 举例 性质
standard [ 1 0 0 0 0 0 0 0 ] [ 0 1 0 0 0 0 0 0 ] .. [ 0 0 0 0 0 0 0 1 ] 效果不好
better 低频(全同): [ 1 1 1 1 1 1 1 1 ] 上下相反: [ 1 1 1 1 1 1 1 1 ] ..高频(棋盘): [ 0 0 0 0 0 0 0 1 ] 不好求逆
Fourier [ 1 1 1 1 1 1 1 1 ] [ 1 ω ω 2 ω 3 ω 4 ω 5 ω 6 ω 7 ] .. [ 1 ω 7 ω 14 ω 21 ω 28 ω 35 ω 42 ω 49 ] 效果好,好求逆
简化小波 [ 1 1 1 1 1 1 1 1 ] [ 1 1 1 1 1 1 1 1 ] [ 1 1 1 1 0 0 0 0 ] [ 0 0 0 0 1 1 1 1 ] [ 1 1 0 0 0 0 0 0 ] [ 0 0 1 1 0 0 0 0 ] [ 0 0 0 0 1 1 0 0 ] [ 0 0 0 0 0 0 1 1 ] 效果好,好求逆

猜你喜欢

转载自blog.csdn.net/mobius_strip/article/details/79313105