1、转置对于二维数组有用,对一位数组无效

2、理解特征值和特征向量的对应关系

a=np.array([[1 ,2, 3],[4, 5, 6],[7, 8, 9]])

a
Out[27]:
array([[1, 2, 3],
[4, 5, 6],
[7, 8, 9]]) w,v = LA.eig(a) w
Out[29]: array([ 1.61168440e+01, -1.11684397e+00, -1.30367773e-15]) v
Out[30]:
array([[-0.23197069, -0.78583024, 0.40824829],
[-0.52532209, -0.08675134, -0.81649658],
[-0.8186735 , 0.61232756, 0.40824829]]) a
Out[31]:
array([[1, 2, 3],
[4, 5, 6],
[7, 8, 9]]) #dot(a[:,:], v[:,i]) = w[i] * v[:,i] w[0]
Out[33]: 16.116843969807043 v[:,0]
Out[34]: array([-0.23197069, -0.52532209, -0.8186735 ]) w[0]*v[:,0]
Out[35]: array([ -3.73863537, -8.46653421, -13.19443305]) np.dot(a[:,:],v[:,0])
Out[37]: array([ -3.73863537, -8.46653421, -13.19443305]) a
Out[38]:
array([[1, 2, 3],
[4, 5, 6],
[7, 8, 9]]) v[:,0]
Out[39]: array([-0.23197069, -0.52532209, -0.8186735 ]) 1*-0.23197069+2*-0.52532209+3*-0.8186735
Out[40]: -3.73863537 4*-0.23197069+5*-0.52532209+6*-0.8186735
Out[41]: -8.46653421 v[:,0]
Out[42]: array([-0.23197069, -0.52532209, -0.8186735 ]) v[:,0].T
Out[43]: array([-0.23197069, -0.52532209, -0.8186735
]) w[0]
Out[44]: 16.116843969807043 w[0]*v[:,0]
Out[45]: array([ -3.73863537, -8.46653421, -13.19443305]) a
Out[
46]:
array([[1, 2, 3],
[4, 5, 6],
[7, 8, 9]]) a.T
Out[47]:
array([[1, 4, 7],
[2, 5, 8],
[3, 6, 9
]])

c = v[0] c
Out[49]: array([-0.23197069, -0.78583024, 0.40824829]) c.T
Out[50]: array([-0.23197069, -0.78583024, 0.40824829])
a
Out[55]:
array([[1, 2, 3],
[4, 5, 6],
[7, 8, 9]]) a*v[:,0]
Out[56]:
array([[-0.23197069, -1.05064419, -2.4560205 ],
[-0.92788275, -2.62661047, -4.912041 ],
[-1.62379481, -4.20257675, -7.36806149]]) v[:,0]
Out[57]: array([-0.23197069, -0.52532209, -0.8186735 ])
y = a*v[:,0]

y
Out[59]:
array([[-0.23197069, -1.05064419, -2.4560205 ],
[-0.92788275, -2.62661047, -4.912041 ],
[-1.62379481, -4.20257675, -7.36806149]]) np.sum(y[0,:])
Out[63]: -3.7386353719172973 np.sum(y[1,:])
Out[64]: -8.4665342116284013 np.sum(y[2,:])
Out[65]: -13.194433051339505

