Gershgorin圆盘定理
创始人
2024-05-14 03:16:12
0

计算方法

  圆盘定理Gershgorin Circle Theorem是用来预估矩阵特征值的范围的。这个定理用起来相当简单。对于一个n×nn\times nn×n的方阵来说,有nnn个行圆盘和列圆盘,它的特征值就位于这些圆盘中。所谓的圆盘Gershgorin disk,就是复平面上的一个圆,它有中心和半径组成,是这样定义的:
D(a,r)={x∈C:∣x−a∣≤r}D(a,r)=\{x \in \mathbb{C}:|x-a|\le r\} D(a,r)={x∈C:∣x−a∣≤r}
  从定义看表示一个点的集合,那么中心和半径怎么确定呢?中心就是矩阵的对角线元素,半径则分两种,行圆盘和列圆盘。行圆盘的半径就是某一行除了对角线元素外,所有元素的模长的和,列圆盘则是某一列除了对角线元素外,所有元素的模长的和。既然每个中心有两个半径,那肯定取最小的那个半径喽,把半径用rrr表示,那么半径就是:
ri=min(∑j=1,j≠in∣aij∣,∑j=1,j≠in∣aji∣)r_i=min(\sum_{j=1,j\ne i}^{n}|a_{ij}|,\sum_{j=1,j\ne i}^{n}|a_{ji}|) ri​=min(j=1,j=i∑n​∣aij​∣,j=1,j=i∑n​∣aji​∣)

举例

  比如计算这个矩阵的圆盘:
A=(30.40.20.20.230.10.10.40.540.5−0.310.51)A=\begin{pmatrix}3 & 0.4 & 0.2 & 0.2\\ 0.2 & 3 & 0.1 & 0.1\\ 0.4 & 0.5 & 4 & 0.5\\ -0.3 & 1 & 0.5 & 1\\ \end{pmatrix} A=​30.20.4−0.3​0.430.51​0.20.140.5​0.20.10.51​
  得到的结果就是:
D(3,0.8),D(3,0.4),D(4,0.8),D(1,0.8)D(3,0.8),D(3,0.4),D(4,0.8),D(1,0.8) D(3,0.8),D(3,0.4),D(4,0.8),D(1,0.8)
  因为D(3,0.4)D(3,0.4)D(3,0.4)比D(3,0.8)D(3,0.8)D(3,0.8)范围更小,所以最终结果是:
D(3,0.4),D(4,0.8),D(1,0.8)D(3,0.4),D(4,0.8),D(1,0.8) D(3,0.4),D(4,0.8),D(1,0.8)
  所以矩阵A的特征值就位于复平面上这些圆盘中。

python实现

  这种实现代码非常简单:

    def gershgorin_disk(self):disks = [[0, 0] for _ in self.__lines]n = len(self.__lines)for i in range(n):d = disks[i]d[0]= self.__lines[i][i]line_radius = 0for j in range(n):if j != i:line_radius += abs(self.__lines[j][i])column_radius = 0for j in range(n):if j != i:column_radius += abs(self.__lines[i][j])d[1] = min(line_radius, column_radius)return disks

相关内容

热门资讯

银河麒麟V10SP1高级服务器... 银河麒麟高级服务器操作系统简介: 银河麒麟高级服务器操作系统V10是针对企业级关键业务...
【NI Multisim 14...   目录 序言 一、工具栏 🍊1.“标准”工具栏 🍊 2.视图工具...
不能访问光猫的的管理页面 光猫是现代家庭宽带网络的重要组成部分,它可以提供高速稳定的网络连接。但是,有时候我们会遇到不能访问光...
Android|无法访问或保存... 这个问题可能是由于权限设置不正确导致的。您需要在应用程序清单文件中添加以下代码来请求适当的权限:此外...
AWSECS:访问外部网络时出... 如果您在AWS ECS中部署了应用程序,并且该应用程序需要访问外部网络,但是无法正常访问,可能是因为...
北信源内网安全管理卸载 北信源内网安全管理是一款网络安全管理软件,主要用于保护内网安全。在日常使用过程中,卸载该软件是一种常...
AWSElasticBeans... 在Dockerfile中手动配置nginx反向代理。例如,在Dockerfile中添加以下代码:FR...
AsusVivobook无法开... 首先,我们可以尝试重置BIOS(Basic Input/Output System)来解决这个问题。...
ASM贪吃蛇游戏-解决错误的问... 要解决ASM贪吃蛇游戏中的错误问题,你可以按照以下步骤进行:首先,确定错误的具体表现和问题所在。在贪...
​ToDesk 远程工具安装及... 目录 前言 ToDesk 优势 ToDesk 下载安装 ToDesk 功能展示 文件传输 设备链接 ...