python最简单怎么弄矩阵?Python求解特征向量和拉普拉斯矩阵
学过线性代数和深度学习先关的一定知道特征向量和拉普拉斯矩阵,这两者是很多模型的基础,有着很重要的地位,那用Python要怎么实现呢?
numpy和scipy两个库中模块中都提供了线性代数的库linalg,scipy更全面些。
特征值和特征向量
import scipy as sc#返回特征值,按照升序排列,num定义返回的个数def eignvalues(matrix, num): return sc.linalg.eigh(matrix, eigvalues(0, num-1))[0]#返回特征向量def eighvectors(matrix): return sc.linalg.eigh(matrix, eigvalues(0, num-1))[1]
调用实例
#创建一个对角矩阵,很容易得知它的特征值是1,2,3matrix = sc.diag([1,2,3])#调用特征值函数,获取最小的特征值minValue = eighvalues(matrix, 1)#调用特征向量函数,获取所有的特征向量vectors = eighvectors(matrix, 3)
拉普拉斯矩阵
很多图模型中都涉及到拉普拉斯矩阵,它有三种形式,这次给出的代码是D-A(度矩阵-邻接矩阵)和第二种标准化的形式:
#laplacian矩阵import numpy as npdef unnormalized_laplacian(adj_matrix): # 先求度矩阵 R = np.sum(adj_matrix, axis=1) degreeMatrix = np.diag(R) return degreeMatrix - adj_matrix def normalized_laplacian(adj_matrix): R = np.sum(adj_matrix, axis=1) R_sqrt = 1/np.sqrt(R) D_sqrt = np.diag(R_sqrt) I = np.eye(adj_matrix.shape[0]) return I - D_sqrt * adj_matrix * D_sqrt
- 03-04户外找蜂视频:寻蜂路上二,真实户外
- 11-17产品设计师要考什么证?产品设计师证怎么考办理什么条件多少钱
- 05-12lng采访doinb,LNG直播惨状,Doinb被喷拉夫你真菜
- 11-0512岁少年离家出走后续:我儿子,15岁留下一封信带16500元离家出走
- 04-29七年级生物期末复习重点归纳,七年级生物期末复习死记硬背勉强及格
- 05-07王俊凯跳舞慢节奏 王俊凯跳舞跳到低血糖
- 05-07早知天命:早知天下事
- 04-07科室特色和优势:完美结合科室特色
- 10-23赛里木湖什么时候最好看?赛里木湖,景超美
- 12-02韩国innisfree悦诗风吟亲肤透润 本站首晒,innisfree
- 03-19关于oppo手机使用小窍门:使用OPPO手机一定要知道这4大功能
- 04-06论语告诉我们的人生道理:每天一句论语,第170天子曰
- 04-04新市区儿童玩具,永嘉教玩具,玩出大市场
- 01-27钱江巡航摩托车哪款好:摩博会钱江新发五款车型
- 05-21佳能相机的四个参数介绍:详解佳能,CanonPowerShot
- 01-29吃什么美白皮肤最快最有效?美白肌肤吃什么效果最快
热门
推荐
- 1党代表候选人评价优秀范文482
- 2古文里面惊艳的句子_古风语录里惊艳人的句子172
- 3关于七夕的300字日记两篇224
- 4写春的现代诗430
- 5小学优秀少先队员的事迹材料322
- 6励志语录的话语146
- 7钢材供货合同270
- 8荔枝酒的酿造方法是什么432