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中is和id的用法

(ob1 is ob2) 等价于 (id(ob1) == id(ob2)) 首先id函数可以获得对象的内存地址,如果两个对象的内存地址是一样的,那么这两个对象肯定是一个对象。和is是等价...

pandas string转dataframe的方法

今天业务上碰到用pandas处理一个大文件的内存不够问题,需要做concat 合并多个文件,每个文件数据在1.4亿行左右。当时第一反应是把dataframe分割成多块小文件处理,后面发现...

Python使用装饰器进行django开发实例代码

本文研究的主要是Python使用装饰器进行django开发的相关内容,具体如下。 装饰器可以给一个函数,方法或类进行加工,添加额外的功能。 在这篇中使用装饰器给页面添加session而不...

Python yield使用方法示例

1. iterator叠代器最简单例子应该是数组下标了,且看下面的c++代码: 复制代码 代码如下:int array[10];for ( int i = 0; i < 10; i...

浅谈pycharm下找不到sqlalchemy的问题

浅谈pycharm下找不到sqlalchemy的问题

平台: windows 10 pycharm 2016.2 python 2.7.12 问题始于我在pycharm下建了一个flask工程,然后导入sqlalchemy的包: fro...