zoukankan      html  css  js  c++  java
  • 线性代数和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 x
    print y

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

    …… 

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

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

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

       列向量线性无关

    最小二乘法——投影解释

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

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

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

        但维数过高容易过拟合。

    行列式

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

    特征值和特征向量

      AX=λX——向量X经过向量A的线性变换,得到的是向量A的倍数。λ称为A的特征值,X称为A的特征向量。

      什么意思呢?X经过A的变换,得到的是X的方向的伸缩。

      在matlab里边有例子,eigshow,然后运行,推重两个向量重合时,重合的向量就是特征向量。两个向量成一条直线的时候,也会有特征向量。然后反方向,所以有四组。从图形中看,是一个椭圆和圆的交合。

      如何求特征值和特征向量?

    如何用特征值和特征向量来求Fibonacci序列

  • 相关阅读:
    zabbix添加自定义监控项
    文件的时间属性(摘自网友的博客)
    svn从Windows服务器上迁移到Linux上
    modelsim 保存波形文件
    Viso Professional 2013版本激活(office 系列产品 -- visio 2013 / project 2013 破解工具
    kindle paperwhite3 连不上WIFI解决方法
    TotoiseSVN的基本使用方法
    SVN 版本服务器搭配全过程详解(含服务端、客户端)
    modelsim读写TXT文件
    如何对多个文件进行MODELSIM仿真? (由于是一个很大的项目,不可能把所有MODULE都放在一个文件里。 如何在ModelSim中对多个.V文件进行仿真?)
  • 原文地址:https://www.cnblogs.com/shixisheng/p/5928215.html
Copyright © 2011-2022 走看看