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绘制矩阵盖尔圆简单实例的全部内容,希望对大家有所帮助。感兴趣的朋友可以继续参阅本站其他相关专题,如有不足之处,欢迎留言指出。感谢朋友们对本站的支持!

相关文章

对numpy中的数组条件筛选功能详解

在程序设计中,时常会遇到数据的唯一化、相同、相异信息的提取等工作,在格式化的向量存储矩阵中南,numpy能够提供比较不错的快速处理功能。 1,唯一化的实现: In [63]: dat...

跟老齐学Python之print详解

eval() 在print干事情之前,先看看这个东东。不是没有用,因为说不定某些时候要用到。 复制代码 代码如下: >>> help(eval)  ...

解决pycharm安装后代码区不能编辑的问题

解决pycharm安装后代码区不能编辑的问题

此问题是由于最新的pycharm在安装时自动装了vimVim插件 你可以在tools Vim emulator将对勾去掉就可以了。 以上这篇解决pycharm安装后代码区不能编辑的问题...

Python通过调用mysql存储过程实现更新数据功能示例

本文实例讲述了Python通过调用mysql存储过程实现更新数据功能。分享给大家供大家参考,具体如下: 一、需求分析 由于管理费率配置错误,生成订单的还本付息表和订单表的各种金额,管理费...

python3 tcp的粘包现象和解决办法解析

这篇文章主要介绍了python3 tcp的粘包现象和解决办法解析,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下 服务器端 impo...