首页 > 其他 > 详细

Numpy库进阶教程(二)

时间:2016-04-04 19:46:34      阅读:243      评论:0      收藏:0      [点我收藏+]

第一篇在这里:Numpy库进阶教程(一)求解线性方程组

求解特征值和特征向量

关于特征值和特征向量的介绍,可以点击这里
首先创建一个矩阵

In [1]: A=mat("3 -2;1 0")

In [2]: A
Out[2]: 
matrix([[ 3, -2],
        [ 1,  0]])

在numpy.linalg模块中,eigvals函数可以计算矩阵的特征值,而eig函数可以返回一个包含特征值和对应特征向量的元组。

使用eigvals函数求解特征值

In [3]: linalg.eigvals(A)
Out[3]: array([ 2.,  1.])

使用eig函数求解特征值和特征向量。该函数将返回一个元祖,按列排放着特征值和对应的特征向量,其中第一列为特征值,第二列为特征向量。

In [5]: B=eig(A)

In [6]: B
Out[6]: 
(array([ 2.,  1.]), matrix([[ 0.89442719,  0.70710678],
         [ 0.4472136 ,  0.70710678]]))

使用dot函数使矩阵相乘验证求得的解是否正确

eigenvalues, eigenvectors = np.linalg.eig(A)
print "First tuple of eig", eigenvalues
print "Second tuple of eig\n", eigenvectors

for i in range(len(eigenvalues)):
      print "Left", np.dot(A, eigenvectors[:,i])
      print "Right", eigenvalues[i] * eigenvectors[:,i]
      print

输出为

Left [[ 1.78885438]
 [ 0.89442719]]
Right [[ 1.78885438]
 [ 0.89442719]]

Left [[ 0.70710678]
 [ 0.70710678]]
Right [[ 0.70710678]
 [ 0.70710678]]

计算矩阵行列式

In [51]: A = np.mat("3 4;5 6")

In [52]: print "A\n", A
A
[[3 4]
 [5 6]]


使用det函数计算行列式

In [53]: print "Determinant", np.linalg.det(A)
Determinant -2.0

Numpy库进阶教程(二)

原文:http://blog.csdn.net/pwiling/article/details/51058028

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!