python求矩阵的谱半径

在学习计算方法的时候，线性方程组的迭代法中的雅可比(Jacobi)迭代法和高斯-塞德尔(Gauss-Seidel)迭代法的收敛条件中需要求矩阵的谱半径，而经过查阅资料，python numpy库中没有直接求谱半径的函数。
谱半径的定义为：
设A是n × n矩阵，λ1, …, λn是其特征值，则A的谱半径ρ(A)定义如下：

即矩阵A的谱半径等于矩阵A的特征值绝对值的最大值。

由于python numpy库中有直接求矩阵特征值的函数，所以我们只要对所得特征值集合求绝对值的最大值即可

import numpy as py
def spectral_radius(M):a,b=np.linalg.eig(M) #a为特征值集合，b为特征值向量return np.max(np.abs(a)) #返回谱半径

python求矩阵的谱半径相关推荐

matlab求矩阵的谱半径,【什么是矩阵的谱半径？怎么求？】百度-谱半径
概述:本道作业题是应镜拇同学的课后练习,分享的知识点是谱半径,指导老师为鲁老师,涉及到的知识点涵盖:[什么是矩阵的谱半径?怎么求?]百度-谱半径,下面是应镜拇作业题的详细. 题目:[什么是矩阵的谱半径 ...
python求矩阵特征值和特征向量_特征值，特征向量，标准正交向量组与numpy
使用python的数值计算库numpy来计算矩阵的特征值,特征向量与标准正交向量组 import numpy as np 1.求矩阵的特征值和各特征值所对应的特征向量 x = np.array([[ ...
Python求矩阵的内积、外积、克罗内克直积、Khatri-Rao积
文章目录矩阵乘法内积和外积直积 Khatri-Rao积矩阵乘法线性代数研究的核心对象是矩阵,所谓矩阵就是由 m m m行 n n n列的数组成的一个举行的数阵,从编程的角度理解,就是二维数组 ...
用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求矩阵的逆
文章目录准备工作 Ready to work 环境 Environment 模块导入 Module import 代数余子式 Algebraic cofactor 定义 Definition 代码实 ...
python 矩阵特征值_用python求矩阵特征值的方法？
我目前正在尝试寻找矩阵H的特征值. 我试过用这两种方法纽比.利纳格.艾格以及西皮.利纳格.艾格,尽管两者显然使用相同的底层方法.在问题是我的矩阵H是纯实的,特征值必须是实的,也是正的.在但是sci ...
python 求矩阵的转置矩阵
1:导入包numpy from numpy import * 2: 创建矩阵 A = mat([[1,2,0],[3,5,9]]) 3: 求逆矩阵 B = A.T
python求矩阵维度必须一致_python数据分析(二)--Numpy
数组运算数组与标量运算:作用于数组的每一个元素a = a / a.mean() 一元函数:np.abs(x).np.fabs(x).np.sqrt(x).np.square(x).np.log(x) ...
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求矩阵的谱半径相关推荐

最新文章

热门文章