numpy.linalg.eig的更多相关文章

  1. numpy 辨异(三)—— hstack/column_stack,linalg.eig/linalg.eigh

    1. np.hstack np.column_stack >>> np.hstack([np.array([1, 2, 3]), np.array([4, 5, 6])]) arra ...

  2. numpy linalg模块

    # 线性代数# numpy.linalg模块包含线性代数的函数.使用这个模块,可以计算逆矩阵.求特征值.解线性方程组以及求解行列式等. import numpy as np # 1. 计算逆矩阵# 创 ...

  3. python中numpy计算数组的行列式numpy.linalg.det()

    numpy.linalg.det numpy.linalg.det(a)[source] 计算任何一个数组a的行列式,但是这里要求数组的最后两个维度必须是方阵. 参数: a : (..., M, M) ...

  4. numpy.linalg.norm(求范数)

    1.linalg=linear(线性)+algebra(代数),norm则表示范数. 2.函数参数 x_norm=np.linalg.norm(x, ord=None, axis=None, keep ...

  5. numpy.linalg.svd函数

    转载自:python之SVD函数介绍 函数:np.linalg.svd(a,full_matrices=1,compute_uv=1) 参数: a是一个形如\((M,N)\)的矩阵 full_matr ...

  6. numpy linalg

    线性代数 np.mat("0 1 0;1 0 0;0 0 1") np.linalg.inv(A)

  7. Python数据分析--Numpy常用函数介绍(9)-- 与线性代数有关的模块linalg

    numpy.linalg 模块包含线性代数的函数.使用这个模块,可以计算逆矩阵.求特征值.解线性方程组以及求解行列式等.一.计算逆矩阵 线性代数中,矩阵A与其逆矩阵A ^(-1)相乘后会得到一个单位矩 ...

  8. [Python学习] python 科学计算库NumPy—矩阵运算

    NumPy库的核心是矩阵及其运算. 使用array()函数可以将python的array_like数据转变成数组形式,使用matrix()函数转变成矩阵形式. 基于习惯,在实际使用中较常用array而 ...

  9. Numpy应用100问

    对于从事机器学习的人,python+numpy+scipy+matplotlib是重要的基础:它们基本与matlab相同,而其中最重要的当属numpy:因此,这里列出100个关于numpy函数的问题, ...

随机推荐

  1. HDU 4123 (2011 Asia FZU contest)(树形DP + 维护最长子序列)(bfs + 尺取法)

    题意:告诉一张带权图,不存在环,存下每个点能够到的最大的距离,就是一个长度为n的序列,然后求出最大值-最小值不大于Q的最长子序列的长度. 做法1:两步,第一步是根据图计算出这个序列,大姐头用了树形DP ...

  2. BackTrack5-r3硬盘安装

    用Universal-USB-Installer-1.9.4.2这个工具,选择系统类型为Backtrack5 R3,选择Backtrack5 R3系统的ISO文件,然后选中U盘盘符,点击”Create ...

  3. c/c++ 函数指针 指针函数 数组的引用 指针数组 数组指针

    1.指针数组数组指针 引用数组 数组的引用 int *a[10] 指针数组 每一个元素都是一个指针 Int (*a)[10] 数组指针 P指向一个含有10个元素的数组 Int (&a)[10] ...

  4. cocos2dx 3.0 之 lua 创建类

    利用3.0里面封装好的模板类来写,不罗嗦,直接上例子(消除星星小游戏中 星星类) require "extern" StarSprite = class("StarSpr ...

  5. Oracle Listener日志位置及压缩转移

    近日由于Oracle Listener异常断开导致应用无法上传数据,需要从listener日志开始分析问题原因.此文介绍如何获取日志位置.由于日志文件大小问题,同时将日志文件进行压缩存放. alert ...

  6. 移动互联网实战--资源类APP的数据存储处理和优化

    前言: 对于资源类的APP, 其音频/图形占据了APP本身很大的比例. 如何存储和管理这些资源文件, 成了一个颇具挑战性的难点. 移动端的碎片化, 高中低端手机的并存, 需要开发者不光是具备基础的存储 ...

  7. Java设计模式(十一) 享元模式

    原创文章,同步发自作者个人博客 http://www.jasongj.com/design_pattern/flyweight/.转载请注明出处 享元模式介绍 享元模式适用场景 面向对象技术可以很好的 ...

  8. java.lang.ClassCastException: java.lang.String cannot be cast to com.jy.hfims.domain 映射实体类型错误

    今天在做 excel导出的时候,出现了一个问题"java.lang.ClassCastException: java.lang.String cannot be cast to com.do ...

  9. 加载json文件没显示图片

    IOS 9新特性要求App内访问的网络必须使用HTTPS协议. 但是现在公司的项目使用的是HTTP协议,使用私有加密方式保证数据安全.现在也不能马上改成HTTPS协议传输. 以下解决办法: 在info ...

  10. Spring触发器配置Quartz

    参考文献:http://blog.csdn.net/liaq325/article/details/8269439 http://sundoctor.iteye.com/blog/441951 业务类 ...