python求矩阵特征值和特征向量_特征值,特征向量,标准正交向量组与numpy
使用python的数值计算库numpy来计算矩阵的特征值,特征向量与标准正交向量组
import numpy as np
1.求矩阵
的特征值和各特征值所对应的特征向量
x = np.array([[-1,0,1],[1,2,0],[-4,0,3]])
a,b=np.linalg.eig(x) ##特征值赋值给a,对应特征向量赋值给b
for i in range(len(a)):
print('特征值',a[i],'对应特征向量为',b[:,i])
特征值 2.0 对应特征向量为 [0. 1. 0.]
特征值 1.0 对应特征向量为 [ 0.40824829 -0.40824829 0.81649658]
特征值 1.0 对应特征向量为 [-0.40824829 0.40824829 -0.81649658]
2.求矩阵
的特征值和各特征值所对应的特征向量
x = np.array([[1,-2,2],[-2,-2,4],[2,4,-2]])
a,b=np.linalg.eig(x) ##特征值赋值给a,对应特征向量赋值给b
for i in range(len(a)):
print('特征值',a[i],'对应特征向量为',b[:,i])
特征值 2.000000000000001 对应特征向量为 [ 0.94280904 -0.23570226 0.23570226]
特征值 -6.999999999999997 对应特征向量为 [-0.33333333 -0.66666667 0.66666667]
特征值 1.9999999999999993 对应特征向量为 [-0.0232036 0.7126935 0.7010917]
3.由向量组
构造一组标准正交向量组
print('循环')
a = np.array([[0,1,0],[0,-1,1],[1,-1,2]])
b = np.zeros(a.shape)
#正交化
for i in range(len(a)):
b[i] = a[i]
for j in range(0,i):
b[i] -= np.dot(a[i],b[j])/np.dot(b[j],b[j])*b[j]
#归一化
for i in range(len(b)):
b[i] = b[i]/np.sqrt(np.dot(b[i],b[i]))
print(b)
循环
[[0. 1. 0.]
[0. 0. 1.]
[1. 0. 0.]]
python求矩阵特征值和特征向量_特征值,特征向量,标准正交向量组与numpy相关推荐
- python求矩阵的谱半径
在学习计算方法的时候,线性方程组的迭代法中的雅可比(Jacobi)迭代法和高斯-塞德尔(Gauss-Seidel)迭代法的收敛条件中需要求矩阵的谱半径,而经过查阅资料,python numpy库中没有 ...
- Python求矩阵的内积、外积、克罗内克直积、Khatri-Rao积
文章目录 矩阵乘法 内积和外积 直积 Khatri-Rao积 矩阵乘法 线性代数研究的核心对象是矩阵,所谓矩阵就是由 m m m行 n n n列的数组成的一个举行的数阵,从编程的角度理解,就是二维数组 ...
- python求矩阵的逆
文章目录 准备工作 Ready to work 环境 Environment 模块导入 Module import 代数余子式 Algebraic cofactor 定义 Definition 代码实 ...
- python中矩阵的表示方法_关于Python表示矩阵的方法详解
这篇文章主要介绍了Python表示矩阵的方法,结合具体实例形式分析了Python表示矩阵的方法与相关操作注意事项,需要的朋友可以参考下 本文实例讲述了Python表示矩阵的方法.分享给大家供大家参考, ...
- 用Python求矩阵的广义逆
文章目录 inv 穆尔-彭罗斯广义逆 pinv 对于两个方阵 A , B A,B A,B,若 A B = E AB=E AB=E,且 E E E为单位阵,则 A , B A,B A,B互逆,可记作 A ...
- python求三个整数最大值_实例讲解Python中整数的最大值输出
实例讲解Python中整数的最大值输出 在Python中可以存储很大的值,如下面的Python示例程序: x = 1000000000000000000000000000000000000000000 ...
- python 矩阵特征值_用python求矩阵特征值的方法?
我目前正在尝试寻找矩阵H的特征值. 我试过用这两种方法纽比.利纳格.艾格以及西皮.利纳格.艾格,尽管两者显然使用相同的底层方法.在 问题是我的矩阵H是纯实的,特征值必须是实的,也是正的.在 但是sci ...
- python求矩阵的最大特征值及对应的特征向量
直接看下面的代码就行啦! import numpy as np#输入矩阵 A = np.array([[1, 1/2, 4, 3, 3],[2, 1, 7, 5, 5],[1/4, 1/7, 1, 1 ...
- python求三个整数最大值_怎么用python比较三个数大小
大部分初学编程的人来说刚开始都会练习判断两个数或者三个数的大小,来熟悉某种语言的特性和最基本的if,else循环,当我们学习了更高级的语法知识后,又会有不同的实现方式,比如依次接收用户输入的3个数,排 ...
最新文章
- R语言使用forestplot包绘制森林图:编码创建森林图仿真数据、汇总线修改、元素位置调整、垂直线、字体、风格、置信区间、线型、图例、刻度、标签等
- C# 数据库连接笔记
- js 获取元素,同级元素下的子元素总结
- android 编辑9图片,Android基础入门教程——1.6 .9(九妹)图片怎么玩
- linux智能电压表设计与实现,STC89C51数字电压表
- 一个数如果刚好与它所有的真因子之和相等,则称该数为一个“完数
- arcgis分隔图层重复出文件_【干货】ArcGIS不可或缺的制图技巧,处理好细节才能让图更专业!...
- autosize px转dp_干货集中营
- 计算机应用202001常规,2001年4月份全国高等教育自学考试计算机应用基础试题
- Java读取Rinex 2.11格式的观测值o文件
- NRF51822基于蓝牙协议栈SDK+FDS-Flash data storage 的使用及应用例程
- HTML生成一维码(带条码
- SQL中DDL、DML、DQL、DCL、TCL是什么意思
- iOS-微信朋友圈图片压缩算法
- 【OpenGL学习笔记⑦】——键盘控制镜头的平移【3D正方体 透视投影 观察矩阵 对LookAt的理解】
- 南阳理工ACM 题4《ASCII码排序》
- Association Class VS Full Class
- 批量爬起百度图片上的数据代码--python
- 关于 Glide 加载图片圆角问题
- 相机camera与小孔成像opengl