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

yipeiwu_com7年前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之list对应元素求和的方法

本次分享将讲述如何在Python中对多个list的对应元素求和,前提是每个list的长度一样。比如:a=[1,2,3], b=[2,3,4], c=[3,4,5], 对a,b,c的对应元...

基于DATAFRAME中元素的读取与修改方法

DATAFRAME中使用iat[1,0]和iloc[0,1]对元素进行修改。 a = [("hahaha",1),("lalala",2),("cacaca",6)] b = pad...

Django获取该数据的上一条和下一条方法

使用Django的ORM操作的时候,想要获取本条,上一条,下一条。 初步的想法是写3个ORM,3个ORM如下: 本条:models.Obj.objects.filter(id=n).fi...

Python字典遍历操作实例小结

本文实例讲述了Python字典遍历操作。分享给大家供大家参考,具体如下: 1 遍历键值对 可以使用一个 for 循环以及方法 items() 来遍历这个字典的键值对。 dict =...

Python骚操作之动态定义函数

在 Python 中,没有可以在运行时简化函数定义的语法糖。然而,这并不意味着它就不可能,或者是难以实现。 from types import FunctionType foo_c...