Python科学计算之NumPy入门教程

前言

NumPy是Python用于处理大型矩阵的一个速度极快的数学库。它允许你在Python中做向量和矩阵的运算，而且很多底层的函数都是用C写的，你将获得在普通Python中无法达到的运行速度。这是由于矩阵中每个元素的数据类型都是一样的，这也就减少了运算过程中的类型检测。

矩阵基础

在 numpy 包中我们用数组来表示向量，矩阵和高阶数据结构。他们就由数组构成，一维就用一个数组表示，二维就是数组中包含数组表示。

创建

# coding: utf-8
import numpy as np

a = np.array([
 [1.73, 1.68, 1.71, 4],
 [1, 2, 3, 4],
 [1, 2, 3, 4]
])
print type(a) # <type 'numpy.ndarray'>

ndarray （N-dimensional array object） 意思就是n维数组。例子中就表示一个3行4列的二维数组。

形状

数组的大小可以通过其 shape 属性获得：

print a.shape # (3L,4L)

数组的元素数量可以通过 ndarray.size 得到：

print a.size # 12

使用 ndarray 的 dtype 属性我们能获得数组元素的类型：

print a.dtype # float64

可以用过 shape 重新设置矩阵的形状或者通过 reshape 方法创建一个改变了尺寸的新数组，原数组的shape保持不变:

a.shape = 4, 3
b = a.reshape((2, 6))
# 尽管b的形状是新的，但是a和b是共享数据存储内存区域的，如果b[0][1] = 8 那么a[0][1] 也会是8

数组生成

可以用过 np.arange 来创建数组,参数与range类似：

x = np.arange(0, 10, 1) # arguments: start, stop, step

也可以用 np.linspace 创建等差数列:

x = np.linspace(1, 10, 5) # arguments: start, stop, num元素个数
# [ 1. 3.25 5.5 7.75 10. ]

# np.logspace 是创建等比数列

矩阵运算

计算将变量直接参与运算符，操作符优先级不变：

a = np.random.rand(5, 5)
b = np.random.rand(5, 5)

print a + b
print a - b
print a * b
print a / b
print a ** 2
print a < b
print a > b

一个数组中除了 dot() 函数，其他这些操作都是单元操作。

np_arr = np.array([2,3,34,5,5])
print np.mean(np_arr) # 平均数
print np.median(np_arr) # 中位数
print np.corrcoef(a[0], a[1]) # 判断两个轴的数据是否有相关性
print np.std(np_arr) # 标准差

数据提取

切片索引语法：M[lower:upper:step]

a = np.array([1,2,3,4,5])
a[1:3] # array([2, 3])

# 进行切片赋值时，原数组会被修改
a[1:3] = [-2, -3] # array([ 1, -2, -3, 4, 5])


b = np.random.rand(5, 5)
b[1:4, 1:4] # 提取 1~4 行，1~4列

b > 0.1 #array([False, False, False, ...])
# 因此要提取可以用, 这是利用了布尔屏蔽这个特性
b[ b > 0.1 ]

# where()函数是另一个有用的方式，当需要以特定条件来检索数组元素的时候。只需要传递给它一个条件，它将返回符合条件的元素列表。
c = np.where(b > 0.1)

矩阵运算

NumPy和Matlab不一样，对于多维数组的运算，缺省情况下并不使用矩阵运算，如果你希望对数组进行矩阵运算的话，可以调用相应的函数。

matrix对象

numpy库提供了matrix类，使用matrix类创建的是矩阵对象，它们的加减乘除运算缺省采用矩阵方式计算，因此用法和matlab十分类似。但是由于NumPy中同时存在ndarray和matrix对象，因此用户很容易将两者弄混。这有违Python的“显式优于隐式”的原则，因此并不推荐在较复杂的程序中使用matrix。

>>> a = np.matrix([[1,2,3],[5,5,6],[7,9,9]])
>>> a*a**-1
matrix([[ 1.00000000e+00, 1.66533454e-16, -8.32667268e-17],
  [ -2.77555756e-16, 1.00000000e+00, -2.77555756e-17],
  [ 1.66533454e-16, 5.55111512e-17, 1.00000000e+00]])

从数组转换为矩阵可以用m = np.matrix(a) 进行转换, 使用 m.T 可以得到m的转置矩阵。

矩阵求逆

m.I * m
=> matrix([[ 1.00000000e+00+0.j, 4.44089210e-16+0.j],
   [ 0.00000000e+00+0.j, 1.00000000e+00+0.j]])

浅拷贝与深拷贝

为了获得高性能，Python 中的赋值常常不拷贝底层对象，这被称作浅拷贝。使用 copy 进行深拷贝：

b = copy(a)

遍历数组元素

