python+pillow绘制矩阵盖尔圆简单实例

yipeiwu_com6年前Python基础

本文主要研究的是使用Python+pillow绘制矩阵盖尔圆的一个实例,具体如下。

盖尔圆是矩阵特征值估计时常用的方法之一,其定义为:

与盖尔圆有关的两个定理为:

定理1:矩阵A的所有特征值均落在它的所有盖尔圆的并集之中。

定理2:将矩阵A的全体盖尔圆的并集按连通部分分成若干个子集,(一个子集由完全连通的盖尔圆组成,不同子集没有相连通的部分),对每个子集,若它恰好由K个盖尔圆组成,则该子集中恰好包含A的K个特征值。

与盖尔圆定理有关的几个推论为:

推论1:孤立盖尔圆中恰好包含一个特征值。

推论2:实矩阵的孤立盖尔圆恰好包含一个实特征值。

推论3:盖尔圆方法中盖尔圆半径可以按列求和。(因为方阵转置后特征值不变)

下面的代码使用Python+pillow绘制给定矩阵的盖尔圆:

当输入矩阵[[10,20,30],[30,40,50],[50,65,70]]时,得到的图形如下:

当输入矩阵[[6,27,33],[39,28,17],[80,60,71]]时,得到的图形如下:

总结

以上就是本文关于python+pillow绘制矩阵盖尔圆简单实例的全部内容,希望对大家有所帮助。感兴趣的朋友可以继续参阅本站其他相关专题,如有不足之处,欢迎留言指出。感谢朋友们对本站的支持!

相关文章

Python实现bilibili时间长度查询的示例代码

Python实现bilibili时间长度查询的示例代码

说明:最近在B站看一些材料力学视频时候,感觉有一些分集狂魔的分集真的很恐怖,有的甚至上百集,因此决定写个小脚本每次分析下到底这个系列视频到底有多长。 好了,下面是分析过程: 第一步当然...

利用python获取当前日期前后N天或N月日期的方法示例

前言 最近因为工作原因,发现一个Python的时间组件,很好用分享出来!(忘记作者名字了,在这里先感谢了),下面话不多说,来一起看看详细的介绍吧。 示例代码: # -*- codi...

详解Python 实现元胞自动机中的生命游戏(Game of life)

详解Python 实现元胞自动机中的生命游戏(Game of life)

简介 细胞自动机(又称元胞自动机),名字虽然很深奥,但是它的行为却是非常美妙的。所有这些怎样实现的呢?我们可以把计算机中的宇宙想象成是一堆方格子构成的封闭空间,尺寸为N的空间就有NN个格...

Python 的类、继承和多态详解

类的定义 假如要定义一个类 Point,表示二维的坐标点: # point.py class Point: def __init__(self, x=0, y=0): s...

Python语言的变量认识及操作方法

今天我给大家介绍的是python中的Number变量,与c++,java有些不同,下面让来为大家介绍: 在python中是不用声明变量类型的,不过在使用变量前需要对其赋值,没有值得变量是...