zoukankan      html  css  js  c++  java
  • 线性代数入门

    前言

    某次模拟赛被矩阵虐哭,补一波线代

    这篇博客偏入门,概念较多,算法相关较少

    大力膜拜(3B1B)线性代数的本质系列

    (参考资料来源,或者干脆叫观影总结吧……)

    完全就是观影总结(qwq)

    记号:不作特殊说明,本文中的大写字母均表示某个矩阵,小写字母均表示某个向量

    顺便一提,如果对标题分布有点混乱,那么建议看一眼博客右侧的分级目录

    向量

    向量的定义

    高考课本讲的向量大概是:具有长度和方向的箭头

    (OI)中向量大部分情况下是:一个(n*1)的矩阵,或者说具有实际含义的一列数字 (egin{bmatrix}1\1\4\5\1\4end{bmatrix})

    然而实际上只要满足两个事物相加或相乘有意义就可以称之为向量,后面再详细介绍这里

    但是这提示了我们:向量相加和相乘在线性代数中起着极为重要的作用

    为了更适合入门,我们下面的讲解都会优先考虑几何意义和二维平面

    在二维坐标系中,我们定义向量为从原点出发有长度和方向的箭头,虽然理论上向量只由长度和方向决定,与起点无关,但在线代中向量通常以原点为起点

    例如:

    这里有一个我们通常意义上的向量

    那么这个箭头和我们的另一种形式:一列数字有什么关联呢?

    我们发现一个向量坐标由一对数字构成,这对数字指引你如何从原点(它的起点)出发,到达它的尖端

    如上图中,我们可以通过尖端(B)的坐标((3,2))了解到,我们可以先由原点向(x)轴正方向走(3)个单位长度,再向(y)轴走(2)个单位长度得到

    如果我们规定原点是一切向量的起点,那么只要一对数字就可以唯一确定一个二维向量

    上图中的箭头可以表示为(egin{bmatrix}3\2end{bmatrix}) 为了与坐标区分,我们采用竖写加方括号

    我们应该可以脑补出三维的场景,这里我就不画图了我实在画不出来了

    向量加法

    如图地我们由两个向量(u)(v)

    对于求出(w=u+v)的结果,只需要将(v)的起点移动到(u)的尖端,然后由(u)的起点指向(v)的尖端

    这是线性代数中为数不多的允许向量的起点离开原点的情况

    为什么这样定义向量加法?

    因为我们可以将一个向量看作一段运动,可以发现在上图中,我们先按照(u)向量所描述的方式运动,再按照(v)向量所描述的方式运动,最后与直接按照(w)向量所描述的方式运动是等价的

    它也刚好对应了向量所代表的数值上的加法运算:

    (egin{bmatrix}3\1end{bmatrix}+egin{bmatrix}2\-2end{bmatrix}=egin{bmatrix}3+2\1-2end{bmatrix}=egin{bmatrix}5\-1end{bmatrix})

    向量数乘

    对于一个向量

    我们将它乘以一个标量(2)得到:

    效果就是将长度伸长为原来的两倍

    如果我们将它乘以(frac{1}{2})呢?

    效果就是先反向,然后长度缩小为原来的(frac{1}{2})

    实际上自始自终,标量在线性代数中的主要作用就是缩放向量

    从数字角度看,每个向量数乘即将每一个分量乘以这个数字

    (2*egin{bmatrix}4\2end{bmatrix}=egin{bmatrix}2*4\2*2end{bmatrix}=egin{bmatrix}8\4end{bmatrix})

    线性组合

    线性组合:

    两个向量通过改变数乘的标量而产生新的向量

    我们可以从另外一个角度去理解一个向量的坐标:对于一个向量的放缩

    对于向量(v)来说,我们可以看作是由(i)伸长为原来的(3)倍,(j)反向伸长为原来的(2)倍再相加得到的

    此时向量(v)的横纵坐标作为一个标量对两个向量进行放缩,放缩后的向量通过加和得到了向量(v)

    其中(i)(j)被称为坐标系的基,或者基向量,它表示它是某个的向量坐标放缩并加和的对象

    则二维空间中所有的向量都是由两个基向量线性组合得到的

    张成的空间

    两个向量线性组合得到的所有向量的集合称为张成的空间

    我们发现,所有的向量都可以由基向量组合表示,也就是说如果我们选取了不同的基向量那么将会产生不同的坐标系

    比如我们选取了下面的两个向量作为基向量

    那么我们仍然可以通过基向量的线性组合产生二维空间中所有的向量,虽然对于同一组向量来说它们的坐标表示也许不太相同

    那么我们称这组新的基向量(u,v)张成的空间为二维平面

    我们考虑一种特殊情况

    两个向量只能组合出与当前向量共线的全部向量,我们称(i,j)张成的空间为一条直线

    对于这种特殊情况,我们概括为:某一个向量的增加对于张成的空间没有添加新的维度,那么我们称这组向量是线性相关的

    或者说其中一个向量是由其他向量线性组合得到的,则这组向量线性相关

    对应地,如果我们每个向量都对张成的空间增加了新的维度,那么我们称这组向量线性无关

    还有更特殊的情况,比如两个向量都是原点,那么张成的空间就是一个点

    然而注意到在特殊情况中我并没有称两个向量为基向量,因为基向量有严格的定义:

    向量空间中的一组基是张成该空间的一个线性无关向量集

    线性变换

    线性变换定义

    变换其实和函数的定义类似,不过变换是输入一个向量,再输出一个向量

    对于一个空间的变换其实由很多,比如拉伸,扭曲,旋转等等,但是我们考虑的只有线性变换

    那么什么才是线性的变换呢?

    1.直线在变换后仍然是直线

    2.原点保持固定

    我相信你们能够脑补出什么是线性变换,什么不是线性变换,因为这个东西我并没有找到合适的画图工具展现出来……

    但是对于一个经典的例子我由必要手绘一下

    对于这样一个变换,它并不是线性的,因为某些我们没有展示在图中的线发生了扭曲

    可以自行脑补一下原图中(F-N-R)这条直线变换之后的形状(草老子画不出来啊)

    我们其实可以认为,一个变换如果保持网格线之间平行并且等距,那么这个变换是线性变换

    那么我们如何通过数值去描绘一个线性变换?

    考虑到基向量的特殊性,我们只需要描述线性变换之后基向量在原空间中的坐标,就可以表示出这个变换的效果

    同时线性变换有重要的推论:对于一个向量(v),它在变换前和变换后都是基向量的特定的线性组合得到的

    即 我们设(d^{'})是向量(d)经过某个线性变换的向量,那么如果

    (v=x*i+y*j)

    (v^{'}=x*i^{'}+y*j^{'})

    其中(x,y)是标量

    举个例子,一个向量((2,4)),如果经过了某个线性变换,基向量由((1,0))((0,1))变为了((2,1))((1,2))

    那么这个向量在变换后的坐标应该为(2*(2,1)+4*(1,2)=(8,10))

    线性变换与矩阵

    由上面可以得出,对于任意一个线性变换,只要记录了变换之后的基向量,和原向量,就可以得到变换之后的新向量

    如果我们的基向量变为了((2,1)和(1,2)),那么对于任意向量(egin{bmatrix}x\yend{bmatrix}),我们可以通过

    (egin{bmatrix}x\yend{bmatrix}^{'}=x*egin{bmatrix}2\1end{bmatrix}+y*egin{bmatrix}1\2end{bmatrix})得到变换后的向量

    我们通常把基向量((a,c),(b,d))封装为一个(2*2)的矩阵,将变换过程称为矩阵乘法

    (egin{bmatrix}a b\c dend{bmatrix} egin{bmatrix}x\yend{bmatrix}=egin{bmatrix}x\yend{bmatrix}^{'})

    在这里一个(2*2)的矩阵即表示一个线性变换,其中每一列是变换后的一个基向量

    我们仔细分解一下矩阵乘法的过程

    (egin{bmatrix}a b\c dend{bmatrix} egin{bmatrix}x\yend{bmatrix}=x*egin{bmatrix}a\cend{bmatrix}+y*egin{bmatrix}b\dend{bmatrix}= egin{bmatrix}x*a+y*b\x*c+y*dend{bmatrix})

    现在也许能更清楚地理解为什么矩阵乘法要这样定义了

    试着用矩阵描述一些变换:

    把空间逆时针旋转(90^{。})

    意味着:(egin{bmatrix} 0 -1\1 0end{bmatrix})(这东西怎么这么丑啊(qwq)

    有一个特殊的概念:如果我们矩阵中的几个向量是线性相关的,会导致你遭受降维打击变换之后空间维度降低

    线性变换复合与矩阵乘法

    有时候我们需要连续描述两个或多个线性变换,但是我们仍然可以追踪基向量的变化来实现对变换的描述

    不过我们可能需要借助一些工具来描述,比如矩阵

    假设我们要对一个空间进行两次变换:先旋转再剪切

    我们先将空间逆时针旋转(90^{。})得到

    (egin{bmatrix} 0 -1\1 0end{bmatrix}egin{bmatrix}x\yend{bmatrix})

    再剪切(想要知道为什么这个矩阵叫做剪切可以自己模拟一下把基向量变成矩阵的两列)

    (egin{bmatrix}1 1\0 1end{bmatrix}igg(egin{bmatrix} 0 -1\1 0end{bmatrix}egin{bmatrix}x\yend{bmatrix}igg))

    那么它将等价于

    (igg(egin{bmatrix}1 1\0 1end{bmatrix}egin{bmatrix} 0 -1\1 0end{bmatrix}igg)egin{bmatrix}x\yend{bmatrix}=egin{bmatrix} 1 -1\1 0end{bmatrix}egin{bmatrix}x\yend{bmatrix})

    我们称左边的矩阵为等价矩阵

    考虑在几何中我们两个变化的描述应该为,基向量先变为了((0,1),(-1,0)),又在新的空间中变为了((1,0),(1,1)),那么以左边两个矩阵的乘法作为两个变换连续的效果是很合理的(矩阵中描述变换从右向左描述)

    如何计算两个矩阵相乘?

    还记得我们矩阵的含义吗,每一列代表了一个向量

    也就是说对于

    (egin{bmatrix}e f\g hend{bmatrix}egin{bmatrix}a b\c dend{bmatrix}),我们可以看作右边的矩阵的每一列代表的向量做左边的变换

    (egin{bmatrix}e f\g hend{bmatrix}egin{bmatrix}a b\c dend{bmatrix}= egin{bmatrix}e f\g hend{bmatrix}egin{bmatrix}a\cend{bmatrix}+egin{bmatrix}e f\g hend{bmatrix}egin{bmatrix}b\dend{bmatrix} =egin{bmatrix}ae+bg af+bh\ce+dg cf+dhend{bmatrix})

    中间应该还有一步向量乘变换,懒得展开了

    当然我们应该可以通过这个过程发现矩阵乘法并不满足交换律

    但是它为什么满足结合律呢?

    考虑矩阵变换是从右向左描述的,所以……矩阵乘法中加不加括号在几何意义上根本没有任何变化

    高维的情况应该可以自己推出来,但是实在是难以表述&作出图,所以就不在这里写了

    行列式

    在线性变换中,也许我们会想要知道一个问题:在经过了某个线性变换之后,原空间受到了多少拉伸或者挤压?

    或者转化成一个在二维平面中等价的问题:一个给定区域的面积大小的变化

    行列式定义

    一个线性变换对“体积”所造成的影响。

    当然由于是在二维平面中的线性变换,所以我们在研究这个问题的时候可以只研究一个特殊部分:两个基向量围成的四边形的面积的变化比例

    比如说:一个(2*2)矩阵的行列式为(3),那么说明(i,j)两个基向量四边形的面积变为原来的三倍

    特殊地:有些矩阵的行列式为(0),例如:(detigg(egin{bmatrix}1 3\0 0end{bmatrix}igg)=0)

    这表示在经过这个线性变换之后空间将降维

    另外还有(detigg(egin{bmatrix}0 1\1 0end{bmatrix}igg)=-1)

    虽然在经过(egin{bmatrix}0 1\1 0end{bmatrix})这个变换之后(i,j)两个基向量之间的四边形面积仍然是(1),但是我们会发现,(i)向量变到了(j)向量的左边,两个向量的相对位置交换了,空间的定向发生了改变,所以得到的行列式值为负

    类比可以推出:在三维空间中,一个线性无关方程组的行列式的值代表的是三个向量确定的平行六面体的体积

    如果在三维情况下行列式的值为负呢?

    这里有一个用来描述三维空间空间定向的方法,叫做 “右手定则”

    在一个三维空间中,我们假设基向量为(i,j,k),分别对应(x,y,z)

    一般情况下,我们右手食指指向(i),中指指向(j),大拇指指向(k)

    如果在某个线性变换之后,你做到这一步需要变成左手,那么说明这个三维空间的定向发生了改变,此时行列式的值应该为负

    对于行列式的计算我打算单独拿出来讲,因为本篇文章偏入门,这里只介绍行列式的定义和概念

    不过由于后面的需要,我需要指出对于(2*2)的矩阵(egin{bmatrix}a b\c dend{bmatrix}),它的行列式为(ad-bc)

    如果你真的理解了行列式的具体含义,你应该可以很轻易的理解下面这个等式

    (M_i)为某个线性变换矩阵

    [det(M_1M_2)=det(M_1)det(M_2) ]

    线性方程组

    在某些情况下,方程组会具有一个非常特殊的形式:

    在每个方程中,未知量只具有常系数,未知量之间只进行加和操作

    [ left{ egin{aligned} 2x+5y+3z=-3\ 4x+0y+8z=0\ 1x+3y+0z=2 end{aligned} ight. ]

    整理这种方程组通常是将所有未知量放在侧坐,将相同的未知量放在同一列对齐(有可能需要补(0)完成这一步),其他常量放在等式右侧

    我们称这种方程组为线性方程组

    其实我们会发现这种方程组对应着一个矩阵的形式

    (egin{bmatrix}2 5 3\4 0 8\1 3 0end{bmatrix}egin{bmatrix}x\y\zend{bmatrix}=egin{bmatrix}-3\0\2end{bmatrix})

    我们称左侧的系数矩阵为(A),包含未知量的矩阵为(x),右侧的结果矩阵为(v)

    则可以简写为(Ax=v),结合前面的知识,我们已经可以发现这其实是一个线性变换了

    那么这个求解线性方程组就具备了它的几何意义:对于一个线性变换后空间中的向量,找到在原空间的对应向量

    那么对于这个问题我们需要稍微分类讨论一下了,因为我们面临两种情况:(det(A)=0)(det(A) e 0)

    先来考虑似乎较为简单的(det(A) e 0)

    那么说明原空间经过(A)之后并没有降低维度,原空间中应该只有一个向量与(v)对应(这个不用解释吧应该?)

    但是我们只有变换之后的向量,如果要追寻原本的向量,我们应该将(v)做一个(A)的逆变换,记作(A^{-1})

    所谓逆的含义应该在学乘法逆元的时候都有所了解了,大概可以总结为(A^{-1}*A=E),其中(E)为恒等变换

    对于(n)维空间,(E)是一个对角全部为(1),其它区域全部为(0)(n*n)的矩阵

    所以可以推出:

    [A^{-1}Ax=A^{-1}v ]

    [x=A^{-1}v ]

    但是特殊的(det(A)=0)怎么办?

    经过(A)之后,空间会直接降维,此时(A)是不存在逆矩阵的,因为没有一种变换能使空间维度提升

    但是解仍然可能存在,如果你的(v)恰好处于降完维之后的空间之中

    至于解是什么,就不打算在本篇说明了(雾)

    经过线性变换之后空间的维度,我们称为秩

    比如一个平面,我们将他旋转,那么旋转之后它的秩为(2)

    但是如果经历了一个(egin{bmatrix}1 3\0 0end{bmatrix}),它的秩为(1)

    可以发现对于每一个降智降秩的变换,矩阵的行列式值都为(0)

    列空间

    对于一个线性变换(A)来说,所有可能的向量(Av)的集合,我们称为(A)的列空间

    当然这个名字的来源很好理解……因为向量(Av)的集合其实就是构成矩阵的每一列的向量的所有线性组合成的向量

    所谓列空间,其实就是列向量张成的空间,所以秩更准确的定义应该是列空间的维数

    当秩达到最大时,与矩阵的列数相等,称之为满秩

    零空间

    值得注意:(0)向量必定在列空间中,因为线性变换必须保持原点不动

    但是(0)向量并不一定只是原点,举个例子,在某个二维变换(A)中,若(det(A)=0),那么某条直线上的所有向量将会全部压缩到原点

    在某个线性变换(A)中,所有变换后落在原点的向量的集合,称为零空间(或者核)

    所以对于(Ax=v)中,如果(v)是零向量,那么所有(x)的解将是整个零空间

    非方阵

    目前为止,我们讨论的变换方阵都是(n*n)的方阵,那么对于非方阵我们应该怎样理解?

    例如一个(3*2)的矩阵(egin{bmatrix}1 3\4 1\2 3end{bmatrix})

    它说明我们在二维空间中的两个基向量(i,j)在经过某个线性变换之后对应到了某个三维空间中的二维平面

    注意它对应到的是某个三维空间中的二维平面,而不是某个三维空间,没有一种变换能让空间升维

    或者令一种(1*2)的矩阵(egin{bmatrix}1 2end{bmatrix})

    意味着某个二维空间经过线性变换之后降低到了一维,但是它仍然有两个基向量

    这对点积有着重要的意义

    点积

    点积定义

    我们来回忆一下高中数学中的点积:

    代数意义:对于两个行数相同的矩阵(egin{bmatrix}a\b\cend{bmatrix}·egin{bmatrix}e\f\gend{bmatrix}=ae+bf+cg),当然不一定是三维的

    几何意义:(u·v)应该是(v)(u)上的投影的长度乘以(u)的模长

    (u·v)=((v)的投影长度)*((u)的模长)

    如果(v)的投影方向与(u)相反那么点积的值应该为负

    现在你有没有一个问题:这两种解释有什么关系??????

    对偶性

    根据线性变换的规定,如果我们有一个高维对一维的线性变换,那么对于原空间中等距分布在一条直线上的一些点,应用完线性变换之后在数轴上仍然等距

    然而对于高维对一维的线性变换,我们所有的基向量都会落在一个数轴上,也就是说会产生一个(1*n)的变换矩阵,例如(egin{bmatrix}2 1end{bmatrix})

    考虑它对原空间某个向量(egin{bmatrix}x\yend{bmatrix})的影响(2*x+1*y)

    等价于矩阵乘法中(egin{bmatrix}2 1end{bmatrix}egin{bmatrix}x\yend{bmatrix})

    然而我们发现这个过程貌似和点积是运算非常相似!

    这个(1*n)的矩阵不就是一个躺倒的向量吗?

    但是线性代数中显然没有躺倒或者直立这种操作……但是这恰好说明一个高维对一维的线性变换和一个高维向量之间存在在某种对应关系

    以二维来举例子

    我们在二维空间中放置一条数轴,其基向量为(u)

    考虑某一种线性变换:对于任意向量(v),求出这条向量在数轴上的投影

    那么我们在描述变换的时候应该考虑的就是,变换后(i,j)的位置

    这里有就非常有趣的东西了!

    我们作(k)是角(GCD)的角平分线,会发现(GC=HC),即(i)向数轴的投影就是(u)(x)轴的投影,而(u)(x)轴的投影就是(u)的横坐标!

    同理(j)向数轴的投影就是(u)的纵坐标

    所以这个线性变换完全可以表示为(egin{bmatrix}u_x u_yend{bmatrix})

    而对于一个数轴非基向量(u),我们可以看作基向量的线性变换再乘以(u)的模长

    这就是为什么对于一个(v)(u)上的投影的长度乘以(u)的模长这样的操作,我们用矩阵来描述就是

    (egin{bmatrix}u_x u_yend{bmatrix}egin{bmatrix}v_x\v_yend{bmatrix})

    而它恰恰与(egin{bmatrix}u_x\u_yend{bmatrix}egin{bmatrix}v_x\v_yend{bmatrix})在数值上是一致的!

    神奇的对偶性

    叉积

    叉积定义

    两个向量的叉积,即两个向量围成的平行四边形的面积(虽然两个向量是两条线段,但是你们都懂我意思对吧~)

    对于(v×u),如果(v)(u)右侧,那么结果为正,否则结果为负

    所以说顺序对叉积的结果是有影响的

    (v×u=-u×v)

    也许你会发现这和行列式有一点相似?这是行列式计算的是基向量围成的面积

    那么我会告诉你,(v×u=detigg(egin{bmatrix}v_x u_x\v_y u_yend{bmatrix}igg))

    但是严格地说它并不是叉积……

    叉积严格的定义是由两个三维向量生成一个新的三维向量

    其长度为面积的值,方向垂直于两个向量所在的平面

    不过是向哪边垂直呢?我们需要右手定则

    右手食指指向(v),中指指向(u),大拇指就是叉积的方向

    我们是怎么求一个叉积的?

    一般计算的时候也许会有个老师告诉你:

    计算两个向量(egin{bmatrix}v_1\v_2\v_3end{bmatrix}×egin{bmatrix}u_1\u_2\u_3end{bmatrix}= detigg(egin{bmatrix}i v_1 u_1\j v_2 u_2\k v_3 u_3end{bmatrix}igg))

    (可能在某些书中或者博客中矩阵是横过来的,不过结果是一样的)

    右边矩阵的第一列居然是一组基向量?

    如果你头铁继续往下算,会得到(i(v_2u_3-v_3u_2)+j(v_1u_3-v_3u_1)+k(v_1u_2-v_2u_1))

    虽然你不知道为什么向量能够放到矩阵中去算行列式,也不知道为什么这样算出来的结果是叉积,但是也许你有很多种方法验证发现得到的结果的确是(v)(u)的叉积

    我们需要一个解释

    对偶性

    这里我们需要用到对偶性

    回忆一下对偶性:

    任何一个高维空间向一维空间的线性变换,都能在高维空间中找到一个向量唯一对应,使得某个向量应用这个线性变换等价于于这个高维向量做点积

    数值上说是因为这个任何一个向量与这个线性变换做矩阵乘法和对应的高维向量做矩阵乘法结果都是一样的

    考虑一下如果我没有告诉你什么是真正的叉积,并且你预先也没有了解,然后在看完二维空间中叉积的计算之后,会怎样考虑三维空间的叉积?

    那么大概率会认为是三个向量的行列式

    然而事实上叉积是接收两个向量并输出一个向量

    那么我们假设有一个函数(figg(egin{bmatrix}x\y\zend{bmatrix}igg)=detigg(egin{bmatrix}x v_1 u_1\y v_2 u_2\z v_3 u_3end{bmatrix}igg))

    这个函数是线性的(线性定义:(f(a+b)=f(a)+f(b),f(ka)=kf(a))其中(k)是标量)

    这说明它对应着一个三维对一维的线性变换

    那么(f)函数其实可以写成

    (egin{bmatrix}p_1 p_2 p_3end{bmatrix}egin{bmatrix}x\y\zend{bmatrix}=detigg(egin{bmatrix}x v_1 u_1\y v_2 u_2\z v_3 u_3end{bmatrix}igg))

    根据对偶性,它等价于

    (egin{bmatrix}p_1\p_2\p_3end{bmatrix}·egin{bmatrix}x\y\zend{bmatrix}=detigg(egin{bmatrix}x v_1 u_1\y v_2 u_2\z v_3 u_3end{bmatrix}igg))

    考虑这个式子的意义:一个向量(p)与一个未知向量(x)的点积等于(x,v)(u)的行列式

    我们已经知道在三维空间中,三个线性无关向量的行列式是一个平行六面形,那我们考虑点积的几何含义:(p)(x)上的投影乘以(x)的模长(交换(x,p)没有关系)

    平行六面体的体积求法应该是底面积乘以第三个向量在垂直于底面上的投影

    那么如果我们把(p)看作第三个向量,那么(x)就应该满足:长度等于平行六面体的底面积,且垂直于底面

    在上面的等式中底面的边是(v,u),这个(x)完美符合(v,u)叉积的定义!

    所以说(v×u=detigg(egin{bmatrix}i v_1 u_1\j v_2 u_2\k v_3 u_3end{bmatrix}igg))

    解出以后的(i(v_2u_3-v_3u_2)+j(v_1u_3-v_3u_1)+k(v_1u_2-v_2u_1))作为一个向量,就是(egin{bmatrix}x\y\zend{bmatrix})

    (upd:)叉积似乎还有更高维的定义,不太了解

    基变换

    不同基向量变换

    在前面我们提到过,对于一个向量的坐标,我们可以解读为对于基向量的放缩

    例如某个向量坐标为((3,2)),我们可以理解为将(i)伸长(3)倍,再将(j)伸长(2)倍,然后通过相加得到向量(egin{bmatrix}3\2end{bmatrix})

    通常意义上我们的(i,j)分别选取((1,0),(0,1))

    但是如果选取不同的基向量呢?

    比如有一个帅比叫动动,它选取了((2,1),(-1,1))两个向量作为基向量,称之为(t1,t2)

    那么他对于(egin{bmatrix}3\2end{bmatrix})应该是(egin{bmatrix}frac{5}{3}\frac{1}{3}end{bmatrix})

    因为在他看来,这个向量是由(frac{5}{3}t1+frac{1}{3}t2)得到的

    但是这个有一个奇特的地方

    在动动的视角中,他认为的基向量是((1,0),(0,1)),所谓的((2,1),(-1,1))只是以我们的视角来命名的

    所以实际上我们的语言是不同的,对于同一向量,我们叫它(egin{bmatrix}3\2end{bmatrix}),而动动叫它(egin{bmatrix}frac{5}{3}\frac{1}{3}end{bmatrix})

    问题就来了

    由于动动是帅比,你非常想和他做朋友,但是你们说着不同的语言,所以我们要学会不同语言之间的转化(不同坐标系之间的转化)

    具体地说,对于动动说出的一个向量(egin{bmatrix}-1\2end{bmatrix}),在我们的语言中它应该如何表示?

    在我们的语言中,(t1=egin{bmatrix}2\1end{bmatrix},t2=egin{bmatrix}-1\1end{bmatrix}),所以向量(v=-1egin{bmatrix}2\1end{bmatrix}+2egin{bmatrix}-1\1end{bmatrix} =egin{bmatrix}-4\1end{bmatrix})

    所以我们所说的(egin{bmatrix}-4\1end{bmatrix})就是动动所说的(egin{bmatrix}-1\2end{bmatrix})

    结合这个形式,很容易发现这其实是一个线性变换的过程

    当然,这只是我们读懂了动动的语言,但是如果要做♂朋♂友,那么显然动动也要读懂我们的语言

    例如我们所说的(egin{bmatrix}3\2end{bmatrix}),动动如何知道在它的语言中是(egin{bmatrix}frac{5}{3}\frac{1}{3}end{bmatrix})表示

    观察刚才我们理解动动语言的过程,实际上动动理解我们的语言应该是一个逆过程

    在刚刚的线性变换中,我们把我们的坐标系变为了动动的坐标系,应用了变换矩阵(egin{bmatrix}2 -1\1 1end{bmatrix})

    那么让动动的坐标系变为我们的坐标系,应该应用这个矩阵的逆矩阵(egin{bmatrix}2 -1\1 1end{bmatrix}^{-1})

    它代表我们所描述的变换逆向进行

    可以算出(egin{bmatrix}2 -1\1 1end{bmatrix}^{-1}=egin{bmatrix}frac{1}{3} frac{1}{3}\frac{-1}{3} frac{2}{3}end{bmatrix})

    动动想要知道我们说的坐标在它的语言中如何表示,应该(egin{bmatrix}frac{1}{3} frac{1}{3}\frac{-1}{3} frac{2}{3}end{bmatrix}egin{bmatrix}3\2end{bmatrix} =egin{bmatrix}frac{5}{3}\frac{1}{3}end{bmatrix})

    其实如果换一个角度:动动是正确的坐标系,我们才是变换过的坐标系,那么这个过程很好理解

    不同基线性变换

    但是实际上动动和我们之间的语言差异并不只在某一向量上

    比如我们描述某一线性变换:逆时针旋转(90^{。}),我们用(egin{bmatrix}0 1\-1 0end{bmatrix})

    但是(egin{bmatrix}0\1end{bmatrix})(egin{bmatrix}-1\0end{bmatrix})是我们语言中变换后的基向量,并不是动动的

    如果我们想对动动表达逆时针旋转(90^{。}),我们该怎么办?

    我们可以先把动动的基向量转化为我们的语言,然后应用线性变换再转化回去!

    比如我们想知道动动语言中的某个向量(v)旋转后在动动语言中如何描述

    (A)为基变换矩阵,(M)为旋转矩阵,则整个过程可以表示为 (A^{-1}MAv),注意矩阵从右向左变换

    其中(A^{-1}MA)就是旋转在动动语言中的矩阵形式

    特征向量与特征值

    特征向量与特征值定义

    考虑某个线性变换(egin{bmatrix}3 0\1 2end{bmatrix})

    在这个线性变换下,大部分向量会离开它张成的空间

    但是有些向量很特殊,它会留在原本张成的空间中

    这个线性变换对它来说只不过是拉伸或者压缩,作用同一个标量

    (egin{bmatrix}1\-1end{bmatrix}),在经过(egin{bmatrix}3 0\1 2end{bmatrix})只是伸长为(2)

    同时说明对于原来(egin{bmatrix}1\-1end{bmatrix})张成空间内的任意向量,都只是伸长为(2)

    对于这类变换之后没有离开张成空间的向量,我们称为变换的特征向量,每一个特征向量都有一个所属的值,称为特征值,即衡量压缩或拉伸比例的值

    你们肯定能猜到特征值可以为负或者为零

    那么它有什么作用呢?

    事实上以前我们研究线性变换过于依赖坐标轴,但是线性变换本身应该和坐标轴是无关的,通过特征向量和特征值也许可以更深刻的理解线性变换

    特征值的计算思想

    对于一个线性变化(A),设(v)是它的某个特征向量,其特征值为(λ)

    那么(Av=λv)

    但是我们根本解不了这么个东西……因为等式左边是一个矩阵乘法而右边是向量数乘

    所以我们应该将(λ)替换为某个矩阵,其效果为将任意向量乘以(λ)

    由于矩阵的每一列代表了变换后的基向量,所以我们只要让主对角线全部为(λ),其它地方全部为(0)

    通常也写作(λE),其中(E)为恒等变换

    我们变形一下式子得到((A-λE)v=0)

    如果(v)本身是零向量那么这个式子恒成立,但是我们一般不需要解出一个零向量作为特征向量

    根据前面的知识,我们知道一个非零向量经过一个线性变换变换为了一个零向量,那么这个线性变换一定降维,或者说行列式的值为(0)

    (det(A-λE)=0)

    但是有些线性变换其实并不存在特征向量,比如旋转

    在这种线性变换中,我们会解出一个虚数作为它的特征值

    特征值优化矩阵乘法

    考虑一个特殊情况:所有的基向量都是特征向量:

    那么我们可以发现,这个线性变换对应的矩阵是一个对角矩阵,即只有主对角线上有值,其它地方都是(0)

    它有什么良好的性质呢?

    性质就是对于一个对角矩阵,在自乘的时候复杂度为(O(n))而非(n^3)

    这对我们有一个启发:如果某个线性变换,它的特征向量很多,多到能够与维数同样多

    那么我们可以变换坐标系,使这些特征向量作为我们的基向量,然后进行快速的矩阵乘法,再将坐标系变换会原本的状态

    我们就可以(O(n^3+nlogk))实现矩阵快速幂

    抽象向量空间

    回到原本的问题:向量是什么?

    回顾行列式和特征向量的定义:

    一个代表某个线性变换对空间的压缩或拉伸程度

    一个代表某个线性变换中没有离开张成空间的向量

    好像和一个箭头或者一列数字都没有直接的联系

    为了解释向量到底是什么,我们引入一个看似与向量毫无关联的东西:函数

    假设我们现在有了函数(f(x))(g(x)),考虑它们的性质

    我们可以把两个函数相加,得到新函数((f+g)(x)=f(x)+g(x)),和向量加法非常类似

    那么对于函数与实数相乘,得到新函数((λf)(x)=λf(x)),也有着和向量数乘极为类似的解释

    其实我们对向量能够进行的操作只有相加和数乘,那么同样具备这些操作的函数理论上来说就可以直接套用所有向量相关的推论

    在这里我们要重新描述一下线性的严格定义

    可加性:(f(x+y)=f(x)+f(y))

    成比例(f(λx)=λf(x))

    如果一个变换在应用之后仍然对于原(x,y)满足这些性质,那么称这个变换是线性变换

    根据这个定义我们可以发现一个也许并没有听说过的定理:求导是线性运算

    为了类比求导与线性变换之间的关系,我们尝试用矩阵描述求导

    假设我们有一个多维空间,其空间组成为:全体多项式

    首先我们要赋予空间坐标的含义,即选取基向量,或者说基函数

    很自然我们应该选取(x)的不用次幂作为基函数,因为多项式次数可以任意高,所以我们的基函数几何应该是无穷大的

    对于某个多项式:(2x^2+3x+1),我们可以表示为(egin{bmatrix}1\3\2\0\vdotsend{bmatrix})

    在这个向量之中,求导可以用一个无限阶的矩阵来描述

    (frac{d}{dx}=egin{bmatrix}0 1 0 0 cdots\0 0 2 0 cdots\0 0 0 3 cdots\vdots vdots vdots ddots\end{bmatrix})

    然后我们对某一个多项式求导(frac{d}{dx}(5x^2+4x+5)=3x^2+10x+4),用矩阵描述为

    (egin{bmatrix}0 1 0 0 cdots\0 0 2 0 cdots\0 0 0 3 cdots\vdots vdots vdots ddots\end{bmatrix}egin{bmatrix}5\4\5\vdotsend{bmatrix} =egin{bmatrix}4\10\3\vdotsend{bmatrix})

    对求导矩阵的构造,只要将基函数每一项进行求导,然后按列排布

    为什么我们要类比求导和矩阵的关系?

    为了说明:只要我们处理的对象的集合满足相加和数乘原则,都可以称之为向量,对应的,向量的所有推论在这些对象中都适用

    我们可以大概概括为:向量是满足了相加和数乘原则的事物的集合

    结尾

    终于肝完了线代入门

    由于技术原因本文难以加入大量形象的图片,还有不了解的可以观看原视频

    再次大力膜拜(3B1B)

  • 相关阅读:
    LINQ 详解
    oracle下查询的sql已经超出IIS响应时间
    IOC应用之 Ninject
    JSONP ---------跨域
    国内各大互联网公司相关技术站点2.0版 (集合腾讯、阿里、百度、搜狐、新浪、360等共49个)
    IO多路复用,以socket为例
    socket机制下实现的多用户与服务器交互
    在一个进程中定义多个线程
    基于tcp的socketserver,即tcp的多线程
    基于upd的socketserver,即udp的多线程
  • 原文地址:https://www.cnblogs.com/knife-rose/p/12107204.html
Copyright © 2011-2022 走看看