通常情况下，我们是希望尽可能避免遍历数组元素的。因为迭代相比向量运算要慢的多。但是有些时候迭代又是不可避免的，这种情况下用 Python 的 for 是最方便的：

v = np.array([1,2,3,4])

for element in v:
 print(element)

M = np.array([[1,2], [3,4]])

for row in M:
 print("row", row)
 for element in row:
  print(element)

总结

以上就是关于Python科学计算之NumPy的全部内容了，希望本文的内容对大家的学习或者工作能带来一定的帮助，如果有疑问大家可以留言交流。

# python  # numpy  # 教程  # python3  # Python中的Numpy入门教程  # Python NumPy库安装使用笔记  # Python Numpy 数组的初始化和基本操作  # Python numpy 提取矩阵的某一行或某一列的实例  # python中找出numpy array数组的最值及其索引方法  # Python编程给numpy矩阵添加一列方法示例  # python中numpy的矩阵、多维数组的用法  # 基于Python Numpy的数组array和矩阵matrix详解  # Python的numpy库中将矩阵转换为列表等函数的方法  # Python中numpy模块常见用法demo实例小结  # 可以用  # 都是  # 这是  # 遍历  # 可以通过  # 情况下  # 用过  # 多维  # 的是  # 组中  # 等比数列  # 等差数列  # 迭代  # 如果你  # 加减乘除  # 也会  # 又是  # 你在  # 很容易  # 这也 

相关文章： 教程网站设计制作软件,怎么创建自己的一个网站？  湖南网站制作公司,湖南上善若水科技有限公司做什么的？  黑客如何利用漏洞与弱口令入侵网站服务器？  网站制作模板下载什么软件,ppt模板免费下载网站？  南阳网站制作公司推荐,小学电子版试卷去哪里找资源好？  专业的网站制作设计是什么,如何制作一个企业网站，建设网站的基本步骤有哪些？  行程制作网站有哪些,第三方机票电子行程单怎么开？  Android滚轮选择时间控件使用详解  如何快速上传建站程序避免常见错误？  建站之星如何助力企业快速打造五合一网站？  如何通过西部数码建站助手快速创建专业网站？  建站主机数据库如何配置才能提升网站性能？  保定网站制作方案定制,保定招聘的渠道有哪些？找工作的人一般都去哪里看招聘信息？  临沂网站制作公司有哪些,临沂第四中学官网？  如何设置并定期更换建站之星安全管理员密码？  学校为何禁止电信移动建设网站？  如何在IIS中新建站点并解决端口绑定冲突？  C#如何使用XPathNavigator高效查询XML  公司网站制作费用多少,为公司建立一个网站需要哪些费用？  建站主机默认首页配置指南：核心功能与访问路径优化  制作网站公司那家好,网络公司是做什么的？  阿里云网站制作公司,阿里云快速搭建网站好用吗？  官网自助建站系统：SEO优化+多语言支持，快速搭建专业网站  在线流程图制作网站手机版,谁能推荐几个好的CG原画资源网站么？  微信小程序制作网站有哪些,微信小程序需要做网站吗？  制作证书网站有哪些,全国城建培训中心证书查询官网？  专业制作网站的公司哪家好,建立一个公司网站的费用.有哪些部分,分别要多少钱？  网站建设制作、微信公众号,公明人民医院怎么在网上预约？  如何在建站主机中优化服务器配置？  电视网站制作tvbox接口,云海电视怎样自定义添加电视源？  如何在局域网内绑定自建网站域名？  建站之星如何防范黑客攻击与数据泄露？  成都网站制作公司哪家好,四川省职工服务网是做什么用？  宁波免费建站如何选择可靠模板与平台？  网站制作大概要多少钱一个,做一个平台网站大概多少钱？  ,网页ppt怎么弄成自己的ppt？  山东云建站价格为何差异显著？  logo在线制作免费网站在线制作好吗,DW网页制作时，如何在网页标题前加上logo？  网站海报制作教学视频教程,有什么免费的高清可商用图片网站，用于海报设计？  临沂网站制作企业,临沂第三中学官方网站？  制作网站的网址是什么,请问后缀为.com和.com.cn还有.cn的这三种网站是分别是什么类型的网站？  如何通过山东自助建站平台快速注册域名？  如何用景安虚拟主机手机版绑定域名建站？  微网站制作教程,我微信里的网站怎么才能复制到浏览器里？  简历在线制作网站免费版,如何创建个人简历？  网站制作壁纸教程视频,电脑壁纸网站？  Python文件管理规范_工程实践说明【指导】  如何配置支付宝与微信支付功能？  魔方云NAT建站如何实现端口转发？  高防服务器：AI智能防御DDoS攻击与数据安全保障