numpy.linalg.eig
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的更多相关文章
- 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 ...
- numpy linalg模块
# 线性代数# numpy.linalg模块包含线性代数的函数.使用这个模块,可以计算逆矩阵.求特征值.解线性方程组以及求解行列式等. import numpy as np # 1. 计算逆矩阵# 创 ...
- python中numpy计算数组的行列式numpy.linalg.det()
numpy.linalg.det numpy.linalg.det(a)[source] 计算任何一个数组a的行列式,但是这里要求数组的最后两个维度必须是方阵. 参数: a : (..., M, M) ...
- numpy.linalg.norm(求范数)
1.linalg=linear(线性)+algebra(代数),norm则表示范数. 2.函数参数 x_norm=np.linalg.norm(x, ord=None, axis=None, keep ...
- numpy.linalg.svd函数
转载自:python之SVD函数介绍 函数:np.linalg.svd(a,full_matrices=1,compute_uv=1) 参数: a是一个形如\((M,N)\)的矩阵 full_matr ...
- numpy linalg
线性代数 np.mat("0 1 0;1 0 0;0 0 1") np.linalg.inv(A)
- Python数据分析--Numpy常用函数介绍(9)-- 与线性代数有关的模块linalg
numpy.linalg 模块包含线性代数的函数.使用这个模块,可以计算逆矩阵.求特征值.解线性方程组以及求解行列式等.一.计算逆矩阵 线性代数中,矩阵A与其逆矩阵A ^(-1)相乘后会得到一个单位矩 ...
- [Python学习] python 科学计算库NumPy—矩阵运算
NumPy库的核心是矩阵及其运算. 使用array()函数可以将python的array_like数据转变成数组形式,使用matrix()函数转变成矩阵形式. 基于习惯,在实际使用中较常用array而 ...
- Numpy应用100问
对于从事机器学习的人,python+numpy+scipy+matplotlib是重要的基础:它们基本与matlab相同,而其中最重要的当属numpy:因此,这里列出100个关于numpy函数的问题, ...
随机推荐
- BZOJ 4034 BIT & Dfs序
调了恒久突然发现输出优化忘记带负号了.. 就是差分树状数组维护Dfs序即可. #include <iostream> #include <cstring> #include & ...
- 进程间通信IPC:消息队列,信号量,共享内存
2015.3.4星期三 阴天 进程间通信:IPC 文件对象:记录文件描述符,文件开关等 IPC标示符:系统全局的流水号两个进程要通信,打开的是唯一的对象进行通讯,通过key操作 XSI IPC:消息队 ...
- 制作.frameWork的最全最真实的解决办法
这个制作流程 本博主 已经完全测试成功 我这边 制作的.frameWork 要接入游戏 我们游戏已经上架 所以这个东西完全可以用 http://www.cocoachina.com/bbs/read. ...
- viewPager动态加载listview数据
废话不多说,先上效果图.(代码见附件) 代码是修改自某大神的,我做了很多修改,之前只能向右滑动,现在可以左右无限滑动,只要数据没加载完就可以一直滑动.过程不算复杂,代码主要的地方都有注释. 附件dem ...
- Node.js高级编程读书笔记 - 6 应用程序构建和调试 - Never
Explanation 现阶段console.log(...),util.inspect(...), JSON.stringify(...)在控制台输出已经够用了[2015/07/19]. 单元测试隶 ...
- ajax上传图片到服务器
通过Ajax方式上传文件,使用FormData进行Ajax请求.上传文件或图片时,要求的格式为enctype ="multipart/form-data"(以二进制格式上传),在使 ...
- 职责链模式(chain of responsibility Pattern)
职责链模式:使多个对象都有机会处理请求,从而避免请求的发送者和接受者之间的耦合关系.将这个对象连成一条链,并沿着这条链传递该请求,直到有一个对象处理他为止. •Handler: 抽象处理者:定义出一个 ...
- Mac下Nginx环境配置
环境信息: Mac OS X 10.11.1 Homebrew 0.9.5 正文 一.安装 Nginx 终端执行: brew search nginx brew install nginx 当前版本 ...
- C++访问托管类(C#类库)
1.新建C# 类库项目Airth,新建类 using System; using System.Collections.Generic; using System.Linq; using System ...
- Unity3D特效-场景淡入淡出
最近公司开始搞Unity3D..整个游戏..特效需求还是比较多的.关于UI部分的特效淡入淡出.看网上用的方法都是用个黑东东遮挡然后设置alpha这么搞....本大神感觉非常的low.而且很渣.故奋笔疾 ...