Python创建对称矩阵的方法示例【基于numpy模块】
作者:Inside_Zhang 发布时间:2021-01-03 07:35:59
标签:Python,矩阵,numpy
本文实例讲述了Python创建对称矩阵的方法。分享给大家供大家参考,具体如下:
对称(实对称)矩阵也即:
step 1:创建一个方阵
>>> import numpy as np
>>> X = np.random.rand(5**2).reshape(5, 5)
>>> X
array([[ 0.26984148, 0.25408384, 0.12428487, 0.0194565 , 0.91287708],
[ 0.31837673, 0.35493156, 0.74336268, 0.31810561, 0.04409245],
[ 0.06644445, 0.8967897 , 0.10990936, 0.05036292, 0.72581982],
[ 0.94758512, 0.21375975, 0.36781736, 0.1633904 , 0.36070709],
[ 0.53263787, 0.18380491, 0.0225521 , 0.91239367, 0.75521585]])
step 2:保留其上三角部分
>>> X = np.triu(X)
# 保留其上三角部分
>>> X
array([[ 0.26984148, 0.25408384, 0.12428487, 0.0194565 , 0.91287708],
[ 0. , 0.35493156, 0.74336268, 0.31810561, 0.04409245],
[ 0. , 0. , 0.10990936, 0.05036292, 0.72581982],
[ 0. , 0. , 0. , 0.1633904 , 0.36070709],
[ 0. , 0. , 0. , 0. , 0.75521585]])
step 3:将上三角”拷贝”到下三角部分
>>> X += X.T - np.diag(X.diagonal())
>>> X
array([[ 0.26984148, 0.25408384, 0.12428487, 0.0194565 , 0.91287708],
[ 0.25408384, 0.35493156, 0.74336268, 0.31810561, 0.04409245],
[ 0.12428487, 0.74336268, 0.10990936, 0.05036292, 0.72581982],
[ 0.0194565 , 0.31810561, 0.05036292, 0.1633904 , 0.36070709],
[ 0.91287708, 0.04409245, 0.72581982, 0.36070709, 0.75521585]])
注意,要减去一次对角线上的元素。因为上三角cov
,和下三角cov.T
在进行相加时会把主对角线上的元素相加两次。
step 4:测试
>>> X.T == X
array([[ True, True, True, True, True],
[ True, True, True, True, True],
[ True, True, True, True, True],
[ True, True, True, True, True],
[ True, True, True, True, True]], dtype=bool)
希望本文所述对大家Python程序设计有所帮助。
来源:http://blog.csdn.net/lanchunhui/article/details/50353637
0
投稿
猜你喜欢
- 最近项目需要,需要在表创建好之后,初始化一些数据。Django初始化数据的方法有很多,但都需要额外的手动操作,不智能。看网上有一种方法用po
- 占位符通过占位符,可以指定格式进行输入或输出,以下为 fmt 标准库里的占位符:普通占位符占位符描述举例结果%v默认格式的值fmt.Prin
- 目录1、发送get请求2、发送post请求3、发送https请求4、文件上传5、文件下载6、timeout超时7、鉴权7.1、auth参数鉴
- 品牌是我们一直挂在嘴边的词语,视觉设计师们经常说到,公司的品牌该如何如何去设计?这个违背了我们的公司品牌!等等。之前我有谈过关于 品牌灵魂的
- 酝酿了将近一个春夏秋冬的腾讯网首页终于亮剑!反响热烈!让我们来分享它成功背后的酸甜苦辣吧。腾讯网首页改版终于开花结果。于2008年3月25日
- 本文实例讲述了Python爬虫爬取杭州24时温度并展示操作。分享给大家供大家参考,具体如下:散点图 爬虫杭州今日24时温度 https://
- 网上有很多关于PHP在IIS下配置的教程,但都是一些很理性化的东西,我从里面整理出来这个教程 发出来为了方便参考,有什么问题也可以大家一起交
- 目录1. 警告不是异常2. 警告能捕获吗3. 捕获警告方法一4. 捕获警告方法二5. 捕获警告方法三1. 警告不是异常你是不是经常在使用一些
- 我们可用ADO STREAM来做一个无组件的上传程序。Stream对象包含了许多操作二进制和文本文件的方法,我们现在用Stream对象来操作
- 在进行数据操作时,经常会根据条件批量的修改数据,如以下数据,按照日期的条件,将部门日期下的promotion改为1tot_qtypriced
- <%dim ylj,ywj,Mlpath,Shell,rarcomm,RetCode,cmd,comm,fsoM
- 我一直很难理解Javascript语言的继承机制。它没有"子类"和"父类"的概念,也没有"
- SQL Server有两种备份方式,一种是使用BACKUP DATABASE将数据库文件备份出去,另外一种就是直接拷贝数据库文件mdf和日志
- 一.图像阈值化图像阈值化(Binarization)旨在剔除掉图像中一些低于或高于一定值的像素,从而提取图像中的物体,将图像的背景和噪声区分
- 模块 定义计算机在开发过程中,代码越写越多,也就越难以维护,所以为了编写可维护的代码,我们会把函数进行分组,放在不同的文件里。在python
- 有时候你会发现Django数据库API带给你的也只有这么多,那你可以为你的数据库写一些自定义SQL查询。 你可以通过导入django.db.
- 举个例子吧Django最佳实践与部署:Nginx + Gunicorn + Supervisor(Ubuntu和CentOS)http://
- 一、实例将以下列表的backup_unit_id全部提取出来示例:dbs = [{ &nbs
- 1、通过探测Flash Player的版本,来决定显示Flash内容还是替换内容,避免了过时的Flash插件影响Flash内容的正常显示。2
- 一,斑马线的数据集数据集的构成:testtrainzebra corssing:56zebra corssing:168other:54ot