首页 > 代码库 > 线性代数和numpy

线性代数和numpy

一、代数是什么

  代数->数的抽象表示->向量空间(线性空间)

  线代->线性代数

 

关系:

  向量空间之间和内部转换是通过线性变换。

实数——一维空间的点

复数——二维空间的点

  如果两个向量的组合可以生成平面,则要求两个向量要线性无关。

  推广一下,N维空间里点可以用N个线性无关的向量来表示。这N个向量就是这个平面的基。

  向量的封闭——对加法和数乘封闭。

    向量V中任意两个向量a,b加法a+b,仍然在V中,实数乘法x*b,仍然也在V中。

 

线性相关——其中的一个向量可以用其他的向量表示出来。

 

矩阵操作在python里编程依赖一个最常用的库——numpy

 

1、矩阵的创建

a=np.arange(1,5)a=np.array([1,2,3,4,5])print a, a.dtype, a.shape, a.size, a.ndim

技术分享

np.arange类似range函数np.array用来生成矩阵dtype是数据类型,有int64, complex, uint16等shape是个元组属性,表示每一维的宽度size是所有元素个数ndim是维数
b=np.array([1,2,3],dtype=float16) # int64, complex, uint16......print b, b.dtype

技术分享

 

m=np.array([np.arange(6),np.arange(6)])print m, m.shape, m.size

技术分享

# 每一个[]代表一维,比如# [[[1,2,3],[4,5,6]],[[7,8,9],[10,11,12]]], 代表矩阵的维度是(2,2,3)# 其中第一个2,代表最外层的两个[],第二个2代表第二层[],第三个3代表最里层的维度。n=np.array([[1,2,3,4],[5,6,7,8]])print n, n[0,2], n[1,1]m=np.array([[[1,2,3],[4,5,6]],[[7,8,9],[10,11,12]]])print m.shape

 

 技术分享
x=m.ravel()y=n.flatten()print xprint y

技术分享

 

 

 

ravel()和flatten()看起来效果一样,都是把矩阵展平了。它们的区别在于ravel()返回的是原有数据的一个映射(view),没有分配新的存储flatten()返回的是新的数据因此如果我们改变它们的值,就可以看出区别numpy还有一些函数有这样的区别,关键在于判断函数返回是原数据的映射还是返回新的数据。

…… 

等等一共三十道题,详情可参阅黑板课老师的notebook。

 

 

正交:A矩阵的转置乘以B等于0。

可逆:方阵,行数等于列数

   列向量线性无关

 

 

最小二乘法——投影解释

    二维平面上,一个点向一条直线投影。

    扩展到N维里边,就是让一个直线拟合若干个点,让这些点到线的距离和最短。y=ax+b

    但有的时候,欠拟合,可以往更高维的空间投影,y=ax^3+bx^2+cx+d,或者更高维的空间。

    但维数过高容易过拟合。

 

行列式

  从平行四边形的面积推出来。

 

线性代数和numpy